Abstract

We investigated influences of risk of predation by mountain lions (Puma concolor), topographic metrics at multiple scales, and vegetation, land, and snow cover on resource selection by Sierra Nevada bighorn sheep (Ovis canadensis sierrae), an endangered taxon, during winters 2002–2007, in the Sierra Nevada, California, USA. We hypothesized that those mountain ungulates would trade off rewards accrued from using critical low-elevation habitat in winter for the safety of areas with reduced risk of predation. Sierra Nevada bighorn sheep did not trade off benefits of forage for reduced risk of predation but selected areas of high solar radiation, a correlate of vegetation productivity, where risk of predation by mountain lions was greatest, while mitigating indirect risk of predation by selecting for steep, rugged terrain. Bighorn sheep selected more strongly for areas where mountain lions were active, than for low-elevation habitat in winter, likely because mountain lions were most active in those areas of bighorn sheep winter ranges overlapping ranges of mule deer (Odocoileus hemionus), where both ungulates accrued forage benefits. We demonstrated reduced benefit of migration to low elevation during drought years, providing an alternative explanation to the predator-induced abandonment hypothesis for the disuse of low-elevation winter range observed during drought years.

1. Introduction

Animals living in temperate or arctic environments, where a seasonal abundance of forage coincides with increased nutrient demands of late gestation and lactation [13], must balance the need to acquire nutrients against constraints from risk of predation [46]. Many populations of ungulates migrate between discrete seasonal ranges [713], with those occupying montane environments obtaining high-quality resources by selecting among elevations that enable exploitation of new growth in forage [11, 14]. Benefits of migration to areas of high-quality forage must outweigh increased risk of predation to comprise an evolutionarily stable strategy [15]. Populations of mountain sheep (Ovis spp.) occupying montane environments may migrate between high-elevation summer ranges and lower-elevation winter ranges, corresponding to the progression of new growth in grasses, forbs, and shrubs [1618]. Variation in temperature, precipitation, and vegetation phenology, however, may alter behavior and habitat selection by mountain sheep [13, 19, 20].

Predator avoidance operates through indirect mechanisms that affect the likelihood of encountering, detecting, or eluding a predator, as well as through direct means by which animals minimize the odds of success of a predator in capturing prey upon detection [21, 22]. Direct mechanisms of predator avoidance include sight [23, 24], as well as olfactory and auditory cues to proximate risk posed by a predator [2527]. Mountain sheep also may assess the likelihood of encountering and escaping a predator indirectly, as a function of habitat, in the absence of overt signs of a predator [13, 28, 29]. Visually unobstructed habitats enable ocular detection of predators [13, 24, 30], whereas proximity to a refuge (e.g., steep and rugged escape terrain [17, 31]) may confer protection from some predators, in part, because coursing predators seldom use those areas (e.g., coyotes, Canis latrans [32]; gray wolves, C. lupus [4, 33]). Differences between the sexes in risk tolerance and predator avoidance strategies [29, 32, 34, 35] must be accounted for in models of resource selection under risk of predation.

Selection of resources is subject to constraints imposed by perceived risk; sublethal effects of risk of predation may have profound influences on behavior and habitat selection [36, 37] and incur physiological costs [38]. There was, however, no association between sublethal response to risk of predation by wolves on pregnancy or body condition in North American elk (Cervus elaphus [39]). Responses to fear of predation have been posited to drive shifts in habitat selection, ultimately resulting in negative demographic effects among large herbivores [40, 41].

Mountain lions (Puma concolor) are the primary predators of mountain sheep in areas of sympatry and have been implicated as a threat to populations of Sierra Nevada bighorn sheep (O. canadensis sierrae [41]), an endangered taxon [42]. These large felids may cause substantial mortality in some populations of mountain sheep [4346], but questions remain about the potential for sublethal effects as a consequence of changes in behavior or habitat selection by bighorn sheep. A precipitous decline in the largest population of Sierra Nevada bighorn sheep at Mount Baxter, California, USA, from 1987 to 1991, was attributed to decreased recruitment after females failed to occupy high-quality, low-elevation winter ranges in response to risk of predation by mountain lions [41]. Wehausen [41] concluded that “range abandonment” by Sierra Nevada bighorn sheep was a response to perceived risk of predation. Alternative explanations, such as the effect of 6 years of drought [47, 48] on forage or the potential demographic consequences of removal of large numbers of bighorn sheep from the winter range for translocation [49], however, have not been adequately investigated [50]. While evaluation of the predator-induced range abandonment hypothesis [41] is limited by the absence of information from telemetered animals, and by a lack of information on habitat conditions, an evaluation of indices of snow cover and forage conditions from the period of reported range abandonment is possible using archived, remotely sensed data. Those habitat factors likely played a principal role in winter habitat selection by bighorn sheep.

Snow cover may be reliably indexed through remote sensing of normalized difference snow index (NDSI) [51, 52] and variation in forage quality and quantity indexed by normalized difference vegetation index (NDVI), calculated from multispectral satellite imagery, at spatial and temporal resolutions dictated by the satellite platform employed [53]. Peak values of NDVI also were correlated with values of fecal crude protein [54], an index of dietary quality [55], thereby validating NDVI as a useful index to quality of forage used by ungulates.

For sublethal effects of predation risk to result in the severe demographic consequences posited by Wehausen [41], bighorn sheep must have made disproportionate tradeoffs by giving up nutrients in exchange for reduced risk of predation. Under that predator-induced abandonment hypothesis, a landscape-scale shift in habitat selection must have occurred, wherein Sierra Nevada bighorn sheep avoided traditional low-elevation winter ranges in response to risk of predation from mountain lions, trading off benefits of forage at low elevation for reduced risk by use of high-elevation habitat. We tested the overarching hypothesis that bighorn sheep subject to substantial risk of predation would select areas of lesser risk, but with lower-quality forage, and concomitantly reduced nutrient availability, in both wet and dry years (Table 1). Alternatively, if landscape-scale movements are not evident, we expected bighorn sheep to remain in areas with high predation risk to access quality forage, while exhibiting antipredator behaviors through fine-scale selection of rugged and steep terrain (Table 1).

2. Methods

The study area was the eastern slope of the Sierra Nevada, a rugged escarpment rising from the Owens Valley and forming the western boundary of the Great Basin in eastern California, USA. The highest peaks and steepest slopes are located in the southern portion of the range in Inyo and Mono counties, California (centroid 37°11′N, 118°23′W), where four distinct populations of Sierra Nevada bighorn sheep occupied rugged granitic slopes (Mono Basin, Wheeler Ridge, Mount Baxter, and Mount Langley [29]; Figure 1). We focused on the two largest populations of Sierra Nevada bighorn sheep: Wheeler Ridge and Mount Baxter. Elevations of winter ranges were 1,500–4,100 m at Wheeler Ridge and 1,300–4,100 m at Mount Baxter. Winter ranges of both populations of bighorn sheep overlap winter ranges of migratory herds of mule deer (Odocoileus hemionus [56]). The predominant vegetation community at low elevation was Great Basin sagebrush steppe, dominated by sagebrush (Artemisia tridentata) and antelope bitterbrush (Purshia tridentata [57]), whereas mountain whitethorn (Ceanothus cordulatus), mountain mahogany (Cercocarpus ledifolius), and horsebrush (Tetradymia canescens) were common at middle elevations. Trees, including Jeffrey pine (Pinus jeffreyi), and juniper (Juniperus occidentalis) were sparsely distributed throughout middle elevations, with grasses (Achnatherum speciosum, Festuca spp.) and buckwheat (Eriogonum spp.) common as understory species. At high elevations (3,300–4,300 m), characteristic alpine vegetation included alpine crested wheatgrass (Agropyron spp.), sedges (Carex spp.), and vetches (Astragalus spp. [29]).

2.1. Animal Capture

During 2002–2007, we captured bighorn sheep with a net gun fired from a helicopter [58]. Bighorn sheep were physically restrained and transported to a central processing area where samples were collected and a GPS (global positioning system) collar (described later) and a small secondary collar containing a VHF (very high frequency) transmitter were fitted to adult (>1 year old) bighorn sheep. Bighorn sheep were returned to the capture location for release.

We captured mountain lions from 2002 to 2009 by pursuit with hounds, employing techniques described by Pierce et al. [59, 60] or by opportunistic darting. Adult mountain lions were immobilized with Telazol (4.4 mg/kg [59]), with an extra-long-range projector (Palmer Cap-Chur, Inc., Powder Springs, Georgia, USA). When proximity to water presented a potential risk of drowning, we used darts equipped with a VHF transmitter (Pneudart Inc., Williamsport, Pennsylvania, USA) [61] to facilitate rapid location of mountain lions after immobilization. Mountain lions were fitted with GPS collars (described later) and monitored at the capture site until they regained mobility. Capture methods for bighorn sheep and mountain lions followed guidelines of the American Society of Mammalogists [62], and were approved by the Institutional Animal Care and Use Committee (protocol #0202) at the University of Alaska Fairbanks.

2.2. Spatial Analysis: Bighorn Sheep

Four models of GPS collars were deployed on bighorn sheep: Advanced Telemetry Systems GPS 2100 (; Isanti, Minnesota, USA), Lotek Wireless 4400s (; Newmarket, Ontario, Canada), and Televilt Simplex () and Televilt Telus (; TVP Positioning AB, Lindesberg, Sweden). Collars used heterogeneous fix schedules of 3–720 GPS locations/day, because some collars were programmed by collaborators addressing different objectives. Data from GPS collars were censored for 2 days after capture to decrease the influence of atypical movements immediately following capture [63], and erroneous GPS locations eliminated by screening in ArcView 3.2 (Environmental Systems Research Institute, Redlands, California, USA; ESRI) with the BadFix extension [64]. We analyzed data for animals with ≥10 days of location data in winter (January–April).

We used Home Range Tools for ArcGIS [65] and ArcGIS 9.2 software (ESRI, Redlands, California, USA) to define the study area as the winter ranges available to the Wheeler Ridge and Mt. Baxter populations from composite home ranges for each population of bighorn sheep in winter. For purposes of defining the study area, data from GPS collars were rarified to one randomly selected location per animal, per week, to minimize autocorrelation and reduce the influence of heterogeneous sample size among individual animals. The 95% adaptive-kernel home range was calculated with a smoothing factor of 80% that prevented polygons within home ranges from fragmenting [66] (Figure 1).

Analysis of resource selection used only high-quality locations, after elimination of 2D locations and those locations where dilution of precision (DOP) was ≥10 [67, 68]. Such screening introduced biases against detecting use of areas with high canopy closure or areas where available sky is obscured by topography [6870]. Closed canopy was rare in our study area and, where present, generally was avoided by mountain sheep [17]. Cain III et al. [69] demonstrated a reduction in 3D fixes of 36% in areas with <33% available sky versus those locations with >66% available sky. We acknowledge a likely bias against detection of use in areas where proximity to steep escape terrain limited available sky necessary for satellite fixes. High-quality GPS data were rarified by selecting one location daily, nearest in time to solar noon, a timing chosen for consistency with historical data based on daytime visual observations [41]. We tested for independence of locations using the program ASSOC1 [71], assuming that animals with a proximity of <200 m for >25% of locations represented animals in the same group.

We calculated mean GPS location error, based on published error estimates [64, 67], and weighted by the proportion of GPS locations from each model of collar in our analysis, to assess appropriate scale for evaluation of third-order habitat selection [72, 73]. Location error of 3D fixes, from previous studies, averaged 6.3 m (i.e., 50% circular error probable; CEP) and 95% CEP = 28.6 m. Ensuring that the accuracy of locations of animals is comparable to the scale of habitat maps is a key consideration when evaluating resource selection [73]. Our estimates of GPS location error were comparable to the accuracy of Landsat data used (approximately 30 m, or 1 pixel; Section 2.6).

2.3. Spatial Analysis: Random Locations

Resources available to bighorn sheep were quantified by casting ten random locations within a 2.4 km aerial radius of each midday bighorn sheep location (Figure 2) with Hawth’s Analysis Tools 2.10 for ArcGIS 9.2 [74]. The 2.4 km radius reflected availability within a distance an animal could readily move in 24 h at the rate of 100 m/h noted for 16.3% of movements during 1–36 h intervals between daily GPS locations, allowing evaluation of fine-scale selection, while also encompassing high and low elevations within the circle of availability [74, 75] (Figure 8). Random locations were further constrained by the 95% composite home range for bighorn sheep; random points falling outside that composite home range for the population were eliminated and new random points cast to achieve 10 random locations per animal location. Animal locations represented used points in the matched-case design [76]. Fewer than 0.01% of random locations were cast within the same pixel as used locations; thus contamination was not a concern [77].

2.4. Spatial Analysis: Mountain Lions

We deployed four models of GPS collars on mountain lions: Lotek Wireless 4400s (), Televilt Simplex (), Televilt Telus (), and Northstar RASSL Tracker (; Northstar Science and Technology, King George, Virginia, USA). Collars were programed to collect 6–10 GPS locations per day; 1–3 locations during midday (10:00–14:00), and 5–7 locations during crepuscular and nighttime hours (18:00–06:00). Risk of predation by mountain lions was indexed by two measures: the distribution of locations of GPS collars from mountain lions active (and likely hunting) within winter range of bighorn sheep and the distribution of locations where bighorn sheep were killed or cached by mountain lions.

We defined locations of “active” mountain lions by selecting locations of GPS collars not associated with cached prey or resting places [78]. Clusters of locations characteristic of cached prey or of resting sites, where mountain lions were not actively hunting, were selected with a C++ program (see Supplement 1; ClusterDetector in Supplementary Material available online at http://dx.doi.org/10.1155/2015/357080), which identified clusters of locations based on temporal and spatial parameters. We set the spatial tolerance to <200 m within ≤6 days [79]. We retained only one GPS location for each cluster, to avoid pseudoreplication of locations where inactive mountain lions were likely resting [78], or feeding at a kill [59] and thus not actively hunting. We retained all GPS locations (i.e., 2D and 3D locations of all DOP values) to increase the sample of locations of active mountain lions for calculation of utilization distributions.

Locations where mountain lions killed and cached bighorn sheep () were identified by evaluating clusters of locations from GPS collars deployed on mountain lions () [79], investigating mortality signals from collared bighorn sheep (), or from kills encountered while following mountain lion tracks or trailing hounds that located a carcass (). Mortalities of bighorn sheep were assigned as “probable” kills by mountain lions when two or more of these characteristics were present: tracks, drag marks, or feces from mountain lions; puncture wounds or hematomas on the neck; presence of an intact rumen; ribs cleanly cut as by carnassial teeth; and caching of the carcass. Where possible, location of a kill site was marked with a hand-held GPS (Model 12 XL, Garmin International Inc., Olathe, Kansas, USA); however, in most instances, location of the cached carcass rather than the kill site was recorded, because tracks, drag marks, or other kinds of evidence (e.g., blood, hair) at the kill location often were not detectable on granitic substrates.

We created one risk surface from all locations of GPS collars of active mountain lions within winter range of bighorn sheep and a second risk surface from all locations where bighorn sheep were killed or cached by mountain lions; these were based on adaptive-kernel utilization distributions (Home Range Tools for ArcGIS; [80]) with a smoothing factor of 80% that prevented polygons within home ranges from fragmenting [66]. Risk values from utilization distributions were scaled from 1 to 99 and areas within the winter range that fell outside of the 1% utilization distribution for each risk surface were assigned the minimum value of 1.

2.5. Spatial Analysis: Topography

We obtained elevation data from 1 : 24,000-scale, 10 m digital elevation models (DEM; US Geological Society Seamless Data Server: http://seamless.usgs.gov, accessed 27 September 2010). Average incident solar radiation was quantified monthly in winter (January–April) by calculating potential relative radiation, a measure incorporating monthly changes in solar elevation, as well as local shading because of topography [81]. We calculated potential relative radiation as the sum of hillshade values for each daylight hour on the 15th day of each month, with solar azimuth and elevation values for that hour (NOAA Solar Calculator: http://www.esrl.noaa.gov/gmd/grad/solcalc/). Terrain ruggedness was calculated with the terrain diversity index of Nicholson et al. [12]. We also calculated a binary value for convexity by comparing the elevation at a location with the average elevation of 10-m pixels within the evaluation radius, differentiating concave drainages (convexity = 0), which are less likely to be used by bighorn sheep, from convex ridgelines (convexity = 1); those features could have an identical terrain-ruggedness value. Positive convexity over 15–30 m likely represents outcroppings or ridges providing enhanced visibility for bighorn sheep [24] or, more generally, escape terrain at the microhabitat scale of the location of the animal, whereas negative convexity at the 100–150 m macroscale probably indicates proximity to steep slopes and rock outcroppings constituting escape terrain. Terrain ruggedness and convexity were calculated over radii of 15, 20, 30, 100, and 150 m to elucidate the scale at which bighorn sheep selected for rugged and convex escape terrain [8284].

2.6. Remote Sensing

We obtained daily data for precipitation from 1949 to 2009 from the Bishop Airport, California, USA, (Western Regional Climate Center: http://wrcc.dri.edu/) located 25 km from winter range at Wheeler Ridge and 42 km from winter range at Mount Baxter. Data were summarized by calculating cumulative precipitation for each date, beginning 1 October, for the water year (1 October–30 September), total precipitation for each water year, and departure from the 60-year average (1949–2009) for each water year.

We used image analysis software (Imagine 2010, ERDAS, Inc., Atlanta, Georgia, USA) to classify vegetation types and calculate vegetation and snow-cover indices from cloud-free Landsat 5 Thematic Mapper (TM; ; 1987–1994 and 2002–2009) and Landsat 7 Enhanced Thematic Mapper Plus (ETM+ 2002-2003; collectively denoted herein as “TM” imagery; USGS; http://earthexplorer.usgs.gov; accessed 2–18 August 2011). Downloaded images were used without postprocessing, with geometric and radiometric correction by the USGS Earth Resources Observation and Science (EROS) data center. Imagery employed Level 1 Terrain Correction, which used both ground control points for geometric correction and 1-arc second DEMs for orthorectification to remove errors introduced by steep slopes (http://landsat.usgs.gov/descriptions_for_the_levels_of_processing.php). We used 1-m resolution, 4-band aerial imagery from National Agriculture Imagery Program (NAIP) collected in 2009 (http://www.atlas.ca.gov/download.html#/casil/imageryBaseMapsLandCover/imagery/naip/naip_2009) to generate the rock cover class at the 1-m scale of those aerial photos.

We calculated NDVI and NDSI [51] from digital number values of 31 TM images collected 2002–2009. To minimize inclusion of misleading values of NDVI that measured snow and not vegetation, we classified pixels having NDSI > 0.2 as snow-covered [52] and nullified the noninformative NDVI values for those pixels. NDVI values were retained for snow-free pixels with NDSI values ≤ 2. We used the more sensitive NDSI threshold of >0.2 to maximize sensitivity to snow, described by Dozier [51] as sensitive to detection of alpine snow cover even in shadowed areas in the southern Sierra Nevada. Pixels covered by snow (NDSI > 0.2), tree-shrub, or rock were assumed to have no forage value and assigned the reference NDVI value of zero. The remaining pixels ( = 10,378) were categorized into three equal bins as “Low” (1), “Medium” (2), and “High” (3).

Sesnie et al. [85] assessed the utility of remote sensing to evaluate bighorn sheep habitat, concluding that MODIS indices were less degraded by effects of low sun angle and topography than were those from Landsat TM. Heterogeneity of cover at the 250-m scale of the MODIS sensor, however, precluded the use of MODIS to assess third-order selection. Further justification for using Landsat is the availability of TM archives dating from 1982, prior to and spanning the period of “range abandonment” by Sierra Nevada bighorn sheep reported by Wehausen [41]. We obtained vegetation classifications from the US Forest Service Existing Vegetation (EVEG) data layers (http://www.fs.usda.gov/detail/r5/landmanagement/resourcemanagement/?cid=stelprdb5347192; image date; 31 August 2001, accessed 8 November 2010). We combined EVEG layers for conifer forest woodland, hardwood forest woodland, and mixed conifer-hardwood woodland into a single tree layer. The EVEG system did not adequately resolve areas with partial canopy closure, often dominated by mountain mahogany or juniper, which EVEG usually classified as “shrub,” a classification that also included much of the open shrub habitat where bighorn sheep foraged. We classified a cover type of open-canopy tree-shrub using Imagine 2010 and a TM image from 1 May 2006, a timing that captured leaf-out of most mountain mahogany, antelope bitterbrush, and mountain whitethorn, characteristic shrubs in this cover class. The resulting classification was validated with 1 m resolution NAIP imagery collected in 2009. Specificity was 81% (79/98), and sensitivity 73% (37/51). We combined our open-canopy, tree-shrub layer and the EVEG tree classification into a composite tree-shrub cover type.

2.7. Forage Tradeoff

We evaluated the potential benefit of migration to low-elevation habitat in winter by assessing differences in the quality and quantity of potential forage, as indexed by NDVI at high elevations (>3,000 m) and low elevations (<2,000 m) used by bighorn sheep in winter (Figure 3). These elevations were chosen as convenient thresholds for delineation of areas encompassing the lowest and highest elevations used by bighorn sheep in winter, in which we also had adequate GPS collar data. Only high-quality locations (3D, DOP < 10) from GPS collars were retained to investigate forage tradeoff. We assumed those locations were representative of the earlier 1987–1995 period, when disuse of low-elevation winter range was reported [41], but for which no telemetry locations were available [50]. Available telemetry data, including an earlier study [86] and recent data from GPS collars, demonstrate use of both high and low elevations in winter, supporting the extrapolation of modern GPS collar data as representative of locations used by bighorn sheep in winter [50].

We used 57 Landsat 5 and Landsat 7 TM scenes (Path 42, Row 34), from 1987 to 1994 and from 2002 to 2009 to assess snow cover (NDSI) and NDVI for pixels used by bighorn sheep in winter. Continuous values for NDVI were calculated after removing TM pixels covered by tree-shrub, rock, or snow (NDSI > 0.2) and, consequently, removed pixels were not indexed for available forage. Within-scene comparisons between high and low elevations were essential to detect the weak signal (≤0.10 NDVI units) among atmospheric variation (e.g., path radiance, atmospheric attenuation, and aerosol-induced distortion) and variation arising from geometric sun angle-illumination (e.g., solar elevation, azimuth angle, and viewing angle), as well as sensor calibration, all of which can alter NDVI values across images [87]. Our use of both Landsat 5 and Landsat 7 sensors added an additional source of between-image measurement error that was mitigated through within-image comparison.

2.8. Statistical Analyses

Resource Selection Modeling.–SAS 9.3 software (SAS Institute, Cary, North Carolina, USA) was used for all statistical analyses. We modeled resource selection functions (RSFs) separately for females and for males [32] but pooled the two study areas to increase inference limits by avoidance of overfitting models to site-specific features. We calculated RSFs using conditional logistic regression (PROC LOGISTIC [88]) in a matched-case design [8, 76, 89, 90], with 10 random locations per used location. Matched-case design is particularly appropriate for evaluation of RSFs when resources vary over time [91]. Logistic models were fit for the populations of males and females, including a random effect conditioned upon the temporal window of each TM image, thus fitting the model to temporal variation in vegetation quality and quantity (indexed by NDVI) snow cover (indexed by NDSI), and solar radiation.

Prior to analysis, we rescaled variables for elevation, potential relative radiation, and terrain ruggedness at all radii (15–150 m), so that a 1-unit change represented 100 units of the change in the raw variable; thus, odds ratios for those variables indicate the relative odds of selection for each 100-unit change in elevation, potential relative radiation, and terrain ruggedness. Mean, SD, and range for those variables in their original scale are provided in Table 5.

2.9. Model Selection and Validation

We retained candidate variables with absolute values of Pearson correlation < 0.60 and screened for multicollinearity using linear regression diagnostics in SAS 9.3 software [92]. Candidate models included all possible combinations of noncorrelated predictor variables with univariate [93]. We calculated Akaike’s information criterion adjusted for small sample size (AICc), , and Akaike weights for candidate models where <4, indicating substantial information-theoretic support [94]. We selected among correlated predictor variables by comparing among models differing only in those correlated predictors. We evaluated the relative importance of variables based on their Akaike importance weights, calculated as the sum of Akaike weight across all models that contained a particular variable [94]. We then calculated model-averaged parameter estimates and unconditional standard errors (SE), thereby incorporating both parameter uncertainty and model-selection uncertainty for each predictor variable [95]. We determined if model-averaged parameter estimates differed from zero by examining whether the unconditional 95% CI overlapped zero [11].

We used -fold cross validation with five partitions [76, 9699] to evaluate predictive strength of the resource selection functions by withholding a randomly selected 20% test set of used and random locations and estimating model parameters with the remaining locations. For each of five iterations, the coefficients calculated from the training set were used to calculate RSF values for random locations in the test set, which were then ranked by RSF and assigned to 1 of 10 equal-area bins (i.e., 10-percentile bins). Training-set coefficients were then used to calculate values of RSF for used locations in the test set. Used locations were placed in the bins according to RSF value, and number of locations summed for each bin. We regressed the count of used locations in each bin versus the median value from random points and used coefficient of determination, slope, and Spearman rank correlations, averaged among the fivefold validation groups, as diagnostic of the predictive success of each model [97, 99].

2.10. Analysis of Forage Tradeoff

We compared NDVI at low elevations (<2,000 m) with NDVI at high elevations (>3,000 m) when ≥10 snow-free pixels used by bighorn sheep were available in each elevation category within a TM scene, to ensure an adequate sample of potential forage. Within each TM scene, mean NDVI was calculated among snow-free, low-elevation pixels used by bighorn sheep and among snow-free, high-elevation pixels used by bighorn sheep and not classified as rock or tree-shrub. We calculated NDVI tradeoff as the difference of mean low-elevation NDVI minus mean high-elevation NDVI for each TM image and assigned precipitation to date, total precipitation by water year, and total precipitation in the prior water year to each TM image. We calculated the mean NDVI across 202–1,396 snow-free, low-elevation pixels, and 10–696 snow-free, high-elevation pixels, and calculated NDVI tradeoff as the difference between mean low-elevation NDVI minus mean high-elevation NDVI for 51 TM scenes having ≥10 snow-free pixels at high and low elevation.

We compared the NDVI tradeoff between years of above-average precipitation, with the NDVI tradeoff in those years of below-average precipitation with the nonparametric Wilcoxon signed-rank test for median difference [100], to account for small samples of snow-free pixels present at high elevation during above-average snow years. We explored the relationship among the NDVI tradeoff and measures of precipitation through linear regression.

3. Results

Precipitation was highly variable (CV = 59%) over the 61 years for which data were available (1949–2009) and exceptionally variable (CV = 73%) during 2002–2007, when GPS data were collected from collared bighorn sheep (Figure 9). The driest (2002; 3.5 cm) and third-wettest (2005; 28.5 cm) years on record occurred during that period, demonstrating extremes in snow cover and resource selection by bighorn sheep during drought and wet years (Figure 3).

3.1. Location of Bighorn Sheep and Mountain Lions

We deployed GPS collars on 28 individual bighorn sheep (19 females, 9 males) and collected a mean of 79 days (SD = 49.2 days, range = 10–172 days) of locations in winter. Data met the assumption of independence among animals; no pairs of GPS collars deployed on bighorn sheep were associated (<200 m) for >25% of locations. Thirty-nine GPS collars (including replacement collars) were fitted to 24 individual mountain lions (12 females, 12 males; mean number of fixes 1,094, SD = 1,046, range = 65–4,304 fixes). Twenty-two mountain lions (12 females and 10 males) used winter ranges of bighorn sheep at Wheeler Ridge () or Mount Baxter (), and two mountain lions overlapped both winter ranges. Active locations comprised 47.9% of 5,453 GPS locations within winter ranges of bighorn sheep.

Risk of predation by mountain lions, quantified by the risk surface of locations of active lions (Figure 4), was negatively correlated with elevation (; Table 6); consequently, candidate models included either one, but not both, predictor variables. The risk surface generated from locations where mountain lions killed and cached bighorn sheep was moderately correlated with the risk surface from active lions (; Figure 4) and elevation (; Table 6); accordingly, both variables measuring risk from mountain lions were included as candidate predictors for modeling RSFs.

3.2. Resource Selection Models

Results of RSF modeling indicated selection for the risk surface quantifying probability of mountain lion activity as one of the strongest predictor variables (i.e., Akaike importance weight = 1.0; Table 2; Figure 5). Although both sexes of bighorn sheep selected for that measure of risk of predation by mountain lions, selection was nearly twice as strong for males () as for females (), with nonoverlapping confidence intervals (Table 2). Selection for areas used by active mountain lions was a better predictor of resource selection during winter than was selection for lower elevations ( for females, for males; Table 3).

Additional predictors having the strongest level of support in all candidate models (i.e., Akaike importance weight = 1.0; Table 2) for the sexes of bighorn sheep included positive selection for areas of high solar radiation, steep slopes, and proximity to rugged terrain, which differed in scale by sex (150 m for females, 100 m for males; Figure 5). Our methods prevented simultaneous evaluation to determine relative strength of the correlated metrics of ruggedness at 100-m and 150-m scales (; Table 6). All candidate models for females included an additional predictor, avoidance of snow (i.e., negative selection for NDSI), whereas candidate models for males indicated avoidance of tree-shrub cover (Akaike importance weight = 0.97; Table 2; Figure 5).

Less influential variables in the best-fit model included the rock cover class, avoided by both sexes of bighorn sheep, whereas females selected for convexity at the 15-m scale, avoided ruggedness at the 30-m scale, and avoided convexity at the 150-m scale; males selected for ruggedness at the 20-m scale but avoided convexity at the 100-m scale (Figure 5). Masking temporally varying snow cover, static tree-shrub, and rock cover classes for noninformative NDVI values collectively eliminated 49.9% of 18,240 random locations in our data set. Vegetation conditions, as indexed by NDVI at the remaining 50.1% of locations, relative to reference areas of rock, snow or tree-shrub, did not enter RSF models. Results of -fold cross validation for RSF models indicated robust fit and substantial accuracy of prediction by models for female and male bighorn sheep (Table 4).

3.3. Evaluation of Forage Tradeoff

Forage tradeoff was evaluated with high-quality locations from GPS collars that identified TM pixels used by bighorn sheep in winter at high elevation ( = 1,272 locations; unique TM pixels) and low elevation ( = 3,830 locations; = 1,396 unique pixels). Values of NDVI pooled across TM images showed high variability and were inconclusive in assessing a potential tradeoff between low-elevation vegetation and high-elevation vegetation condition in drought years compared with wet years (Figure 6(a)). When data were analyzed by comparison of the median of differences between low-elevation habitat and those at high elevations within TM images, however, a significant tradeoff in NDVI was demonstrated for years of above-average precipitation compared with years of below average precipitation (Figure 6(b); Wilcoxon statistic = 589, ). The best single predictor of the magnitude of an NDVI tradeoff was precipitation to date (Figure 7). Additional variation was explained ( adjusted = 0.51) by addition of precipitation in the prior water year, yielding the linear equation: NDVI tradeoff = 0.017 + 0.0030 × (precipitation to date) + 0.0013 (precipitation in prior water year).

4. Conclusions

4.1. Hypothesis I-A: Resource Selection Was Influenced by Direct Risk of Predation

The hypothesis predicting an influence by direct risk of predation on resource selection by male and female Sierra Nevada bighorn sheep (Table 1) was supported, but in a direction opposite to that predicted under the predation risk hypothesis. One of the strongest positive predictors of resource selection by bighorn sheep was a measure of direct risk of predation by mountain lions; all candidate models (i.e., Akaike importance weight = 1.0) for female and for male bighorn sheep included selection for proximity to locations where mountain lions were active (Figure 5; Table 2), and likely hunting [78]. This outcome may have resulted from bighorn sheep selecting areas with suitable forage and escape terrain and mountain lions hunting in areas where mule deer overlapped winter range of bighorn sheep [50, 56]. This outcome also was consistent with the determination by Pierce et al. [101] that the distribution of mountain lions was predicated upon the distribution of mule deer, their primary prey, available to hunt. Positive selection for areas of high use by mountain lions was also reported by Greene [102], in her assessment of a cohort of Sierra Nevada bighorn sheep collared in 2007-2008 on the Mount Baxter winter range.

Challenges to assessing the influence of direct risk of predation included the choice of appropriate spatial and temporal scales. We quantified direct risk of predation with a static probability density function generated from the aggregate of locations of active mountain lions over eight winters. Such a metric allowed evaluation of resource selection by bighorn sheep relative to generalized predation risk but did not attempt to quantify short-term responses to transient risk, which are evident in correlation of locations of bighorn sheep and mountain lions, when viewed in a temporal window of days (J. T. Villepique, unpublished data). Although mountain lions are protected in California [103], a total of three mountain lions were removed from the two winter ranges, which were managed under identical protocols, after those animals were determined to have killed endangered bighorn sheep. Removal of three mountain lions identified as a threat to endangered bighorn sheep may have reduced the influence of those mountain lions more likely to hunt bighorn sheep on the probability distribution generated from GPS-collared mountain lions.

4.2. Hypothesis I-B: Resource Selection Was Influenced by Indirect Risk of Predation

The hypothesis positing influence of indicators of indirect risk of predation was supported by selection of rugged topographic features conferring protection from some predators. All models of RSF for the sexes indicated significant selection for steep slopes and selection for proximity to escape terrain (Table 2; Figure 5). Differences in scale among the sexes of bighorn sheep likely were related to group size [23, 29, 104]. The indirect risk hypothesis received additional support from improved model fit when measures of ruggedness and convexity in immediate proximity to the animal were added to RSFs (Table 2, Figure 5). An unexpected result was that females selected against ruggedness at the intermediate 30 m scale, albeit with weak evidence (Akaike importance weight = 0.09). Rachlow and Bowyer [13] similarly documented selection of less-rugged microhabitats, proximate to rugged escape terrain, by female Dall’s sheep (O. dalli). Proximity to escape terrain (represented by the strongly selected ruggedness within 150 m) likely mitigated indirect predation risk for female bighorn sheep in our study (Table 2, Figure 5). A remaining corollary supporting the indirect risk hypothesis was the avoidance of tree-shrub cover by males although, surprisingly, not by females, possibly because availability of tree-shrub cover was lower at midelevations occupied by females compared with riparian and shrub areas proximate to the lowest elevations, used exclusively by males. Another unexpected outcome was avoidance of rock cover by both sexes. This outcome was best explained by inclusion of large areas of steep cliffs inaccessible to bighorn sheep in that cover class.

Evaluation of resource selection by bighorn sheep, relative to habitat-mediated, indirect risk of predation, presented challenges related to biases inherent in GPS collars, which are less likely to obtain high-quality fixes when proximity to escape terrain blocks available sky. Difficulties also arose in differentiating qualitative differences within the rock cover class, present at only 1.7% of random points (Table 5). Rock cover was associated with steep escape terrain used by bighorn sheep in our study area but also included massive granitic plutons, favored by human technical rock climbers, which were too steep for bighorn sheep to traverse. Our analysis of resource selection also did not account for effects of sociality, specifically group size, an important determinant of foraging decisions under risk of predation in bighorn sheep [23, 29, 104].

4.3. Hypothesis II: Habitat Selection Is Influenced by Forage Availability

The forage-availability hypothesis was supported by significant selection in all candidate models for potential relative radiation, a proxy for insolation that is correlated with patterns of vegetation growth [81]. Similarly, females, but not males, avoided snow cover (i.e., selected for negative NDSI), as predicted by the forage-availability hypothesis. Models of resource selection did not, however, include NDVI as an informative predictor. Despite shortcomings in the application of TM imagery to evaluate forage in a use-versus-availability model at the scale of third-order selection, we successfully used radiance values from TM sensors to elucidate strong avoidance of snow cover by females (i.e., females selected for negative NDSI). Males, however, sometimes remained in snow-covered, high-elevation areas rather than moving to lower elevations (J. T. Villepique, personal observations), as was typical for females. This divergent strategy relative to snow was likely related to the larger body size of males [17], allowing movement and foraging in snow at a lower cost, relative to that for smaller-bodied females. Male bighorn sheep in the Sierra Nevada selected areas with significantly higher biomass of shrubs and consumed significantly greater proportions of sagebrush than did females [29]. Those results are consistent with size differences between male and female digestive systems [2, 3], allowing males to use lower-quality woody forage emerging from snow, whereas females would have more difficulty doing so.

Application of satellite NDVI as a metric for spatial and temporal patterns of forage condition within the home range (i.e., third-order selection) [72] was problematic in our study area. A matrix of rock, trees, and snow, where NDVI is noninformative for indexing forage conditions, is common at the spatial scale of available satellite sensors [53]. Those cover types composed 49.9% of random locations in the TM images we used. Consequently, one-half of pixels were assigned the reference NDVI value of zero, reducing the power to detect selection among the remaining categories of NDVI.

We also did not postprocess TM images to convert raw radiance to reflectance, decreasing the power to detect variation in NDVI. Application of those computationally intensive radiometric corrections for each image entails an added layer of complexity when different sensors are used (e.g., Landsat 5 TM and Landsat 7 ETM+) but is essential for temporal comparisons among images [85]. Our analyses, however, did not make between-image comparisons, because RSFs were calculated with “used” versus “available” locations conditioned within the common temporal window of the TM scene, thereby reducing the need for corrections.

4.4. Hypothesis III: Tradeoff between Forage versus Risk of Predation Was Reduced during Drought

The hypothesis that drought reduced the forage benefit of migration to low elevations was supported; significantly greater forage benefits resulted from migration to low elevations, where predation risk was highest, in above-average precipitation years compared with the benefits accrued in lower than average precipitation years. Differences in vegetation condition, as indexed by NDVI at low-elevation locations used by bighorn sheep minus NDVI at high-elevation locations used by bighorn sheep, were greatest in years with above-average precipitation (Figure 6(b)). The NDVI tradeoff had a positive linear relationship with quantity of precipitation, with variation largely explained ( adjusted = 0.51) by cumulative precipitation to date and precipitation the prior water year (Figure 7). The greatest forage benefit of migration, therefore, accrued in sequential years having above-average precipitation, while the lowest benefit was present in consecutive drought years. Conversely, Wehausen [41] concluded that it was the timing, rather than quantity, of precipitation that best correlated with crude protein content in bighorn sheep feces, a proxy for dietary quality [55]. His investigation at Mount Baxter, however, identified a measure of Julian date itself as explanatory of most variation in fecal crude protein [41] () and went on to conclude that the best single predictor () of fecal crude protein was the number of days elapsed after a substantial precipitation event, termed, “the initiating storm” [105]. The conclusion of the import of timing is suspect, however, as this measure was simultaneously the predictor and predicted, because it was predicated on the occurrence of a storm event of sufficient magnitude (i.e., quantity of precipitation), defined by Wehausen [105] as having, “initiated perennial grass growth.”

Locations used by bighorn sheep at high elevation (>3,000 m) during winter often were covered by snow during years of above-average precipitation, leaving a limited number of snow-free pixels in wind-blown or south-facing areas. Consequently, forage conditions could be indexed by NDVI only in small samples of pixels within TM scenes, yielding additional variation to inherent variability among TM scenes (Figure 6). Application of within-scene comparisons of NDVI, calculated from uncorrected TM radiance, however, nullified temporal variation among NDVI values and allowed comparison of forage tradeoff at static high and low-elevation locations used by bighorn sheep in winter.

4.5. How Results Challenge Predator-Induced Range Abandonment

Our models of RSF by bighorn sheep were robust because of pooling of data among diverse winter ranges at Mount Baxter and Wheeler Ridge and because the years 2002–2007, when data were collected from GPS-collars, were characterized by highly variable precipitation. The driest (2002) and wettest (2005) years since 1987 (Figure 9), when Wehausen [41] reported the disuse of low-elevation winter range, occurred during this 6-year span. Wehausen [41] asserted that that use of low-elevation habitat in winter was the rule, unaffected by sequential years of drought, but this is not consistent with our data (Figure 3) or past observations [86].

In our study, one of the strongest positive predictors of resource selection by Sierra Nevada bighorn sheep was a measure of direct risk of predation by mountain lions, implicating the activity of mountain lions as a strong predictor of habitat selection by Sierra Nevada bighorn sheep—as posited by Wehausen [41]—but in the opposite direction he suggested. Our analysis demonstrated that rather than abandon preferred habitat, bighorn sheep continued to select the same areas where mountain lions actively hunted in bighorn sheep winter range (Figure 4). We established an alternative to the notion of predator-induced range abandonment [41], by contrasting forage tradeoff in wet versus dry years, revealing significantly reduced benefits of migration to low elevation in drought years (Figure 6). The period of “range abandonment” [41] occurred during 6 consecutive years of drought (1987–1992 [48]) when snow cover, which was avoided by female bighorn sheep, was sparse at high elevations, and forage benefits of migration to low elevation were significantly lower than in years of above-average precipitation.

Top-down effects, in the form of removal by helicopter capture for translocation, however, may have played a part in the proposed range abandonment in addition to the role of drought as a bottom-up driver. Intensive helicopter capture and removal of 103 bighorn sheep on low elevation winter range at Mount Baxter [41, 106] preceded the 1987–1995 period of proposed range abandonment [41]. Those captures not only preceded but indeed continued into March 1988, when eight female and three male bighorn sheep were removed [106109] from the low-elevation range that Wehausen [41] inexplicably defined as then “abandoned.” Helicopter capture and translocation was analogous to “direct predation,” in that the result of removing animals continuing to use low-elevation winter range was to leave only those animals using a parallel strategy of overwintering at high elevation. Furthermore, while short-term disturbance and flight responses by mountain sheep are well documented in response to low-level helicopter flights [107110], cumulative effects of repeated helicopter captures are unknown. Repeated helicopter captures that remove conspecifics may alter behavior and elicit responses analogous to those behaviors that mitigate indirect risk of predation, such as avoidance of areas lacking escape terrain, such as those low-elevation areas subject to recurrent helicopter captures. More research is warranted to elucidate behavioral responses to repeated helicopter capture efforts, especially for this endangered taxon.

Our alternative bottom-up hypothesis is supported by data demonstrating a reduced tradeoff by bighorn sheep remaining in snow-free, high-elevation winter habitat during drought years, providing an empirically and theoretically consistent explanation free from the inconsistencies of the predator-induced hypothesis of range abandonment. The bottom-up explanation resolves the inherent conflict between the range abandonment hypothesis and the theory of natural selection: maladaptive behaviors cannot persist in natural systems. Mountain lions share millennia of sympatry with bighorn sheep in North America [111, 112], while hydroclimatic reconstruction demonstrates the nonuniqueness of severe drought [113] over the same period. Thus, evidence does not support a novel set of conditions during the 1987–1995 period of range abandonment. The proposition that predator-induced range abandonment led to demographic catastrophe fundamentally conflicts with the concept of an evolutionarily stable strategy [15] for resource selection by bighorn sheep under risk of predation by a non-novel predator, the mountain lion. If avoidance of low-elevation winter range was, on average, more costly in terms of nutritional tradeoff, energetic demands, or risks of stochastic mortality inherent to life at high elevation in winter [41], compared with use of low-elevation winter range, such a strategy would experience strong negative selection. Natural selection would favor adoption of a “cheater” strategy, wherein some animals make use of low-elevation winter range and benefit from nutritious forage, while avoiding the perils of life at high elevation. Survivors among these “cheaters” would be strongly selected for. To support the claim that Sierra Nevada bighorn sheep responded to risk of predation by avoiding resources critical to survival and reproduction (e.g., low-elevation habitat in winter) [41] in a manner incompatible with an ESS, extraordinary evidence is necessary but remains lacking.

The bottom-up hypothesis explaining avoidance of low-elevation winter range by bighorn sheep is consistent with the findings of Kucera [114] and Pierce et al. [47], which documented drought-induced reduction in forage available to mule deer wintering in Round Valley, concurrent with the hypothetical predator-induced range abandonment by bighorn sheep [41]. The winter range of the Round Valley mule deer herd overlapped one-third of the winter range of bighorn sheep at Wheeler Ridge [50]. Kucera [114] concluded that drought and limited availability of forage in winter, and not predation by mountain lions, caused the 50% decline in the number of mule deer wintering adjacent to Wheeler Ridge during 1985–1988. Pierce et al. [47] analyzed long-term demographic and physiological conditions in this mule deer population during 1984–1998 and concluded that bottom-up processes, including prolonged drought, regulated mule deer populations through 1990, when the mule deer population began to recover. Diets of mule deer were dominated by bitterbrush [114, 115], also a major forage species for bighorn sheep in winter [18]. Kucera noted an absence of leader growth in bitterbrush during the 1986-1987 drought, as well as heavily browsed shrubs and frequent observations of mule deer feeding while standing on hind legs to reach bitterbrush leaves [114]. In the absence of data refuting drought-induced limits to forage as the major factor responsible for elevational shifts in use of winter range [41], bottom-up processes offer a stronger alternative explanation for declines in Sierra Nevada bighorn sheep.

Investigations into effects of behavioral responses to risk of predation in large herbivores have yielded diverse results, leaving many unanswered questions regarding the applicability and magnitude of effects from trait-mediated, indirect interactions in terrestrial systems [116]. Ripple et al. [40] reported a behaviorally mediated trophic cascade (BMTC) in elk in the Greater Yellowstone Ecosystem (hereafter, “Yellowstone”), wherein elk responded to risk of predation by wolves by avoiding risky stands of quaking aspen (Populus tremuloides), thereby reducing herbivory and allowing increased recruitment of aspens. Kauffman et al. [117] evaluated evidence for a BMTC of sufficient magnitude to structure aspen communities in Yellowstone, by assessing foraging by elk under risk of predation by wolves and by comparing historical aspen recruitment inside and outside of exclosures, and concluded that historic precipitation and number of elk, and not a BMTC, explained recruitment of aspen.

Similar conflicting results characterize studies of physiological and demographic costs accrued from energetic demands of antipredator behavior. Creel et al. [38] investigated the correlation between risk of predation by wolves and recruitment in elk in the Yellowstone system, with results indicating a negative relationship between predation risk and reproduction. Middleton et al. [39], conversely, detected no significant effects on pregnancy or body condition in elk subject to frequent proximity and risk of predation by wolves in the Yellowstone system. Those examples illustrate the possibility that conclusions about prey behavior and demography attributed to indirect risk of predation may, upon further examination, be better explained by bottom-up processes and direct effects of predation. We find such a bottom-up explanation of resource selection by Sierra Nevada bighorn sheep far more plausible, owing to empirical and theoretical support, than the predator-induced range abandonment hypothesis proffered by Wehausen [41].

Appendix

See Tables 5 and 6 and Figures 8 and 9.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank helicopter pilot S. DeJesus and California Department of Fish and Wildlife capture specialists R. Teagle and T. Glenner, for safe capture of bighorn sheep, and houndsmen J. Davis and J. Ostergard and their assistants V. Davis, M. Kiner, M. Leonard-Cahn, and J. Schmalenberger for their work capturing mountain lions. Additional capture support was provided by personnel from the CDFW Wildlife Investigations Laboratory including B. Gonzales, A. Hunter, K. Jones, L. Konde, J. Schultz., P. Swift, and S. Torres. They thank the many employees and volunteers who helped collect data, including C. Baker, L. Bowermaster, K. Ellis, J. Erlenbach, A. Feinberg, J. Fusaro, D. German, L. Greene, D. Jensen, R. Noles, N. Partridge, P. Partridge, C. Schroeder, A. Stephenson, and T. Taylor. Discussions with D. Delehanty, J. Kie, R. Long, C. Peterson, T. Stephenson, and J. Wehausen were valuable in guiding this investigation. They thank three anonymous reviewers for helpful comments that improved this paper. This study was funded by California Department of Fish and Wildlife Sierra Nevada Bighorn Sheep Recovery Program, California Deer Association, Desert Bighorn Council, California Chapter of the Wild Sheep Foundation, Department of Biological Sciences at Idaho State University, and the Institute of Arctic Biology at the University of Alaska Fairbanks. This is a contribution from the Sierra Nevada Bighorn Sheep Recovery Program and is Professional Paper 103 from the Eastern Sierra Center for Applied Population Ecology.

Supplementary Materials

ClusterDetector is a compressed file containing a software program used to identify spatial and temporal clusters of locations, as from GPS collars deployed on mountain lions. The archive contains a file, “readme.txt”, containing a detailed description of the program, the software, as both a Microsoft Windows executable and C++ source code, and a sample data file.

  1. Supplementary Material