Using empirical location data from individuals to model habitat quality and species distributions is valuable towards understanding habitat use of wildlife, especially for conservation and management planning. Incorporating measures of reproductive success or survival into these models helps address the role of vital rates (a surrogate of fitness) in affecting a species' distribution. We used 24-year datasets of Arctic peregrine falcon Falco peregrinus tundrius nest-site locations and productivity from the Colville River Special Area, Alaska, USA to model suitability of breeding habitat and the relative quality of used and potential nest sites. We used zero-inflated negative binomial regression models and covariates describing nest-site productivity, area of surrounding prey habitat, geology, topography and land-cover type to model and predict intensity of Arctic peregrine falcon nest-site use along the Colville River, and developed a predictive map of intensity of nest-site use. Regions of higher predicted intensity of use were characterized by steeper slopes, greater area of prey habitat, and higher average productivity, which are likely attributed to minimizing predation risk, gaining advantages for hunting, having sufficient prey resources, site quality, and overall fitness. Including productivity in intensity of nest-site use models improved the models, supporting our supposition that adding a fitness parameter enhanced the predictive capability of the species distribution model. Areas predicted to have higher intensity of use by our model can be used to focus efforts of continued protection of areas with frequently occupied and productive nest sites, and conversely, identify areas where protection of nest sites is likely to have few conservation benefits.
Species distribution modeling based on habitat use has become a valuable tool for wildlife management and conservation. Through identification of areas of higher use and their associated attributes, species distribution models have helped guide management efforts through better understanding of wildlife habitat-use patterns and leading to protection of important habitats (Dellinger et al. 2013, Squires et al. 2013). Such modeling may be used to examine covariates affecting spatial use of a variety of resources key to a species, including breeding, foraging, and wintering areas; movement corridors; and areas providing reduced predation risk (Bergman et al. 2006, Bruggeman et al. 2007, Dzialak et al. 2012).
Data input into species distribution models, and associated limitations, must be considered when evaluating their utility for management and conservation strategies. Empirical location data of individuals are used to parameterize species distribution models and can be of many forms, such as presence–absence, use–availability, presence-only and count data (Buckland and Elston 1993, Manly et al. 2002, Pearce and Boyce 2006). These models are a function of covariates describing attributes of each individual or location, which then provide information about variables related to a species' distribution, resource use or behavior (Manly et al. 2002). However, interpreting high quality habitats as those with high likelihood of use may be problematic because individual fitness is related to a complex set of factors, including habitat selection patterns that may sometimes be misleading (Gaillard et al. 2010, Nielsen et al. 2010). For example, maladaptive habitat selection occurs when individuals choose habitats that reduce reproductive success or survival (Misenhelter and Rotenberry 2000, Delibes et al. 2001, Streby et al. 2014b). Incorporating measures of fecundity and/or survival into species distribution models may provide a means of alleviating issues associated with maladaptive habitat selection (Mosser et al. 2009, DeCesare et al. 2014, Peterson et al. 2016). Limitations often exist on obtaining sufficient spatially-explicit data for quantifying measures of fitness; however, surrogates of fitness also may be useful for informing species distribution models. Roever et al. (2013) and McGreer et al. (2015), for example, quantified mortality and predation risk using location data of mortality events and predators, respectively, to further inform their resource selection function models and identify areas of higher conservation importance.
Species that breed in Arctic regions and at high latitudes have been the focus of recent conservation concern owing to potential impacts of climate change (Post et al. 2009). In part because of these potential climate-related impacts, species distribution models have been used to assess a variety of species of previous and current conservation concern in the Arctic (Ferguson et al. 2000, McLoughlin et al. 2002, Booms et al. 2010, Wilson et al. 2012). Many of these studies used combinations of landscape, climate and seasonal variables, along with breeding or parental state (e.g. pregnant or not pregnant; with or without young) to assess factors related to distribution, but to our knowledge none have incorporated fitness parameters into their models.
Arctic peregrine falcons Falco peregrinus tundrius (hereafter Arctic peregrine), which were a species of conservation concern during the 1950s–1990s, breed at high latitudes in portions of Alaska, Canada and Greenland (White 1968, US Dept of the Interior [USDOI] 2008). Arctic peregrine and many peregrine falcon F. peregrinus populations suffered drastic declines during the 1950s–1970s due primarily to exposure to DDT and other organochlorine pesticides that affected reproduction and survival (Ratcliffe 1970). Arctic peregrines were listed as endangered in 1973 under the US Endangered Species Act (ESA) and recovered sufficiently enough to be removed from the ESA in 1994 (US Fish and Wildlife Service 1994). Alaska's Colville River Special Area (CRSA) was established in 1977, during the period of population recovery, to conserve nesting and foraging habitat of Arctic peregrines (USDOI 2008). The CRSA is located within the National Petroleum Reserve in Alaska (NPR-A), which allows for oil and gas mining and exploration, and protective regulations for Arctic peregrines still exist under the CRSA Management Plan to minimize disturbance and preserve nesting and foraging habitat (USDOI 2008). However, additional information needs were identified to improve knowledge of Arctic peregrine ecology in the CRSA, better inform management decisions, and evaluate effectiveness of existing regulations (USDOI 2008).
We used 24 years of Arctic peregrine nest-site location and productivity data collected during the population's recovery (Cade et al. 2003) to develop a species distribution model of intensity of nest-site use with a surrogate fitness covariate for breeding Arctic peregrines throughout a portion of the CRSA. We incorporated productivity data from nest sites into models as a surrogate measure of fitness for two reasons. First, we used productivity data to identify areas with higher predicted intensity of nest-site use that were not confounded by potential maladaptive habitat selection occurring during the recovery. Second, during the recovery the Arctic peregrine population exhibited density dependence (Bruggeman et al. 2015, Swem and Matz 2018), which may negatively affect productivity. We also used knowledge of factors related to Arctic peregrine occupancy of nest sites and abundance on nesting cliffs in the CRSA (Bruggeman et al. 2015, 2016) to select covariates describing nest-site attributes (e.g. habitat; topography; prey habitat availability) to evaluate in models. The goals of our work were to: 1) model and predict Arctic peregrine productivity throughout a portion of our study area; 2) use results from 1) to help predict and map the intensity of Arctic peregrine nest-site use for the purpose of identifying areas in our study area along the Colville River predicted to have higher intensity of use for nesting; 3) evaluate the predictive capability of our model; and 4) assess how our findings related to current CRSA protective regulations.
Our study area consisted of a 347-km stretch of the Colville River and its surrounding landscape < 3 km from the river located in the 1 000 000-ha CRSA in Alaska, USA (centroid 69°18′43″N, 155°22′90″W, Fig. 1A; Bruggeman et al. 2016). Oil and gas exploration, fieldwork associated with monitoring Arctic peregrines and other natural resources, and recreation were primary activities in the CRSA during our study (USDOI 2008). The CRSA contained numerous wetlands and vegetation was characterized by tundra plant communities except for the Colville River floodplain, where willow Salix spp. and alder Alnus spp. communities coincided with perennial herb pioneer communities (Bliss and Cantlon 1957).
Migratory Arctic peregrines began arriving to the CRSA in late April, nested May–August on cliffs, bluffs and escarpments along the Colville River, and returned to wintering areas after young fledged in August and September (Ambrose and Riddle 1988). Ted Swem led two surveys per year for Arctic peregrines by boat along the Colville River during 1981, 1982, 1985, 1987–2002, 2005 and 2011. B. Dittrick, P. Schempf, and J. Silva led surveys in 1983, 1984 and 1986, respectively. Surveys were standardized among all years; no surveys were conducted in 2003, 2004 and 2006–2010. The first survey occurred during egg-laying and incubation in June; the second survey occurred during the nestling period in late July–early August. At each nest site encountered during each survey, observers counted numbers of adults and young (second survey only), mapped the nest-site location, and recorded location by GPS when feasible. We digitized nest-site locations into a GIS layer and assigned a measure of precision based on an assessment of location certainty.
We obtained GIS layers of elevation (US Geological Survey 2017), land cover (Homer et al. 2004), sub surficial geology (Beikman 1980), surficial geology (Karlstrom 1964), and streams in the CRSA. We used the elevation layer to generate aspect and slope layers at 10-m resolution in ArcGIS 9.2 (ESRI, CA, USA). We used the land-cover layer at 30-m resolution to classify each raster cell as open water, wetlands with woody vegetation, wetlands with emergent herbaceous vegetation, barren, dwarf scrub, shrub or ‘other’ categories (Homer et al. 2004).
In GIS, we delineated a 1152-ha sampling universe within the CRSA study area that: encompassed potential nesting cliffs, escarpments, and bluffs on both sides of the Colville River; spanned the length of the river repeatedly surveyed over 24 years; and excluded the Colville River and its major tributary streams (Fig. 1B). We generated a lattice grid of 12 748 765 points spaced 10-m apart spanning the sampling universe to use for model predictions. We used a grid with 10-m resolution to correspond to the resolution of the elevation, aspect and slope layers.
Using the annual productivity for each nest-site location (i.e. number of young enumerated in the nest during the second survey) recorded during 24 years of surveys, we calculated the average annual productivity for each nest-site location based on the number of years the nest site was occupied for surveys having occupied–unoccupied, unoccupied–occupied, or occupied–occupied patterns during the first and second surveys. We defined a response variable as the average annual productivity for each nest-site location, defined eight covariates (Table 1), extracted values for each covariate from each location, and developed linear regression models consisting of all possible covariate combinations. We fitted models in R (< www.r-project.org >) and selected the model with the highest adjusted-R 2 value as the best model (Neter et al. 1996). We used adjusted-R 2 values to rank and select models of average productivity because we were using linear regression models and adjusted-R 2 provides a measure of variability explained while accounting for the number of covariates in the model (Neter et al. 1996). We used the best model to estimate average productivity for each of the 10-m lattice points for the purpose of deriving a productivity covariate for our intensity of nest-site use modeling. To estimate average productivity, we extracted values from GIS layers for covariates included in the best model at each point, used the values to estimate average productivity at each point, and developed a raster map of hypothetical average productivity across the sampling universe.
Covariates used in analyses modeling factors related to average productivity of nest sites and intensity of nest-site use of Arctic peregrine falcons along the Colville River, AK, USA.
Based on guidance provided by Aarts et al. (2012) for modeling species distribution, we developed a grid of 123 173 100-m cells in GIS that spanned our sampling universe. The 100-m resolution represented a tradeoff between larger cells that may have contained more than one nest site and smaller cells that increased computational complexity and added little additional information to the model. We defined a response variable as the total number of years each cell was occupied by an Arctic peregrine nest site during the 24 years of surveys (i.e. intensity of nest-site use). Aarts et al. (2012) demonstrated that as the number of cells, or availability points, increases in a given area, the likelihood for the resulting Poisson generalized linear model becomes a discrete approximation of the inhomogeneous Poisson point process (IPP) likelihood. The IPP likelihood is a function of the intensity of presences (Warton and Shepherd 2010) and can be used to model count, presence–absence or use-availability data (Aarts et al. 2012).
We considered eight covariates (Table 1) and extracted values for each from each nest-site location or the centroid of each 100-m cell if no nest site was located in the cell. We also defined a covariate for average productivity (productivity, Table 1) estimated from the raster map (as described above) and extracted values for cell centroids. We used the observed average productivity value determined from surveys for the 108 cells containing a nest site and predicted average productivity for cells not containing a nest site. We developed 255 Poisson regression models consisting of all possible covariate combinations of the original eight covariates and with each model containing productivity. We included productivity in all models to incorporate a surrogate of fitness that accounted for both consistency and numbers of young produced. We fitted Poisson models in R and calculated an AIC value for each model (Burnham and Anderson 2002). We also examined the suitability of negative binomial (NB), zero-inflated Poisson (ZIP), and zero-inflated negative binomial (ZINB) distributions to our data. For zero-inflated models, we used the same model structure for both the count (Poisson) and zero-inflated (binomial) components of the model. We separately fitted models in R using package mass (Venables and Ripley 2002) for NB models and package pscl (Zeileis et al. 2008, Jackman 2015) for ZIP and ZINB models, calculated AIC values for each model, examined residual plots for best-approximating models with ΔAIC < 2 for each distribution, and compared AIC values among the four distributions. Based on the distribution that had best-approximating models with the lowest AIC values, we ranked and selected the best-approximating models using ΔAIC values and calculated an Akaike weight (w) for each model (Burnham and Anderson 2002). We used AIC values to rank and select models of intensity of nest-site use because a comparable R 2 statistic does not exist for zero-inflated models. We used the best-approximating model to estimate the intensity of Arctic peregrine nest-site use throughout the sampling universe by extracting values from GIS layers for covariates included in the best model at each 10-m lattice grid point and using those values to estimate intensity of nest-site use at each point. We then developed a raster map of intensity of use for the CRSA in ArcGIS. Again, we note that predicted intensity of use includes a productivity covariate as either average productivity for nest sites or predicted average productivity for locations with no observed nest sites. Because there are no formal validation measures for the type of model that was the best-approximating model, we assessed the predictive capability of our model using GIS-derived statistics summarizing measures of predicted intensity of use for each nest-site location as detailed in Supplementary material Appendix 1. To examine whether inclusion of productivity improved model fit, we fitted ZINB regression models without productivity, calculated an AIC value for each model, and compared AIC values to AIC values of our model results with productivity.
Arctic peregrine surveys
We detected Arctic peregrines at 108 unique nest-site locations during 24 years of surveys. Estimates of detection probability, as determined from previous work, were > 0.8 for the first survey in all years except 1982, when it was 0.7, and slightly lower for the second survey (Bruggeman et al. 2016). Total number of times nest sites were occupied over 24 years ranged from 1–24 (mean = 12.0; SE = 0.655). The number of nest sites at which we detected Arctic peregrines ranged from 28 in 1981, 1982, and 1983, to 69 in 2001 (mean = 52.5; SE = 2.89, n = 24). Total maximum number of adult Arctic peregrines enumerated during surveys increased during the 24-year survey period, ranging from 27 birds in 1982 to 121 birds in 1998 (mean = 84.2; SE = 5.72, n = 24). Productivity of individual nest sites ranged from 0–4 young (mean = 0.604; SE = 0.023; n = 2592 nest sites). Productivity was negatively correlated with year (p < 0.001) when not accounting for repeated territory observations, and autocorrelation function (ACF) values ranged from 0.010–0.126 over 24 years with minimum and maximum ACF values occurring at time lags of 11 years and 2 years, respectively (Supplementary material Appendix 1 Fig. A1). We summarize values of covariates used in modeling intensity of nest-site use in Table 2.
Predicting average productivity
The model of average productivity with highest adjusted-R 2 value (0.204) included aspect, elevgain, geology and subgeology covariates (Table 3) with R2 = 0.286 (F11,96 = 3.5, p < 0.001). The model with the second highest adjusted-R2 value (0.203) also contained aspect, geology and subgeology covariates, but included height instead of elevgain with R2 = 0.285 (F11,96 = 3.5, p < 0.001). Correlation between model-predicted average productivity (productivitymodel) and observed average productivity (productivitysurveys) recorded in 108 nest sites during surveys was R2 = 0.246 (F1,106 = 34.6, p < 0.001) with a line-of-best-fit productivitymodel = 0.799 + 0.277 × productivitysurveys. Model-predicted average productivity ranged from 0.077–3.05 (mean = 1.32; SD = 0.612; n = 12 748 765) across the sampling universe (Fig. 2). Annual average productivity was positively correlated with the number of years the nest site was occupied (Fig. 3A; p = 0.019, F1,106 = 5.71). However, there was no correlation between the variance of annual average productivity and the number of years the nest site was occupied (Fig. 3B; p = 0.758, F1,96 = 0.095).
Ranges, means, and standard errors for continuous covariates used to model intensity of Arctic peregrine falcon nest-site use along the Colville River, AK, USA during 1981–2002, 2005 and 2011. Summary statistics are separated into locations used by Arctic peregrine falcons as nest sites (i.e. nest sites observed during surveys) and locations classified as available for nesting. Covariates are defined in Table 1.
Covariate coefficient estimates and 95% confidence intervals from the best regression model of factors related to average productivity of Arctic peregrine falcon nest sites along the Colville River, AK, USA during 1981–2002, 2005 and 2011. Covariates are defined in Table 1.
Predicting and mapping intensity of nest-site use
We found ZINB models had the lowest AIC values among best-approximating models of intensity of nest-site use for all four distributions evaluated. There was one best-approximating model (AIC = 1329, w = 0.791) that included aspect, elevgain, habitat, productivity, slope, subgeology and waterarea covariates for both the Poisson and binomial model components (Table 4). For the Poisson component, productivity, slope, and north and south aspect had 95% confidence intervals (CI) that did not include 0, whereas waterarea, west aspect, and dwarf scrub habitat had 95% CI that contained, but were not centered on, 0 (Table 4). All seven covariates for the binomial component had 95% CI that did not include 0 (Table 4). There were no other competing models of intensity of nest-site use with ΔAIC < 2. The minimum AIC among models without productivity was 1341, indicating inclusion of the covariate improved the models.
Model-predicted intensity of Arctic peregrine nest-site use ranged from 0–17.5 (mean = 0.005; SD = 0.111, n = 12 748 765) across the sampling universe (Fig. 4,; see Supplementary material Appendix 2 for more detailed maps). Along the upriver portion of the Colville River in our study area, most areas with highest predicted intensity of use were isolated cliffs and bluffs that were previously used by Arctic peregrines as nest sites (Fig. 4). However, there were several upriver areas not previously occupied by Arctic peregrines, including some farther from the river, which had a higher predicted intensity of use (Fig. 4). Downriver, areas with highest predicted intensity of use were more extensive and located along larger cliffs and bluffs, including those previously occupied by Arctic peregrines (Fig. 5). Of the 108 nest-site locations, between 17.6–21.3% were classified as high observed use/high predicted use and 33.3% were classified as low observed use/low predicted use (Supplementary material Appendix 1 Table A1).
We developed an empirical model of Arctic peregrine nest-site distribution by incorporating a measure of productivity (a surrogate of fitness) to identify areas of higher predicted intensity of nesting use in our study area along the Colville River, which provides nesting habitat for ~25% of Alaska's migratory Arctic peregrine population. Our work provides an example of using an alternative approach with count data to model and map a species' distribution, which is often done using logistic regression analyses based on point location data (Aarts et al. 2012). Depending on the extent of logistic regression analyses and whether they used repeated measures techniques to account for multiple observations at the same location, they may not impart extra importance to locations used more than once or, in our case, that had relatively high productivity. Our intensity of nest-site use model incorporated productivity data, which resulted in model predictions of areas both likely to be occupied by Arctic peregrines and where there was higher nest site productivity to address the role of fitness in resource selection. Incorporating this fitness metric improved the predictive capability of the model, supporting our original supposition.
Covariate coefficient estimates and 95% confidence intervals (CI) from the best-approximating zero-inflated negative binomial regression model of factors related to intensity of Arctic peregrine falcon nest-site use along the Colville River, AK, USA during 1981–2002, 2005 and 2011. Estimates for the Poisson (count) component and binomial (zero-inflated) component are provided for each covariate. Covariates are defined in Table 1.
Selection of habitat and resources by animals should maximize fitness, which is generally measured empirically through reproductive success or adult survival (Gaillard et al. 2010), but can also be assessed through mortality (predation) risk (Nielsen et al. 2010), juvenile survival (Streby et al. 2014a, 2014b), litter (clutch) size (Forsman et al. 2007), or daily nest survival (Chalfoun and Martin 2007). The value of incorporating fitness parameters (or surrogates of fitness) into species distribution models addresses the possibility there is not a positive relationship between an individual's selection of resources and their associated benefits (i.e. maladaptive selection). Relying on animal density alone to provide a measure of habitat quality or the habitats that confer the highest fitness may be problematic because density, resource availability, and individual fitness interact to determine an animal's choice of habitats and associated habitat quality (Van Horne 1983, Mosser et al. 2009). We used nest-site productivity as a surrogate of fitness because we did not have site-specific data on either juvenile survival (i.e. recruitment) or number of young fledged per nest (i.e. nest success) owing to logistical challenges due to the large length of river surveyed and timing of surveys to determine fledgling survival or success of all nest sites. Because nest-site productivity provided only a measure of the number of young in the nest at the time of the survey and not the number of young raised to independence in a breeding season, fledgling survival would have been a more accurate surrogate of individual fitness (Streby et al. 2014b).
Individual fitness components may be highly important to nest-site selection with the selection process influenced by abiotic and biotic factors acting across multiple spatial scales (Newton 1985, Clark and Shutler 1999, Martínez et al. 2003). Minimizing predation risk, protecting the nest from inclement weather, and having sufficient nearby prey resources are factors that are related to site quality and may affect an individual's choice of nest-site location (Møller 1988, Martin 1995, Kolbe and Janzen 2002). Occupancy of nest sites has been associated with site quality, which may be related to individual fitness (Sergio and Newton 2003), and it is likely the Arctic peregrine nest sites used most often during our study and areas predicted to have higher use by our model are indicative of higher-quality habitat, as evidenced by the positive relationship between average productivity and number of years the site was occupied (Fig. 3A). However, the implications of nest-site choice on fitness are complicated and may involve tradeoffs to balance opposing selection pressures (Streby et al. 2014b), resulting in choices among individual fitness components to maximize total fitness (Schluter et al. 1991, Chalfoun and Schmidt 2012). Total reproductive fitness for an individual or breeding pair is dependent on fledgling survival after a nest is successful and, depending on the species, a female may make multiple attempts at nesting to produce a successful nest. Arctic peregrines are not known to make multiple nesting attempts within one breeding season in the CRSA, so the results of our study and previous work (Bruggeman et al. 2015, 2016, Swem and Matz 2018) lend insight into relationships between abiotic and biotic factors and nest-site decisions of breeding pairs to maximize total fitness. In addition to steeper slopes and greater area of surrounding prey habitat, these factors include a lack of snow for establishing a nest site early in the season to increase the likelihood of nest success, locating nests on cliffs situated higher above the Colville River to reduce predation risk on juveniles and increase advantage in capturing prey, and establishing nest sites in locations where warmer temperatures may minimize negative effects of cold rain and wind on juvenile survival (Franke et al. 2010, Anctil et al. 2014, Bruggeman et al. 2015).
Spatial patterns of higher predicted Arctic peregrine use in the CRSA were generally consistent with observations of nesting habitat used most frequently during surveys. The western, upriver portion of the Colville River in our study area was characterized by smaller, discrete, separated cliffs that often only provided sufficient space for one nesting pair (Swem and Matz 2018). Upriver, surveyors occasionally found nest sites on small outcroppings expected to be unsuitable for nesting. In contrast, larger more extensive cliffs and escarpments supporting multiple nesting pairs typified the eastern downriver segment of the Colville River. Nesting density of Arctic peregrines was also greater downriver compared to upriver (Fig. 4, 5), suggesting downriver cliffs offered some combination of greater resource availability, better access to nest sites, higher quality habitat, or greater fitness benefits. Our models of predicted average productivity (Fig. 2) and intensity of nest-site use (Fig. 4, 5) support these observations with areas of higher predicted use being fewer in number, smaller in size, and more spatially separated upriver relative to downriver. There were also more nest-site locations upriver than downriver for which the model did not predict use and it is possible that factors other than those related to productivity and nest-site use (e.g. factors related to fledgling survival) influenced nest-site selection at these locations.
With inclusion of productivity data into analyses to account for fitness, our model correctly predicted between 51–55% of nest-site locations that had high and low observed use combined. There are limitations of our study that may explain why our best-approximating model did not accurately predict ~45% of nest-site use. First, factors that we could not quantify, such as competition for nest sites, territory size, biotic characteristics (e.g. prey abundance), or presence of snow early during the nesting period, may have influenced nesting decisions. Arctic peregrines require a snow-free substrate at nest scrapes and it is possible observed nest-site locations not predicted to have use were selected because other potential sites were covered with snow when nest-site selection occurred. Previous work documented Arctic peregrine apparent survival and arrival rates on cliffs were higher during years with earlier snowmelt and milder winters based on broad-scale climate patterns (Bruggeman et al. 2015). Second, the CRSA on the North Slope of Alaska is remote and GIS layers for the region are generally updated less frequently and available at lower resolution than other areas in the US. Because Arctic peregrines and other birds may select nest sites based on fine-scale attributes (Brambilla et al. 2006), we were not able to fully incorporate these factors in our modeling. Finally, we were not able to incorporate detection probability into our average productivity and intensity of nest-site use models. Although detection probability was high, especially during the first survey (Bruggeman et al. 2016), our estimates of average productivity and intensity of nest-site use could be biased low, resulting in underestimation of both variables.
Improving understanding of factors affecting spatial dynamics of a species and visualizing this distribution using GIS-based mapping has become a valuable tool to assist with conservation for wildlife managers. Incorporating fitness (or surrogates of fitness) parameters has increasingly become an important part of species distribution modeling to better understand the roles of reproductive success, adult and juvenile survival, and mortality risk in affecting habitat selection (Gaillard et al. 2010). However, acquiring the data needed to include fitness parameters in models may be resource intensive. Therefore, there are tradeoffs between the benefit of using fitness-related data in species distribution analyses and the added cost for data collection. Our work provides an example of applying a long-term dataset of occupancy and productivity to model, predict, and map nest-site distribution of a once endangered population during its recovery. Our models provide a means of identifying factors associated with occupancy and productivity, which is information that can help target conservation efforts where they are likely to be most effective.
The CRSA Management Plan provides guidance for the protection of Arctic peregrine nest sites, nesting cliffs, and foraging habitat that governs all nest sites and habitat under the same regulations (USDOI 2008). For example, a 1600-m buffer exists along the Colville River in which no permanent oil and gas facilities and related development are allowed (USDOI 2008). Our results identified areas along the Colville River predicted to have higher intensity of use for Arctic peregrine nesting to help focus efforts of continued protection and preservation of these cliffs, especially those with frequently occupied and highly productive nest sites. For cliffs with no predicted use, particularly where no observed nest sites have been located during surveys, considerations could be given for relaxing some protective measures after evaluating historical and current Arctic peregrine use of the cliff and associated nest-site productivity. In areas where observed Arctic peregrine nest-site locations overlapped those of higher predicted use, management actions are likely to have the most impact on population-level processes. In contrast, areas where observed nest-site locations overlapped those of low predicted use may have lower conservation priority and be where management actions may have less effect on population dynamics.
Arctic peregrines failed to produce any young at 11 of the 108 nest-site locations during the study, suggesting these sites had minimal or no contribution at the population level. These locations were distributed along the entire Colville River in our study area and not necessarily concentrated downriver where nesting density was higher and Arctic peregrines may have been more likely to establish alternative nest sites during years of greater competition. The broad distribution of unproductive nest-site locations, as opposed to being concentrated along one stretch of river, may make relaxing protective measures challenging because of potential negative impacts to nearby, productive nest sites. This is especially relevant downriver where some unproductive nest sites occurred among higher densities of high use, productive nest sites. Because a substantial number of observed nest sites were not correctly classified in our analysis, we do not recommend changing regulations for cliffs with observed, productive nest sites that were predicted as having low intensity of use. We classified low intensity of use for 49 observed nest-site locations, of which the model predicted 36 to have low use. The majority of the 36 low-use nest sites had low total productivity (< 10 young) and average productivity (< 2 young) even if occupied several years, suggesting these locations or the birds that occupied them were of poor quality.
Regulation changes around consistently unproductive nest-site locations and those not frequently occupied may have minimal population-level impacts, but additional assessments of disturbance on Arctic peregrines in the CRSA may be warranted. It is also possible unproductive nest sites could be productive in future years under different climate scenarios or Arctic peregrine population sizes. Our models were based on data collected during a period when fossil fuel exploration was minimal and human activity and potential related disturbance were limited to occasional recreational boaters and natural resources field workers. The effects of human activity on birds vary widely depending on the species, a particular population's extent of habituation to the activity, type of human activity (e.g. vehicle or foot travel), timing of activity during the phase of the breeding season (i.e. courtship, incubation, fledging), and degree of cover provided by vegetation or topography (Boyle and Samson 1985, Richardson and Miller 1997, Brambilla et al. 2004). Our work provides information related to spatial variation in nesting for the CRSA Arctic peregrine population, and our model provides a tool to help minimize potential negative impacts of human activity on a breeding Arctic peregrine population while also identifying areas where relaxing restrictions on human activity will potentially have minimal negative impacts. There is potential to use this tool for other populations with long-term datasets that include fitness or surrogate fitness parameters.
We thank J. Fieberg for advice on modeling; C. Hamfler for providing GIS data and support; B. Dittrick, P. Schempf, and J. Silva for their lead in surveys in 1983, 1984 and 1986; and 33 colleagues who provided insight and helped collect data. We also thank K. Proffitt for useful comments on a previous draft of this manuscript.
Funding – Provided by the Bureau of Land Management and U.S. Fish and Wildlife Service in Fairbanks, AK. Any use of trade, firm or product names is for descriptive purposes only and does not imply endorsement by the US Government, University of Minnesota, or Oregon State University.
- Aarts, G. et al. 2012. Comparative interpretation of count, presence-absence and point methods for species distribution models. – Methods Ecol. Evol. 3: 177–187. Google Scholar
- Ambrose, R. E. and Riddle, K. E. 1988. Population dispersal, turnover, and migration of Alaska Peregrines. – In: Cade, T. J. et al. (eds), Peregrine falcon populations: their management and recovery. The Peregrine Fund, pp. 677–684. Google Scholar
- Anctil, A. et al. 2014. Heavy rainfall increases nestling mortality of an Arctic top predator: experimental evidence and long-term trend in peregrine falcons. – Oecologia 174: 1033–1043. Google Scholar
- Beikman, H. M. 1980. Geological map of Alaska. – US Geological Survey. Google Scholar
- Bergman, E. J. et al. 2006. Assessment of prey vulnerability through analysis of wolf movements and kill sites. – Ecol. Appl. 16: 273–284. Google Scholar
- Bliss, L. C. and Cantlon, J. E. 1957. Succession on river alluvium in northern Alaska. – Am. Midl. Nat. 58: 452–469. Google Scholar
- Booms, T. L. et al. 2010. Gyrfalcon nest distribution in Alaska based on a predictive GIS model. – Polar Biol. 33: 347–358. Google Scholar
- Boyle, S. A. and Samson, F. B. 1985. Effects of nonconsumptive recreation on wildlife: a review. – Wildl. Soc. Bull. 13: 110–116. Google Scholar
- Brambilla, M. et al. 2004. Rock climbing and raven Corvus corax occurrence depress breeding success of cliff-nesting peregrines Falco peregrinus. – Ardeola 51: 425–430. Google Scholar
- Brambilla, M. et al. 2006. Factors affecting breeding habitat selection in a cliff-nesting peregrine Falco peregrinus population. – J. Ornithol. 147: 428–435. Google Scholar
- Bruggeman, J. E. et al. 2007. Covariates affecting spatial variability in bison travel behavior in Yellowstone National Park. – Ecol. Appl. 17: 1411–1423. Google Scholar
- Bruggeman, J. E. et al. 2015. Dynamics of a recovering Arctic bird population: the importance of climate, density dependence and site quality. – Ecol. Appl. 25: 1932–1943. Google Scholar
- Bruggeman, J. E. et al. 2016. Multi-season occupancy models identify biotic and abiotic factors influencing a recovering Arctic peregrine falcon Falco peregrinus tundrius population. – Ibis 158: 61–74. Google Scholar
- Buckland, S. T. E. and Elston, D. A. 1993. Empirical models for the distribution of wildlife. – J. Appl. Ecol. 30: 478–495. Google Scholar
- Burnham, K. P. and Anderson, D. R. 2002. Model selection and multi-model inference. – Springer. Google Scholar
- Cade, T. J. et al. 2003. Return of the peregrine: a saga of North American tenacity and teamwork. – The Peregrine Fund. Google Scholar
- Chalfoun, A. D. and Martin, T. E. 2007. Assessments of habitat preferences and quality depend on spatial scale and metrics of fitness. – J. Appl. Ecol. 44: 983–992. Google Scholar
- Chalfoun, A. D. and Schmidt, K. A. 2012. Adaptive breeding-habitat selection: is it for the birds? – Auk 129: 589–599. Google Scholar
- Clark, R. G. and Shutler, D. 1999. Avian habitat selection: pattern from process in nest-site use by ducks? – Ecology 80: 272–287. Google Scholar
- DeCesare, N. J. et al. 2014. Linking habitat selection and predation risk to spatial variation in survival. – J. Anim. Ecol. 83: 343–352. Google Scholar
- Delibes, M. et al. 2001. Effects of an attractive sink leading into maladaptive habitat selection. – Am. Nat. 158: 277–285. Google Scholar
- Dellinger, J. A. et al. 2013. Habitat selection of a large carnivore, the red wolf, in a human-altered landscape. – Biol. Conserv. 157: 324–330. Google Scholar
- Dzialak, M. R. et al. 2012. Temporal and hierarchical spatial components of animal occurrence: conserving seasonal habitat for greater sage-grouse. – Ecosphere 3(4): 30. Google Scholar
- Ferguson, S. H. et al. 2000. Influence of sea ice dynamics on habitat selection by polar bears. – Ecology 81: 761–772. Google Scholar
- Forsman, J. T. et al. 2007. Mechanisms and fitness effects of interspecific information use between migrant and resident birds. – Behav. Ecol. 18: 888–894. Google Scholar
- Franke, A. et al. 2010. Long-term trends of persistent organochlorine pollutants, occupancy and reproductive success in peregrine falcons (Falco peregrinus tundrius) breeding near Rankin Inlet, Nunavut, Canada. – Arct. 63: 442–450. Google Scholar
- Gaillard, J.-M. et al. 2010. Habitat–performance relationships: finding the right metric at a given spatial scale. – Phil. Trans. R. Soc. B 365: 2255–2265. Google Scholar
- Homer, C. et al. 2004. Development of a 2001 national land cover database for the United States. – Photogramm. Eng. Remote Sens. 70: 829–840. Google Scholar
- Jackman, S. 2015. pscl: classes and methods for R developed in the political science computational laboratory, Stanford University. – <pscl.stanford.edu/> accessed 10 September 2016. Google Scholar
- Karlstrom, T. N. V. 1964. Surficial geology map of Alaska. – US Geological Survey. Google Scholar
- Kolbe, J. J. and Janzen, F. J. 2002. Impact of nest-site selection on nest success and nest temperature in natural and disturbed habitats. – Ecology 83: 269–281. Google Scholar
- Manly, B. F. J. et al. 2002. Resource selection by animals: statistical design and analysis for field studies. – Kluwer Academic Publishers. Google Scholar
- Martin, T. E. 1995. Avian life history evolution in relation to nest sites, nest predation and food. – Ecol. Monogr. 65: 101–127. Google Scholar
- Martínez, J. A. et al. 2003. Predictive model for habitat preferences of the Eurasian eagle owl Bubo bubo: a multiscale approach. – Ecography 26: 21–28. Google Scholar
- McGreer, M. T. et al. 2015. Selection for forage and avoidance of risk by woodland caribou (Rangifer tarandus caribou) at coarse and local scales. – Ecosphere 6(12): 288. Google Scholar
- McLoughlin, P. D. et al. 2002. Hierarchical habitat selection by barren-ground grizzly bears in the central Canadian Arctic. – Oecologia 132: 102–108. Google Scholar
- Misenhelter, M. D. and Rotenberry, J. T. 2000. Choices and consequences of habitat occupancy and nest site selection in sage sparrows. – Ecology 81: 2892–2901. Google Scholar
- Møller, A. P. 1988. Nest predation and nest site choice in passerine birds in habitat patches of different size: a study of magpies and blackbirds. – Oikos 53: 215–221. Google Scholar
- Mosser, A. et al. 2009. Serengeti real estate: density vs fitness-based indicators of lion habitat quality. – Ecol. Lett. 12: 1050–1060. Google Scholar
- Neter, J. et al. 1996. Applied linear statistical models. – McGraw-Hill. Google Scholar
- Newton, I. 1985. Lifetime reproductive output of female sparrowhawks. – J. Anim. Ecol. 54: 241–253. Google Scholar
- Nielsen, S. E. et al. 2010. Dynamic wildlife habitat models: seasonal foods and mortality risk predict occupancy–abundance and habitat selection in grizzly bears. – Biol. Conserv. 143: 1623–1634. Google Scholar
- Pearce, J. L. and Boyce, M. S. 2006. Modelling distribution and abundance with presence-only data. – J. Appl. Ecol. 43: 405–412. Google Scholar
- Peterson, S. M. et al. 2016. Spatially-explicit models of full-season productivity and implications for landscape management of golden-winged warblers in the western Great Lakes region. – In: Streby, H. M. et al. (eds), Studies in avian biology. CRC Press, pp. 141–160. Google Scholar
- Post, E. et al. 2009. Ecological dynamics across the Arctic associated with recent climate change. – Science 325: 1355–1358. Google Scholar
- Ratcliffe, D. A. 1970. Changes attributable to pesticides in egg breakage frequency and eggshell thickness in some British birds. – J. Appl. Ecol. 7: 67–115. Google Scholar
- Richardson, C. T. and Miller, C. K. 1997. Recommendations for protecting raptors from human disturbance: a review. – Wildl. Soc. Bull. 25: 634–638. Google Scholar
- Roever, C. L. et al. 2013. Incorporating mortality into habitat selection to identify secure and risky habitats for savannah elephants. – Biol. Conserv. 164: 98–106. Google Scholar
- Schluter, D. et al. 1991. Conflicting selection pressures and life history tradeoffs. – Phil. Trans. R. Soc. B 246: 11–17. Google Scholar
- Sergio, F. and Newton, I. 2003. Occupancy as a measure of territory quality. – J. Anim. Ecol. 72: 857–865. Google Scholar
- Squires, J. R. et al. 2013. Combining resource selection and movement behavior to predict corridors for Canada lynx at their southern range periphery. – Biol. Conserv. 157: 187–195. Google Scholar
- Streby, H. M. et al. 2014a. Redefining reproductive success in songbirds: moving beyond the nest success paradigm. – Auk 131: 718–726. Google Scholar
- Streby, H. M. et al. 2014b. Retirement investment theory explains patterns in songbird nest-site choice. – Phil. Trans. R. Soc. B 281: 20131834. Google Scholar
- Swem, T. and Matz, A. C. 2018. Natural history and recovery of Arctic peregrine Falcons along the Colville River, Alaska, 1981–2011. – In: Shuford, W. D. et al. (eds), Avifaunal change in western North America, studies in western birds, No. 3. Western Field Ornithologists, pp. 132–147. Google Scholar
- US Dept of the Interior, USDOI 2008. Colville River Special Area management plan and environmental assessment. – Bureau of Land Management Arctic Field Office. Google Scholar
- US Fish and Wildlife Service 1994. Endangered and threatened wildlife and plants; removal of the Arctic peregrine falcon from the list of endangered and threatened wildlife. – Fed. Regist. 59: 50796–50805. Google Scholar
- US Geological Survey 2017. National Elevation Dataset. – <nationalmap.gov/elevation.html>. Accessed 15 Sep 2017. Google Scholar
- Van Horne, B. 1983. Density as a misleading indicator of habitat quality. – J. Wildl. Manage. 47: 893–901. Google Scholar
- Venables, W. N. and Ripley, B. D. 2002. Modern applied statistics with S. – Springer. Google Scholar
- Warton, D. I. and Shepherd, L. C. 2010. Poisson point process models solve the “pseudo-absence problem” for presence-only data in ecology. – Ann. Appl. Stat. 4: 1383–1402. Google Scholar
- White, C. M. 1968. Diagnosis and relationships of the North American tundra-inhabiting peregrine falcons. – Auk 2: 179–191. Google Scholar
- Wilson, R. R. et al. 2012. Summer resource selection and identification of important habitat prior to industrial development for the Teshekpuk caribou herd in northern Alaska. – PLoS One 7(11): e48697. Google Scholar
- Zeileis, A. et al. 2008. Regression models for count data in R. – J. Stat. Softw. 27: 1–25. Google Scholar
Supplementary material (available online as Appendix wlb.00475 at < www.wildlifebiology.org/appendix/wlb-00475 >). Appendix 1–2.