127 ESTIMATION OF MOOSE PARTURITION DATES IN COLORADO: INCORPORATING IMPERFECT DETECTIONS Eric J. Bergman1, Forest P. Hayes2, and Kevin Aagaard1 1Colorado Parks and Wildlife, 317 W. Prospect Ave., Fort Collins, CO 80521, USA; 2Wildlife Biology Program, University of Montana, Forestry 108, 32 Campus Drive, Missoula, MT 59812, USA. ABSTRACT: Researchers and managers use productivity surveys to evaluate moose populations for harvest and population management purposes, yet such surveys are prone to bias. We incorporated detection probability estimates (p) into spring and summer ground surveys to reduce the influence of observer bias on the estimation of moose parturition dates in Colorado. In our study, the cumulative parturition probability for moose was 0.50 by May 19, and the probability of parturition exceeded 0.9 by May 27. Timing of moose calf parturition in Colorado appears synchronous with parturition in more northern latitudes. Our results can be used to plan ground surveys in a manner that will reduce bias stemming from unobservable and yet-born calves. ALCES VOL. 56: 127–135 (2020) Key Words: Alces, calf-at-heel, detection probability (p), ground surveys, parturition, recruitment Throughout North America and Europe, researchers and managers use surveys of moose productivity to evaluate populations for harvest management purposes (Boertje et al. 2007, Grøtan et al. 2009, Milner et al. 2013); however, surveys are prone to bias (Williams et al. 2001, White 2005). When surveying for newborn moose calves, one source of bias is associated with the detection probability (p) of moose calves-at-heel (Bergman et al. 2020). More specifically, if a calf is observed, then p is conceptually 1 for that individual during that occasion. However, if a calf is not observed, then uncertainty about its presence exists (i.e., was the calf simply not observed, or was there no calf to be observed). If surveys are conducted near the peak time of parturition, this bias is confounded by the possibility that cows may have not yet given birth. Calf-at-heel estimates are also prone to bias as calf mortality occurs. However, multiplying monthly (or daily) calf survival rates by calf-at-heel ratios provides a numerical correction for bias stemming from calf mortality (Bergman et al. 2020). No simple multiplicative, numerical correction exists for pre-parturition observations. Fortunately, accounting and accommo- dating for many types of bias is possible in both modelling and survey design. First, estimates of p can be modelled from repeated observations (Bergman et al. 2020). Once estimated, p is used to inflate calf-at-heel or calf:cow ratios to reduce bias in estimates. An example of such an approach was com- pleted by Bergman et al. (2020) who used occupancy modelling and 3 years of ground observation data from radio-collared cow moose to generate a summertime estimate of p = 0.80. We suggest that a supplemental approach to reducing bias stemming from unborn calves is to quickly and efficiently conduct calf-at-heel surveys after the bulk of parturition has occurred. PARTURITION DATES IN COLORADO – BERGMAN ET AL. ALCES VOL. 56, 2020 128 Under ideal conditions and with modern technology, timing of moose parturition can be estimated with minimal uncertainty. For instance, Vaginal Implant Transmitters (VIT) are used to alert researchers to the timing and location of a birthing event (Patterson et al. 2013, 2016, McLaren et al. 2017). This approach is often used when the objective is to capture and collar newborn calves. However, it requires capturing adult females to assess pregnancy status and to deploy a VIT. The recent development of satel- lite-based VITs minimizes the previous need for daily ground or aerial monitoring, but the technology remains cost prohibitive for most routine management purposes. A second and increasingly tractable approach to estimate parturition dates and locations for large herbivores is also tied to satellite technology. Satellite collars now allow researchers to shorten the duration between sequential locations of animals and achieve nearly real-time transmission of data. Movement algorithms, or even close scrutiny of sequential data points can be used to identify clustered locations that are often indicative of birthing events (Severud et al. 2015, McLaren et al. 2017, Cameron et al. 2018). However, neither traditional VHF radio-collars, store-on-board GPS collars, nor early generation satellite collars provide the frequency of locations and the nearly real-time transmission of data necessary to identify birth sites. Our objectives for this research were twofold. First, using ground observation data, we estimated a parturition date curve for moose in Colorado. Managers in Colorado and elsewhere will benefit from estimates of the timing of parturition made more precise by incorporation of p, such that they can implement recruitment surveys when a threshold (such as >90%) of birthing events is predicted. Our second objective was to correct estimates of parturition timing for p in this modelling process, thereby improving the precision of parturition date estimates. We hypothesized that accounting for p would shift the date of cumulative births to an earlier date, thereby allowing managers to initiate calf surveys at an earlier date without pre-parturition bias. STUDY AREA We conducted this research across 3 study areas in Colorado. The 2 most north- erly were located in Jackson (North Park) and Larimer (Laramie River) Counties, with the southern study area (San Juan Mountains) in Hinsdale and Mineral Counties (Fig. 1). North Park was a high elevation (2,400–2,750 m), wide (14–46 km) moun- tain valley surrounded on the west by the Park Range mountains, on the south by the Rabbit Ears mountain range, and on the east by the Rawah and Never Summer mountain ranges. To the north of the study area was a mix of private and public lands managed pri- marily for agricultural and open rangeland purposes. Moose habitat in North Park fol- lowed small rivers and creeks comprised of a diversity of willow (Salix spp.) communi- ties. Moose also used aspen (Populus tremu- loides), lodgepole pine (Pinus contorta), and Englemann spruce (Picea engelmannii) for- ests. Much of the pine and spruce forests in North Park and throughout Colorado experi- enced mountain pine beetle (Dendroctonus ponderosae) outbreaks during the latter part of the 20th and first decade of the 21st cen- tury, placing these forests into an array of successional stages (Hayes 2020). The Laramie River study area was located ~ 40 km northeast of North Park with the Rawah mountain range (3,200–3,840 m) separating them. It was also a high elevation mountain valley (2,470–2,800 m), although the valley floor was not as wide as North Park (3–9 km). Diverse willow stands located along the rivers and creek corridors gave way ALCES VOL. 56, 2020 PARTURITION DATES IN COLORADO – BERGMAN ET AL. 129 to more upland aspen, lodgepole pine, and Englemann spruce forests. The San Juan Mountains study area in southern Colorado at 2,750–3,130 m ele- vation was higher than the North Park and Laramie River study areas. It was comprised of narrow valleys (0.5–1.5 km wide) with vegetation communities similar to those in the northern study areas. Management authority of moose belonged to Colorado Parks and Wildlife (CPW) in all 3 study areas, and each sus- tained limited cow and bull harvest. Predator assemblages were consistent across study areas with black bears (Ursus americanus) and mountain lions (Puma concolor) the pri- mary predators of moose, although coyotes (Canis latrans) could potentially kill new- born moose calves; wolves (Canis lupus) and grizzly bears (Ursus arctos horribilis) were absent. Predation pressure was consid- ered low in each study area. METHODS Field methods We captured cow moose (>1 year-old) via helicopter darting for 4 winters between mid-December and the end of January, 2015–2018. We sedated moose using one of three different drug combinations: 1) BAM (54.6 mg of butorphanol, 18.2 mg of azaper- one, and 21.8 mg of medetomidine) in combination with ketamine (200 mg), 2) carfentanil (3 mg) in combination with xyla- zine (100 mg), or 3) thiafentanil (10 mg) in combination with xylazine (25 mg). Once sedated, we blindfolded each animal and administered oxygen (via nasal canula) to minimize the risk of adult and fetal hypoxia. We fitted moose with either a VHF radio- collar (Advanced Telemetry Systems, Isanti, Minnesota, USA; USA model: M2520B), a store-on-board GPS/VHF collar ([Advanced Telemetry Systems; USA model: G2110D], or a satellite/VHF telemetry collar Fig. 1. Map of Colorado, USA (black rectangular perimeter) depicting 3 study areas in relation to nearby cities and communities. Study units are depicted by gray filled polygons. PARTURITION DATES IN COLORADO – BERGMAN ET AL. ALCES VOL. 56, 2020 130 [Vectronics Aerospace GmbH, Berlin, Germany; model: Vertex Plus, and Advanced Telemetry Systems; USA model: G5-2D]). Blood samples were taken to determine pregnancy status using Pregnancy Specific Protein B (PSPB, Wood et al. 1986). After handling, capture drugs were antagonized with naltrexone (100 mg, antagonist for carfentanil and thiafentanil), tolazoline (500 mg, antagonist for azaperone and xyla- zine), and atipamezole (100–150 mg, antag- onist for medetomidine and xylazine). All capture and handling methods were approved by the Institutional Animal Care and Use Committees at Colorado Parks and Wildlife (#08-2013) and the University of Montana (#032-17CBWB-060517). During the first year, no moose had been previously captured or collared. In subse- quent years, previously collared moose were neither targeted nor avoided by the capture crew. As a result, this random process meant that the pregnancy status of some collared moose was unknown. Each spring and summer we conducted ground surveys to evaluate the calf-at-heel status of each collared cow. We began obser- vations in early May and continued through August. Pregnant moose, based on PSPB results at the time of capture, were priori- tized for observation. Once these animals were observed, we completed observations of radio-collared animals with unknown pregnancy status. Typically, ground observa- tions were completed by a single observer by relying on previous known locations and using VHF signals for ground tracking. A second observer was used for individual moose that consistently evaded observation by a single observer. In cases with two observers, one homed in on the moose using the described techniques, with the second observer stationed along the expected exit route with the goal of observing the moose as it passed by. Moose observations typically fell into 2 categories: stationary or moving. Stationary observations were made of moose that were either bedded or standing idly while they foraged. Stationary observations often lasted from 5 to 20 min and ended when a moose stood and moved or foraged out of sight. Moving observations were those of moose displaced by an observer. To be considered a completed observation, observers needed to see the entire moose and the surrounding 2 m of space. Repeat obser- vations were made on animals throughout the summer to increase the detection proba- bility of calves, and to help determine the fate of calves. We recorded date, time, and location of each observation. Analytical methods Our objective was to estimate the partu- rition date for moose in Colorado. Thus, only cow moose that were eventually observed with calves-at-heel were included in analy- ses. Cow moose that were never observed with a calf helped inform calf-at-heel and calf:cow ratios (Bergman et al. 2020), but did not inform estimates of birth dates. The date of each observation was standardized against the date of the earliest survey (27 April), allowing for simple numerical progression throughout the survey period. We used a hierarchical Bayesian model to evaluate the probability of parturition during the study period (McCarthy 2007, Gelman et al. 2009). The hierarchical com- ponent refers to the multiple levels included in the model, which are ultimately integrated to estimate posterior estimates for each parameter of interest (Gelman et al. 2009). The base model for parturition probabil- ity included an estimate of the probability of detection and followed a logistic regression with a “logit” link. The model had the following form: logit (ρ i ) ~ α + (β × ϑ i ) + γδi + γτi, ALCES VOL. 56, 2020 PARTURITION DATES IN COLORADO – BERGMAN ET AL. 131 where α is the global intercept, ϑ is the date of observation (with corresponding coeffi- cient β), and random effects (γ) of year (δ) and cow (τ) for each cow, i. The predicted calf presence, ρ, was influenced by an esti- mated probability of detection 0.8 (Bergman et al. 2020) and followed a Bernoulli distri- bution, modeled as: Y i ~ Bern (φ i ); φ i = μ × ω i ; ω i ~ Bern (ρ i ), where the observed calf detections (Y) fol- low a Bernoulli distribution with probability (φ) informed by the product of the detection probability (μ) and estimated calf detection (ω, per cow i). This is standard practice for including detection probability in Bayesian models (McCarthy 2007). The random effects were given vague normal priors with uniform precision (inverse of variance): θ θ σ σ  =  N U 0, ; 1 ; ~ 0, 10 . 2 2 The global intercept and β coefficient were given vague normal priors: N [0, 1.0 × 10−6]. We ran the model using the “runjags” pack- age (Denwood 2016) in R (R Core Team 2019) including 3 chains, with 10,000 iter- ations per sample and a burn-in of 5,000 iterations and a thinning parameter of 5. We determined convergence when R-hat < ~1.1 for monitored parameters (Gelman and Rubin 1992). RESULTS We captured 46 individual cow moose that were observed with spring or summer calves-at-heel, providing for 86 unique animal-by-year observations (i.e., some cows were observed multiple years). We made a total of 213 unique observations of these individuals. Within a single year, the minimum number of observations of an indi- vidual moose was 1 (when a cow was observed with a calf during the first observa- tion and subsequent observations were not made), and the maximum was 5. We made an average of 1.72 (SD = 0.96) observations of each cow. Our earliest survey was on 27 April 2016 (day 1), and our earliest observa- tion of a calf was on 17 May 2016. Based on raw observation data, the median annual parturition dates ranged from 8 June through 18 July; however, these dates reflect uncorrected adjustments. As expected, accounting for p shifted dates for the predicted probability of parturition to an earlier period. The cumulative parturition probability for moose was 0.50 by day 22 (19 May, with a 95% credible interval [CI] range from 19 to 20 May). In consideration of cumulative parturition among all moose, without correcting for detection probability, the probability of parturition exceeded 0.90 on day 64 (June 30, 95% CI = 20 June to 19 July; Fig. 2). When detection probability was incorporated into the model, the cumu- lative probability of parturition exceeded 0.90 by 30 days after 27 April (27 May; 95% CI = 26 May to 28 May; Fig. 2). DISCUSSION Evolutionary theory suggests that for many large, northern ungulates, the peak and duration of parturition periods are shaped to occur within a narrow window of time (Rutberg 1987). But because our ground-based field methods to estimate parturition dates were laborious and obser- vation rates low (0–5/day), the date of first observation for many cows extended well into summer. However, analytical adjust- ments to the estimation of parturition dates PARTURITION DATES IN COLORADO – BERGMAN ET AL. ALCES VOL. 56, 2020 132 (i.e., accounting for p) buffered the bias associated with our slower field methods and led to the prediction that 50% of partu- rition events had occurred by 19–20 May. This 2-day window aligned very closely with the range of median parturition date of 19–22 May reported by Gasaway et al. (1983) and Keech et al. (2000) for the inte- rior of Alaska. Similarly, Bowyer et al. (1998) reported a mean parturition date of 25 May for moose in Denali National Park, and concluded that 95% of births occurred during a 16-day window. The median partu- rition date for moose calves in southwest Yukon was also 25 May (Larsen et al. 1989). Our results also aligned with parturi- tion dates for moose in the eastern United States and Scandinavia. In New Hampshire, Musante et al. (2010) and Jones et al. (2017) reported a median date of 19 May with 78 and 90% of births occurring between 13 and 27 May, respectively. Parturition dates in Norway were also simi- lar (23 May), but dates were sensitive to the number of mature bulls in the population (Sæther et al. 2003). Finally, Severud et al. (2015) reported a slightly earlier mean par- turition date (14 May) for moose in Minnesota, but a 1-month range of parturi- tion (2 May–2 June). This earlier mean par- turition date aligned with that reported in Ontario (13 May; Patterson et al. 2016). We estimated that the cumulative proba- bility of parturition increased from 0.50 to Fig. 2. Predicted probability of parturition by date (shown as days since 27 April), modeled with (solid) and without (dashed) including the probability of detection, Colorado, USA. The horizontal dotted-dashed line indicates a 90% parturition probability. The vertical lines indicate the days on which that 90% parturition probability was estimated to have been achieved for each model (dotted lines represent 95% credible intervals). Raw observational data are depicted as black dots. ALCES VOL. 56, 2020 PARTURITION DATES IN COLORADO – BERGMAN ET AL. 133 0.90 between 19 and 27 May, indicating that Colorado has a similarly narrow parturition period as reported across much of moose range. This narrow window may be shaped by the interaction of habitat and season (Rutberg 1987, Bowyer et al. 1998), as well as predation (Bergerud 1975, Testa 2002). Perhaps less intuitive was that Colorado’s moose appear to calve in synchrony with moose at more northern latitudes. In com- parison, the onset of spring and summer is generally earlier and winter later in Colorado. This variation in seasonality could poten- tially afford moose in Colorado and other southern populations flexibility from tight parturition patterns identified in northern populations; however, no shifts in parturi- tion date are apparent. While the seasonality of Colorado’s southern latitude may be mediated by high elevations, the parturition synchrony within the species across latitudes may prove to be relevant and informative in the face of a generally warming environ- ment. More specifically, moose occupy a wide geographical and latitudinal range, over which seasons are not perfectly syn- chronous. Yet, they apparently maintain tight synchrony in the timing of parturition across this range. Thus, concerns over the shifting of seasonality due to global warming (i.e., earlier spring and delayed winter) may not intrinsically, or negatively impact the timing of parturition. From a management perspective, the estimation of parturition dates in Colorado was particularly useful to design field surveys. One goal of refining productivity surveys is to reduce bias, and as noted, one source of bias is p and its confounding effects when moose are transitioning between pregnancy and calf-at-heel. Ideally, surveys should be implemented post-parturition, but early enough that neonatal mortality is minimal. After applying the probability of detection at which the probability of parturition reached 90%, we recommend that surveys in Colorado be initiated on 27 May, or 34 days prior to the date predicted without considering probability of detection. In addition, incorporating probability of detection decreased the credible interval (by about 28 days) associated with the predicted date at which 90% of parturition events occurred. While ground surveys cannot fully replicate the results of aerial surveys, managers can use our results to improve and facilitate the timing of ground surveys. For example, a concerted ground survey effort at the end of May, conducted typically on foot or horseback with a large number of volunteers and field personnel, should produce a productivity estimate with minimal bias from pregnant females and calf mortality. Importantly, the narrower credible interval indicates that earlier surveys need not compromise confidence, and applying our refined, optimal date to initiate earlier surveys will produce measurable savings in effort and agency resources. ACKNOWLEDGENTS This research was funded in part by a United States Fish and Wildlife Service Federal Aid Research Grant, CPW game cash fund, and auction and raffle grants administered by CPW. We are also indebted to the efforts of S. Boyle, B. Cimpher, R. Cordova, A. Howell, A. Maclean, S. Peterson, B. Smith, and K. Yeager who spent many hours tracking and observing moose. Early drafts of this manuscript were improved by the comments provided by M.Alldredge, D. Tripp, M. Lankester, and 2 anonymous reviewers. REFERENCES Bergerud, A. T. 1975. The reproductive sea- son of Newfoundland caribou. Canadian Journal of Zoology 53: 1213–1221. doi: 10.1139/z75-145 PARTURITION DATES IN COLORADO – BERGMAN ET AL. ALCES VOL. 56, 2020 134 BergmAn, e. J., F. P. HAyes, P. m. LukAcs, and c. J. BisHoP. 2020. Moose calf detection probabilities: quantification and evaluation of a ground-based sur- vey technique. Wildlife Biology. doi: 10.2981/wlb.00599 BoerTJe, r. d., k. A. keLLie, c. T. seATon, m. A. keecH, d. d. young, B. W. dALe, L. g. AdAms, and A. r. AdermAn. 2007. Ranking Alaska moose nutrition: signals to begin liberal antlerless harvests. Journal of Wildlife Management 71: 1494–1506. doi: 10.2193/2006-159 BoWyer, r. T., V. VAn BALLenBergHe, and J. g. kie. 1998. Timing and synchrony of parturition in Alaskan moose: long-term versus proximal effects of climate. Journal of Mammalogy 79: 1332–1344. doi: 10.2307/1383025 cAmeron, m. d., k. JoLy, g. A. Breed, L. s. PArreTT, and k. kieLLAnd. 2018. Movement-based methods to infer par- turition events in migratory ungulates. Canadian Journal of Zoology 96: 1187–1195. doi: 10.1139/cjz-2017-0314 denWood, m. J. 2016. Runjags: an R Package providing interface utilities, model tem- plates, parallel computing methods and additional distributions for MCMC models in JAGS. Journal of Statistical Software 71: 1–25. doi: 10.18637/jss. v071.i09 gAsAWAy, W. c., r. o. sTePHenson, J. L. dAVis, P. e. k. sHePHerd, and o. e. Burris. 1983. Interrelationships of wolves, prey, and man in interior Alaska. Wildlife Monographs 84: 1–50. geLmAn, A., J. B. cArLin, H. s. sTern, and d. B. ruBin. 2009. Bayesian Data Analysis, 2nd Edition. Chapman & Hall Press, Boca Raton, Florida, USA. _____, and d. B. ruBin. 1992. Inference from iterative simulation using multiple sequences. Statistical Science 7: 457–511. doi: 10.1214/ss/1177011136 grøTAn, V., B.-e. sæTHer, m. LiLLegård, e. J. soLBerg, and s. engen. 2009. Geographical variation in the influence of density dependence and climate on the recruitment of Norwegian moose. Oecologia 161: 685–695. doi: 10.1007/ s00442-009-1419-5 HAyes, F. P. 2020. Resource selection and calving success of moose in Colorado. M. S. Thesis, University of Montana, Missoula, Montana, USA. Jones, H., P. J. Pekins, L. e. kAnTAr, m. o’neiL, and d. eLLingWood. 2017. Fecundity and summer calf survival of moose during 3 successive years of win- ter tick epizootics. Alces 53: 85–98. keecH, m. A., r. T. BoWyer, J. m. Ver HoeF, r. d. BoerTJe, B. W. dALe, and T. r. sTePHenson. 2000. Life-history conse- quences of maternal condition in Alaskan moose. Journal of Wildlife Management 64: 450–462. doi: 10.2307/3803243 LArsen, d. g., d. A. gAuTHier, and r. L. mArkeL. 1989. Causes and rate of moose mortality in southwest Yukon. Journal of Wildlife Management 53: 548–557. doi: 10.2307/3809175 mccArTHy, m. A. 2007. Bayesian Methods for Ecology. Cambridge University Press, Cambridge, United Kingdom. mcLAren, A. A. d., J. F. Benson, and B. r. PATTerson. 2017. Multiscale habitat selection by cow moose (Alces alces) at calving sites in central Ontario. Canadian Journal of Zoology 95: 891–899. doi: 10.1139/cjz-2016-0290 miLner, J. m., F. m. VAn BeesT, e. J . soLBerg, and T. sTorAAs. 2013. Reproductive success and failure: the role of win- ter body mass in reproductive allocation in Norwegian moose. Oecologia 172: 995–1005. doi: 10.1007/s00442-012- 2547-x musAnTe, A. r., P. J. Pekins, and d. L. scArPiTTi. 2010. Characteristics and dynamics of a regional moose Alces alces population in northeastern United States. Wildlife Biology 16: 185–204. doi: 10.2981/09-014 PATTerson, B. r., J. F. Benson, k. r. middeL, k. J. miLLs, A. siLVer, and m. e. oBBArd. ALCES VOL. 56, 2020 PARTURITION DATES IN COLORADO – BERGMAN ET AL. 135 2013. Moose calf mortality in central Ontario, Canada. Journal of Wildlife Management 77: 832–841. doi: 10.1002/ jwmg.516 _____, k. J. miLLs, k. r. middeL, J. F. Benson, and m. e. oBBArd. 2016. Does predation influence the seasonal and diel timing of moose calving in central Ontario, Canada? PLoS One 11: e0150730. doi: 10.1371/ journal.pone.0150730 r core TeAm. 2019. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project. org (accessed April 2020). ruTBerg, A. T. 1987. Adaptive hypotheses of birth synchrony in ruminants: an inter- specific test. The American Naturalist 130: 692–710. doi: 10.1086/284739 sæTHer, B.–e., e. J. soLBerg, and m. Heim. 2003. Effects of altering sex ratio struc- ture on the demography of an isolated moose population. Journal of Wildlife Management 67: 455–466. doi: 10.2307/ 3802703 seVerud, W. J., g. deLgiudice, T. r. oBermoLLer, T. A. enrigHT, r. g. WrigHT, and J. d. ForesTer. 2015. Using GPS collars to determine parturition and cause-specific mortality of moose calves. Wildlife Society Bulletin 39: 616–625. doi: 10.1002/wsb.558 TesTA, J. W. 2002. Does predation on neo- nates inherently select for earlier births? Journal of Mammalogy 79: 1332–1344. WHiTe, G. C. 2005. Correcting wildlife counts using detection probabilities. Wildlife Research 32: 211–216. doi: 10.1071/WR03123 WiLLiAms, B. k., J. d. nicHoLs, and m. J. conroy. 2001. Analysis and Management of Animal Populations. Academic Press, Inc. San Diego, California, USA. Wood, A. k., r. e. sHorT, A. dArLing, g. L. dusek, r. g. sAsser, and c. A. ruder. 1986. Serum assays for detecting preg- nancy in mule and white-tailed deer. Journal of Wildlife Management 50: 684–687. doi: 10.2307/3800981 https://www.R-project.org https://www.R-project.org