Pregnancy status in harbour seals can be estimated from concentrations of progesterone in blubber as well as in blood samples, which are significantly higher in pregnant than non-pregnant animals. This study investigated the accuracy of estimating pregnancy rates using samples from live-captured and released harbour seals from three regions around Scotland, coupled with observed pregnancy outcomes. Concentrations of progesterone in blood (plasma) and blubber were obtained during the capture of animals early in the year (February to May). Individual animals were identified from the unique markings on their pelage, with a proportion (n = 51) of females re-sighted during the subsequent breeding season and the reproductive outcomes determined (pregnant or possibly non-pregnant) during observations from long-term photo-identification studies. Generalised linear models with a binomial link function were fitted to training (60% of the data) and test datasets (40% of the data) to estimate pregnancy status from progesterone concentrations in blubber, plasma or both, and a received operating curves (ROC) approach was used to evaluate the performance of each classifier. The accuracy for the plasma concentrations was 85% with a high classification performance (as estimated from an area under the curve (AUC) of 0.82). The Youden method to determine the cut-point (threshold) and bootstrapping the training dataset resulted in a cut-point of 58 ng ml−1 (95th percentiles, 25–102 ng ml−1). For blubber, the accuracy was 77% (AUC = 0.86) with an optimal cut-point of 56 ng g−1 (95th percentiles, 26–223 ng g−1). In the combined analysis (both blubber and plasma), the accuracy was 87.5% (AUC 0.81) with the cut-points of 72 ng ml−1 (95th percentiles, 25–103 ng ml−1) in plasma and 56 ng g−1 (95th percentiles, 26–223 ng g−1) in blubber. These thresholds were then used to estimate the pregnancy proportions among adult females at the three study sites, including those that were not included in the photo-id studies. Proportions were high at all sites, (63%–100%) regardless of which matrices were used and were not statistically significantly different from each other but suggested that analysing concentrations in both sample matrices would minimise the uncertainty.