Determining Timing of Births and Habitat Selection to Identify Lambing Period Habitat for Bighorn Sheep

Timing of births in ungulates is influenced by ecological factors, and differences in seasonality of births have evolutionary implications for these mammals. Birthing habitat is one of the most important home-range attributes for ungulates, and disturbances during this time can decrease survival of young and population growth. We calculated timing of births and quantified habitat selected by nursery groups (i.e., females with young) of desert bighorn sheep (Ovis canadensis nelsoni) to produce and project a geographic information system (GIS) model of lambing period habitat (i.e., birthing locations and areas used up to 6 weeks post-parturition) in southeastern Utah, United States. We then applied that model to identify suitable lambing period habitat in an adjacent area for a population of Rocky Mountain bighorn sheep (O. c. canadensis). We monitored 19 global positioning system (GPS) collared females from 2012 to 2013 to delineate the birthing season. We used GPS locations during that time and logistic regression within a model-selection framework to differentiate between lambing period habitat and random locations based on habitat and anthropogenic covariates. We used model coefficients to produce and project a GIS model of lambing period habitat. Across both years, 89% of 45 births occurred in May (range = April 27–June 9). We quantified covariate values at 750 lambing period and 750 random locations during that time in the San Rafael study area. The top model included elevation, slope, terrain ruggedness, aspect, existing vegetation type, and distance to trails and roads. Those variables predicted lambing period habitat for desert bighorn sheep (rho = 0.99, P = 0.02). We then overlaid that model on 1,130 GPS locations from 17 female Rocky Mountain bighorns in the Green River study area. In both areas, a mean of 71% of locations were in the top two (of five) probability categories (San Rafael high = 37%, medium high = 31%; Green River high = 50%, medium high = 23%). Females in lambing period habitat preferred steep, north-facing slopes, rugged terrain, low elevation, and avoided roads. Our GIS model projected areas of lambing period habitat for adjacent desert and Rocky Mountain bighorns and provided land managers with a map of habitat in areas where resource extraction and recreation are increasing. Identifying timing of parturition, and high-probability lambing period habitat, will help managers mitigate temporal and spatial conflicts between bighorn sheep and anthropogenic activities. Similar data regarding timing of births and a GIS model can be applied to conserve habitat for other ungulates also.

Timing of births in ungulates is influenced by ecological factors, and differences in seasonality of births have evolutionary implications for these mammals. Birthing habitat is one of the most important home-range attributes for ungulates, and disturbances during this time can decrease survival of young and population growth. We calculated timing of births and quantified habitat selected by nursery groups (i.e., females with young) of desert bighorn sheep (Ovis canadensis nelsoni) to produce and project a geographic information system (GIS) model of lambing period habitat (i.e., birthing locations and areas used up to 6 weeks post-parturition) in southeastern Utah, United States. We then applied that model to identify suitable lambing period habitat in an adjacent area for a population of Rocky Mountain bighorn sheep (O. c. canadensis). We monitored 19 global positioning system (GPS) collared females from 2012 to 2013 to delineate the birthing season. We used GPS locations during that time and logistic regression within a model-selection framework to differentiate between lambing period habitat and random locations based on habitat and anthropogenic covariates. We used model coefficients to produce and project a GIS model of lambing period habitat. Across both years, 89% of 45 births occurred in May (range = April 27-June 9). We quantified covariate values at 750 lambing period and 750 random locations during that time in the San Rafael study area. The top model included elevation, slope, terrain ruggedness, aspect, existing vegetation type, and distance to trails and roads. Those variables predicted lambing period habitat for desert bighorn sheep (rho = 0.99, P = 0.02).
We then overlaid that model on 1,130 GPS locations from 17 female Rocky Mountain bighorns in the Green River study area. In both areas, a mean of 71% of locations were in the top two (of five) probability categories (San Rafael high = 37%, medium high = 31%; Green River high = 50%, medium high = 23%). Females in lambing period habitat preferred steep, north-facing slopes, rugged terrain, low elevation, and avoided roads. Our GIS model projected areas of lambing period habitat for adjacent desert and Rocky Mountain bighorns and provided land managers with a map of habitat in areas where resource extraction and recreation are increasing. Identifying timing of

INTRODUCTION
Timing of births in ungulates is influenced by many ecological factors (Sadleir, 1969;Bowyer, 1991;Loe et al., 2005), and differences in seasonality of births for these mammals can have evolutionary implications (Marshall and Cambridge, 1937;Rutberg, 1987;Ims, 1990), such as reduced survival of young and future reproductive potential (Clutton-Brock et al., 1987;Festa-Bianchet et al., 2000;Keech et al., 2000). For bighorn sheep (Ovis canadensis), latitude, elevation, growing season length, climate, nutrition, and photoperiod are factors influencing reproductive seasonality (Bunnell, 1982;Thompson and Turner, 1982;Whiting et al., 2012). Generally, bighorn sheep occupying northern latitudes and higher elevations give birth late in spring, because of a constricted birthing period and shortened growing season (Bunnell, 1982;Thompson and Turner, 1982). Conversely, bighorn sheep occupying southern latitudes can give birth during most months, likely because growing seasons are much less predictable (Lenarz, 1979;Thompson and Turner, 1982;Rubin et al., 2000). Understanding the timing of births for ungulates can guide conserving habitat by mitigating temporal conflicts between these mammals and anthropogenic activities (e.g., recreation and mining) during that time (Singh et al., 2010;Kaze et al., 2016).
Female ungulates select birth sites and habitat for nursery groups (i.e., females with young) based on a variety of trade-offs (Bowyer, 1991;Bowyer, 1991, 1994). These tradeoffs include predation risk, exposure to extreme weather, and forage quality and availability for the mother to fulfill nutritional requirements during late gestation and lactation (Festa-Bianchet, 1988b;Berger, 1991;Severud et al., 2019). For example, female bighorn sheep often return to the same general area each year to give birth (Geist, 1971;Etchberger and Krausman, 1999;Whiting et al., 2012). Young are usually born in relatively flat areas of high elevation in rugged, steep terrain close to perennial water, on south and west facing slopes, and away from anthropogenic disturbance (Geist, 1971;Bangs et al., 2005;Smith et al., 2015). After giving birth, female bighorns generally form nursery groups with other females and young in rugged areas of high elevation and steep slopes (Geist, 1971;Bangs et al., 2005;Karsch et al., 2016). These groups then move cohesively among patches of nursery habitat for the subsequent months (Bangs et al., 2005;Whiting et al., 2012;Wiedmann and Bleich, 2014). In general, Rocky Mountain bighorn sheep (O. c. canadensis) give birth later in spring, and the birthing period is constricted allowing females to exploit the shortened growing season during summer (Bunnell, 1982;Thompson and Turner, 1982). Conversely, desert bighorn sheep (O. c. nelsoni) give birth in almost any month, ostensibly because of a less predictable growing season (Bunnell, 1982;Thompson and Turner, 1982;Rubin et al., 2000). Delineation of birthing habitat has important implications for the conservation, reproductive biology, and perpetuation of bighorn populations (Etchberger and Krausman, 1999;Wiedmann and Sargeant, 2014;Smith et al., 2015), and is a topic that is receiving increased attention for land-use planning (Severud et al., 2019).
Anthropogenic activities can influence habitat selection by female ungulates during birthing and when animals congregate in nursery groups (Stankowich, 2008;Dzialak et al., 2011;Kaze et al., 2016), which can reduce forage intake and suppress population growth (Ciuti et al., 2012b;Sproat et al., 2019). For example, Saiga antelope (Saiga tatarica) avoided areas of human use while selecting birthing sites (Singh et al., 2010). Parturient bison (Bison bison) selected birthing areas that were away from recreational trails, roads, and buildings (Kaze et al., 2016). In an area that was intentionally disturbed by humans during parturition, counts of young/female elk (Cervus canadensis) were 22.5% lower than in areas with undisturbed herds (Phillips and Alldredge, 2000). Populations of bighorn sheep can be negatively influenced by anthropogenic disturbances, especially during the birthing season (Papouchis et al., 2001;Wiedmann and Bleich, 2014;Smith et al., 2015). These disturbances can cause females to abandon previously used nursery habitat. Thereafter, females use less-suitable habitat (Longshore et al., 2013), which may increase predation risk for lambs (Papouchis et al., 2001), and lower recruitment of young, resulting in a declining population (Papouchis et al., 2001;Wiedmann and Bleich, 2014). Delineating birthing habitat and habitat for nursery groups of ungulates can reduce spatial conflict between these mammals and anthropogenic activities (Kaze et al., 2016;McLaren et al., 2017).
Conservation and land-use planning will increasingly rely on identifying preferred wildlife habitats in areas of varying levels of human use (Margules and Pressey, 2000;Dzialak et al., 2011;Harju et al., 2011). Biologists can reduce and mitigate disturbance from anthropogenic activities to female ungulates during birthing by identifying when and where animals give birth Kaze et al., 2016;McLaren et al., 2017). Determining timing of parturition in ungulates to identify high probability nursery areas can help managers mitigate temporal and spatial conflicts (Kaze et al., 2016;McLaren et al., 2017), and guide habitat management (Severud et al., 2019). Global positioning system (GPS) data can effectively document such patterns at fine scales (Longshore et al., 2013), especially for parturient females (Smith et al., 2015;McLaren et al., 2017;Severud et al., 2019). We calculated timing of births to delineate habitat for nursery groups of GPS-collared desert bighorn sheep and Rocky Mountain bighorn sheep in two of the largest populations of those subspecies in Utah, United States. We used GPS locations during that time and logistic regression within a model-selection framework to differentiate between lambing period habitat (i.e., birthing locations and areas used up to 6 weeks post-parturition; Zeigenfuss et al., 2000) and random locations based on habitat and anthropogenic covariates. We then constructed a GIS model that projected lambing period habitat in both areas. That model provided biologists and managers with a map of high probability lambing period habitat. By identifying timing of parturition, and high-probability lambing period habitat, managers will be able to mitigate temporal and spatial conflicts between bighorn sheep and recreation, mining, and domestic livestock grazing. Data regarding timing of births and a similar GIS model can be applied to manage and conserve habitat for ungulate species in other areas as well.

Study Areas
The population of desert bighorn sheep we studied occupied the North San Rafael Swell (hereafter, San Rafael). The San Rafael is located in Emery County, Utah (Figure 1; 38 • 58 N, 110 • 37 W). Bighorns were native and abundant in the San Rafael, but were likely extirpated from that area; the last confirmed sighting occurred in 1964 (Dalton and Spillett, 1971). That population was re-established in the 1970s and 1980s with five translocations of 57 bighorns from Canyonlands National Park and the San Juan Population, Utah (Utah Division of Wildlife Resources [UDWR], 2018). Those translocation efforts were successful, and in 2001 the UDWR estimated 543 animals in that population. However, bighorn sheep started to decline in the 2000s, and 1 month prior to our study, biologists estimated 143 individuals in that area.
The San Rafael is characterized by steep canyons in the Wingate Formation with broad mesa tops in Navajo and Entrada Sandstone formations (Rigby and Beus, 1987). Desert bighorn sheep habitat in that area ranged in elevation from 1,700 to 2,100 m. Vegetation consisted of species typical of salt desert shrub environments. Common shrubs included blackbrush (Coleogyne ramosissima) and fourwing saltbush (Atriplex canescens). Pinyon pine (Pinus edulis) and juniper (Juniperus spp.) were predominate on mesa tops and on north facing canyon slopes. The San Rafael is dry with annual precipitation averaging <20 cm per year. Daily high temperatures during summer (June-September) averaged 31 • C and often exceeded 35 • C. Winters (November-February) were typified with daily low temperatures averaging −12 • C. Native populations of pronghorn (Antilocapra americana) and mule deer (Odocoileus hemionus) inhabited the study area, although mule deer persisted at low densities. Wild burros (Equus asinus) and domestic cattle also occupied portions of the San Rafael. Mammalian predators included mountain lions (Puma concolor), bobcat (Lynx rufus), and coyotes (Canis latrans). The San Rafael study area was predominantly managed by the Bureau of Land Management (BLM). Recreation in the San Rafael included on and off-road travel, rock climbing, biking, hiking, camping, and hunting. Peak recreation occurred in spring during mild weather (Wade Paskett, pers. comm., UDWR).
We also studied a population of Rocky Mountain bighorn sheep in the Green River corridor (hereafter, Green River) located in Emery and Carbon counties, Utah (Figure 1). Bighorns were native to that area, but were extirpated by 1960 (Dalton and Spillett, 1971). Currently, bighorns in that area are a result of nine transplants beginning in 1970 with animals from Utah, Montana, Wyoming, and Colorado, United States and Alberta and British Columbia, Canada. Rocky Mountain bighorn in Green River form the largest metapopulation of this subspecies in Utah (Figure 1), and that area contains one of the largest contiguous sections of Rocky Mountain bighorn sheep habitat in Utah. An estimated 1,000 individuals occupied that area at the beginning of our study.
Elevation used by bighorn sheep in the Green River varied from 1,250 to 3,000 m. Riparian areas in that area were dominated by willow (Salix spp.), cottonwood (Populus fremontii), and box elder (Acer negundo) trees. General vegetation of the Green River was a desert environment dominated by salt shrubs, bunch grasses, pinyon, and juniper. Populations of pronghorn, mule deer, bison, elk, and domestic cattle occupied that area; and mammalian predators were the same as those that occurred in San Rafael, except the Green River has a large black bear (Ursus americanus) population. The Green River is of mixed land ownership, most of which is managed by the BLM. Recreation in the Green River included on and off-road travel, rock climbing, biking, hiking, camping, fishing, and hunting. The Green River corridor is also governed by the 1979 River Management Plan (Barry, 1979), which allows for private and commercial float permits and establishes a carrying capacity of 35,000 user days per year. The recreation season is divided into high and low use periods. High use occurred from May 15 through August 15. Total user days on the river have increased by 43% from 2003 to 2018 (Jaydon Mead, pers. comm., US BLM). Despite our two study populations being separated by <20 km, no genetic mixing has been documented between these populations.

Methods
In January 2012, 30 adult, female bighorn sheep were captured by aerial net gunning (Krausman et al., 1985) in the San Rafael, and 17 adult females were captured in the Green River. Efforts were made to sample a wide distribution of bighorns across the units by distributing collars proportionally to aerial count numbers and locations of bighorns established 2 months prior. Animals were equipped with Lotek 6000SD GPS/VHF collars (Lotek Wireless Inc., Ontario, Canada) with mortality transmitters and preprogrammable drop off mechanisms. Collars deployed in the San Rafael were programmed to acquire a GPS fix every 8 h, and collars deployed in the Green River were programmed to acquire a GPS fix on either a 2.5, 4, or 6 h schedule. In January 2013, eight additional females in the San Rafael and five in the Green River were captured and collared to replace mortalities, bringing the total number of collared females to 38 in the San Rafael and 22 in the Green River. Due to collar malfunctions and mortalities, we used data from 36 GPS collars (San Rafael = 19 and Green River = 17) for our GIS analyses; however, all females in the San Rafael with active collars at the time of lambing were observed for collection of data for parturition timing. We retrieved collars after mortalities or at the end of the study in January 2014, and GPS locations were uploaded from the onboard data storage of the collars. Wildlife biologists from the UDWR used care in capturing, handling, and attaching GPS collars to bighorns (Sikes et al., 2016). Using GPS data is an extremely effective way to document parturition sites in bighorn sheep, and these data should be used over other sources (i.e., VHF collars) when delineating critical birthing habitat for this species (Smith et al., 2015), especially when considering anthropogenic influences on habitat use (Longshore et al., 2013).
To determine timing of births, we monitored collared bighorn sheep in the San Rafael weekly using radio telemetry for 2 years until January 2014. We relocated collared and uncollared females with binoculars and spotting scopes to record birthdates from 25 April to 25 June during 2012 and 2013. We searched the San Rafael a mean (±SD) of every 2 days ± 1.6 days in 2012, and a mean of every 2 days ± 1.7 days in 2013. To estimate parturition dates, we observed the behavior of marked females before, during, and after parturition; as well as first sighting, motor skills, size, and behavior of neonates (Festa-Bianchet, 1988a;Whiting et al., 2008Whiting et al., , 2011. To determine birthdates for neonates of uncollared females, we compared their young with neonates of estimated ages of collared females when all females congregated in nursery groups after parturition (Côté and Festa-Bianchet, 2001;Whiting et al., 2008Whiting et al., , 2012. When female and lamb pairings were questionable, we waited until the lamb nursed to identify its mother (Festa-Bianchet, 1988a). We exercised care not to disturb females with young (Sikes et al., 2016).
We estimated birthdates of young, pooled them into sampling intervals and calculated corrected means (timing of births) and SD values (synchrony of births) for the San Rafael population in each year (Johnson et al., 2004;Whiting et al., 2011). This technique allowed robust calculations of unequal sampling intervals (bin sizes) in determining timing and synchrony of births (Johnson et al., 2004). We then calculated a date range using ±2 SDs from the mean for lambing period habitat for our spatial analyses using data from GPS collars in the San Rafael population. A mean of 88 (±3.6) adult females occupied the San Rafael during our study.
We did not relocate collared and uncollared females with binoculars and spotting scopes in the Green River during our study. The general estimate of the peak birthing period, however, for that population was from 25 May to 5 June (Whiting et al., 2011). We therefore used the date range of May 1-June 10 for our analyses of lambing period habitat in that population. A mean (±SD) of 462 (±69.3) adult females occupied the Green River during our study.

Resource Selection Functions
Using the appropriate spatial scale when defining habitat availability to animals is critical when making inferences about habitat selection at the population level, making it important to define an area biologically relevant to the species of interest (Johnson, 1980;Boyce, 2006). We used the reproducible home range (rhr) package in Program R 3.1 to delineate a 95% minimum convex polygon for each study area (Figure 1). We created those polygons using bighorn birthing and lambing period locations during the date range calculated for both study populations.
We evaluated bighorn lambing period habitat selection in those 95% minimum convex polygons using a resource selection function within a use-availability study design (Manly et al., 2002) where the response variable was 1 for a use location or 0 for a random site. Random points were extracted with the random points tool in ArcMap 10.3 (Esri, Redlands, CA, United States) and R (R Core Team, 2015). We used a mixed-effects, logistic regression with a random intercept for individual bighorn (radio collar ID) and analyzed covariates at use versus random locations. We evaluated lambing period habitat using locations only from the San Rafael study area. We then produced a model from those variables and applied it to the Green River to verify if our GIS model captured lambing period habitat for that population.
We performed a database query in ArcGIS 10.3 (Redlands, CA, United States) to eliminate locations with <3D accuracy for both study areas. Because GPS collars across both study areas were programmed at different sampling intervals, we used only one randomly selected location per collar per day both for resource selection function analysis (San Rafael) and model testing (Green River). After removing <3D locations and locations that occurred outside the 95% minimum convex polygon, use locations totaled 750 for the San Rafael, and testing locations totaled 1,130 for the Green River. We generated 750 random locations within the San Rafael to be equal to the number of use locations and assigned them equal weight. Because random locations were cast within the boundary of the study area and not associated with individual home ranges, our modeling of resource selection generally corresponded to Johnson's second order of selection (Johnson, 1980). To ensure that 750 random locations adequately characterized our study area, we calculated the true mean values (i.e., mean of all pixels within the study area) for continuous variables and compared our sample means with 95% CIs to those values (Long et al., 2014). In every case, the confidence intervals of our samples overlapped the true mean values indicating that the 750 random locations were adequate to characterize our study area. We selected the following landscape level features potentially influencing bighorn sheep habitat selection: slope, ruggedness, elevation, aspect, Landfire existing vegetation type (LANDFIRE 1 ) (Smith et al., 1991;Bleich et al., 1997;Sappington et al., 2007;Rollins, 2009), and distance to trails and roads (Smith et al., 1991;Longshore et al., 2013;Wiedmann and Bleich, 2014). We extracted topographic features using a United States Geological Survey 30 m digital elevation model. We calculated slope using the slope tool in the Spatial Analyst Tools extension. We calculated ruggedness using the Vector Ruggedness Measure Tool in the Terrain Tools extension in ArcGIS (Sappington et al., 2007;Lowrey and Longshore, 2017). That tool measures terrain ruggedness as the variation in three-dimensional orientation of grid cells within a neighborhood. Vector Ruggedness Measure values can range from 0 (no terrain variation) to 1 (complete terrain variation). We calculated aspect using the aspect tool in the Spatial Analyst Tools extension and that variable was divided into the four cardinal directions (north, east, south, west). The LANDFIRE existing vegetation type layer consisted of five types (barren, sparse, herb, shrub, and tree). We calculated distance to roads and trails using the Generate Near Table Tool in the Analysis Tools extension. Distance to water sources was not included in our analysis. Because our study occurred in spring and due to the extensive nature of ephemeral water sources formed in holes in sandstone rock after rain or from seasonal seeps, it was not feasible to map water source locations accurately or comprehensively. Additionally, major rivers were in each study area. We evaluated explanatory variables for multicollinearity and did not include any variables with a correlation coefficient |r| > 0.6.
We developed models using an information theoretic approach (Burnham and Anderson, 2002) within a mixedeffects logistic regression (Hosmer et al., 2013) and used a random intercept to account for individual heterogeneity. We used R package lme4 for mixed-effect modeling analysis (Bates et al., 2013). We used combinations of covariates (all permutations) to build models. All continuous covariates (slope, ruggedness, elevation, distance to trails, and distance to roads) were standardized before model development (x i −X/s (Baxter et al., 2017). We used AIC c values to identify the most supported models, and to identify which variables were informative (Burnham and Anderson, 2002;Arnold, 2010). We evaluated models based on AIC c , AIC c values, and AIC c model weight.
To assess predictive ability of our top model, we performed k-folds cross validation with k = 5 (Long et al., 2009;Villepique et al., 2015;Kaze et al., 2016) in the San Rafael. We randomly sorted observations into five partitions, with an equal number of locations in each partition. During each iteration of this procedure, we used four partitions (80% of the data) as the training set to estimate model coefficients and the remaining partition (20% of the data) to test model predictions. We repeated this procedure until all observations were used for the test set and part of the training set. We then applied the coefficients from the predictive model to calculate relative probability of use in the San Rafael and Green River study areas. We generated two predictive maps, one for the San Rafael and one for the Green River, by applying this procedure to each raster pixel in each study area. We then used five equal-area bins to categorize the relative probabilities of use for each pixel from low to high (Sawyer et al., 2007;Kaze et al., 2016).

RESULTS
In 2012, we estimated birthdates for 29 young. Mean (±2 SDs) birthdate for bighorn sheep in the San Rafael during that year was 21 May (±19 days). In 2013, we estimated birthdates for 16 young, and mean birthdate for bighorn sheep in that area during that year was 20 May (±21 days). In both years, 89% of 45 births occurred in May (range = April 27-June 9). The date range we used for our analyses of lambing period habitat from GPS collar data for 2012 was 2 May-9 June. Whereas, for 2013 that range of dates was 29 April-10 June.
Our habitat analyses resulted in two top models that accounted for 99% of the AICc weight (Table 1). Because the top model accounted for 86% of AIC c weight and included all covariates from the second model, we did not need to model average and only report parameter estimates from the top model ( Table 2). That model included measures of elevation, slope, ruggedness, aspect, existing vegetation type, distance to trails, and distance to roads. Estimates for variables with evidence of selection (p < 0.05) were positive for barren vegetation, slope, ruggedness, north facing slopes, and distance to roads (indicating a preference for areas farther from roads). Estimates for variables with significant negative correlation were elevation (indicating a preference for lower elevations), south facing slopes, herb, shrub, and tree vegetation types, and distance to trails (indicating a preference for areas closer to trails). Five-fold cross validation (rho = 0.99, P = 0.02) indicated that these variables predicted lambing period habitat of bighorns in our study area. In the San Rafael, predicted high-probability habitat was clumped around canyon rims surrounding the mesa top. Large areas of low-use predicted habitat occurred in that area, and locations of females during birthing were concentrated in small canyons (Figure 2). Of our 750 GPS locations, 68% were in the top two probability categories (high = 37%, medium high = 31%), 17% were within the medium category, 11% in the medium-low category, and 4% were within the low category (Figure 2). In the Green River, predicted high-probability habitat occurred along the Green River; however, some locations of females during birthing were concentrated in small canyon draws near areas of low habitat suitability (Figure 3). In the Green River, of our 1,130 GPS locations, 73% were in the top two categories of probability (high = 50%, medium high = 23%), 14% fell into the medium category, 10% fell within the medium-low category, and 3% fell within the low category (Figure 3).

DISCUSSION
Management and conservation decisions for ungulates and their habitat are increasingly being made that emphasize birthing areas (Singh et al., 2010;Kaze et al., 2016;Severud et al., 2019). For bighorn sheep, additional research is needed to understand and predict habitat use during birthing and nursery periods (Bangs et al., 2005). Our description of timing of births and GIS model predicted lambing period habitat for adjacent populations of desert and Rocky Mountain bighorn sheep. In similar, adjacent habitat these methods could be applied to conserve and manage habitat for these ungulates and to identify lambing period habitat for reintroduction and translocation sites. Reintroduction and translocation continue to be useful management tools for bighorn conservation (Wiedmann and Sargeant, 2014;Boyce and Krausman, 2018;Robinson et al., 2019). In Utah >1,000 Rocky Mountain bighorn sheep and >850 desert bighorn sheep have been released in areas of historical habitat since 1966 Shannon et al., 2008; Utah Division of Wildlife Resources [UDWR], 2018). Currently in Utah, 26 units/subunits are identified as potential sites for augmentation or reintroduction, and 27% of those units/subunits are in eastern Utah near our study area (Utah Division of Wildlife Resources [UDWR], 2018). Our results can provide a model that can identify potential lambing period habitat for both subspecies in these areas. Similar site-specific data could be used in other areas to identify birthing period habitat for other populations of adjacent ungulates.
Our GIS model showed where bighorn sheep gave birth and identified lambing period habitat; that model also indicated areas of low-probability lambing period habitat. These results can be important for land-use planning in our study areas. In these areas, outdoor recreation and energy development have increased in the last 40 years and are predicted to continually increase (Smith and Burr, 2011;Sproat et al., 2019). Several of the forms of recreation We report the number of parameters (K), Akaike's Information Criterion adjusted for small sample sizes (AICc), difference in AICc value from top model ( AICc), AICc model weight (ωi), and log likelihood (LL).
that can potentially effect bighorn sheep are off-highway vehicles, mountain biking, hiking, rafting, rock climbing, canyoneering, and camping near rivers (Papouchis et al., 2001;Longshore et al., 2013;Sproat et al., 2019). Hiking can have a pronounced effect because of the unpredictable locations of that activity and people approaching bighorn sheep, especially in spring when females are giving birth (Macarthur et al., 1979;Papouchis et al., 2001;Wiedmann and Bleich, 2014). Researchers recommend that known habitat for nursery groups of bighorn sheep be closed to hiking to benefit females and lambs (Papouchis et al., 2001).
Our results can be used by biologists and land managers in land-use planning to identify areas of low-probability lambing period habitat that can be designated for increased recreation and energy development. There is an increase in demand for special recreation or special-use permits granted by land management agencies for commercial, competitive, or organized purposes, often for large groups. Given the concentrated level of anthropogenic activity generally associated with these permits, they can be especially disruptive to bighorn sheep. Our results can be helpful for land managers when reviewing permit applications for uses that overlap timing of births and lambing period habitat. Expectedly, slope and ruggedness were significant variables in our models. Many other researchers indicate that these habitat components are important to minimize risk of predation for bighorn sheep with young (Smith et al., 1991;Bleich et al., 1997;Zeigenfuss et al., 2000;Bangs et al., 2005). Smith et al. (2015) also documented these components selected for at parturition sites and south and west facing slopes. Our model indicated, however, a preference for north facing slopes and an avoidance of south facing slopes. That response was due to a difference in latitude and temperature relative to different study areas, because bighorns are ostensibly selecting north-facing slopes in our study area to avoid warm, south-facing slopes during lactation. Smith et al. (2015) were also modeling parturition sites specifically, while we analyzed lambing period habitat. Lactating females have high water demands (Geist, 1971;Robbins, 2001), and travel of young lambs is restricted, necessitating the use of areas near water (Zeigenfuss et al., 2000). This could explain the preference for north facing slopes and lower elevation where water collects and is retained in our study area. While our GIS model indicated that bighorn sheep avoided roads, largely corroborating past work Smith et al., 1991;Papouchis et al., 2001), that model also indicated a preference for areas near trails. That result was unexpected and likely occurred because most trails in the San Rafael traverse up the bottoms and around the side rims of canyons. The locations of those trails allow recreationists easy access to unobstructed views of the scenery and landscape. Many of those areas were also used by nursery groups of bighorn sheep. Trail traffic, however, was relatively low in the San Rafael. For example, a motion-sensor camera placed on a canyon of rim trail within high-use lambing period habitat recorded 316 spring visitor days by hikers and bikers. Previous researchers have reported an avoidance of areas by bighorn sheep with more intensive human recreation (X = 15,925 visitor days, SD = 6,038 visitor days) (Wiedmann and Bleich, 2014). Therefore, in the San Rafael the number of hikers might not be affecting bighorn sheep like other areas with high-use trails. We hypothesize that that effect will increase with increased recreation in these areas. Our methods can be applied to other areas of desert and Rocky Mountain bighorn sheep habitat in Utah. For example, eastern Utah is the largest, contiguous area in the state consisting of important habitat for both subspecies of bighorn sheep (Utah Division of Wildlife Resources [UDWR], 2018). This area has been historically important bighorn sheep habitat (Flinders et al., 2002;Shannon et al., 2008), and some areas have experienced five times the growth in outdoor recreation since 1979 (Sproat et al., 2019), which often peaks in May (Papouchis et al., 2001). Oil and gas exploration and mining are also increasing industries in this area. Human disturbance has caused the abandonment of habitat by desert bighorn sheep (Etchberger et al., 1989;Lowrey and Longshore, 2017) and Rocky Mountain bighorn sheep (Wiedmann and Bleich, 2014), decreased population performance (Wiedmann and Bleich, 2014), and interrupted metapopulation dynamics for these animals (Epps et al., 2005;Bleich et al., 2016). The UDWR now has over 300 GPS collared bighorn sheep in 15 herds across the state. Our methods can be used in these areas to document lambing period habitat. Another important issue in southeastern Utah for which our GIS model can be used is in Grand Staircase-Escalante and Bears Ears national monuments. These areas are now being substantially reduced and opened for potential mining, livestock grazing, public access, and recreation (McBrayer and Roberts-Cady, 2018;Wilson et al., 2018;Reese, 2019). Timing of births data and a GIS model similar to ours could help inform natural resource managers about bighorn lambing period habitat in those areas facing immediate conservation and management changes.
Management of bighorn sheep habitat remains an important issue across much of western North America (Gutierrez-Espeleta et al., 2001;Lowrey and Longshore, 2017;Bleich, 2018;Boyce and Krausman, 2018). Proximity of bighorn sheep to domestic livestock allotments on public land (Cahn et al., 2011;Carpenter et al., 2014;O'Brien et al., 2014), and proximity of those wild ungulates to areas occupied by domestic sheep (O. aries) and goats (Capra hircus) on private land (Turner et al., 2004;Shannon et al., 2014), are controversial issues in the western United States. Domestic sheep and goats can transmit diseases to bighorn sheep (McClintock and White, 2007;Wehausen et al., 2011;Besser et al., 2012), and these issues are critical for the conservation and management of bighorn sheep populations (Cahn et al., 2011;Carpenter et al., 2014;O'Brien et al., 2014). Additionally, domestic cattle can affect habitat use by bighorns (Garrison et al., 2016). With female bighorn sheep showing general fidelity to birth sites and habitat for nursery groups (Etchberger and Krausman, 1999;Whiting et al., 2012), our data can be used to guide land-management decisions when assessing domestic livestock grazing allotments on public land, and when working with land owners on private land, which will help minimize conflict between these domestic and native ungulates.
Timing and synchrony of births in ungulates can have evolutionary implications (Marshall and Cambridge, 1937;Rutberg, 1987;Ims, 1990), such as reduced survival of young and reduced future reproductive potential (Clutton-Brock et al., 1987;Festa-Bianchet et al., 2000;Keech et al., 2000). Bighorn sheep from differing source populations that are released into adjacent areas can take up to 5 years to adapt timing of births Frontiers in Ecology and Evolution | www.frontiersin.org 9 April 2020 | Volume 8 | Article 97 to local environments, which can hinder population growth (Whiting et al., 2011). Further, females in those populations do not compensate for late births by increasing maternal care, which possibly reduces survival of young (Whiting et al., 2010). After released bighorn sheep have adapted to local environments, differences of up to 29 days can still exist in timing of births across adjacent populations, because of discrepancies in peak green-up of vegetation (Whiting et al., 2012). All that information underscores the importance of collecting site-specific birthing data to produce GIS models of lambing period habitat for bighorn sheep. Documenting the effects of human disturbance on ungulates is important for the ecology, evolution, and conservation of these animals, especially in an increasingly human-dominated world Ciuti et al., 2012a,b). As anthropogenic activity increases in these areas, conservation planning and habitat management will increasingly rely on identifying habitats used by wildlife in areas of varying levels of human use (Goodson et al., 1999;Margules and Pressey, 2000;Harju et al., 2011). Biologists can substantially reduce effects from human activities to female ungulates by identifying where and when animals give birth and rear young (Singh et al., 2010;Kaze et al., 2016;McLaren et al., 2017), and this is a current, pressing conservation issue . Identification of lambing period habitat is becoming increasingly more important for bighorn conservation (Smith et al., 2015). Additionally, with GPS data being more readily accessible in wildlife studies, our methods can be used with site-specific timing of births and abiotic data and applied to other areas occupied by Rocky Mountain and desert bighorn sheep -and other ungulates -to identify important lambing period areas, which will aid in conservation and landuse planning.

DATA AVAILABILITY STATEMENT
The datasets analyzed for this study can be found by contacting the Utah Division of Wildlife Resources.

ETHICS STATEMENT
Ethical review and approval was not required for the animal study because all animal handling was done by the Utah Division of Wildlife Resources who held the pertinent permits to do such.

AUTHOR CONTRIBUTIONS
RR, JW, and JS designed the study. RR performed field work. JW, RR, and RL analyzed the data. RR, TS, and JW wrote an initial draft of the manuscript and all authors finalized it.