Next Article in Journal
Reducing the Extinction Risk of Populations Threatened by Infectious Diseases
Next Article in Special Issue
Rusty Blackbird (Euphagus carolinus) Foraging Habitat and Prey Availability in New England: Implications for Conservation of a Declining Boreal Bird Species
Previous Article in Journal
An Overview of the Helminths of Moor Frog Rana arvalis Nilsson, 1842 (Amphibia: Anura) in the Volga Basin
Previous Article in Special Issue
Occupancy of the American Three-Toed Woodpecker in a Heavily-Managed Boreal Forest of Eastern Canada
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Flock Size Predicts Niche Breadth and Focal Wintering Regions for a Rapidly Declining Boreal-Breeding Passerine, the Rusty Blackbird

by
Brian S. Evans
1,2,
Luke L. Powell
1,3,4,*,
Dean W. Demarest
5,
Sinéad M. Borchert
6 and
Russell S. Greenberg
1,†
1
Migratory Bird Center, Smithsonian Conservation Biology Institute, National Zoological Park, P.O. Box 37012, Washington, DC 20013, USA
2
Department of Biology, Georgetown University, 37th and O Streets NW, Washington, DC 20057, USA
3
Institute of Animal Health and Comparative Medicine, Graham Kerr Building, University of Glasgow, Glasgow G128QQ, UK
4
Department of Biosciences, Durham University, Stockton Road, Durham DH13LE, UK
5
U.S. Fish & Wildlife Service, Division of Migratory Birds, Atlanta, GA 30345, USA
6
U.S. Fish & Wildlife Service, Panama City Field Office, Panama City, FL 32405, USA
*
Author to whom correspondence should be addressed.
Deceased.
Submission received: 30 November 2020 / Revised: 28 January 2021 / Accepted: 29 January 2021 / Published: 4 February 2021
(This article belongs to the Special Issue Boreal Bird Ecology, Management and Conservation)

Abstract

:
Once exceptionally abundant, the Rusty Blackbird (Euphagus carolinus) has declined precipitously over at least the last century. The species breeds across the Boreal forest, where it is so thinly distributed across such remote areas that it is extremely challenging to monitor or research, hindering informed conservation. As such, we employed a targeted citizen science effort on the species’ wintering grounds in the more (human) populated southeast United States: the Rusty Blackbird Winter Blitz. Using a MaxEnt machine learning framework, we modeled patterns of occurrence of small, medium, and large flocks (<20, 20–99, and >99 individuals, respectively) in environmental space using both Blitz and eBird data. Our primary objective was to determine environmental variables that best predict Rusty Blackbird occurrence, with emphasis on (1) examining differences in key environmental predictors across flock sizes, (2) testing whether environmental niche breadth decreased with flock size, and (3) identifying regions with higher predicted occurrence (hotspots). The distribution of flocks varied across environmental predictors, with average minimum temperature (~2 °C for medium and large flocks) and proportional coverage of floodplain forest having the largest influence on occurrence. Environmental niche breadth decreased with increasing flock size, suggesting an increasingly restrictive range of environmental conditions capable of supporting larger flocks. We identified large hotspots in floodplain forests in the Lower Mississippi Alluvial Valley, the South Atlantic Coastal Plain, and the Black Belt Prairie.

1. Introduction

The Rusty Blackbird (Euphagus carolinus) is a Nearctic icterid that exhibits near exclusive reliance on forested wetland habitats in boreal, transitional and deciduous biomes throughout its annual cycle [1,2,3]. The species is among the most steadily and precipitously declining temperate North American landbirds, [4,5,6] and is recognized as “Vulnerable” by the International Union for the Conservation of Nature. There is some anecdotal evidence that factors in boreal breeding regions of northeast North America may be limiting the population (timber management and nest success: [7]; mercury: [8]). However, demographic rates for Rusty Blackbirds across the Boreal are similar to those of other songbirds [1]), suggesting that the overall population is likely limited by factors outside the breeding grounds.
Specifically, habitat modification on the species’ wintering grounds has been implicated as a likely driver of the decline [5,9], stressing the importance of studying this period of the species’ life cycle. Further, population monitoring and broad assessments of habitat selection across the boreal breeding range are extremely challenging given the general remoteness and the very low density of roads, humans, and Rusty Blackbirds themselves. The Breeding Bird Survey, for example, detects so few Rusty Blackbirds that with the exception of more developed parts of Alaska it is not possible to use these data to monitor population trends. All of these factors coupled with the species’ considerably more social (flocking) nature on the wintering grounds and the high densities of citizen scientists (birders) in the southeast United States point to the advantages of studying the species during winter. The wintering grounds are not without their challenges, however; diffuse and localized distribution, nomadism, and difficulty in accessing often ephemeral wetland habitats constrain our ability to dependably study wintering Rusty Blackbirds [9]. Conventional means of investigation, such as use of mark-recapture (e.g., to evaluate survival as a function of winter habitat metrics), may not be viable due to high vagility and low site fidelity between years (e.g., [10,11].
Once reported in wintering flocks numbering tens of thousands [4], Rusty Blackbirds are now occasionally observed in flocks approaching a few thousand individuals, with much smaller groups being more common [12,13]. Still, flocking behavior in Rusty Blackbirds may provide a proximate cue allowing researchers to identify and characterize winter habitat (e.g., key sites, condition, relative quality or suitability for foraging). Foraging theory predicts that abundance and distribution of food resources influence group size, with more resource rich environments favorable to flocks of larger size (sensu [14,15,16,17,18]). Relationships between group size and resource availability have been shown in a wide variety of taxa, including primates (e.g., Lemur catta [19]), fish (e.g., Fundulus diaphanous [20]), and birds (e.g., Junco phaeonotus [21]). Where food resources are less available, interference competition can moderate flock size as individuals seek to optimize foraging efficiency beyond the confines of competitive groups [15]. Demographic studies in a number of bird species support that the size of conspecific aggregations is positively associated with higher quality habitats [22,23].
Because flock size is expected to vary as a function of habitat quality or compatibility, evaluating environmental features associated with flocks of different sizes may provide evidence for conditions and locations most favorable for greater numbers of wintering Rusty Blackbirds. Assessing distributions in environmental space, as opposed to geographic space, is necessary for the species as wintering individuals often travel long distances across the landscape, and thus site level habitat characteristics have not consistently predicted occupancy [2]. An understanding of such relationships could offer immediate relevance in informing conservation programs and facilitating further research and monitoring supporting full annual cycle stewardship of this species, as well as yield new insights regarding niche dynamics in flocks of different sizes.
Here, we evaluate the relationships of climatic and landcover attributes in predicting the occurrence of Rusty Blackbirds across their wintering range and ask whether these relationships vary with flock size. Previous studies have observed that, at local scales, large flocks of Rusty Blackbirds are constrained to a narrow set of environmental conditions and geographic areas, whereas smaller flocks are comparatively widespread (e.g., [12,24,25]). We sought to quantify these observations at the spatial extent of winter range of the species, predicting that (1) large flocks of Rusty Blackbirds (>100 observed individuals) occupy a different environmental niche than small or medium flocks (<20, 20–99, respectively) and (2) large flocks occupy a narrower environmental niche breadth than small or medium flocks [26]. We define niche breadth as the multidimensional range of specific biotic and abiotic conditions that a group of organisms occupies [27,28,29]. We used occurrence data from two citizen science programs, the Rusty Blackbird Winter Blitz (Blitz; http://rustyblackbird.org/outreach/migration-blitz/) and eBird [30], to test these predictions by modeling the distributions of Rusty Blackbird flock size classes in environmental space. In a closely related project, we performed a three-year “Spring Blitz” for migrating Rusty Blackbirds, but those results are addressed elsewhere [31]. We compared biotic and abiotic attributes that best predicted occurrences of flock size classes and evaluated differences in environmental niche breadth among classes. We suggest focal regions for wintering Rusty Blackbirds (hotspots), defined as regions with higher predicted occurrence, and discuss implications to future research and conservation. Finally, given the proliferation of citizen science monitoring projects, we evaluated whether the Blitz, a targeted volunteer effort requiring considerable investment and coordination, yielded observational data that improved predictive strength of relationships relative to use of comparable eBird data alone.

2. Methods

We used a maximum entropy (MaxEnt) modeling approach to examine habitat suitability across flock size classes and observational methods. In species distribution modeling, habitat suitability is defined as the likelihood of occurrence of a species in association with environmental variables ([32], but see [33]). MaxEnt is a machine-learning method that compares occurrence data with that of background samples in environmental space (see Supplementary File S1 1.1; [34]). In this context, training samples represent subsets of locations in which Rusty Blackbirds were observed. Background samples represent subsets of locations that were sampled but no Rusty Blackbirds were observed.
We used Rusty Blackbird observations submitted from the Blitz, as well as those collected through traditional eBird [30] protocols. The Blitz was a coordinated effort among citizen scientists to search for Rusty Blackbirds across the species’ wintering range by conducting traveling counts of 8 km or less during target dates of 1 January–28 February 2009–2011. Blitz observations were submitted using a special eBird portal (https://ebird.org/home) and recorded date, location, species, numbers, effort and other attributes (see [35]). To supplement Blitz data, we used traditional eBird checklists that reported Rusty Blackbirds on traveling counts conducted in 2009–2011 during corresponding Blitz periods [30]. We considered each eBird and Blitz checklist as an independent observation, omitting duplicates (i.e., observations submitted by multiple persons). We georeferenced observation locations with a 4-km resolution raster grid (see below). To avoid pseudoreplication, if a grid cell contained >1 independent observation we retained only the observation recording the highest Rusty Blackbird count.
We classified Rusty Blackbird observations by flock size (small: <20; medium: 20–99; large: >99 individuals), and sampling protocol (Blitz or eBird). Flock size classes were determined by the frequency of observations of a given number of individuals and approximate the lower (small flocks) and upper (large flocks) quantiles of observations within the combined pool of eBird and Blitz samples. While we use the term “flock” to designate the number of birds observed, both Blitz and traditional eBird observations were reported as traveling counts. Thus, it is possible that in some cases, aggregates designated as “large flocks” may be comprised of several small or medium-sized true flocks.
Given that flocks aggregate in common roosts at night and may be observed traveling in smaller groups during the day, we examined whether the time of observation biased our results. We determined the time of observation for each checklist used in our analysis (i.e., the maximum observed Rusty Blackbirds within a given raster cell). All time values were standardized to Coordinated Universal Time as a function of their geographic coordinates. Given that day length varies by geographic location and over the two-month period of this study, we used the R package activity [36] to transform the times of observation to solar times, which are defined here as radian time values relative to civil sunrise times for a given date and location [37,38]. We then fit Von Mises circular kernel density distributions to detections within each flock size class [36]. Following Ridout and Linkie [39], we calculated the degree of overlap ( Δ ^ 4 ) between fitted distributions and used randomization (with replacement; n = 10,000 iterations) to test the null hypothesis that detections come for the same distribution (α = 0.05). For each of the flock size class pairs, we failed to reject the null hypothesis (File S1), providing supportive evidence that the time of day in which an observation occurred was unlikely to introduce bias into our results.
A limitation to occurrence-only modeling methods such as MaxEnt is that occurrence data, especially data collected opportunistically, are often biased toward areas of higher sampling effort (e.g., as a function of accessibility, convenience, or availability of human observers; see [40]. This confounds inferences regarding distribution, habitat use and environmental niche, and may cause inflated measures of model performance [41,42]. To mitigate this potential bias in our models, we generated background data by summarizing eBird checklists (2009–2011; n = 13,218 checklists) across each 4 km grid cell. In doing so, this method assumes that the sampling bias associated with checklists observing Rusty Blackbirds is similar to that of eBird checklists where Rusty Blackbirds were not observed [43,44]. By using random eBird sampling locations as background points, rather than random locations from the entire study extent, we were able to construct models that evaluate the occurrence of Rusty Blackbird observations in environmental space relative to points on the landscape presumed to be representative of any geographic biases in our Rusty Blackbird observations (see [45]). Additionally, we limited background point selection, and thus the extent of our distribution models [46], to the geographic extent that contained 99% of Rusty Blackbird observations, which corresponded to the land area of the Eastern United States, bounded in the west at a longitude of −100°.
We constructed models with up to 14 environmental covariates to predict Rusty Blackbird habitat suitability, including two climatic and 12 land cover covariates (see Table S1). We used the raster package in Program R [47,48] for all environmental layer processing. We obtained mean precipitation and minimum temperature data for the months of January and February 2009–2011 (4-km resolution, [49]). Though there was considerable annual variation in average minimum temperature and precipitation within our study area, small annual sample sizes of Rusty Blackbird observations necessitated averaging these covariates across the three winters. We obtained 30-m-resolution land cover data from the US Gap Analysis project [50] and reclassified Gap classes into the following 12 categories: floodplain forest, hardwood forest, plantation hardwood forest, upland forest, mixed forest, woodland, shrub, wetland, grassland, pasture, row crops, and developed land (Figure 1, Table S1). We calculated the mean proportional cover of each category within a 4 km resolution raster grid. This spatial scale is roughly equivalent to winter home range sizes observed for Rusty Blackbirds [12]. The values of all environmental variables were extracted to the spatial locations of Rusty Blackbird observations and eBird background samples.
Models were constructed using MaxEnt 3.3.3 [34] and implemented in the R package dismo [51]. We randomly partitioned observations into five replicates (k-fold partitioning with cross-validation), with 80% of observations for each replicate used to fit the MaxEnt model (i.e., training samples) and the remaining 20% of the observations held aside to evaluate model performance (i.e., test samples). To minimize the number of features used in model construction and maximize the interpretability of individual covariate effects, we used linear feature constraints of environmental covariates in model construction [52]. However, because Rusty Blackbirds may exhibit a non-linear distribution in response to minimum winter temperatures, we included a quadratic form of the minimum temperature model covariate. To limit model over-fitting, we calibrated models by selecting the most parsimonious models for each flock size class (see Supplementary File S1 1.2).

3. Model Evaluation

We evaluated model performance by comparing model sensitivity (the true positive rate—the proportion of correctly identified samples at a given threshold of habitat suitability) and specificity (the false positive rate—the proportional predicted area for the model) for each flock size class and observational method. To do so, we assessed the area under the receiver operator curve (AUC) for test samples. AUC describes the sensitivity and specificity of a model at a given threshold of suitability and represents how well the model predicts the Rusty Blackbird observations (see Figure 2). AUC values of 0.5 represent equivalent model performance relative to random, 0.5–0.6 “poor” performance, 0.6–0.7 “fair” performance, 0.7–0.8 “good” performance, 0.8–0.9 “very good” performance, and 0.9–1.0 “excellent” model performance [53,54]. We tested whether the predictive capacity of models varied by flock size by comparing training AUC across folds for a given flock size class against a null distribution developed by permuting two flock size classes (i.e., suitability models developed by shuffling large and small flock observations). To determine whether Blitz data improved model performance, we compared observed AUC against a null distribution developed by randomizing eBird and Blitz samples.

4. Influence of Environmental Variables

To determine which environmental variables best predict Rusty Blackbird occurrence for each flock size class, we compared the influence of variable inclusion or removal on model performance. The coefficient for each variable (λ) describes the influence of that variable on suitability estimates, with positive values representing enhanced suitability and negative values representing lower suitability. To determine contribution of a variable to model fit, we evaluated jackknife estimates of the increase and decrease in AUC with the inclusion and removal of a given variable, averaged across replicate runs [55]. Additionally, we examined the difference in the occurrence of environmental covariates between flock sizes, observation protocols, and background data by permuting covariate values between classes. We compared empirical and null distributions of each environmental variable with a one-tailed test of the null hypothesis that the distributions are not statistically different (α = 0.05) using the R package permute (n = 10,000, [56]).

5. Difference in Environmental Niche by Flock Size

To assess differences in Rusty Blackbird environmental niche breadth among flock size classes, we evaluated the predicted niche breadth using ENMTools [57]. This metric follows Levins’ [29] definition of niche breadth, which describes the degree of uniformity of resource states within a distribution of individuals. Niche breadth is standardized to range between 0 and 1, where 0 represents a habitat specialist and 1 represents a habitat generalist. To evaluate our prediction that niche breadth will vary by flock size, we compared the empirical niche breadth metric for each size class with species distribution maps developed by permuting flock size class assignments. We further explored Rusty Blackbird use of niche space by assessing model prevalence, which is the MaxEnt estimate of the average probability of presence at background sites [58], and fractional predicted area of suitable habitat. We calculated fractional predicted area at the logistic threshold of equal model sensitivity and specificity. We statistically evaluated differences in prevalence and fractional predicted area between flock size classes by comparing the observed values for a given flock with a null distribution developed by permuting flock size class labels (α = 0.05).
To test whether realized ecological niches differed among flock size classes, we evaluated the suitability distributions using niche equivalency analysis and assessed the distribution of samples across individual environmental variables. Niche equivalency between flock size classes was determined by calculating the degree of similarity (modified Hellinger distances, I) between habitat suitability maps [59] and implemented in the R package phyloclim [60]. I ranges in value from 0, in which there is no overlap in the environmental niche between classes, to 1, in which the two classes have identical distributions [59]. We compared observed similarity against a null distribution developed by randomly permuting flock class labels (n = 99 permutations). We assessed statistical significance (α = 0.05) by comparing the observed and null distributions, with a one-tailed test of the null hypothesis that environmental niche models are equivalent. We compared the observed distributions of flock size classes of Rusty Blackbirds for each environmental variable using two-tail Kolmogorov–Smirnov tests of the null hypothesis that samples were drawn from the same distribution.

6. Results

The Blitz produced 678 independent checklists east of −100° longitude. eBird provided an additional 1429 traditional checklists with Rusty Blackbird observations corresponding to the same area and time periods. Limiting checklists to one per 4-km raster grid cell resulted in 495 Blitz and 714 traditional eBird checklists (Table 1). Predicted suitability maps (Figure 2) across all flock size classes suggested hotspots within the Lower Mississippi Alluvial Valley (LMAV), Black Belt Prairie (BBP), and South Atlantic Coastal Plain (SACP). For small flocks, there were additional broad areas of moderately high to high predicted suitability, whereas for large flocks the extent of high predicted suitability was more clearly limited to the aforementioned areas albeit much more restricted within the LMAV.

7. Model Performance and Environmental Correlates by Flock Size

Model performance increased with flock size, with fair performance for small flocks (AUC: 73.7 across samples) and good performance for medium and large flocks (AUCs 83.2 and 88.0 across samples, respectively; Figure 3B). Compared to models using eBird data alone, models augmented with Blitz data showed improved performance for medium and large flocks as suggested by higher AUC values (0.73 ± 0.01 and 0.72 ± 0.01 vs. 0.83 ± 0.02 and 0.88 ± 0.02, respectively; Figure 3B).
Among the environmental covariates, average minimum temperature most strongly predicted habitat suitability, contributing >60% of the models’ predictive capacity for each flock size class (Figure 4 and Figure 5A). Larger flock sizes were associated with higher average minimum temperatures (small: −1.74 °C, 95% CI [−1.74, −1.56]; medium: −0.40 °C, CI [−0.52, −0.27]; large: 0.05 °C, CI [−0.09, 0.19]), with peak densities of small, medium and large flocks occurring at −3.81, 1.91, and 1.92 °C, respectively (Figure 5A). Average minimum temperatures associated with the species’ occurrence were considerably higher than that of background points (−3.59 °C, CI [−3.62, −3.56]).
Among land cover covariates (Figure 4 and Figure 5, Table 2), proportion of floodplain forest had the greatest influence on predicted habitat suitability. The variable contribution of floodplain forest to large flock distributions was nearly twice as that of small and medium flocks. Other landcover variables that contributed to predicted habitat suitability across all flock size classes were proportion of row crops and pasture, both of which were positive influences (Table 2). Proportion of highly developed land was negatively associated with large flock distributions but was not predictive of small or medium-sized flocks. Likewise, woodland and shrub land cover classes were negatively associated with small and medium-sized flocks but were not predictive of large flocks. Surprisingly, the data provided only limited support for a positive association between habitat suitability and non-floodplain wetlands, with emergent wetlands only associated with medium-sized flock distributions and wooded wetlands for only small flocks.

8. Environmental Niche

Niche breadth decreased with increasing flock size. Predicted niche breadth for small flocks was 0.67, medium flocks were 0.47, and large flocks 0.37. Prevalence (i.e., average probability of presence at background sites) was inversely associated with flock size, with small, medium, and large flocks having predicted prevalence values of, 0.44, 0.28, and 0.20, respectively. The observed prevalence for large flocks was significantly lower than that of small flocks (p < 0.001) but statistically indistinguishable from medium-sized flocks (p = 0.773). Observed prevalence for medium flocks was significantly lower than that of small flocks (p = 0.005). Likewise, the fractional predicted area of suitable habitat increased with flock size, ranging from 33 percent for small flocks, 24 percent for medium-sized flocks, and 20 percent for large flocks. Fractional predicted area of suitable habitat was significantly lower for large and medium-sized flocks relative to small flocks (p < 0.001) but was statistically indistinguishable for large and medium flocks (p = 0.72).
The three flock size classes occupied similar but statistically distinct realized environmental niches. Across niche axes, we found significant differences between the small flock size class and the medium (i = 0.939, p < 0.001) and large flock size classes (i = 0.906, p < 0.001) but only marginal differences between the medium and large flock size classes (i = 0.979, p = 0.102). Among individual environmental variables, we found strong evidence of differences between the distribution of samples across flock size classes, most notably regarding the distribution about minimum temperatures and floodplain forests (Table 3). Whereas the modes of the temperature distributions for all three flock size classes were similar, medium and large flocks were observed in a narrower range of temperatures than small flocks (Figure 5B). Likewise, while all flock size classes were positively associated with floodplain forest (Figure 5B), the degree of association to floodplain forest increased with flock size (Figure 5C).

9. Discussion

We used novel targeted citizen science data from the Blitz and eBird to identify and describe suitable habitat for a widespread, highly vagile, and declining species, the Rusty Blackbird. Across two winters, occurrence was most strongly linked to minimum temperatures and the proportion of floodplain forest in the surrounding landscape. Although all flock sizes were associated with proportion of floodplain forest, pasture, and row crops, and the convex quadratic effect of average minimum temperature, we found considerable differences among flock size distributions. Predictably, large and medium-sized flocks were more similar in their environmental distributions relative to small flocks. Large and medium flocks had narrower environmental niches than small flocks, showing a greater preference for floodplain forest. We identified Rusty Blackbird hotspots in the LMAV, the Black Belt region of Alabama and Mississippi and the Southeast Coastal Plain. The Blitz was ultimately successful: adding Blitz data to eBird data significantly increased our predictive power.

10. Environmental Predictors of Occurrence

Rusty Blackbirds exhibited a strong association with floodplain forest, with the degree of influence of this variable positively associated with flock size (Table 3, Figure 5). Given that floodplain forest is expected to represent optimal foraging habitat for Rusty Blackbirds, this provides supportive evidence that large flocks are representative of higher quality habitats as a consequence of coarse-level local enhancement (e.g., [22,23]). Winter floodplains in the southeastern United States historically provided vast areas of shallow water crucial for migratory birds [61]. These environments, which represent seasonally flooded forested habitats, provide an abundance of invertebrates for Rusty Blackbirds, which spend much of the winter foraging in, or adjacent to, shallow water [1,62]. In the southeastern United States, floodplain forests have experienced extraordinary rates of direct landscape modification, historically due primarily to agricultural conversion and more recently as a result of urban development [9,61,63,64]. Moreover, hydrologic alteration has greatly reduced the extent and integrity of floodplain wetlands throughout the region [65]. Combined, these pressures have resulted in a 75–85% loss of the historic extent of floodplain forests in the Southeast [66,67] and the remnant patches of this habitat are highly fragmented [68,69]. Within our study extent, the LMAV hosts the largest area of floodplain forest and this region yielded the highest predicted habitat suitability for Rusty Blackbirds across flock size classes (Figure 2). Christmas Bird Count data suggest that the LMAV, considered to be the core of the species’ range, has experienced the steepest population declines and thus, maintenance of floodplain forests in this region is likely critical for Rusty Blackbird conservation [70].
Surprisingly, neither emergent nor woody wetlands had substantial influence on habitat suitability for Rusty Blackbirds. Although these land cover classes appear to provide important habitat features for the species [71], they may have been underrepresented by Blitz and eBird samples. Our ability to detect an influence of these land cover classes may have also been limited by the coarse spatial grain of our analysis. For example, though Rusty Blackbirds often forage on wet lawns in suburban areas [25] it would not be feasible to detect these habitat features at our scale of analysis, especially when using traveling count data. Sampling constraints also limited our ability to assess the effects of pecan orchards on Rusty Blackbird occurrences. Lipid-rich tree mast from pecan orchards may help the birds prepare for cold weather events [25], and wintering individuals that utilized pecan groves in Mississippi (mostly adult males), were in better body condition than those in wet forest or along creeks [24]. Despite the expected importance of this anthropogenic resource, there was not adequate representation of this land cover class in background nor Rusty Blackbird samples to include the variable in our analysis.
The influence of human land use on habitat suitability varied with flock size. Large flocks, unlike medium or small flocks, were negatively associated with high intensity development; i.e., when the proportional cover of high intensity development exceeded 70%, the probability of encountering a large flock approached zero. These results are consistent with those of Mettke-Hofmann et al. [24], Newell [10], and Newell Wohner et al. [25] who observed that smaller groups of Rusty Blackbirds can use certain agricultural and suburban areas, taking advantage of abundant pecans and earthworms, respectively. The species’ positive (though relatively weak) association with pasture and row crops is more puzzling, as it is not generally associated with these habitats (pecans are not classified as a row crop; [1]). It is possible that the low-lying, nutrient-rich floodplain landscapes that the species favors also make good croplands and pasture for farmers, and as such, both floodplain forest that Rusty Blackbirds presumably prefer, and that these agricultural habitats were often included in the 4-km grid cells we employed. Combined, our findings suggest that, though suburban and agricultural areas appear to offer resources for smaller groups of Rusty Blackbirds, larger tracts of floodplain forest will have to be maintained to support the persistence of large flocks.
Average minimum temperature was the most important variable driving habitat suitability for all flock sizes, with medium and large flocks predicted to occur at higher temperatures than small flocks. Habitat suitability estimates for medium and large flocks peaked where minimum mean temperatures occurred at means slightly higher than the freezing point of water: 1.91 °C and 1.92 °C, respectively. Suitability estimates for small flocks peaked at −3.8 °C—nearly 6 °C less than medium and large flocks. The warmer temperatures associated with medium and large flocks likely reflect the foraging behavior of Rusty Blackbirds, which spend much of their time wading in shallow water searching for aquatic prey [1]. Shallow water tends to freeze before deeper water, which would quickly prohibit birds from foraging in the substrate, making aquatic prey unavailable. Moreover, wetland invertebrate abundance is positively associated with water temperature [72,73,74]. Conversely, the colder temperatures associated with small flocks may represent marginal habitats in which birds utilize terrestrial resources in freezing temperatures (see [24,25]).
Because temperatures were averaged across years of the study, we were unable to assess the annual occurrences of Rusty Blackbirds. The species varies widely in its wintering distribution [9] and site fidelity [10], most likely because of differences in winter temperatures and water levels across years. Rusty Blackbirds’ association with temperatures just above freezing supports our theory that the species, which is a facultative migrant not tightly linked to day length [1], migrates south from the Boreal forest to avoid freezing water. Birds may then spend many weeks in the northern United States during stopover periods in fall [75,76], only to move south again when frozen foraging habitat pushes them further. Future work should assess whether carryover effects due to variable wintering conditions impact the birds across their full annual cycle [77,78], including on their breeding grounds in the Boreal forest.

11. Larger Flocks Have a Narrower Niche Breadth

We found strong evidence that large and medium flocks have narrower environmental niches than small flocks, as measured by Levins’ [29] niche breadth metric and, consequently, prevalence and fractional predicted area of suitable habitat was higher for small flocks than medium and large flock size classes. Whereas such differences in niche breadth have been observed across a variety of taxa using environmental niche modeling (e.g., [79]), this represents among precious few examples of the relationship between flock size and niche [26]. This pattern is abundantly clear across large scales, as small flocks have a much broader geographic range of suitable habitat than large and medium flocks. Conversely, the area of suitable habitat (Figure 2) for large flocks is much smaller and is largely constrained within the center of the species range. Again, this highlights the relative vulnerability of larger flocks of Rusty Blackbirds and the relative flexibility of smaller flocks. Small flocks may be able to persist in heterogenous landscapes with small patches of ephemeral resources (e.g., fruiting pecan orchards, lawns with surfacing earthworms), whereas larger flocks may be limited to the highest quality habitats. Future studies (e.g., using tracking technology) should examine space use of large vs. small flocks to understand how medium (i.e., <4 km) and small-scale habitat features affect fitness and niche space.
The use of flock size herein to approximate habitat quality represents a limitation to the inference of our results. Blitz and traditional eBird observations were traveling counts, thus in some instances, birds observed on a given count could have comprised one or several true flocks. As such, it was not technically possible to distinguish between a single large flock or a higher density of smaller flocks; however, given our own observations of the species (e.g., [12]), we believe the vast majority of the counts probably represented single flocks, so we are comfortable interpreting results based on this assumption. Moreover, traveling counts limit the spatial grain of the analysis, as counts may occur at any point along the observer’s route. To address this, we highly recommend the use of stationary count data to allow for assessing the distribution of flocks as it relates to land cover data (e.g., evaluating resource dispersion hypothesis as described by [80]). Even without these limitations, density itself is not a clear indicator of habitat quality [11,81]. Further research (e.g., survival) is therefore needed to confirm that large flocks are representative of environments of higher quality habitat [33]. Sex or age-related segregation of the observed flocks may provide additional evidence of habitat quality [82]. Dominance hierarchies can indicate habitat quality in a number of bird species (e.g., the American Redstart, Setophaga ruticilla, [83]). For example, in Rusty Blackbirds, Mettke-Hofmann et al. [24] observed that adult males maintained a higher body condition and occupied habitats with higher pecan mast production than females and immature birds. Unfortunately, age and sex data from the Blitz were not adequate to address this issue in the current study.
Ultimately, a demographic response, most appropriately overwinter survival, is necessary to truly ascertain whether flock size is an appropriate indication of winter habitat quality for Rusty Blackbirds [84]. In particular, despite some location-specific evidence hinting at population limitation on the breeding grounds in the Boreal forest [7,8], a full-annual cycle population model suggests that wintering juvenile survival is the demographic parameter most tied to the species’ rate of population change (Clark Rushing, Steve Matsuoka, Luke L. Powell et al. unpublished analysis). Whereas estimating overwinter survival may be very challenging given the species’ low site fidelity, alternative proximate habitat quality metrics that can be utilized in combination for Rusty Blackbird flocks include but are not limited to: departure time (e.g., [85]), body condition (e.g., [86,87], but see [84,88]), and telomere length (e.g., [89]).

12. Regional Hotspots and Conservation

The LMAV was a particularly suitable area for all flock sizes (Figure 2, [5]). Though we also found a broad swath across the southeast coastal plain that was generally suitable, with small hotspots in north coastal South Carolina [25], the majority of the East Coast of the United States appears less suitable than the LMAV. Stable isotope data [90] and light level geolocators [75] indicate that Rusty Blackbirds wintering in the LMAV primarily migrate to Alaska and western Canada—the Western Boreal—whereas those from the southeast coastal plain breed in northeast North America—the Southeastern Boreal. It remains unclear whether flock sizes in the southeast coastal plain are constrained by anthropogenic disturbance or if it is the quality and extent of floodplain forests of the LMAV that leads to larger flock sizes in that region.
A novel finding from our study was the emergence of the Black Belt region (Ecoregion 65, [91,92,93]) as a hotspot across all Rusty Blackbird flock size classes. The Black Belt, which arcs along the boundary of the piedmont and coastal plain in Mississippi and Alabama, was historically a mosaic of prairie and forested habitat [94]. Whereas most of the region’s native vegetation was cleared for cotton-based agriculture [94,95,96,97], the Black Belt remains embedded within a vast matrix of floodplain forest. To date, we know of no Rusty Blackbird research projects within 500 km of the Black Belt—thus, our results highlight a need for targeted field studies to further qualify the use and importance of this region to Rusty Blackbirds.
The emergence of the Black Belt region as a hotspot underscores the benefits and limitations of species distribution modeling. Whereas comparably few eBird and Blitz samples were recorded within the region, the MaxEnt model, as implemented here, is carried out in environmental space [52], which allowed us to generate predictions for the Black Belt as a function of its environmental characteristics [98]. As a cautionary note, however, this approach is limited in its ability to directly inform conservation efforts, as it yields information about the type of environment an organism uses rather than occurrences per se [99]. Moreover, our models are undoubtedly imperfect—for example, some areas of high predicted suitability for large flocks were of low predicted suitability for medium-sized flocks (e.g., north-central LMAV). This was likely driven by the difference in habitat preferences of large vs. medium flocks (e.g., large flocks had higher preference for floodplain forest, Table 3) as well as by stochasticity in where birds were detected. As such, we suggest that systematic sampling of predicted Rusty Blackbird hotspots—particularly in the Black Belt region—is a crucial step towards linking our results with management efforts.

13. Did the Blitz Help Relative to eBird Alone?

The Blitz significantly improved the predictive power of spatial models compared to using eBird alone (Figure 3), highlighting the value of this novel citizen science approach to improve understanding of a wide-ranging species of conservation concern. Further, the addition of the Blitz especially improved our predictive power for large flocks (Figure 3B), as we effectively doubled the number of large flocks detected (Table 1). This is important, as it allows us to precisely concentrate targeted studies on sites at which we predict only the very highest flock sizes of Rusty Blackbirds. We also utilized a somewhat novel way to address biases associated with modeling eBird data by using background samples from eBird lists across taxa: our background data were generated with equivalent bias to our occurrence data. Our comparison of background and occurrence data therefore provides conservative estimates of Rusty Blackbirds in environmental space. An important caveat to our Blitz effort is that participants were searching for Rusty Blackbirds in habitats that were expected to be suitable for the species–this has the potential to bias suitability predictions towards environments that are pre-determined to be suitable. A preliminary analysis of our results, however, produced estimates from eBird and Blitz data that were similar, suggesting that any such bias was nominal. In addition to the quantitative benefits of the Blitz, it was clear the Blitz raised considerable awareness for the plight of the Rusty Blackbird (e.g., Audubon magazine, eBird website, conservation partners participating in the Blitz, etc.). Given that citizen scientists submitted 678 checklists specifically for the Blitz, the program was clearly successful in engaging birders with ornithological research and conservation concerns to which they may not previously have been exposed.
Taken together, the data provided by citizen science participants of the Blitz and eBird program have enhanced our understanding of the occurrence of this Boreal-breeding species throughout its wintering areas. Whereas wintering Rusty Blackbirds can be found across most of the southeastern United States, only a narrow subset of the region was found to be suitable for large flocks of the species. By targeting these regions for future research, we can maximize sampling efficiency by searching areas predicted to provide high quality habitat. Targeted research is crucial given the limitations associated with presence-only observations. Using the MaxEnt framework, with presence-only data, locations where birds were not observed (i.e., background data) cannot be treated as true absences (see [100,101]). This is especially problematic when detection probability varies and cannot be incorporated into the model [102] and imparts a limitation to how citizen science data can be used to inform conservation in this context [103]. The use of presence-absence data to model occupancy (e.g., [2]), can also allow researchers to address temporal variation in the presence of Rusty Blackbird flocks, which may yield further insight into habitat quality (e.g., abundance-occupancy relationships, [104]). Moreover, identifying and counting individual flocks from fixed points in the landscape and evaluating these observations at multiple spatial scales would provide researchers with the ability to assess alternative explanations for flock size distributions, such as the resource dispersion hypothesis (reviewed in [80]). Future work that utilizes systematic sampling to estimate the geographic distribution of Rusty Blackbird flocks and assess the relationship between flock size and habitat quality is therefore a critical next step towards conserving this imperiled species.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1424-2818/13/2/62/s1: Table S1: rubl_gap_reclassification.csv, Supplementary File S1: Additional methods and results, Supplementary File S2: rubl_distribution_maps.zip.

Author Contributions

Conceptualization, R.S.G.; methodology, B.S.E.; formal analysis, B.S.E.; data curation, B.S.E.; writing—original draft preparation, B.S.E., L.L.P., D.W.D., and S.M.B.; supervision, R.S.G.; project administration, R.S.G.; funding acquisition, D.W.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the associated APC were supported by the U.S. Fish & Wildlife Service/USFWS.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: https://github.com/SMBC-NZP/rubl.

Acknowledgments

Russ Greenberg conceived of, advertised, and led the execution of the Blitz; our hope is that this manuscript honors the memory of a brilliant ornithologist. This paper was written in association with the International Rusty Blackbird Working Group (RustyBlackbird.org; founded and led by Russ). We thank the United States Fish and Wildlife Service for funding. We deeply thank the regional coordinators and many citizen scientists that contributed to eBird and to the Blitz, and we thank the Cornell Lab of Ornithology and eBird for creating a portal for Blitz data entry. The findings and conclusions in this article are those of the author(s) and do not necessarily represent the views of the U.S. Fish and Wildlife Service. Any mention of trade names is purely coincidental and does not represent endorsement by the author(s) and / or their respective organization.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Avery, M.L. Rusty Blackbird (Euphagus Carolinus), version 1.0. In Birds of the World; Poole, A.F., Ed.; Cornell Lab of Ornithology: Ithaca, NY, USA, 2020; Available online: https://birdsoftheworld.org/bow/species/rusbla/cur/introduction (accessed on 10 February 2018).
  2. Luscier, J.D.; Lehnen, S.E.; Smith, K.G. Habitat Occupancy by Rusty Blackbirds Wintering in the Lower Mississippi Alluvial Valley. Condor 2010, 112, 841–848. [Google Scholar] [CrossRef]
  3. Powell, L.L.; Hodgman, T.P.; Fiske, I.J.; Glanz, W.E. Habitat occupancy of Rusty Blackbirds (Euphagus carolinus) breeding in northern New England, USA. Condor 2014, 116, 122–133. [Google Scholar] [CrossRef]
  4. Greenberg, R.; Droege, S. On the Decline of the Rusty Blackbird and the Use of Ornithological Literature to Document Long-Term Population Trends. Conserv. Biol. 1999, 13, 553–559. [Google Scholar] [CrossRef] [Green Version]
  5. Greenberg, R.; Demarest, D.W.; Matsuoka, S.M.; Mettke-Hofmann, C.; Evers, D.; Hamel, P.B.; Luscier, J.; Powell, L.L.; Shaw, D.; Avery, M.L.; et al. Understanding Declines in Rusty Blackbirds. Boreal Birds North Am. 2011, 41, 107–126. [Google Scholar] [CrossRef]
  6. Sauer, J.R. The North American Breeding Bird Survey, Results and Analysis 1966–2015. Version 1.30.2015. USGS Patuxent Wildlife Research Center: Laurel, MD 2015. Available online: http://www.mbr-pwrc.usgs.gov.bbs/ (accessed on 15 August 2016).
  7. Powell, L.L.; Hodgman, T.P.; Glanz, W.E.; Osenton, J.D.; Fisher, C.M. Nest-Site Selection and Nest Survival of the Rusty Blackbird: Does Timber Management Adjacent to Wetlands Create Ecological Traps? Condor 2010, 112, 800–809. [Google Scholar] [CrossRef]
  8. Edmonds, S.T.; Evers, D.C.; Cristol, D.A.; Mettke-Hofmann, C.; Powell, L.L.; McGann, A.J.; Armiger, J.W.; Lane, O.P.; Tessler, D.F.; Newell, P.; et al. Geographic and Seasonal Variation in Mercury Exposure of the Declining Rusty Blackbird. Condor 2010, 112, 789–799. [Google Scholar] [CrossRef]
  9. Greenberg, R.; Matsuoka, S.M. Special section: Rangewide ecology of the declining Rusty Blackbird. Rusty Blackbird: Mysteries of a species in decline. Condor 2010, 112, 770–777. [Google Scholar] [CrossRef]
  10. Newell, P.J. Winter Ecology of the Rusty Blackbird (Euphagus Carolinus). Ph.D. Thesis, University of Georgia, Athens, GA, USA, 2013. [Google Scholar]
  11. Van Horne, B. Density as a Misleading Indicator of Habitat Quality. J. Wildl. Manag. 1983, 47, 893–901. [Google Scholar] [CrossRef]
  12. Borchert, S.M. Site-Specific Habitat and Landscape Associations of Rusty Blackbirds Wintering in Louisiana. Master’s Thesis, Louisiana State University, Baton Rouge, LA, USA, 2015. [Google Scholar]
  13. eBird Basic Dataset. Version: EBD_relNov-2017; Cornell Lab of Ornithology: Ithaca, NY, USA, 2017; p. 8. [Google Scholar]
  14. Beauchamp, G.; Belisle, M.; Giraldeau, L.A. Influence of conspecific attraction on the spatial distribution of learning foragers in a patchy habitat. J. Anim. Ecology 1997, 66, 671–682. [Google Scholar] [CrossRef]
  15. Pöysä, H. Group Foraging in Patchy Environments: The Importance of Coarse-Level Local Enhancement. Ornis Scand. 1992, 23, 159. [Google Scholar] [CrossRef]
  16. Stolen, E.D.; Collazo, J.A.; Percival, H.F. Group-Foraging Effects on Capture Rate in Wading Birds. Condor 2012, 114, 744–754. [Google Scholar] [CrossRef]
  17. Thiebault, A.; Mullers, R.H.; Pistorius, P.A.; Tremblay, Y. Local enhancement in a seabird: Reaction distances and foraging consequence of predator aggregations. Behav. Ecology 2014, 25, 1302–1310. [Google Scholar] [CrossRef]
  18. Ward, P.; Zahavi, A. The importance of certain assemblages of birds as “information-centres” for food-finding. Ibis 1973, 115, 517–534. [Google Scholar] [CrossRef]
  19. Pride, R.E. Optimal group size and seasonal stress in ring-tailed lemurs (Lemur catta). Behav. Ecol. 2005, 16, 550–560. [Google Scholar] [CrossRef] [Green Version]
  20. Hoare, D.; Couzin, I.; Godin, J.-G.; Krause, J. Context-dependent group size choice in fish. Anim. Behav. 2004, 67, 155–164. [Google Scholar] [CrossRef]
  21. Caraco, T. Time Budgeting and Group Size: A Theory. Ecology 1979, 60, 611–617. [Google Scholar] [CrossRef]
  22. Beerens, J.M.; Noonburg, E.G.; Gawlik, D.E. Linking dynamic habitat selection with wading bird foraging distributions across resource gradients. PLoS ONE 2015, 10, e0128182. [Google Scholar] [CrossRef]
  23. Gopi Sundar, K.S. Flock Size, Density and Habitat Selection of Four Large Waterbirds Species in an Agricultural Landscape in Uttar Pradesh, India: Implications for Management. Waterbirds 2006, 29, 365–374. [Google Scholar]
  24. Mettke-Hofmann, C.; Hamel, P.B.; Hofmann, G.; Zenzal, T.J., Jr.; Pellegrini, A.; Malpass, J. Competition and Habitat Quality Influence Age and Sex Distribution in Wintering Rusty Blackbirds. PLoS ONE 2015, 10, e0123775. [Google Scholar] [CrossRef]
  25. Wohner, P.J.; Cooper, R.J.; Greenberg, R.S.; Schweitzer, S.H. Weather affects diet composition of rusty blackbirds wintering in suburban landscapes. J. Wildl. Manag. 2015, 80, 91–100. [Google Scholar] [CrossRef]
  26. Arevalo, J.E.; Gosler, A.G. The behaviour of Treecreepers Certhia familiaris in mixed-species flocks in winter. Bird Study 1994, 41, 1–6. [Google Scholar] [CrossRef]
  27. Colwell, R.K.; Futuyma, D.J. On the Measurement of Niche Breadth and Overlap. Ecology 1971, 52, 567–576. [Google Scholar] [CrossRef]
  28. Hutchinson, G.E. Concluding remarks. Cold Spring Harb. Symp. Quant. Biol. 1957, 22, 415–427. [Google Scholar] [CrossRef]
  29. Levins, R. Evolution in Changing Environments; Princeton University Press: Princeton, NJ, USA, 1968. [Google Scholar]
  30. Sullivan, B.L.; Wood, C.L.; Iliff, M.J.; Bonney, R.E.; Fink, D.; Kelling, S. eBird: A citizen-based bird observation network in the biological sciences. Biol. Conserv. 2009, 142, 2282–2292. [Google Scholar] [CrossRef]
  31. Evans, B.S. Summary Report: Assessing Rusty Blackbird Habitat Suitability on Wintering Grounds and during Spring Migration Using A Large Citizen-Science Dataset; Smithsonian Migratory Bird Center: Washington, DC, USA, 2016; Available online: http://rustyblackbird.org/wp-content/uploads/RUBL_Blitz10pgReportEvans23jan2017.pdf (accessed on 10 February 2018).
  32. Hirzel, A.H.; Le Lay, G. Habitat suitability modelling and niche theory. J. Appl. Ecology 2008, 45, 1372–1381. [Google Scholar] [CrossRef]
  33. Kellner, C.J.; Brawn, J.D.; Karr, J.R. What Is Habitat Suitability and How Should It Be Measured? In Wildlife 2001: Populations; McCullough, D.R., Barrett, R.H., Eds.; Elsevier Applied Science: New York, NY, USA, 1992; pp. 476–488. [Google Scholar]
  34. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. Ecology Model. 2006, 190, 231–259. [Google Scholar] [CrossRef] [Green Version]
  35. Sullivan, B.L.; Aycrigg, J.L.; Barry, J.H.; Bonney, R.E.; Bruns, N.; Cooper, C.B.; Damoulas, T.; Dhondt, A.A.; Dietterich, T.; Farnsworth, A.; et al. The eBird enterprise: An integrated approach to development and application of citizen science. Biol. Conserv. 2014, 169, 31–40. [Google Scholar] [CrossRef]
  36. Rowcliffe, M. Activity: Animal Activity Statistics. R Package Version 1.3. 2019. Available online: https://CRAN.R-project.org/package=activity (accessed on 15 July 2019).
  37. Nouvellet, P.; Rasmussen, G.S.A.; Macdonald, D.W.; Courchamp, F. Noisy clocks and silent sunrises: Measurement methods of daily activity pattern. J. Zool. 2011, 286, 179–184. [Google Scholar] [CrossRef] [Green Version]
  38. Vazquez, C.; Rowcliffe, J.M.; Spoelstra, K.; Jansen, P.A. Comparing diel activity patterns of wildlife across latitudes and seasons: Time transformations using day length. Methods Ecology Evol. 2019, 10, 2057–2066. [Google Scholar] [CrossRef] [Green Version]
  39. Ridout, M.S.; Linkie, M. Estimating overlap of daily activity patterns from camera trap data. J. Agric. Biol. Env. Stat. 2009, 14, 322–337. [Google Scholar] [CrossRef]
  40. Dickinson, J.L.; Zuckerberg, B.; Bonter, D.N. Citizen science as an ecological research tool: Challenges and benefits. Annu. Rev. Ecology Evol. Syst. 2010, 41, 149–172. [Google Scholar] [CrossRef] [Green Version]
  41. Hortal, J.; Jiménez-Valverde, A.; Gómez, J.F.; Lobo, J.M.; Baselga, A. Historical bias in biodiversity inventories affects the observed environmental niche of the species. Oikos 2008, 117, 847–858. [Google Scholar] [CrossRef]
  42. Newbold, T.; Reader, T.; El-Gabbas, A.; Berg, W.; Shohdi, W.M.; Zalat, S.; El Din, S.B.; Gilbert, F. Testing the accuracy of species distribution models using species records from a new field survey. Oikos 2010, 119, 1326–1334. [Google Scholar] [CrossRef]
  43. Dudík, M.; Phillips, S.J.; Schapire, R.E. Correcting sample selection bias in maximum entropy density estimation. Adv. Neural Inf. Process. Syst. 2006, 18, 323–330. [Google Scholar]
  44. Elith, J.; Phillips, S.J.; Hastie, T.; Dudík, M.; Chee, Y.E.; Yates, C.J. A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 2010, 17, 43–57. [Google Scholar] [CrossRef]
  45. Phillips, S.J.; Dudík, M.; Elith, J.; Graham, C.H.; Lehmann, A.; Leathwick, J.; Ferrier, S. Sample selection bias and occurrence-only distribution models: Implications for background and pseudo-absence data. Ecology Appl. 2009, 19, 181–197. [Google Scholar] [CrossRef] [Green Version]
  46. Anderson, R.P.; Raza, A. The effect of the extent of the study region on GIS models of species geographic distributions and estimates of niche evolution: Preliminary tests with montane rodents (genus Nephelomys) in Venezuela. J. Biogeogr. 2010, 37, 1378–1393. [Google Scholar] [CrossRef]
  47. Hijmans, R.J.; van Etten, J.; Cheng, J.; Mattiuzzi, M.; Sumner, M.; Greenberg, J.A.; Lamigueiro, O.P.; Bevan, A.; Racine, E.B.; Shortridge, A.; et al. Raster: Geographic Data Analysis and Modeling. In R Package Version 2.8-19. 2015. Available online: https://CRAN.R-project.org/package=raster (accessed on 27 September 2020).
  48. R Core Team 2013. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013; Available online: http://www.r-project.org/ (accessed on 10 March 2012).
  49. PRISM Climate Group. Oregon State University. Available online: http://prism.oregonstate.edu (accessed on 1 June 2013).
  50. Scott, J.M.; Davis, F.; Csuti, B.; Noss, R.; Butterfield, B.; Groves, C.; Anderson, H.; Caicco, S.; D’Erchia, F.; Edwards, T.C., Jr.; et al. Gap Analysis: A Geographic Approach to Protection of Biological Diversity. In Wildlife Monographs; John Wiley & Sons: Hoboken, NJ, USA, 1993; pp. 3–41. [Google Scholar]
  51. Hijmans, R.J.; Phillips, S.; Leathwick, J.; Elith, J.; Hijmans, M.R.J. Package ‘Dismo’. Circles 2017, 9, 1–68. [Google Scholar]
  52. Merow, C.; Smith, M.J.; Silander, J.A. A practical guide to MaxEnt for modeling species’ distributions: What it does, and why inputs and settings matter. Ecography 2013, 36, 1058–1069. [Google Scholar] [CrossRef]
  53. Duan, R.Y.; Kong, X.Q.; Huang, M.Y.; Fan, W.Y.; Wang, Z.G. The predictive performance and stability of six species distribution models. PLoS ONE 2014, 9, e112764. [Google Scholar] [CrossRef] [Green Version]
  54. Swets, J.A. Measuring the accuracy of diagnostic systems. Science 1988, 240, 1285–1293. [Google Scholar] [CrossRef] [Green Version]
  55. Shcheglovitova, M.; Anderson, R.P. Estimating optimal complexity for ecological niche models: A jackknife approach for species with small sample sizes. Ecology Model. 2013, 269, 9–17. [Google Scholar] [CrossRef]
  56. Simpson, G.L. Permute: Functions for Generating Restricted Permutations of Data. R Package Version 0.9-5. 2019. Available online: https://CRAN.R-project.org/package=permute (accessed on 10 March 2012).
  57. Warren, D.L.; Glor, R.E.; Turelli, M. ENMTools: A toolbox for comparative studies of environmental niche models. Ecography 2010, 33, 607–611. [Google Scholar] [CrossRef]
  58. Phillips, S.J.; Elith, J. On estimating probability of presence from use-availability or presence-background data. Ecology 2013, 94, 1409–1419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Warren, D.L.; Glor, R.E.; Turelli, M. Environmental niche equivalency versus conservatism: Quantitative approaches to niche evolution. Evolution 2008, 62, 2868–2883. [Google Scholar] [CrossRef] [PubMed]
  60. Heibl, C.; Calenge, C. Phyloclim: Integrating Phylogenetics and Climatic Niche Modeling. In R Package Version 0.9-4. 2013. Available online: https://CRAN.R-project.org/package=phyloclim (accessed on 27 September 2020).
  61. Brown, C.R.; Baxter, C.; Pashley, D.N. The Ecological Basis for the Conservation of Migratory Birds in the Mississippi Alluvial Valley. In Strategies for Bird Conservation: The Partners in Flight Planning Process; Bonney, R., Pashley, D.N., Cooper, R.J., Niles, L., Eds.; Cornell Lab of Ornithology: Ithaca, NY, USA, 1999; pp. 3–6. [Google Scholar]
  62. Wehrle, B.W.; Kaminski, R.M.; Leopold, B.D.; Smith, W.P. Aquatic invertebrate resources in Mississippi forested wetlands during winter. Wildlife Society Bulletin 1995, 23, 774–783. [Google Scholar]
  63. Faulkner, S.P. Urbanization impacts on the structure and function of forested wetlands. Urban Ecosyst. 2004, 7, 89–106. [Google Scholar] [CrossRef]
  64. Lockaby, B.G. Floodplain ecosystems of the Southeast: Linkages between forests and people. Wetl 2009, 29, 407–412. [Google Scholar] [CrossRef]
  65. Harmar, O.P.; Clifford, N.J.; Thorne, C.R.; Biedenharn, D.S. Morphological changes of the Lower Mississippi River: Geomorphological response to engineering intervention. River Res. Appl. 2005, 21, 1107–1131. [Google Scholar] [CrossRef]
  66. Hefner, J.M.; Brown, J.P. Wetland trends in the southeastern U.S. Wetlands 1984, 4, 1–11. [Google Scholar] [CrossRef]
  67. Hefner, J.M.; Wilen, B.O.; Dahl, T.E.; Frayer, W.E. Southeastern Wetlands: Status and Trends, Mid-1970s to Mid-1980s; U.S. Fish and Wildlife Service and U.S. Environmental Protection Agency: Atlanta, GA, USA, 1994.
  68. Rudis, V.A. Regional forest fragmentation effects on bottomland hardwood community types and resource values. Landsc. Ecology 1995, 10, 291–307. [Google Scholar] [CrossRef]
  69. Twedt, D.J.; Loesch, C.R. Forest area and distribution in the Mississippi alluvial valley: Implications for breeding bird conservation. J. Biogeogr. 1999, 26, 1215–1224. [Google Scholar] [CrossRef]
  70. Niven, D.K.; Sauer, J.R.; Butcher, G.S. Christmas bird count provides insights into population change in land birds that breed in the boreal forest. Am. Birds 2004, 58, 10–20. [Google Scholar]
  71. DeLeon, E.E. Ecology of Rusty Blackbirds Wintering in Louisiana: Seasonal Trends, Flock Compositions and Habitat Associations. Master’s Thesis, Louisiana State University, Baton Rouge, LA, USA, 2012. [Google Scholar]
  72. Batema, D.L.; Kaminski, R.M.; Magee, P.A. Wetland invertebrate communities and management of hardwood bottomlands in the Mississippi Alluvial Valley. In Ecology and Management of bottomland Hardwood Systems: The State of our Understanding; Fredrickson, L.H., King, S.L., Kaminski, R.M., Eds.; Gaylord Memorial Laboratory Special Publication 10, University of Missouri-Columbia: Puxico, MO, USA, 2005; pp. 173–190. [Google Scholar]
  73. Fredrickson, L.H.; Batema, D.L. Greentree Reservoir Management Handbook; Gaylord Memorial Lab, University of Missouri: Puxico, MO, USA, 1992. [Google Scholar]
  74. White, D.C. Lowland hardwood wetland invertebrate community and production in Missouri. Arch. Hydrobiol. 1985, 103, 509–533. [Google Scholar]
  75. Johnson, J.A.; Matsuoka, S.M.; Tessler, D.F.; Greenberg, R.; Fox, J.W. Identifying Migratory Pathways Used by Rusty Blackbirds Breeding in Southcentral Alaska. Wilson J. Ornithol. 2012, 124, 698–703. [Google Scholar] [CrossRef]
  76. Wright, J.R.; Powell, L.L.; Tonra, C.M. Automated telemetry reveals staging behavior in a declining migratory passerine. Auk 2018, 135, 461–476. [Google Scholar] [CrossRef] [Green Version]
  77. Harrison, X.A.; Blount, J.D.; Inger, R.; Norris, D.R.; Bearhop, S. Carry-over effects as drivers of fitness differences in animals. J. Anim. Ecology 2010, 80, 4–18. [Google Scholar] [CrossRef]
  78. Marra, P.P.; Cohen, E.B.; Loss, S.R.; Rutter, J.E.; Tonra, C.M. A call for full annual cycle research in animal ecology. Biol. Lett. 2015, 11, 20150552. [Google Scholar] [CrossRef]
  79. Saupe, E.E.; Qiao, J.; Hendricks, J.R.; Portell, R.W.; Junter, S.J.; Sobseron, J.; Lieberman, B.S. Niche breadth and geographic range size as determinants of species survival on geological time scales. Glob. Ecology Biogeogr. 2015, 24, 1159–1169. [Google Scholar] [CrossRef]
  80. Johnson, D.D.; Kays, R.; Blackwell, P.G.; Macdonald, D.W. Does the resource dispersion hypothesis explain group living? Trends Ecology Evol. 2002, 17, 563–570. [Google Scholar] [CrossRef]
  81. Battin, J. When Good Animals Love Bad Habitats: Ecological Traps and the Conservation of Animal Populations. Conserv. Biol. 2004, 18, 1482–1491. [Google Scholar] [CrossRef]
  82. Kawecki, T.J. Adaptation to Marginal Habitats. Annu. Rev. Ecology Evol. Syst. 2008, 39, 321–342. [Google Scholar] [CrossRef] [Green Version]
  83. Powell, L.L.; Ames, E.M.; Wright, J.R.; Matthiopoulos, J.; Marra, P.P. Interspecific competition between resident and wintering birds: Experimental evidence and consequences of coexistence. Ecology 2020, in press. [Google Scholar] [CrossRef]
  84. Johnson, M.D. Measuring habitat quality: A review. Condor 2007, 109, 489–504. [Google Scholar] [CrossRef]
  85. Studds, C.E.; Marra, P.P. Rainfall-induced changes in food availability modify the spring departure programme of a migratory bird. Proc. R. Soc. B Biol. Sci. 2011, 278, 3437–3443. [Google Scholar] [CrossRef]
  86. Bearhop, S.; Hilton, G.M.; Votier, S.C.; Waldron, S. Stable isotope ratios indicate that body condition in migrating passerines is influenced by winter habitat. Proc. R. Soc. Lond. Ser. B 2004, 271, 215–218. [Google Scholar] [CrossRef] [Green Version]
  87. Latta, S.C.; Faaborg, J. Demographic and population responses of Cape May Warblers wintering in multiple habitats. Ecology 2002, 83, 2502–2515. [Google Scholar] [CrossRef]
  88. Rintamäki, P.T.; Stone, J.R.; Lundberg, A.; Moore, F. Seasonal and diurnal body-mass fluctuations for two nonhoarding species of Parus in Sweden modeled using path analysis. Auk 2003, 120, 658–668. [Google Scholar] [CrossRef]
  89. Young, R.C.; Kitaysky, A.S.; Barger, C.P.; Dorresteijn, I.; Ito, M.; Watanuki, Y. Telomere length is a strong predictor of foraging behavior in a long-lived seabird. Ecosphere 2015, 6, 1–26. [Google Scholar] [CrossRef]
  90. Hobson, K.A.; Greenberg, R.; Van Wilgenburg, S.L.; Mettke-Hofmann, C. Migratory Connectivity in the Rusty Blackbird: Isotopic Evidence From Feathers of Historical and Contemporary Specimens. Condor 2010, 112, 778–788. [Google Scholar] [CrossRef] [Green Version]
  91. Omernik, J.M. Ecoregions: A spatial framework for environmental management. In Biological Assessment and Criteria: Tools for Water Resource Planning and Decision Making; Davis, W.S., Simon, T.P., Eds.; Lewis Publishers: Boca Raton, FL, USA, 1995; pp. 49–62. [Google Scholar]
  92. Omernik, J.M. Perspectives on the Nature and Definition of Ecological Regions. Environ. Manag. 2004, 34, S27–S38. [Google Scholar] [CrossRef]
  93. Omernik, J.M.; Griffith, G.E. Ecoregions of the Conterminous United States: Evolution of a Hierarchical Spatial Framework. Env. Manag. 2014, 54, 1249–1266. [Google Scholar] [CrossRef]
  94. Webster, G.R.; Samson, S.A. On Defining the Alabama Black Belt: Historical Changes and Variations. Southeast. Geogr. 1992, 32, 163–172. [Google Scholar] [CrossRef]
  95. Barone, J.A. Historical occurrence and distribution of prairies in the Black Belt of Mississippi and Alabama. Castanea 2005, 70, 170–183. [Google Scholar] [CrossRef]
  96. Cleland, H.F. Black Belt of Alabama. Geogr. Rev. 1920, 10, 375–387. [Google Scholar] [CrossRef]
  97. Wilson, T.H. Natural history of the black belt prairie. J. Ala. Acad. Sci. 1981, 52, 10–19. [Google Scholar]
  98. Phillips, S.J.; Dudík, M.; Schapire, R.E. A Maximum Entropy Approach to Species Distribution Modeling. In Proceedings of the Twenty-First International Conference on Machine Learning, Banff, AB, Canada, 4–8 July 2004; ACM: New York, NY, USA, 2004. [Google Scholar]
  99. Elith, J.; Graham, C.H.; Anderson, R.P.; Dud’ik, M.; Ferrier, S.; Guisan, A.; Hijmans, R.J.; Huettmann, F.; Leathwick, J.R.; Lehmann, A. Novel methods improve prediction of species’ distributions from occurrence data. Ecography 2006, 29, 129–151. [Google Scholar] [CrossRef] [Green Version]
  100. Guillera-Arroita, G.; Lahoz-Monfort, J.J.; Elith, J. Maxent is not a presence–absence method: A comment on Thibaud et al. Methods Ecology Evol. 2014, 11, 1192–1197. [Google Scholar] [CrossRef]
  101. Li, W.; Guo, Q.; Elkan, C. Can we model the probability of presence of species without absence data? Ecography 2011, 34, 1096–1105. [Google Scholar] [CrossRef]
  102. Lahoz-Monfort, J.J.; Guillera-Arroita, G.; Wintle, B.A. Imperfect detection impacts the performance of species distribution models. Glob. Ecology Biogeogr. 2014, 23, 504–515. [Google Scholar] [CrossRef]
  103. Guillera-Arroita, G.; Lahoz-Monfort, J.J.; Elith, J.; Gordon, A.; Kujala, H.; Lentini, P.E.; McCarthy, M.A.; Tingley, R.; Wintle, B.A. Is my species distribution model fit for purpose? Matching data and models to applications. Glob. Ecology Biogeogr. 2015, 24, 276–292. [Google Scholar] [CrossRef]
  104. Gaston, K.J.; Blackburn, T.M.; Greenwood, J.J.D.; Gregory, R.D.; Quinn, R.M.; Lawton, J.H. Abundance-occupancy relationships. J. Appl. Ecology 2000, 37, 39–59. [Google Scholar] [CrossRef]
Figure 1. Processing steps from classified land cover layer to distribution modeling.
Figure 1. Processing steps from classified land cover layer to distribution modeling.
Diversity 13 00062 g001
Figure 2. Habitat suitability estimates for small, medium, and large flocks of Rusty Blackbirds (<20, 20–99, and >99 individuals, respectively). Across all flock size classes, maps suggest focal wintering regions within the Lower Mississippi Alluvial Valley, Black Belt Prairie, and South Atlantic Coastal Plain. These are also available in Supplementary File S2 as spatially-referenced kmz files suitable for Google Earth.
Figure 2. Habitat suitability estimates for small, medium, and large flocks of Rusty Blackbirds (<20, 20–99, and >99 individuals, respectively). Across all flock size classes, maps suggest focal wintering regions within the Lower Mississippi Alluvial Valley, Black Belt Prairie, and South Atlantic Coastal Plain. These are also available in Supplementary File S2 as spatially-referenced kmz files suitable for Google Earth.
Diversity 13 00062 g002
Figure 3. Receiver operator curves (A) and Area Under the curve (AUC; B) for small, medium, and large flocks of Rusty Blackbirds using Blitz data, eBird data alone, or Blitz and eBird data (“combined”). Solid lines (A) represent receiver operator curves for replicate runs across sampling methods and flock size classes. AUC values of 0.5 represent equivalent MaxEnt model performance relative to random, 0.5–0.6 “poor” performance, 0.6–0.7 “fair”, 0.7–0.8 “good”, 0.8–0.9 “very good”, and 0.9–1.0 “excellent”.
Figure 3. Receiver operator curves (A) and Area Under the curve (AUC; B) for small, medium, and large flocks of Rusty Blackbirds using Blitz data, eBird data alone, or Blitz and eBird data (“combined”). Solid lines (A) represent receiver operator curves for replicate runs across sampling methods and flock size classes. AUC values of 0.5 represent equivalent MaxEnt model performance relative to random, 0.5–0.6 “poor” performance, 0.6–0.7 “fair”, 0.7–0.8 “good”, 0.8–0.9 “very good”, and 0.9–1.0 “excellent”.
Diversity 13 00062 g003
Figure 4. The relative contribution of environmental covariates to predicted Rusty Blackbird occurrence. Variables are subset to those with the greatest explanatory power, with the bar graph on the right excluding minimum temperature (i.e., showing landcover variables only). Legend symbols (+,−) represent whether Rusty Blackbirds were positively or negatively associated with the variable. Rusty Blackbird observations peaked at intermediate minimum temperatures, and is thus denoted with +/−. For a complete list of variable contributions, see Table S1.
Figure 4. The relative contribution of environmental covariates to predicted Rusty Blackbird occurrence. Variables are subset to those with the greatest explanatory power, with the bar graph on the right excluding minimum temperature (i.e., showing landcover variables only). Legend symbols (+,−) represent whether Rusty Blackbirds were positively or negatively associated with the variable. Rusty Blackbird observations peaked at intermediate minimum temperatures, and is thus denoted with +/−. For a complete list of variable contributions, see Table S1.
Diversity 13 00062 g004
Figure 5. The gaussian kernel density distributions of average minimum temperature and floodplain forest (A,B) at occurrence points for each flock size class. Plot (C) describes the proportional density distribution of floodplain forest, which is the ratio of the kernel density of observations relative to that of the background points.
Figure 5. The gaussian kernel density distributions of average minimum temperature and floodplain forest (A,B) at occurrence points for each flock size class. Plot (C) describes the proportional density distribution of floodplain forest, which is the ratio of the kernel density of observations relative to that of the background points.
Diversity 13 00062 g005
Table 1. Number of checklists used in this analysis by flock size class and observation method. The sample size reflects the number of samples used after subsetting samples to one observation per raster grid cell.
Table 1. Number of checklists used in this analysis by flock size class and observation method. The sample size reflects the number of samples used after subsetting samples to one observation per raster grid cell.
Sampling MethodFlock Size Class
Small (1–19)Medium (20–99)Large (>99)Total
Blitz28112886495
eBird38723493714
Total6683621791209
Table 2. Land cover variable coefficients (λ) and contribution of each variable to model performance and small, medium, and large flocks. Variables represent the proportional cover of the land use category within a 4 km grid cell. Empty cells represent variables that were excluded from the candidate model sets due to lack of statistical support, as evaluated with AICc.
Table 2. Land cover variable coefficients (λ) and contribution of each variable to model performance and small, medium, and large flocks. Variables represent the proportional cover of the land use category within a 4 km grid cell. Empty cells represent variables that were excluded from the candidate model sets due to lack of statistical support, as evaluated with AICc.
Land Cover VariableVariable Coefficient (λ)Variable Contribution (%)
Flock SizeFlock Size
SmallMediumLargeSmallMediumLarge
Highly developed −0.83 1.2
Low-intensity development 1.171.41 1.31.1
Floodplain1.842.443.3314.914.126.4
Mixed forest 2.33 2.9
Pasture0.701.581.472.13.31.6
Row crop1.091.932.085.94.93.3
Shrub−3.26−2.29 3.71.1
Emergent wetand 2.50 1.5
Woody Wetland1.63 2.9
Woodland−1.73−1.64 1.91.3
Table 3. The difference (D) in the distribution of each environmental covariate between Rusty Blackbird flock size classes. p-values below 0.05, in bold, suggest that the flocks occupy a different environmental space.
Table 3. The difference (D) in the distribution of each environmental covariate between Rusty Blackbird flock size classes. p-values below 0.05, in bold, suggest that the flocks occupy a different environmental space.
Environmental VariableSmall vs. Medium FlocksSmall vs. Large FlocksMedium vs. Large Flocks
Dp-ValueDp-ValueDp-Value
High-intensity dev.0.0580.1970.0840.1640.0660.573
Low-intensity dev.0.0530.3240.0570.6230.3600.993
Floodplain forest0.112<0.0010.208<0.0010.1390.001
Hardwood forest0.0960.0050.0450.8660.1130.055
Mixed forest0.0240.9940.0990.0640.0960.153
Grassland0.117<0.0010.0840.1640.0790.338
Pasture0.0460.4890.0550.6630.0620.651
Row crop0.0360.7880.0980.0660.0960.146
Shrub0.0280.9560.0140.0020.1410.007
Upland forest0.0460.4890.0221.0000.0500.873
Emergent wetland0.0340.8360.0990.0650.06960.516
Woody wetland0.0610.1830.1160.0180.0870.240
Woodland0.0100.0030.1390.0020.0400.979
Precipitation0.0880.0130.1530.0010.0930.179
Minimum temperature0.170<0.0010.271<0.0010.1120.030
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Evans, B.S.; Powell, L.L.; Demarest, D.W.; Borchert, S.M.; Greenberg, R.S. Flock Size Predicts Niche Breadth and Focal Wintering Regions for a Rapidly Declining Boreal-Breeding Passerine, the Rusty Blackbird. Diversity 2021, 13, 62. https://0-doi-org.brum.beds.ac.uk/10.3390/d13020062

AMA Style

Evans BS, Powell LL, Demarest DW, Borchert SM, Greenberg RS. Flock Size Predicts Niche Breadth and Focal Wintering Regions for a Rapidly Declining Boreal-Breeding Passerine, the Rusty Blackbird. Diversity. 2021; 13(2):62. https://0-doi-org.brum.beds.ac.uk/10.3390/d13020062

Chicago/Turabian Style

Evans, Brian S., Luke L. Powell, Dean W. Demarest, Sinéad M. Borchert, and Russell S. Greenberg. 2021. "Flock Size Predicts Niche Breadth and Focal Wintering Regions for a Rapidly Declining Boreal-Breeding Passerine, the Rusty Blackbird" Diversity 13, no. 2: 62. https://0-doi-org.brum.beds.ac.uk/10.3390/d13020062

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop