Assessing Allowable Take of Migratory Birds

ABSTRACT Legal removal of migratory birds from the wild occurs for several reasons, including subsistence, sport harvest, damage control, and the pet trade. We argue that harvest theory provides the basis for assessing the impact of authorized take, advance a simplified rendering of harvest theory known as potential biological removal as a useful starting point for assessing take, and demonstrate this approach with a case study of depredation control of black vultures (Coragyps atratus) in Virginia, USA. Based on data from the North American Breeding Bird Survey and other sources, we estimated that the black vulture population in Virginia was 91,190 (95% credible interval = 44,520-212,100) in 2006. Using a simple population model and available estimates of life-history parameters, we estimated the intrinsic rate of growth (r_max) to be in the range 7–14%, with 10.6% a plausible point estimate. For a take program to seek an equilibrium population size on the conservative side of the yield curve, the rate of take needs to be less than that which achieves a maximum sustained yield (0.5 × r_max). Based on the point estimate for r_max and using the lower 60% credible interval for population size to account for uncertainty, these conditions would be met if the take of black vultures in Virginia in 2006 was <3,533 birds. Based on regular monitoring data, allowable harvest should be adjusted annually to reflect changes in population size. To initiate discussion about how this assessment framework could be related to the laws and regulations that govern authorization of such take, we suggest that the Migratory Bird Treaty Act requires only that take of native migratory birds be sustainable in the long-term, that is, sustained harvest rate should be <r_max. Further, the ratio of desired harvest rate to 0.5 X r_max may be a useful metric for ascertaining the applicability of specific requirements of the National Environmental Protection Act.     

Spatial risk modeling of cattle depredation by black vultures in the midwestern United States

Negative economic impacts resulting from wildlife disrupting livestock operations through depredation of stock are a cause of human-wildlife conflict. Management of such conflict requires identifying environmental and non-environmental factors specific to a wildlife species' biology and ecology that influence the potential for livestock depredation to occur. Identification of such factors can improve understanding of the conditions placing livestock at risk. Black vultures (Coragyps atratus) have expanded their historical range northward into the midwestern United States. Concomitantly, an increase in concern among agricultural producers regarding potential black vulture attacks on livestock has occurred. We estimated area with greater or lesser potential for depredation of domestic cattle by black vultures across a 6-state region in the midwestern United States using an ensemble of small models (ESM). Specifically, we identified landscape-scale spatial factors, at a zip code resolution, associated with reported black vulture depredation on cattle in midwestern landscapes to predict future potential livestock depredation. We hypothesized that livestock depredation would be greatest in areas with intensive beef cattle production close to preferred black vulture habitat (e.g., areas with fewer old fields and early successional vegetation paired with more direct edge between older forest and agricultural lands). We predicted that the density of cattle within the county, habitat structure, and proximity to anthropogenic landscape features would be the strongest predictors of black vulture livestock-depredation risk. Our ESM estimated the relative risk of black vulture-cattle depredation to be between 0.154–0.631 across our entire study area. Consistent with our hypothesis, areas of greatest predicted risk of depredation correspond with locations that are favorable to vulture life-history requirements and increased potential to encounter livestock. Our results allow wildlife managers the ability to predict where black vulture depredation of cattle is more likely to occur in the future. It is in these areas where extension and outreach efforts aimed at mitigating this conflict should be focused. Researchers and wildlife managers interested in developing or employing tools aimed at mitigating livestock-vulture conflicts can also leverage our results to select areas where depredation is most likely to occur.     

Estimating allowable take for an increasing bald eagle population in the United States

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.     

Prescribed take levels of downy and hairy woodpeckers in the eastern United States

Hairy (Dryobates villosus) and downy (Dryobates pubescens) woodpeckers occur in high densities in residential areas of the eastern United States. In many areas of their range, they cause damage to wooden structures through foraging, excavation of nesting cavities, and drumming behaviors, causing requests for allowable take permits. Both species hold year-round territories, which could make them vulnerable to local extirpation with excess take. To meet the requirements of the Migratory Bird Treaty act, the United States Fish and Wildlife Service (USFWS) requested scientifically informed evaluation of take to minimize population effects as part of its approach to reduce human–wildlife conflict. We used a prescribed take approach, which uses data from population, demographic, and management parameter estimates to determine the allowable take from Louisiana to Minnesota and all states east. Furthermore, we used 2 different methods of estimating growth rates to control for demographic uncertainties. The resulting estimates provide take at the state and USFWS regional scales to improve stakeholder choices when setting allowable take. Current authorized take (2016–2018) is below the take that could be sustained by current populations, and current rates of take are not likely to cause population-level effects. These results were largely consistent across methodologies for calculating the rate of growth for both species. Take still needs to be managed to prevent local extirpation of these resident species. Allowable take estimates should be periodically updated to reflect changing management and population needs for both species.     

Demographics of Black Vultures in North Carolina

Abstract: Understanding the contributions of vital rates to species population growth is critical to developing new management protocols. We constructed a model population for black vultures (Coragyps atratus) in North Carolina, USA, based on demographic data from a 14-year study. The model population was similar in stage structure to the reference population, and adult survival was the primary contributor to the annual rate of increase (10.6%). We suggest that the North Carolina black vulture population is experiencing high rates of survival and fertility, potentially breeding at an age younger than previously assumed, and growing rapidly.     

Comparative Analysis of Mourning Dove Population Change in North America

ABSTRACT Mourning doves (Zenaida macroura) are surveyed in North America with a Call-Count Survey (CCS) and the North American Breeding Bird Survey (BBS). Analyses in recent years have identified inconsistencies in results between surveys, and a need exists to analyze the surveys using modern methods and examine possible causes of differences in survey results. Call-Count Survey observers collect separate information on number of doves heard and number of doves seen during counting, whereas BBS observers record one index containing all doves observed. We used hierarchical log-linear models to estimate trend and annual indices of abundance for 1966–2007 from BBS data, CCS-heard data, and CCS-seen data. Trend estimates from analyses provided inconsistent results for several states and for eastern and central dove-management units. We examined differential effects of change in land use and noise-related disturbance on the CCS indices. Changes in noise-related disturbance along CCS routes had a larger influence on the heard index than on the seen index, but association analyses among states of changes in temperature and of amounts of developed land suggest that CCS indices are differentially influenced by changes in these environmental features. Our hierarchical model should be used to estimate population change from dove surveys, because it provides an efficient framework for estimating population trends from dove indices while controlling for environmental features that differentially influence the indices.     

Composite analysis of black duck breeding population surveys in eastern North America

Waterfowl are monitored in eastern Canada and the northeastern United States with 2 surveys: a transect survey from fixed-wing aircraft and a plot survey conducted from helicopters. The surveys vary in extent, but overlap exists in a core area of 9 strata covering portions of all provinces from Ontario east to Newfoundland. We estimated population change for American black ducks (Anas rubripes) and mallards (Anas platyrhynchos) from these surveys using a log-linear hierarchical model that accommodates differences in sample design and visibility associated with these survey methods. Using a combined analysis of the surveys based on total indicated birds, we estimate the American black duck population to be 901,700 (95% CI: 715,200–1,274,000) in 2011, with 526,900 (95% CI: 357,500–852,300) mallards in the surveyed area. Precision of estimates varies widely by species and region, with transect surveys providing less precise results than plot surveys for black ducks in areas of overlap. The combined survey analysis for black ducks in the eastern survey region produced estimates with an average yearly coefficient of variation (CV) of 12.1% for the entire area and an average CV of 6.9% in the plot survey area. Mallards, which had a more limited distribution in the region, had an average yearly CV of 22.1% over the entire region, and an average CV of 27.7% in the plot survey area. Hierarchical models provide a rich framework for analyzing and combining results from complex survey designs, providing useful spatial and temporal information on population size and change in these economically important species. © 2012 The Wildlife Society.     

Evidence for the buffer effect operating in multiple species at a national scale

A long-standing aim of ecologists is to understand the processes involved in regulating populations. One such mechanism is the buffer effect, where lower quality habitats are increasingly used as a species reaches higher population densities, with a resultant average reduction in fecundity and survival limiting population growth. Although the buffer effect has been demonstrated in populations of a number of species, a test of its importance in influencing population growth rates of multiple species across large spatial scales is lacking. Here, we use habitat-specific population trends for 85 bird species from long-term national monitoring data (the UK Breeding Bird Survey) to examine its generality. We find that both patterns of population change and changes in habitat preference are consistent with the predictions of the buffer effect, providing support for its widespread operation.     

Combining Breeding Bird Survey and Christmas Bird Count Data to Evaluate Seasonal Components of Population Change in Northern Bobwhite

Abstract Annual surveys of wildlife populations provide information about annual rates of change in populations but provide no information about when such changes occur. However, by combining data from 2 annual surveys, conducted in different parts of the year, seasonal components of population change can be estimated. We describe a hierarchical model for simultaneous analysis of 2 continent-scale monitoring programs. The Christmas Bird Count is an early winter survey, whereas the North American Breeding Bird Survey is conducted in June. Combining information from these surveys permits estimation of seasonal population variance components and improves estimation of long-term population trends. The composite analysis also controls for survey-specific sampling effects. We applied the model to estimation of population change in northern bobwhites (Colinus virginianus). Over the interval 1969–2004, bobwhite populations declined, with trend estimate of −3.56% per year (95% CI = [−3.80%, −3.32%]) in the surveyed portion of their range. Our analysis of seasonal population variance components indicated that northern bobwhite populations changed more in the winter and spring portion of the year than in the summer and fall portion of the year. (JOURNAL OF WILDLIFE MANAGEMENT 72(1):44–51; 2008)     

Diversity of birds in eastern North America shifts north with global warming

The distribution of diversity along latitudinal and elevation gradients, and the coupling of this phenomenon with climate, is a pattern long recognized in ecology. Hypothesizing that climate change may have altered this pattern over time, we investigated whether the aggregate of reported northward shifts of bird ranges in North America is now detectable in community‐level indices such as richness and diversity. Here, we report that bird diversity in North America increased and shifted northward between 1966 and 2010. This change in the relationship of diversity to the latitudinal gradient is primarily influenced by range expansions of species that winter in the eastern United States as opposed to species which migrate to this area from wintering grounds in the tropics. This increase in diversity and its northward expansion is best explained by an increase in regional prebreeding season temperature over the past 44 years.     

Bird diversity: a predictable function of satellite-derived estimates of seasonal variation in canopy light absorbance across the United States

Aim  To investigate the relationships between bird species richness derived from the North American Breeding Bird Survey and estimates of the average, minimum, and the seasonal variation in canopy light absorbance (the fraction of absorbed photosynthetically active radiation, fPAR) derived from NASA's Moderate Resolution Imaging Spectroradiometer (MODIS).
Location  Continental USA.
Methods  We describe and apply a 'dynamic habitat index' (DHI), which incorporates three components based on monthly measures of canopy light absorbance through the year. The three components are the annual sum, the minimum, and the seasonal variation in monthly fPAR, acquired at a spatial resolution of 1 km, over a 6-year period (2000–05). The capacity of these three DHI components to predict bird species richness across 84 defined ecoregions was assessed using regression models.
Results  Total bird species richness showed the highest correlation with the composite DHI [ R ² = 0.88, P  < 0.001, standard error of estimate (SE) = 8 species], followed by canopy nesters ( R ² = 0.79, P  < 0.001, SE = 3 species) and grassland species ( R ² = 0.74, P  < 0.001, SE = 1 species). Overall, the seasonal variation in fPAR, compared with the annual average fPAR, and its spatial variation across the landscape, were the components that accounted for most ( R ² = 0.55–0.88) of the observed variation in bird species richness.
Main conclusions  The strong relationship between the DHI and observed avian biodiversity suggests that seasonal and interannual variation in remotely sensed fPAR can provide an effective tool for predicting patterns of avian species richness at regional and broader scales, across the conterminous USA.     

Synergistic effects of climate and land‐use change influence broad‐scale avian population declines

Climate and land‐use changes are expected to be the primary drivers of future global biodiversity loss. Although theory suggests that these factors impact species synergistically, past studies have either focused on only one in isolation or have substituted space for time, which often results in confounding between drivers. Tests of synergistic effects require congruent time series on animal populations, climate change and land‐use change replicated across landscapes that span the gradient of correlations between the drivers of change. Using a unique time series of high‐resolution climate (measured as temperature and precipitation) and land‐use change (measured as forest change) data, we show that these drivers of global change act synergistically to influence forest bird population declines over 29 years in the Pacific Northwest of the United States. Nearly half of the species examined had declined over this time. Populations declined most in response to loss of early seral and mature forest, with responses to loss of early seral forest amplified in landscapes that had warmed over time. In addition, birds declined more in response to loss of mature forest in areas that had dried over time. Climate change did not appear to impact populations in landscapes with limited habitat loss, except when those landscapes were initially warmer than the average landscape. Our results provide some of the first empirical evidence of synergistic effects of climate and land‐use change on animal population dynamics, suggesting accelerated loss of biodiversity in areas under pressure from multiple global change drivers. Furthermore, our findings suggest strong spatial variability in the impacts of climate change and highlight the need for future studies to evaluate multiple drivers simultaneously to avoid potential misattribution of effects.     

Road exposure and the detectability of birds in field surveys

Road ecology, the study of the impacts of roads and their traffic on wildlife, including birds, is a rapidly growing field, with research showing effects on local avian population densities up to several kilometres from a road. However, in most studies, the effects of roads on the detectability of birds by surveyors are not accounted for. This could be a significant source of error in estimates of the impacts of roads on birds and could also affect other studies of bird populations. Using road density, traffic volume and bird count data from across Great Britain, we assess the relationships between roads and detectability of a range of bird species. Of 51 species analysed, the detectability of 36 was significantly associated with road exposure, in most cases inversely. Across the range of road exposure recorded for each species, the mean positive change in detectability was 52% and the mean negative change was 36%, with the strongest negative associations found in smaller-bodied species and those for which aural cues are more important in detection. These associations between road exposure and detectability could be caused by a reduction in surveyors’ abilities to hear birds or by changes in birds’ behaviour, making them harder or easier to detect. We suggest that future studies of the impacts of roads on populations of birds or other taxa, and other studies using survey data from road-exposed areas, should account for the potential impacts of roads on detectability.     

Changes in the number and distribution of Greater Sandhill Cranes in the Eastern Population

Once nearly extirpated, the Eastern Population (EP) of Greater Sandhill Cranes (Grus canadensis tabida) has increased in number and expanded its range in breeding and wintering areas. Data from Christmas Bird Counts (CBCs) and Breeding Bird Surveys (BBSs) were used to delineate changes in the wintering and breeding area distributions during the period from 1966 to 2013. Crane densities were plotted to the centroid of CBC circles or BBS routes, and the Geographic Mean Centers (GMCs) for wintering and breeding populations were calculated. The number of Greater Sandhill Cranes detected during the breeding season has steadily increased since 1966, with just six birds observed in 1966 and 1046 observed in 2013. The GMC of the Sandhill Crane breeding population has remained in Wisconsin during the 47‐yr time frame. The total number of Sandhill Cranes counted in the eastern United States during CBCs grew from 423 in 1965–1966 to 46,194 in 2012–2013, with a peak number of 55,826 in 2011–2012. The GMC of wintering Greater Sandhill Cranes was located in Florida during the periods from 1966 to 1977 and 1978 to 1989, but shifted north‐northwest by nearly 4° of latitude (into Georgia) by 1990–2001. By 2002–2013, the GMC had shifted an additional degree north as well as almost a degree west in longitude. Greater Sandhill Cranes in the EP may continue to winter further north and remain in more northerly areas later in the fall before migrating further south. Factors such as annual weather, long‐term climate change, and changes in land use may influence future population trends and changes in both the breeding and wintering ranges of the EP of Sandhill Cranes.     

Filling in the GAPS: evaluating completeness and coverage of open‐access biodiversity databases in the United States

Primary biodiversity data constitute observations of particular species at given points in time and space. Open‐access electronic databases provide unprecedented access to these data, but their usefulness in characterizing species distributions and patterns in biodiversity depend on how complete species inventories are at a given survey location and how uniformly distributed survey locations are along dimensions of time, space, and environment. Our aim was to compare completeness and coverage among three open‐access databases representing ten taxonomic groups (amphibians, birds, freshwater bivalves, crayfish, freshwater fish, fungi, insects, mammals, plants, and reptiles) in the contiguous United States. We compiled occurrence records from the Global Biodiversity Information Facility (GBIF), the North American Breeding Bird Survey (BBS), and federally administered fish surveys (FFS). We aggregated occurrence records by 0.1° × 0.1° grid cells and computed three completeness metrics to classify each grid cell as well‐surveyed or not. Next, we compared frequency distributions of surveyed grid cells to background environmental conditions in a GIS and performed Kolmogorov–Smirnov tests to quantify coverage through time, along two spatial gradients, and along eight environmental gradients. The three databases contributed >13.6 million reliable occurrence records distributed among >190,000 grid cells. The percent of well‐surveyed grid cells was substantially lower for GBIF (5.2%) than for systematic surveys (BBS and FFS; 82.5%). Still, the large number of GBIF occurrence records produced at least 250 well‐surveyed grid cells for six of nine taxonomic groups. Coverages of systematic surveys were less biased across spatial and environmental dimensions but were more biased in temporal coverage compared to GBIF data. GBIF coverages also varied among taxonomic groups, consistent with commonly recognized geographic, environmental, and institutional sampling biases. This comprehensive assessment of biodiversity data across the contiguous United States provides a prioritization scheme to fill in the gaps by contributing existing occurrence records to the public domain and planning future surveys.     

Representation of Landcover Along Breeding Bird Survey Routes in the Northern Plains

ABSTRACT The North American Breeding Bird Survey (BBS) is used extensively to make inferences about populations of many North American bird species and is increasingly being used for avian conservation planning. How well BBS routes represent the landscape is poorly known, even though accuracy of representation could significantly affect inferences made from BBS data. We used digital landcover data to examine how well landcover within 400-m buffers around BBS routes represented the surrounding landscape (the route neighborhood) for 52 routes in the Prairie Pothole Region of North Dakota and South Dakota. Differences in composition between landcover along BBS routes and the route neighborhood were not statistically significant for upland cover classes. The area of temporary and seasonal wetland basins was accurately represented by BBS routes in our study area, but the area of semipermanent and permanent wetland basins was significantly underrepresented along BBS routes. Number of wetland basins and upland patches was higher along routes. Area of urban, forest, and hay landcover classes was higher along routes, although differences were not statistically significant. Amount of bias in landcover representation was negatively correlated with the proportion of each landcover type in the study area, but bias was not correlated with area of the route neighborhoods. Differences between landcover along BBS routes and the route neighborhood were primarily attributable to increased anthropogenic activity along roads and siting of roads away from relatively large, deep water bodies. Our results suggest that inferences made from BBS data in our study region are likely biased for species that are associated with deeper-water habitats or are strongly influenced by landscape fragmentation. Inferences made from BBS data for species associated with uplands or shallow wetlands are less likely to be biased because of differences in landcover composition.     

Evaluating the ability of regional models to predict local avian abundance

Spatial modeling over broad scales can potentially direct conservation efforts to areas with high species-specific abundances. We examined the performance of regional models for predicting bird abundance at spatial scales typically addressed in conservation planning. Specifically, we used point count data on wood thrush (Hylocichla mustelina) and blue-winged warbler (Vermivora cyanoptera) from 2 time periods (1995–1998 and 2006–2007) to evaluate the ability of regional models derived via Bayesian hierarchical techniques to predict bird abundance. We developed models for each species within Bird Conservation Region (BCR) 23 in the upper midwestern United States at 800-ha, 8,000-ha, and approximately 80,000-ha scales. We obtained count data from the Breeding Bird Survey and land cover data from the National Land Cover Dataset (1992). We evaluated predictions from the best models, as defined by an information-theoretic criterion, using point count data collected within an ecological subregion of BCR 23 at 131 count stations in the 1990s and again in 2006–2007. Competing (Deviance Information Criteria <5) blue-winged warbler models accounted for 67% of the variability and suggested positive associations with forest edge and proportion of forest at the 8,000-ha scale, and negative associations with forest patch area (800 ha) and wetness (800 ha and 80,000 ha). The regional model performed best for blue-winged warbler predicted abundances from point counts conducted in Iowa during 1995–1996 (r_s = 0.57; P = 0.14), the survey period that most closely aligned with the time period of data used for regional model construction. Wood thrush models exhibited positive correlations with point count data for all survey areas and years combined (r_s = 0.58, P ≤ 0.001). In comparison, blue-winged warbler models performed worse as time increased between the point count surveys and vintage of the model building data (r_s = 0.03, P = 0.92 for Iowa and r_s = 0.13, P = 0.51 for all areas, 2006–2007), likely related to the ephemeral nature of their preferred early successional habitat. Species abundance and sensitivity to changing habitat conditions seems to be an important factor in determining the predictive ability of regional models. Hierarchical models can be a useful tool for concentrating efforts at the scale of management units and should be one of many tools used by land managers, but we caution that the utility of such models may decrease over time for species preferring relatively ephemeral habitats if model inputs are not updated accordingly. © 2012 The Wildlife Society.