As climate change alters the global environment, it is critical to understand the relationship between shifting climate suitability and species distributions. Key questions include whether observed changes in population abundance are aligned with the velocity and direction of shifts predicted by climate suitability models and if the responses are consistent among species with similar ecological traits. We examined the direction and velocity of the observed abundance-based distribution centroids compared with the model-predicted bioclimatic distribution centroids of 250 bird species across the United States from 1969 to 2011. We hypothesized that there is a significant positive correlation in both direction and velocity between the observed and the modeled shifts. We then tested five additional hypotheses that predicted differential shifting velocity based on ecological adaptability and climate change exposure. Contrary to our hypotheses, we found large differences between the observed and modeled shifts among all studied bird species and within specific ecological guilds. However, temperate migrants and habitat generalist species tended to have higher velocity of observed shifts than other species. Neotropical migratory and wetland birds also had significantly different observed velocities than their counterparts, which may be due to their climate change exposure. The velocity based on modeled bioclimatic suitability did not exhibit significant differences among most guilds. Boreal forest birds were the only guild with significantly faster modeled-shifts than the other groups, suggesting an elevated conservation risk for high latitude and altitude species. The highly idiosyncratic species responses to climate and the mismatch between shifts in modeled and observed distribution centroids highlight the challenge of predicting species distribution change based solely on climate suitability and the importance of non-climatic factors traits in shaping species distributions.

Effectively managing take of wildlife resulting from human activities poses a major challenge for applied conservation. Demographic data essential to decisions regarding take are often expensive to collect and are either not available or based on limited studies for many species. Therefore, modeling approaches that efficiently integrate available information are important to improving the scientific basis for sustainable take thresholds. We used the prescribed take level (PTL) framework to estimate allowable take for bald eagles (Haliaeetus leucocephalus) in the conterminous United States. We developed an integrated population model (IPM) that incorporates multiple sources of information and then use the model output as the scientific basis for components of the PTL framework. Our IPM is structured to identify key parameters needed for the PTL and to quantify uncertainties in those parameters at the scale at which the United States Fish and Wildlife Service manages take. Our IPM indicated that mean survival of birds >1 year old was high and precise (0.91, 95% CI = 0.90-0.92), whereas mean survival of first-year eagles was lower and more variable (0.69, 95% CI = 0.62-0.78). We assumed that density dependence influenced recruitment by affecting the probability of breeding, which was highly imprecise and estimated to have declined from approximately 0.988 (95% CI = 0.985-0.993) to 0.66 (95% CI = 0.34-0.99) between 1994 and 2018. We sampled values from the posterior distributions of the IPM for use in the PTL and estimated that allowable take (e.g., permitted take for energy development, incidental collisions with human made structures, or removal of nests for development) ranged from approximately 12,000 to 20,000 individual eagles depending on risk tolerance and form of density dependence at the scale of the conterminous United States excluding the Southwest. Model-based thresholds for allowable take can be inaccurate if the assumptions of the underlying framework are not met, if the influence of permitted take is under-estimated, or if undetected population declines occur from other sources. Continued monitoring and use of the IPM and PTL frameworks to identify key uncertainties in bald eagle population dynamics and management of allowable take can mitigate this potential bias, especially where improved information could reduce the risk of permitting non-sustainable take.

Wildlife biologists monitor the status and trends of American woodcock Scolopax minor populations in the eastern and central United States and Canada via a singing-ground survey, conducted just after sunset along roadsides in spring. Annual analyses of the survey produce estimates of trend and annual indexes of abundance for 25 states and provinces, management regions, and survey-wide. In recent years, researchers have used a log-linear hierarchical model that defines year effects as random effects in the context of a slope parameter (the S model) to model population change. Recently, researchers have proposed alternative models suitable for analysis of singing-ground survey data. Analysis of a similar roadside survey, the North American Breeding Bird Survey, has indicated that alternative models are preferable for almost all species analyzed in the Breeding Bird Survey. Here, we use leave-one-out cross-validation to compare model fit for the present singing-ground survey model to fits of three alternative models, including a model that describes population change as the difference in expected counts between successive years (the D model) and two models that include t-distributed extra-Poisson overdispersion effects (H models) as opposed to normally distributed extra-Poisson overdispersion. Leave-one-out cross-validation results indicate that the Bayesian predictive information criterion favored the D model, but a pairwise t-test indicated that the D model was not significantly better-fitting to singing-ground survey data than the S model. The H models are not preferable to the alternatives with normally distributed overdispersion. All models provided generally similar estimates of trend and annual indexes suggesting that, within this model set, choice of model will not lead to alternative conclusions regarding population change. However, as in Breeding Bird Survey analyses, we note a tendency for S model results to provide slightly more extreme estimates of trend relative to D models. We recommend use of the D model for future singing-ground survey analyses.

The North American Breeding Bird Survey (BBS) provides data that can be used in complex, multiscale analyses of population change, while controlling for scale-specific nuisance factors. Many alternative models can be fit to the data, but most model selection procedures are not appropriate for hierarchical models. Leave-one-out cross-validation (LOOCV), in which relative model fit is assessed by omitting an observation and assessing the prediction of a model fit using the remainder of the data, provides a reasonable approach for assessing models, but is time consuming and not feasible to apply for all observations in large data sets. We report the first large-scale formal model selection for BBS data, applying LOOCV to stratified random samples of observations from BBS data. Our results are for 548 species of North American birds, comparing the fit of four alternative models that differ in year effect structures and in descriptions of extra-Poisson overdispersion. We use a hierarchical model among species to evaluate posterior probabilities that models are best for individual species. Models in which differences in year effects are conditionally independent (D models) were generally favored over models in which year effects are modeled by a slope parameter and a random year effect (S models), and models in which extra-Poisson overdispersion effects are independent andt-distributed (H models) tended to be favored over models where overdispersion was independent and normally distributed. Our conclusions lead us to recommend a change from the conventional S model to D and H models for the vast majority of species (544/548). Comparison of estimated population trends based on the favored model relative to the S model currently used for BBS summaries indicates no consistent differences in estimated trends. Of the 18 species that showed large differences in estimated trends between models, estimated trends from the default S model were more extreme, reflecting the influence of the slope parameter in that model for species that are undergoing large population changes. WAIC, a computationally simpler alternative to LOOCV, does not appear to be a reliable alternative to LOOCV.

Setting land aside has long been a primary approach for protecting biodiversity; however, the efficacy of this approach has been questioned. We examined whether protecting lands positively influences bird species in the U.S., and thus overall biodiversity. We used the North American Breeding Bird Survey and Protected Areas Database of the U.S. to assess effects of protected and multiple-use lands on the prevalence and long-term population trends of imperiled and non-imperiled bird species. We evaluated whether both presence and proportional area of protected and multiple-use lands surrounding survey routes affected prevalence and population trends for imperiled and non-imperiled species. Regarding presence of these lands surrounding these survey routes, our results suggest that imperiled and non-imperiled species are using the combination of protected and multiple-use lands more than undesignated lands. We found no difference between protected and multiple-use lands. Mean population trends were negative for imperiled species in all land categories and did not differ between the land categories. Regarding proportion of protected lands surrounding the survey routes, we found that neither the prevalence nor population trends of imperiled or non-imperiled species was positively associated with any land category. We conclude that, although many species (in both groups) tend to be using these protected and multiple-use lands more frequently than undesignated lands, this protection does not appear to improve population trends. Our results may be influenced by external pressures (e.g., habitat fragmentation), the size of protected lands, the high mobility of birds that allows them to use a combination of all land categories, and management strategies that result in similar habitat between protected and multiple-use lands, or our approach to detect limited relationships. Overall, our results suggest that the combination of protected and multiple-use lands is insufficient, alone, to prevent declines in avian biodiversity at a national scale.

Species extinctions have defined the global biodiversity crisis, but extinction begins with loss in abundance of individuals that can result in compositional and functional changes of ecosystems. Using multiple and independent monitoring networks, we report population losses across much of the North American avifauna over 48 years, including once-common species and from most biomes. Integration of range-wide population trajectories and size estimates indicates a net loss approaching 3 billion birds, or 29% of 1970 abundance. A continent-wide weather radar network also reveals a similarly steep decline in biomass passage of migrating birds over a recent 10-year period. This loss of bird abundance signals an urgent need to address threats to avert future avifaunal collapse and associated loss of ecosystem integrity, function, and services.

ABSTRACT Analysis of North American Breeding Bird Survey (BBS) data requires controls for factors that influence detectability of birds along survey routes. Identifying factors that influence the counting process and incorporating them into analyses is a primary means of limiting bias in estimates of population change. Twedt (2015) implemented an alternative counting protocol on operational and nonrandom BBS survey routes in the southeastern United States. Observers on selected routes employed a time–distance protocol in which they recorded birds in 1–min intervals and in 2 distance categories. We hypothesized that processing and recording observations using this time–distance protocol could cause observers to count fewer birds relative to observers using the standard protocol. We used a hierarchical log–linear model with a categorical covariate associated with protocol (standard vs. time–distance) to assess whether use of the time–distance protocol had a measurable effect on counting birds along BBS routes. We applied this model to BBS data from portions of 8 states in which the time–distance protocol was implemented and estimated a protocol effect for 167 bird species. We documented a significant overall effect of the time–distance protocol on observers'counts of birds. On average, the effect of the time–distance protocol was a 10% decline in counted birds; 80% of species had lower counts when the time–distance protocol was used on a survey route. However, because the time–distance protocol was only used on a small portion of the operational BBS routes and for a limited time, including the covariate for the time–distance protocol data had insignificant effects on analysis of population change. Although the covariate controlled for the effects of the time–distance protocol in BBS data, the results emphasize the importance of standardization as well as a need to track and, if necessary, control in analyses for changes in counting procedures along BBS routes.

N-mixture models provide an appealing alternative to mark-recapture models, in that they allow for estimation of detection probability and population size from count data, without requiring that individual animals be identified. There is, however, a cost to using the N-mixture models: inference is very sensitive to the model's assumptions. We consider the effects of three violations of assumptions that might reasonably be expected in practice: double counting, unmodeled variation in population size over time, and unmodeled variation in detection probability over time. These three examples show that small violations of assumptions can lead to large biases in estimation. The violations of assumptions we consider are not only small qualitatively, but are also small in the sense that they are unlikely to be detected using goodness-of-fit tests. In cases where reliable estimates of population size are needed, we encourage investigators to allocate resources to acquiring additional data, such as recaptures of marked individuals, for estimation of detection probabilities.

ABSTRACTBlack vultures (Coragyps atratus) have been increasing in density and expanding their range in the eastern United States since at least the 1960s. In many areas, their densities have increased to the level where they are causing damage to property and livestock and the number of requests for allowable take permits has increased throughout these areas. The United States Fish and Wildlife Service (USFWS) requires updated information to help inform the number of take permits that could reduce conflicts while meeting obligations under the Migratory Bird Treaty Act. We expanded analyses used to estimate allowable take in Virginia to cover the range of black vultures in the eastern United States. We used the prescribed take level approach, which integrates demographic rates, population size estimates, and management objectives into an estimate of allowable take. We provide estimates of allowable take at 4 different scales: individual states, Bird Conservation Regions, USFWS administrative regions, and flyways. Our updated population time series provides evidence of rapidly increasing black vulture populations in many regions of the eastern United States, with an overall population estimate of approximately 4.26 million in 2015 in the Atlantic and Mississippi Flyways. Estimated allowable take ranged from a few hundred individuals per year in states at the northern end of the species range to approximately 287,000/year over the entire eastern United States. The USFWS has no legal mandate regarding the spatial scale at which take should be managed and we found little biological evidence of subpopulation structure for black vultures in the eastern United States. We suggest that allowable take for the species be implemented at a scale that meets stakeholder objectives (e.g., reducing conflict, and ensuring that black vultures are not extirpated from local areas) and is efficient for administrative and monitoring purposes. Published 2018. This article is a U.S. Government work and is in the public domain in the USA.

Shifts in species distributions are major fingerprint of climate change. Examining changes in species abundance structures at a continental scale enables robust evaluation of climate change influences, but few studies have conducted these evaluations due to limited data and methodological constraints. In this study, we estimate temporal changes in abundance from North American Breeding Bird Survey data at the scale of physiographic strata to examine the relative influence of different components of climatic factors and evaluate the hypothesis that shifting species distributions are multidirectional in resident bird species in North America. We quantify the direction and velocity of the abundance shifts of 57 permanent resident birds over 44 years using a centroid analysis. For species with significant abundance shifts in the centroid analysis, we conduct a more intensive correlative analysis to identify climate components most strongly associated with composite change of abundance within strata. Our analysis focus on two contrasts: the relative importance of climate extremes vs. averages, and of temperature vs. precipitation in strength of association with abundance change. Our study shows that 36 species had significant abundance shifts over the study period. The average velocity of the centroid is 5.89 km·yr −1 . The shifted distance on average covers 259 km, 9% of range extent. Our results strongly suggest that the climate change fingerprint in studied avian distributions is multidirectional. Among 6 directions with significant abundance shifts, the northwestward shift was observed in the largest number of species ( n = 13). The temperature/average climate model consistently has greater predictive ability than the precipitation/extreme climate model in explaining strata‐level abundance change. Our study shows heterogeneous avian responses to recent environmental changes. It highlights needs for more species‐specific approaches to examine contributing factors to recent distributional changes and for comprehensive conservation planning for climate change adaptation.

ABSTRACT The U.S. Fish and Wildlife Service (USFWS) uses data from the North American Breeding Bird Survey (BBS) to assist in monitoring and management of some migratory birds. However, BBS analyses provide indices of population change rather than estimates of population size, precluding their use in developing abundance-based objectives and limiting applicability to harvest management. Wood Ducks (Aix sponsa) are important harvested birds in the Atlantic Flyway (AF) that are difficult to detect during aerial surveys because they prefer forested habitat. We integrated Wood Duck count data from a ground-plot survey in the northeastern U.S. with AF-wide BBS, banding, parts collection, and harvest data to derive estimates of population size for the AF. Overlapping results between the smaller-scale intensive ground-plot survey and the BBS in the northeastern U.S. provided a means for scaling BBS indices to the breeding population size estimates. We applied these scaling factors to BBS results for portions of the AF lacking intensive surveys. Banding data provided estimates of annual survival and harvest rates; the latter, when combined with parts-collection data, provided estimates of recruitment. We used the harvest data to estimate fall population size. Our estimates of breeding population size and variability from the integrated population model (N̄ = 0.99 million, SD = 0.04) were similar to estimates of breeding population size based solely on data from the AF ground-plot surveys and the BBS (N̄ = 1.01 million, SD = 0.04) from 1998 to 2015. Integrating BBS data with other data provided reliable population size estimates for Wood Ducks at a scale useful for harvest and habitat management in the AF, and allowed us to derive estimates of important demographic parameters (e.g., seasonal survival rates, sex ratio) that were not directly informed by data.

N-mixture models describe count data replicated in time and across sites in terms of abundance N and detectability p. They are popular because they allow inference about N while controlling for factors that influence p without the need for marking animals. Using a capture-recapture perspective, we show that the loss of information that results from not marking animals is critical, making reliable statistical modeling of N and p problematic using just count data. One cannot reliably fit a model in which the detection probabilities are distinct among repeat visits as this model is overspecified. This makes uncontrolled variation in p problematic. By counter example, we show that even if p is constant after adjusting for covariate effects (the constant p assumption) scientifically plausible alternative models in which N (or its expectation) is non-identifiable or does not even exist as a parameter, lead to data that are practically indistinguishable from data generated under an N-mixture model. This is particularly the case for sparse data as is commonly seen in applications. We conclude that under the constant p assumption reliable inference is only possible for relative abundance in the absence of questionable and/or untestable assumptions or with better quality data than seen in typical applications. Relative abundance models for counts can be readily fitted using Poisson regression in standard software such as R and are sufficiently flexible to allow controlling for p through the use covariates while simultaneously modeling variation in relative abundance. If users require estimates of absolute abundance, they should collect auxiliary data that help with estimation of p.

The North American Breeding Bird Survey (BBS) contains data for >700 bird species, but analyses often focus on a core group of similar to 420 species. We analyzed data for 122 species of North American birds for which data exist in the North American Breeding Bird Survey (BBS) database but are not routinely analyzed on the BBS Summary and Analysis Website. Many of these species occur in the northern part of the continent, on routes that fall outside the core survey area presently analyzed in the United States and southern Canada. Other species not historically analyzed occur in the core survey area with very limited data but have large portions of their ranges in Mexico and south. A third group of species not historically analyzed included species thought to be poorly surveyed by the BBS, such as rare, coastal, or nocturnal species. For 56 species found primarily in regions north of the core survey area, we expanded the scope of the analysis, using data from 1993 to 2014 during which >= 3 survey routes had been sampled in 6 northern strata (Bird Conservation regions in Alaska, Yukon, and Newfoundland and Labrador) and fitting log-linear hierarchical models for an augmented BBS survey area that included both the new northern strata and the core survey area. We also applied this model to 168 species historically analyzed in the BBS that had data from these additional northern strata. For both groups of species we calculated survey-wide trends for the both core and augmented survey areas from 1993 to 2014; for species that did not occur in the newly defined strata, we computed trends from 1966 to 2014. We evaluated trend estimates in terms of established credibility criteria for BBS results, screening for imprecise trends, small samples, and low relative abundance. Inclusion of data from the northern strata permitted estimation of trend for 56 species not historically analyzed, but only 4 of these were reasonably monitored and an additional 13 were questionably monitored; 39 of these species were likely poorly monitored because of small numbers of samples or very imprecisely estimated trends. Only 4 of 66 "new'' species found in the core survey area were reasonably monitored by the BBS; 20 were questionably monitored; and 42 were likely poorly monitored by the BBS because of inefficiency in precision, abundance, or sample size. The hierarchical analyses we present provide a means for reasonable inclusion of the additional species and strata in a common analysis with data from the core area, a critical step in the evolution of the BBS as a continent-scale survey. We recommend that results be presented both 1) from 1993 to the present using the expanded survey area, and 2) from 1966 to the present for the core survey area. Although most of the "new'' species we analyzed were poorly monitored by the BBS during 1993-2014, continued expansion of the BBS will improve the quality of information in future analyses for these species and for the many other species presently monitored by the BBS.

The first day of spring 2017 marked the passing, at the age of 98, of Chandler S. Robbins, one of the great American ornithologists.Chan, as he was called, was a man of boundless energy, unflagging positivity, and exceptional intellectual capabilities.He spent his long career entirely with the U.S. government, which provided him with a forum that perfectly complemented his skills.As a field ornithologist arriving

We used North American Breeding Bird Survey data (1966-91) to estimate distribution, relative abundance, and population trends of quails. Population trends in grassland/shrub birds sympatric with northern bobwhite (Colinus virginianus) were also examined. Northern bobwhite and scaled quail (Callipepl,a squamata) populations have declined since 1966. Rates of decline for these quails have increased during the past decade. California quail (C. califomica), Gambel's quail (C. gambelii), and mountain quail (Oreortyx pictus) populations have been stable over the long-term (1966-91). However, the short-term (1982-91) trend for California quail is positive, whereas Gambel's quail appear to be declining. Patterns in trends indicate similar factors may be negatively affecting breeding populations of grassland/shrub birds throughout the bobwhite's range. We discuss plausible hypotheses to explain population trends and recommend future action.

Management interest in North American birds has increasingly focused on species that breed in Alaska, USA, and Canada, where habitats are changing rapidly in response to climatic and anthropogenic factors. We used a series of hierarchical models to estimate rates of population change in 2 forested Bird Conservation Regions (BCRs) in Alaska based on data from the roadside North American Breeding Bird Survey (BBS) and the Alaska Landbird Monitoring Survey, which samples off-road areas on public resource lands. We estimated long-term (1993-2015) population trends for 84 bird species from the BBS and short-term (2003-2015) trends for 31 species from both surveys. Among the 84 species with long-term estimates, 11 had positive trends and 17 had negative trends in 1 or both BCRs; negative trends were primarily found among aerial insectivores and wetland-associated species, confirming range-wide negative continental trends for many of these birds. Three species with negative trends in the contiguous United States and southern Canada had positive trends in Alaska, suggesting different population dynamics at the northern edges of their ranges. Regional population trends within Alaska differed for several species, particularly those represented by different subspecies in the 2 BCRs, which are separated by rugged, glaciated mountain ranges. Analysis of the roadside and off-road data in a joint hierarchical model with shared parameters resulted in improved precision of trend estimates and suggested a roadside-related difference in underlying population trends for several species, particularly within the Northwestern Interior Forest BCR. The combined analysis highlights the importance of considering population structure, physiographic barriers, and spatial heterogeneity in habitat change when assessing patterns of population change across a landscape as broad as Alaska. Combined analysis of roadside and off-road survey data in a hierarchical framework may be particularly useful for evaluating patterns of population change in relatively undeveloped regions with sparse roadside BBS coverage.

Integrated population models (IPMs) provide a unified framework for simultaneously analyzing data sets of different types to estimate vital rates, population size, and dynamics; assess contributions of demographic parameters to population changes; and assess population viability. Strengths of an IPM include the ability to estimate latent parameters and improve the precision of parameter estimates. We present a hierarchical IPM that combines two broad-scale avian monitoring data sets: count data from the North American Breeding Bird Survey (BBS) and capture-recapture data from the Monitoring Avian Productivity and Survivorship (MAPS) program. These data sets are characterized by large numbers of sample sites and observers, factors capable of inducing error in the sampling and observation processes. The IPM integrates the data sets by modeling the population abundance as a first-order autoregressive function of the previous year's population abundance and vital rates. BBS counts were modeled as a log-linear function of the annual index of population abundance, observation effects (observer identity and first survey year), and overdispersion. Vital rates modeled included adult apparent survival, estimated from a transient Cormack-Jolly-Seber model using MAPS data, and recruitment (surviving hatched birds from the previous season + dispersing adults) estimated as a latent parameter. An assessment of the IPM demonstrated it could recover true parameter values from 200 simulated data sets. The IPM was applied to data sets (1992-2008) of two bird species, Gray Catbird (Dumetella carolinensis) and Wood Thrush (Hylocichla mustelina) in the New England/Mid-Atlantic coastal Bird Conservation Region of the United States. The Gray Catbird population was relatively stable (trend +0.4% per yr), while the Wood Thrush population nearly halved (trend -4.5% per yr) over the 17-yr study period. IPM estimates of population growth rates, adult survival, and detection and residency probabilities were similar and as precise as estimates from the stand-alone BBS and CJS models. A benefit of using the IPM was its ability to estimate the latent recruitment parameter. Annual growth rates for both species correlated more with recruitment than survival, and the relationship for Wood Thrush was stronger than for Gray Catbird. The IPM's unified modeling framework facilitates integration of these important data sets.