Noise, avian abundance, and productivity at banding stations across the Continental United States

Noise is an increasingly common component of the natural world, due in large part to human activity. Anthropogenic noise negatively impacts abundance, health, and reproduction in many songbird populations. A few studies have reported altered abundance at larger scales. But whether continental trends are being detected at banding stations, which also offer data on productivity and survivorship, is unknown. Further, it is not known whether localized trends correlate with population trends observed at larger scales. We used breeding season data from 1160 constant-effort banding stations (Monitoring Avian Productivity and Survivorship; MAPS) and a spatially explicit noise model to determine whether abundance and productivity were related to mean noise level or spatial heterogeneity (SD) in noise within a 1-km station radius for 72 passerine species. We also determined whether particular life history traits were predictive of noise responses, and compared continental results to those from local studies. Increasing mean noise level was associated with declines in abundance for 27.1% of species and productivity in 22.1% species. Increasing heterogeneity was associated with declines in abundance for 14.3% species and productivity in 14.7% species. The relationship between noise and abundance was not correlated with the relationship between noise and productivity, and acoustic and life history traits were not related to noise responses. Continental results were similar to localized data in 43.1% and 21.4% of species for abundance and productivity, respectively. Although some patterns differed between the local and continental scale, our results indicate that the MAPS banding dataset is capable of detecting noise-associated impacts on abundance and productivity. This is currently the only large-scale dataset capable of quantifying the relationship between noise and productivity in the continental USA, although other datasets exist elsewhere that may also contribute to our understanding of noise impacts at the larger scale. Bruit, abondance et productivité des oiseaux aux stations de baguage des États-Unis continentaux RÉSUMÉ. Le bruit est une composante de plus en plus courante du monde naturel et causé en grande partie par l'activité humaine. Le bruit d'origine anthropique a un impact négatif sur l'abondance, la santé et la reproduction de nombreuses populations d'oiseaux chanteurs. Quelques études ont fait état d'une diminution d'abondance à plus grande échelle. On ignore si des tendances continentales sont détectées aux stations de baguage, lesquelles fournissent également des données sur la productivité et la survie. De plus, on ne sait pas si les tendances locales sont corrélées aux tendances de populations observées à plus grande échelle. Nous avons utilisé les données en saison de reproduction issues de 1160 stations de baguage à effort constant (programme « Monitoring Avian Productivity and Survivorship »; MAPS) et un modèle de bruit spatialement explicite pour déterminer si l'abondance et la productivité étaient liées au niveau de bruit moyen ou à l'hétérogénéité spatiale (SD) du bruit dans un rayon d'1 km de la station pour 72 espèces de passereaux. Nous avons également déterminé si des traits particuliers du cycle de vie étaient prédictifs des réactions au bruit, et nous avons comparé les résultats continentaux à ceux d'études locales. L'augmentation du niveau de bruit moyen a été associée à une diminution de l'abondance chez 27,1 % des espèces et de la productivité chez 22,1 % des espèces. Une hétérogénéité croissante a été associée à des baisses d'abondance chez 14,3 % des espèces et des baisses de productivité chez 14,7 % des espèces. La relation entre le bruit et l'abondance n'a pas été corrélée avec la relation entre le bruit et la productivité, et les caractéristiques acoustiques et du cycle de vie n'ont pas été liées aux réactions au bruit. Les résultats à l'échelle continentale étaient similaires aux données à l'échelle locale chez 43,1 % et 21,4 % des espèces pour l'abondance et la productivité, respectivement. Bien que certaines tendances diffèrent aux échelles locale et continentale, nos résultats indiquent que l'ensemble de données de baguage MAPS a la capacité de détecter les impacts du bruit sur l'abondance et la productivité. Il s'agit actuellement du seul ensemble de données à grande échelle permettant de quantifier la relation entre le bruit et la productivité aux États-Unis continentaux, bien qu'il existe d'autres jeux de données qui pourraient également contribuer à notre compréhension des impacts du bruit à plus grande échelle.


INTRODUCTION
Human population growth and development have resulted in a rapid and widespread increase in noise levels (Buxton et al. 2017). In many animals, noise pollution alters acoustic communication, abundance, and migratory behaviors (Reijnen and Foppen 1995, Slabbekoorn and Peet 2003, Shannon et al. 2016. Songbirds are especially sensitive to anthropogenic noise because their mating and territory defense rely on acoustic communication (Francis et al. 2011). In particular, occupancy and abundance often decline when noise levels increase (Bayne et al. 2008, McClure et al. 2013. Noise aversion may even outweigh the effects of land cover when songbirds are selecting places to breed (Kleist et al. 2017). A handful of studies also suggest that ambient noise impacts fitness through reduced pairing (Habib et al. 2007) and reproductive success (Halfwerk et al. 2011). Consequently, increasing anthropogenic noise may threaten songbird persistence (Barber et al. 2010). Throughout this paper we use the term "noise" to refer specifically to ambient sound levels at a location, whether environmental or anthropogenic in origin.
Most noise studies on abundance and productivity occur at small spatial extents (< 50 km²; e.g., Francis et al. 2009, McClure et al. 2013, Proppe et al. 2013. Intensive local studies provide fine-scale detail that can account for potentially confounding localized processes such as interspecific competition (Schoener 1983), legacy effects (Cuddington 2011), density-dependent processes (Haldane 1956) and source-sink dynamics (Pulliam 1988). However, these processes often differ geographically, which limits the generalizability of results gained from local studies. For example, negative impacts of noise on reproductive success have been documented in Great Tits (Parus major;Halfwerk et al. 2011) and House Sparrows (Passer domesticus;Schroeder et al. 2012), but recent studies did not detect negative impacts on productivity in either species (Meillère et al. 2015, Halfwerk et al. 2016. Increasing the spatial extent of a study often comes at a cost to fine-scale detail, but access to larger datasets can reveal trends not immediately visible at smaller scales (Rahbek 2005). For example, juvenile and adult male American Redstart (Setophaga ruticilla) distributions are negatively correlated at small scales because of territoriality, but positively correlated at larger scales because of similar habitat requirements (Sherry and Holmes 1988). Four studies have examined the effect of noise on songbird abundance and occurrence at the continental scale, although all were using large-scale datasets to determine whether vocal pitch or particular life history traits were related to abundance and presence in noisy and urbanized areas (Hu and Cardoso 2009, Cardoso 2014, Francis 2015, Moiron et al. 2015. These studies relied on many sources of localized abundance data or published reports describing urban tolerance. Although all four studies documented declines in abundance or occurrence of some species near noise sources and in cities, the importance of vocal pitch differed between studies (discussed in Cardoso et al. 2018). Clearly, a combination of data sources, from both localized and large-scales studies, is needed to best understand the impact of anthropogenic noise on songbird persistence.
The Institute of Bird Populations has been developing a largescale dataset that could provide additional information on noise impacts at the continental level. The MAPS (Monitoring Avian Productivity and Survivorship) program has been tracking avian demographics, including abundance, productivity, and survivorship during the breeding season at banding stations since 1989. This ongoing program has expanded to over 1200 stations, more than 2 million bird captures, and over 230 publications on avian demographics and vital rates. Effort is standardized across all MAPS stations so that data can be compared across time and space. To our knowledge, this is the only large-scale dataset that also provides information on productivity in North America, although similar datasets exist elsewhere (e.g., Eglington et al. 2015). Quantifying productivity provides a measure of birth rates in a population, and is often strongly tied to population trends (Sallabanks et al. 2000).
Two additional datasets provide information about noise and urbanization at each MAPS banding station. First, the United States National Park Service (NPS) has developed a geospatial sound model that uses 109 natural and anthropogenic factors associated with ambient noise levels (Mennitt et al. 2014) to predict the L 50 dBA at a 270 m resolution for the a typical summer day from 2003 to 2013 across the entire continental United States (CONUS). Second, the United States Geological Survey (USGS) produced CONUS land cover maps in 2001,2006, and 2011 that include habitat classifications and a measure of development at 30 m resolution (Homer et al. 2015).
We assessed the relationship between model-predicted noise levels and the abundance and productivity of 72 passerine species by computing the weighted mean of noise within 1 km of each banding station (Fig. 1). Visual investigation of noise data revealed that a number of MAPS stations had high levels of noise heterogeneity within the 1 km radius. To determine whether this noise pattern impacted species differently than uniformly high or low noise levels, we assessed the impact of noise heterogeneity on abundance and productivity. Species were selected because they were among the 30 most commonly captured at MAPS stations or because localized data on noise impacts were available. Local studies included 49 species, seven which were also among the 30 most common species. After controlling for the effects of location, time, development, and habitat, we determined whether noise was an important predictor of abundance or productivity for each species. We also calculated whether the directional associations of noise with abundance and productivity were correlated, such that knowing the relationship of noise with one measure facilitated a predictive understanding of the other. Finally, we determined whether the relationship with noise was predicted by particular life history traits, and we compared continental data to local data to determine the level of agreement between these two scales.

Datasets
Data on avian demographics from 1989 to 2016 was provided by The Institute of Bird Populations MAPS banding program. Each participating station consists of ten 12-m mist nets located on an area approximately 8 ha in size that is contained within a contiguous landscape of > 20 ha. Mist nets are opened for six hours (beginning at dawn) on one day in a 10-day period. Based on the latitude of the station, six to ten 10-day periods are run consecutively during the breeding season (May-August) each year. Age, sex, and reproductive status are recorded for each bird, and a metal USGS band is attached to the tarsus before release to facilitate the collection of recapture data.
The United States National Park Service used 109 natural and anthropogenic factors associated with ambient noise levels (Mennitt et al. 2014) to predict the L 50 dBA at a 270 m resolution for the entire continental United States. Noise predictions were tested at 319 acoustic monitoring sites within 216 national parks and 34 urban locations in 28 unique cities spanning the continental United States. The model explained 71.74% of the noise variance across all sites. Anthropogenic sources increased the noise level predicted by biophysical processes alone in 78% of the landmass. The model is limited to a single noise measure per pixel, calculated for a typical summer day from 2003 to 2013.
The USGS periodically produces a 30-m resolution layer of habitat classifications and development levels as part of the National Land Cover Database (NLCD). We used 2001 data because this was the year closest to the median year of bird captures available in the MAPS dataset. For our analysis we incorporated two measures of development and seven measures of habitat to control for their effect on bird demographics. Development variables in the NLCD were derived from the percent cover of impervious surfaces, which is a strong predictor of the impact urbanization has on songbirds (Melles et al. 2003, McClure et al. 2015, and the type and quantity of human structures (Homer et al. 2012). We combined open (< 20% anthropogenic cover) and low (20-49%) development into a single category and mid (50-79%) and high (80-100%) development into a second category. Habitat variables included the proportion of land covered by water, wetlands, forest, herbaceous, shrubland, barren, and planted/cultivated.
To determine whether life history traits were predictive of noiseassociated changes in abundance and productivity, we categorized bird species into several guilds based on (1) song frequency, (2) breeding habitat, (3) diet, (4) nest location, and (5) migratory behavior. Vocal frequencies were derived from two to eight highquality songs for each species housed on Thayer Birding Software (Birds of North America, Gold Edition). Because multiple elements of song frequency have been implicated in noise response and can be modulated independently, we measured peak frequency, bandwidth (25%-75% summed energy), lowest frequency (5% of summed energy), lowest quartile (25% of summed energy), and highest frequency (95% of summed energy) in Raven Pro (V1.4,Ithaca,NY). For species lacking a song, the most common call was utilized. All other life history traits were collected from the Cornell Lab of Ornithology's website (https:// www.allaboutbirds.org/news/). Breeding habitat classifications included: grassland, forest, open woodland, scrub, and town. Birds were classified as insectivores or noninsectivores based on the primary components of their diet. Nest location was classified as cavity, ground, shrub, or tree. For migration, species were categorized as nonmigrants, migratory, or partial migrants. Species were categorized as partial migrants when range maps indicated that populations were resident in 25-50% of their range, but otherwise migratory.
To eliminate the potential for implicit bias, we evaluated regional studies for inclusion in our study during September of 2016, prior to initiating the processing of banding data. We established two a priori criteria for inclusion. First, each study must have occurred within the North American region where both MAPS and noise data were available, or examined species commonly found within this region. Second, each study must have included directional results for abundance or productivity in relation to noise or proximity to noise-producing sources, e.g., roads or oil wells. Although the majority of the literature on anthropogenic noise addresses behavioral responses in songbirds, we identified 13 suitable studies using multiple subsets and variants of the search terms; songbird, noise, abundance, reproduction, and productivity in Web of Science and Google Scholar (Appendix 1). Many of these studies examined multiple species, resulting in local demographic information for 49 species. Abundance data was available for 44 species, and productivity was available for 16 species. Of these species, 11 had data on both abundance and productivity.

Data processing
For each species, abundance was calculated as the number of adults captured per 600 net hours, and productivity was the ratio of juveniles (AGE code 4 or 2) to adults (all other AGE codes, except 0 and 9 where age is unknown; DeSante and Kaschube 2009). Although effort was standardized, we cannot assume that capture rates were equal across all sites. However, there are no indications that noise alters capture rates in a systematic fashion when birds are on the breeding grounds (although it may impact site selection during migration, when birds are actively selecting stopover sites; McClure et al. 2017). This information, along with our large sample size, should overcome any potential spurious associations based on differences in site-based capture rates. Survivorship was not analyzed because birds are often recaptured long after they were initially banded, and at sites with noise levels and surface characteristics that are quite different from the original banding location. Thus, a myriad of factors beyond the particular characteristics at the site of origin are likely to impact survivorship. For any given species, only stations within that species' range (breeding status code B, U, O, and T) were included (DeSante et al. 2015). A minority of MAPS stations are operated more than once per 10-day period. For such stations we only included captures during the first sampling subperiod (code A) to standardize effort across all stations. Because MAPS regulations suggest that data from locations < 1 km from the banding station's center are not independent (DeSante et al. 2015) and breeding birds are unlikely to venture > 1 km from their breeding territory, we calculated the mean noise level within a 1 km radius of each MAPS station. The mean was linearly weighted to increase the influence of more central pixels by dividing the noise level in each pixel by the distance between that pixel and the station center, and then averaging the results for all pixels. To evaluate spatial heterogeneity we calculated the standard deviation of noise levels in all pixels contained within the 1 km radius. Although heterogeneity can be measured in many ways, we chose standard deviation because it is used regularly in the literature (e.g., Yang et al. 2015, Khan et al. 2019, is easy to interpret, and effectively creates a standard measure of variability comparable across sites with high and low mean noise levels. Habitat and development variables were calculated as the proportion of all pixels within a 1 km radius surrounding each station. We used a principal component analysis on the development and habitat variables to obtain a set of uncorrelated variables for inclusion as predictors in regression models. We retained eight principal components that accounted for > 80% of the variation (Appendix 2). Each component was given a relevant name based on the top one to three contributing variables. In the order of explanatory power, these components increased with "development" (PC1; 20.08% variance explained), "wetland" (PC2; 11.75%), "shrub" (PC3; 10.78%), "barren" (PC4; 9.22%), "marsh" (PC5; 8.15%), "grassland" (PC6; 7.49%), "water" (PC7; 7.46%), and "urban" (PC8; 6.84%).

Statistical analysis
For each of 72 species of interest, we fitted two generalized additive mixed-effects models (GAMMs): one with abundance as the response variable, and one with productivity. Abundance models were Poisson or quasi-Poisson (quasi-Poisson, in cases where the overdispersion parameter estimate exceeded 2) GAMMs with log link functions (Hastie and Tibshirani 1990, Zuur et al. 2009, Wood 2017. Predictors included the weighted mean noise level, the standard deviation of the noise level, the eight habitat and development principal components, and location (2-dimensional smooth term; thin-plate regression spline with basis of dimension 16; Wood 2017). An offset (logarithm of net-hours) accounted for differential survey effort between stations and years (Zuur et al. 2009). Productivity models were Poisson or quasi-Poisson GAMMs with log link functions (Zuur et al. 2009). The response variable in these models was the count of juvenile birds, with an offset (logarithm of the count of adult birds), so we effectively modeled the number of juveniles per adult (Zuur et al. 2009). Candidate predictors were the same as for abundance models. Year (smooth term; thin-plate regression spline with basis of dimension 5) was included as an additional predictor in all models to account for nonmonotonic patterns over time in a variety of factors, including land use changes. All models also included a random effect of station. Although these GAMM models detected counts, rather than yielding absolute bird abundance or productivity estimates, species-specific modeling and inclusion of covariates related to location, habitat, and time should still allow reliable assessment of noise effects on indices of abundance and productivity (Dunn and Ralph 2004). We used variance inflation factors to confirm the absence of problematic multicollinearity between predictors (particularly noise and development/habitat PCs; Dormann et al. 2013). Four candidate models were considered for abundance and productivity of each species, including (a) both noise variables, (b) mean noise only, (c) heterogeneity of noise only, and (d) no noise variables. Habitat, development, and spatial predictors were retained in all candidate models to control for their effects on abundance or productivity. Best models for each species and response variable were selected using AIC (or QIC, for quasi-Poisson models). Where best models included a noise term, but the AIC/QIC improvement relative to the next-best model without the noise term was below 2 (ΔAIC/QIC < 2), it was noted by labeling the effect as "weak"; (Appendix 3; Burnham andAnderson 2002, Burnham andAnderson 2004). Model fitting was carried out in R version 3.6.1 (R Core Team 2019) using package mgcv (Wood 2011(Wood , 2017. To investigate the relationship between noise-associated changes in abundance and productivity, we plotted and modeled the productivity regression coefficients as a function of the abundance coefficients. Separate models were fitted for weighted mean noise levels and noise heterogeneity coefficients. Although top abundance and productivity models excluded noise terms for some species, those terms were included here to produce effect sizes for all species. To account for uncertainty in the GAMM coefficient estimates used as input to these models relating abundance and productivity results, we used measurement error models, using the standard error of the estimates as the measurement error standard deviations (McElreath 2020). Models were fitted in R using package brms (Bürkner 2018). To assess the relationship between the abundance and productivity noise-effect coefficients, we examined 95% credible intervals for the abundance coefficients in the measurement error models, and additionally compared each model to its intercept-only version via leave-one-out cross-validation (brms function loo(); Vehtari et al. 2017).
To investigate relationships between noise responses and life history traits, we fitted a measurement error model (using standard errors in the GAMM coefficient estimates as the measurement errors, as above) for each noise coefficient (abundance and productivity, weighted mean noise and noise standard deviation) as a function of the suite of vocal and life history traits described above. These included peak frequency, bandwidth, lowest frequency, lowest quartile, highest frequency, habitat, insectivory, nesting, and migration. Models were fitted with horseshoe shrinkage priors (Piironen and Vehtari 2017) in R package brms (Bürkner 2018), and the best model for each noise coefficient was selected after examining model output (estimates and credible intervals), considering output from projection predictive forward variable selection (using R package projpred; Piironen et al. 2019), and as verification, leave-one-out cross-validation scores for the top models identified by the variable selection.
To determine whether regional noise impacts were similar to our landscape results, we compared the directional effect (+, 0, -) of increased noise on abundance and productivity across scales and noted the percentage of cases demonstrating agreement, and vice versa.

Continental models for abundance and productivity
The distribution of model-predicted noise levels at MAPS stations (mean = 37.7, SD = 5.1, range = 23.9-54.8, n = 47,048) was similar to that of the full NPS dataset (mean = 35.6, SD = 5.3, range = 20.1-67.0, n = 183,300,975). Our best species-specific landscape models retained a negative relationship between the weighted mean of noise and abundance in 19 species (27.1%), a positive relationship in 19 species (27.1%), and no clear statistical relationship in 32 species (45.8%; Table 1, species-specific effect sizes and model rankings can be found in Appendix 3 and 4, respectively). As spatial heterogeneity in noise increased, abundance decreased in 10 species (14.3%), increased in 27 species (38.6%), and had no clear statistical relationship in 33 species (47.1%). MAPS data was insufficient to build continental-scale abundance models for two species (Golden-cheeked Warbler [Setophaga chrysoparia] and Tennessee Warbler [Leiothlypis peregrina]) and these species were not used in the calculation of percentages above (n = 70). Productivity was negatively related to increasing mean noise levels in 15 species (22.1%), positively related in 15 species (22.1%), and not statistically related in 38 species (55.8%). Productivity and noise heterogeneity were negatively related in 10 species (14.7%), positively related in 22 species (32.4%), and not statistically related in 36 species (52.9%). MAPS data was insufficient to build continental-scale productivity models for four species (Golden-cheeked Warbler, Tennessee Warbler, Eastern Meadowlark [Sturnella magna], and White-crowned Sparrow [Zonotrichia leucophrys]) and these species were not used in the calculation of percentages above (n = 68). Model failures in these species resulted from low capture numbers, which was primarily due to a narrow geographic range, breeding primarily outside the continental United States, or the inability to capture juvenile grassland species with MAPS protocols.
In 16 species, no changes in abundance and productivity attributable to the weighted mean of noise were identified. Noise heterogeneity had no detectable association in 19 species. Measurement error models indicated no correlation between coefficients for effects of weighted mean noise on abundance and productivity (estimated slope coefficient (with shrinkage) was 0.00, and the intercept-only model was the best model, with a leave-one-out information criteria (LOOIC) difference of just 0.02; Figure 2). The corresponding results for the effects of noise spatial heterogeneity were similar (estimated slope coefficient 0.00 with shrinkage; full model LOOIC was lower than the interceptonly model, but only by 0.03; Figure 3). Neither acoustic nor life history traits improved on the intercept-only model for mean or heterogeneity of noise in abundance or productivity (estimated slope coefficient 0.00 with shrinkage; LOOIC intercept-only model lowest or within 0.13 of all other models; Appendix 5).
Principal components related to habitat and development were retained in all candidate models to control for the expected effects of these conditions rather than undertaking an explicit assessment of their contribution. However, 44 of 70 abundance models contained significant principal component predictors (p < 0.05) and 25 of 68 productivity models also contained principal component predictors (Table 1). In contrast, one or more noise terms were retained in 58 abundance models and 38 productivity models. Although beyond the scope of this report, information about the relative impact of habitat and development on particular species can be gleaned from these results (Appendix 3).

Fig. 2.
Model coefficients for the effects of weighted mean noise level on the abundance and productivity of individual species (Δ Abundance or Productivity/Decibel[dBA SPL]). Coefficients, and their associated standard errors, were derived from speciesspecific models that predicted the relationship between mean noise level and the demographic trait. Error bars indicate 95% confidence intervals. Inset provides an expanded view of species near the origin.

Fig. 3.
Model coefficients for the effects of noise standard deviation (SD) on the abundance and productivity of individual species (Δ Abundance or Productivity/SD Unit). Coefficients, and their associated standard errors, were derived from species-specific models that predicted the relationship between noise heterogeneity and the demographic trait. Error bars indicate 95% confidence intervals. Inset provides an expanded view of species near the origin.
Erratum: Insets were added to Figures 2 and 3 (PDF version) to improve comprehensibility and replaced those in the original publication on 4 August 2020. Table 1. Results from top models depicting the impact of noise level (mean) and heterogeneity (sd) on abundance and productivity for each species. Where available, local results and agreement with Monitoring Avian Productivity and Survivorship (MAPS) continental scale results are presented. For ease of reading, effects are presented as positive (+), negative (-), or no effect (0). Principal components significant at p < 0.05 are included as positive (no sign) or negative (-) effects. Numerical effect sizes can be found in Appendix 4. See Appendix 3 for scientific species names. Where data from multiple local studies is presented, numbered superscripts refer to study superscripts found in Appendix 1.

Comparison of continental and local studies
Results from local studies were only compared to our continental results for the weighted mean of noise because each local study used the mean to quantify noise, rather than a measure of spatial heterogeneity. For abundance, data at both scales was available for 51 cases, which included 43 of the 44 available unique species (one continental model failed) and eight cases where multiple local studies were available on a species (two studies in six species, three in one species; Table 1). The directional effect of noise across multiple local studies was congruent for three of these species (42.9%). The directional effects of noise on abundance were congruent between scales in 22 of 51 comparisons (43.1%), although the trend was reversed (i.e., + to -or vice versa) in only seven studies (13.7%). One local, multispecies study assessed noise impacts during migration (McClure et al. 2013), a time period that may exert different pressures than during the breeding season when MAPS data was collected. However, removing this study did not improve agreement (10 of 29 cases; 34.5%).
For productivity, data at both scales was available for 14 cases, which included 13 of the 16 available unique species (three continental models failed) and one case where an additional local study was available on a species (Table 1). The directional effect of noise between studies was not congruent in this species. The directional effects of noise on productivity between spatial scales was congruent in only 3 of 14 cases (21.4%), although the trend was reversed for only two species (14.2%). . This likely stems from the fact that productivity is impacted by multiple reproductive variables, e.g., fecundity, number of broods, or postfledging survival. Publication skew could also contribute, because local studies are more likely to be conducted on populations where a trend is suspected a priori (Møller and Jennions 2001). It is also plausible that noise is not related to productivity in the majority of species at the continental scale because lowered abundance in noisy areas reduces competition for habitat and foraging resources (Francis et al. 2011). As a result, smaller populations may continue to be proportionally as productive as larger populations. However, if productivity does not increase when abundance decreases, it is likely that populations in noisy areas will remain small or decrease in size over time. Of the 20 species that became less common as mean noise levels increased, only two also had positive productivity (Gray Flycatcher [Empidonax wrightii] and Red-eyed Vireo [Vireo olivaceus]).

DISCUSSION
Perhaps one of the primary benefits of MAPS banding data is the ability to compare noise-associated changes in abundance and productivity simultaneously across a large number of sites and years. In our results, the relationship with noise in one population parameter did not consistently predict the direction of the other. In general, ecologists do not expect abundance and productivity to be correlated, largely because of density-dependent processes (Van Horne 1983, Murdoch 1994. However, in the presence of an external stressor, such as noise or urbanization, it is plausible that both processes could be reduced in parallel by a reduction in fitness that directly impacts habitat selection and the production of offspring simultaneously (Bock and Jones 2004). Our results indicate that the interaction between noise and population demographics is more complex, and that monitoring long-term changes using both measures in tandem might be the most useful for predicting population persistence in noisy areas. In addition to providing multiple demographic measures, data from banding stations is unaffected by some of the confounding variables that typically plague noise studies. Specifically, banding data is robust to the problem of observer detectability. Variation in aural abilities between observers can limit comparisons from aural surveys, a primary method for assessing abundance in avian studies (Faanes and Bystrak 1981). Ambient noise level is also known to impact detectability (Simons et al. 2009), making it difficult to compare data across studies unless similar methodology allows for statistical corrections (e.g., Buckland et al. 1993). Barring the presence of a point source of noise physically located near a mist net, there is little reason to believe that ambient noise has any impact on the likelihood of capturing resident birds that are present on the breeding or winter grounds. Further, having a bird in hand, where photographs, previous bands, and expert help can be utilized, reduces error and variability in identification and counts between observers.
Several studies have investigated whether particular life history traits underlie differential species responses to noise. In our analysis, neither acoustic traits, nor a suite of guilds defined by life history traits were predictive. The lack of any effect was somewhat surprising given recent literature on the subject that has implicated vocal frequency, diet, and nesting as traits related to noise responses (Hu and Cardoso 2009, Cardoso 2014, Francis 2015. However, in each of these published studies, life history traits were compiled for several hundred species, often selected to cover a range of variation in vocal features and life history traits. Despite their extent, some results from these large-scale studies are still conflicting (Moiron et al. 2015). Our study included 72 species that were selected because they were commonly captured in mist nests at MAPS stations or because continental data could be compared to data from local studies. Because our primary objective was to evaluate the viability of detecting noiseassociated changes in MAPS data, we made no attempt to capture a range of vocal or life history traits. Thus, variation in some traits was low. For example, only four species of grassland birds were included in our analysis, and very few of our species were not insectivores. Further, noise in our study was not constricted to low-frequency sources, which may reduce the role of vocal frequency often seen in urban adaptation (Halfwerk and Slabbekoorn 2009). Although a more explicit study of life history traits across species captured at MAPS stations could be valuable, we suggest that our results should be limited to evaluating the contribution of MAPS data for detecting noise-associated patterns in bird abundance and productivity.
Although mean noise level is the most common unit used to assess the impacts of noise, spatially heterogeneous acoustic habitats may impact birds differently (Gill et al. 2015). Acoustic niches, for example, could be more diverse when noise levels vary across space. As a result, birds with both philic and phobic reactions to noise might inhabit the same general area. In keeping with this hypothesis, 85.7% of species increased or remained stable as noise heterogeneity increased, while only 72.9% did so as mean noise level increased. Although our effects sizes are small, the skew toward more positive relationships with noise heterogeneity (Figs. 2 and 3) may provide an acoustic example of the intermediate disturbance hypothesis, which predicts that the highest number of species will thrive when habitats are neither uniformly pristine nor highly degraded (Connell 1978). At minimum, the negative relationship with spatially heterogeneous noise is less prominent than that of a uniform increase in noise. This suggests that banding stations with elevated noise levels on only a portion of their localized habitats may support more diverse songbird communities than those completely inundated by elevated noise levels. Perhaps this is not surprising given that localized noise sources, such as rivers, have been on the landscape for much of geologic history. Another possibility is that these communities are experiencing a release from predation by noise averse predators (Francis et al. 2009, Barber et al. 2010, without incurring the full costs associated with uniformly increased noise levels. Clearly, a better understanding of the spatial extent at which noise impacts songbirds is needed. More generally, studies of animal ecology must continue to move away from interpretations based solely on averages. Variability in many sensory cues, habitat characteristics, temporal patterns, and even life-history traits may impact populations differently. Combining results from the continental-scale MAPS banding database with localized studies may also be important for interpreting noise impacts. Congruent results from the two spatial extents might provide confidence for widespread species-specific management techniques. One example is the Red-breasted Nuthatch (Sitta canadensis), where abundance was negatively associated with noise in two localized studies and in our continental assessment. Taken together, these results indicate that noise mitigation will likely be needed to maintain healthy population levels in this species. Conversely, the abundance of the Western Tanager (Piranga ludoviciana) was negatively impacted in multiple regional studies, but not in our continental model. In this case, the fine-grain detail of localized studies may have detected patterns not present or detectable in MAPS data, and localized or general noise-mitigation might still be recommended for species management. Finally, the MAPS dataset detected no relationship between noise and productivity in House Sparrows, a pattern that is supported in one localized study (Meillère et al. 2015), but not in another (Schroeder et al. 2012). The lack of an outcome at the larger spatial extent might reduce concern that noise is impacting this species, despite the more variable localized results.
Determining the impact of noise on populations dynamics in particular species is critical as we prioritize mitigation objectives for songbirds. Extrapolating solely from local studies may lead to misguided management because of the presence of differential local and landscape level processes. Given that broad metaanalyses to date have been unable to consistently predict noise impacts based on a single set of life-history traits, large scale datasets capable of determining patterns of change directly for individual species may be a valuable tool for assessing noise impacts. Our results suggest that MAPS data can be used to detect noise-related changes at the continental level. The value of this dataset is enhanced by the ability to incorporate habitat and development variables alongside noise terms to assess patterns in multiple population parameters, including, abundance, productivity, and perhaps eventually survivorship. In our models, over a quarter of the species surveyed experienced declines in abundance or productivity beyond what could be attributed to other variables, highlighting the importance of noise mitigation for maintaining songbird diversity. Because the direction of noiseassociated changes was not correlated between abundance and productivity, and continental and localized results were not consistently congruent, best management practices should strive to incorporate multiple demographic measurements from multiple scales. Agreement across scales and population parameters may reveal which species are most and least sensitive to noise. Future attempts to model noise over time, or take standardized noise measurements at MAPS stations, might increase the sensitivity of this dataset. However, the current MAPS dataset can be used to derive continental, noise-associated trends for any bird species that is regularly captured in mist nets, adding a new tool for biologists seeking to understand and manage songbird populations in an age of increasing noise levels.

Data Availability Statement
Data (allspecies_final.csv) that supports the findings of this study is available online in the OSF public repository (https://osf.io/ cka9z/). R code for single species model analysis is available in Appendix 6 and code for the measurement error models is in Appendix 7.  Table A4.138: Model selection information for the productivity model for the Yellow Warbler (Setophaga petechia). The model was a quasi-Poisson GAMM with log link function. Based on QIC, support for including the noise level as a predictor of productivity was: none; support for including the standard deviation of noise was: none.