Breeding Bird Survey statistical methods

The North American Breeding Bird Survey (BBS) is conducted annually in the United States, Canada, and, as of 2008, northern Mexico. The survey is designed to monitor trends in relative abundance of North American breeding birds at continental, national, and regional scales. Environment and Climate Change Canada's Canadian Wildlife Service and the United States Geological Survey's (USGS) Patuxent Wildlife Research Center jointly coordinate the BBS program.

Data collection

BBS data are collected along randomly established roadside routes by skilled volunteer observers. Routes consist of 50 stops spaced 0.8 km apart along a 39.2 km length of all-weather secondary roads. The starting point and direction of routes are selected randomly within a given degree block (degree blocks are geographic areas 1 degree longitude by 1 degree latitude in size and are the basic sample design units for the BBS) in order to sample a range of habitats representative of the region. In Canada, the acceptable dates for running a BBS route are between 28 May and 7 July although most observations are made during the first three weeks of June. Observers are encouraged to run their routes for as many consecutive years as possible in order to reduce the effects of observer variability on trend estimates. Because the BBS is designed to monitor long-term changes in bird populations, changes in habitat along routes are expected. However, routes may be altered or abandoned if the road system is altered, increasing ambient noise interferes with an observer's ability to detect birds, or if traffic conditions are unsafe for the observer.

At each of the 50 stops, observers conduct a three-minute point count observation, during which they record the total number and species of all birds heard, and all birds seen within 0.4 km of their location. Start and finish times for the routes are recorded. Sky and wind conditions are recorded at the beginning of each 10-stop block and at the end of the route. The number of vehicles passing and other excess noise are recorded during each stop. Data are entered into the online database by observers, or are submitted to one of the BBS national offices for entry. All BBS data are stored in the North American BBS database housed at the USGS Patuxent Wildlife Research Center and are available on the USGS BBS website. Original data sheets for Canada are stored by Environment and Climate Change Canada's Canadian Wildlife Service. For a comprehensive description of the BBS field methods, see Robbins et al. (1986).

Data screening

Several factors, in addition to changes in bird populations, contribute to variation within BBS data. These include changes in weather, date, starting and finishing times, as well as differences among observers. To help reduce the extraneous variability in the data, the data are screened to remove surveys run under unacceptable conditions (e.g., rain and/or high winds, outside of the range of acceptable dates, or outside of the range of acceptable times of day). All criteria are available on the USGS BBS website.

Geographic analysis strata

The data are analyzed within geographic strata defined by the intersection of Bird Conservation Regions (BCRs) and the Provinces/Territories. For example, British Columbia is divided into four strata by BCRs 4, 5, 9, and 10. For a given species, only strata with sufficient data are included in the analysis. To be included, a stratum must contain at least 3 routes on which the species was observed, at least one of those routes was surveyed in 5 or more years, and at least one route on which the species has been observed in 3 or more years. Trends and estimates of relative abundance are modeled separately within each of these analytical strata to account for differences in a species' population status and trends among different regions, and to allow the stratum-level estimates to be combined into composite regional estimates based on BCRs or political units.

The geographic nature of these strata means that users of these results should carefully consider the spatial scope of the national and regional estimates of trends and annual indices from the BBS data. The estimates apply to the specific strata in which the data were collected and these strata may not completely cover a given region of interest. For example, we report "national" trends; however, the BBS does not cover Canada in its entirety, so some portion of a species' population may not be sampled by the BBS. To this end, every trend and annual index estimate is presented with an associated map showing the regions to which the estimates apply. The maps provide an at-a-glance assessment of the spatial scope of each trend and annual index estimate.

Additionally, these maps now include information on the trends in each portion of the species' range. Regions are coloured based on the direction (increases in blues or decreases in reds) and the intensity of the colour indicates the magnitude of the trend (darker colours for steeper rates of change).

Analytical methods

We analyze the BBS using a hierarchical Bayesian model that accounts for the effects of variation among observers and routes, first-year observer effects, variations in abundance and the non-linear pattern of population change among strata, and annual population fluctuations. The observed counts are assumed to have an over-dispersed Poisson distribution. The model is fit using standard approaches to Bayesian modeling: Hamilton Monte Carlo (HMC) methods with Stan software and the R-package "bbsBayes" (Edwards and Smith 2020). The model used here estimates the pattern of population change through time using a non-linear smoothing component (hierarchical generalized additive model) with additional yearly fluctuations. The model is described in a scientific paper, available as an open access article, published in the scientific journal Ornithological Applications (Smith and Edwards 2020).

The hierarchical structure of the non-linear model makes efficient use of the data and also means estimates are much less sensitive to random sampling error (e.g., variation in the number of routes run in a given year and region). The population trajectories in each geographic stratum are estimated in a way that shares some information across the species' range. Therefore, estimates of population trajectories in strata with relatively few data are more similar to the species' range-wide average trajectory. As a result, estimated rates of change in regions with relatively sparse data are generally more precise and less extreme (relative to the species' range-wide average slope) than estimates from models that do not share information across the species' range.

The relative influence of local data and data from other parts of the species' range on the estimates of the population trajectory depends on two factors: 1) the quality of the data in the region and 2) the relative shape of the local population trajectory, compared to the range-wide mean trajectory (i.e., range-wide means have greater influence when local estimates are very imprecise and/or are very different from the range-wide estimate). We believe this statistical assumption is reasonable because local populations of a given species are linked across the species' range, due to mixing of individuals through dispersal and migration. Also, many users of BBS trend estimates are likely comfortable with the idea of considering trend estimates from other regions to inform their understanding of a species' status in areas where local trend estimates are particularly poor (i.e., low precision) or are of an unexpected magnitude (e.g., increasing steeply when other regional trends suggest the species is stable or in decline). One important consequence is that BBS data from outside of Canada are now incorporated into the Canadian trend analyses. This is a justifiable feature of this new approach because: 1) populations of bird species are linked across national borders; 2) it improves the precision of regional and national trend estimates; 3) for regions with high-quality data, it has little or no influence; 4) it reduces the likelihood of estimating anomalously extreme trends of very low quality (e.g., steeply increasing trends with exceptionally low precision); and 5) it allows us to estimate trends for more species, particularly those with limited distributions within Canada.

Reliability measures

To quantify the reliability of the BBS trend estimates in describing overall population trends, three measures of reliability are calculated: precision, coverage and local data weight. In addition to these three credibility measures, one may also want to consider the following measures reported in the final columns of the trend data tables: the number of routes used to generate the trends, the estimates of the probabilities that a given population has decreased or increased, or the probabilities that the magnitude of the population change falls within one of the six population change categories. These probabilities provide a concise and intuitive assessment of the uncertainty around some commonly desired interpretations. For example, a relatively strong positive trend that is imprecisely estimated may in fact have a greater probability of being a decline (i.e., the trend was estimated in the wrong direction), than a stable or slightly positive, but precisely estimated trend.

Precision: This measure reflects the width of the 95% credible interval (CI). Smaller values indicate a more precise trend estimate. The 95% CI defines a range of trend values that have a 95% probability of including the true population trend value, given the data and the model. For example, the following statement is a valid interpretation of a hypothetical 95% CI with lower limit = 1.2 and upper limit = 2.3: given the BBS data and the accuracy of the model, there is a 95% probability that the average annual increase in the population is somewhere between 1.2%/year and 2.3%/year.

Geographic coverage: This measure indicates what proportion of the regional population had some probability of influencing the estimated trend. It is calculated as the estimated proportion of the regional population inside degree-blocks with BBS data. The reported value is a geometric average annual proportion covered during the time period of the trend, where degree blocks are considered covered for all years between the first and last years that contributing routes were run in the degree-block. Estimates of regional and degree-block specific population proportions were derived using maps of a species' relative density across its Canadian range. These maps were created by overlaying existing species range maps with estimates of relative population density within the species' range, which were derived using the average number of individuals on BBS routes run from 1997 through 2007, and comparable counts from other monitoring programs.

Local data weight: Indicates the proportion of the region and time-series for which local data contributed to estimating the local population trajectory. That is, the proportion of the region with some data in each year, averaged across all years included in the trend. A value of 1.0 indicates that there are data from every year and stratum that makes up the region. Values less than 1.0 indicate that some portion of the region (i.e., one or more of the strata included in the region) does not include data from the full range of the time series, and therefore relies on the data from other regions of the species range to estimate those years in the local trajectory.

Precision, coverage, and local data weight are combined into an overall reliability score for each trend. These overall reliability scores provide an overview that alerts the user to potentially unreliable trends. However, specific measures are also available to make an informed decision about the uncertainty associated with particular conclusions based on the trends. For example, a given trend may have a low reliability score because it is relatively imprecise (i.e., large credible interval), but, if for example the magnitude of the trend estimate is sufficiently large, this imprecision may not affect an overall conclusion that the species' population has increased.

Literature cited

Back to (Breeding Bird Survey results)
Date modified: