- Open Access
Temporal migration patterns between natal locations of ruby-throated hummingbirds (Archilochus colubris) and their Gulf Coast stopover site
© The Author(s). 2018
- Received: 8 November 2017
- Accepted: 22 December 2017
- Published: 10 January 2018
Autumn latitudinal migrations generally exhibit one of two different temporal migration patterns: type 1 where southern populations migrate south before northern populations, or type 2 where northern populations overtake southern populations en route. The ruby-throated hummingbird (Archilochus colubris) is a species with an expansive breeding range, which allows opportunities to examine variation in the timing of migration. Our objective was to determine a relationship between natal origin of ruby-throated hummingbirds and arrival at a Gulf coast stopover site; and if so, what factors, such as differences in body size across the range as well as the cost of migration, might drive such a pattern. To carry out our objectives, we captured hummingbirds at a coastal stopover site during autumn migration, at which time we collected feathers from juveniles for analysis of hydrogen stable isotopes. Using the hydrogen stable isotope gradient of precipitation across North America and published hydrogen isotope values of feathers from populations of breeding ruby-throated hummingbirds, we assigned migrants to probable natal latitudes.
Our results confirm that individuals from across the range (30–50° N) stopover along the Gulf of Mexico and there is a positive relationship between arrival day and latitude, suggesting a type 1 migration pattern. We also found no relationship between fuel load (proxy for migration cost) or fat-free body mass (proxy for body size) and natal latitude.
Our results, coupled with previous work on the spatial migration patterns of hummingbirds, show a type 1 chain migration pattern. While the mechanisms we tested do not seem to influence the evolution of migratory patterns, other factors such as resource availability may play a prominent role in the evolution of this migration system.
- Stable-hydrogen isotope ratio
- Spatial patterns
- Temporal patterns
- Gulf of Mexico
- Ruby-throated hummingbirds
The majority of forest-dwelling, avian species that breed in eastern North America migrate between temperate breeding areas and tropical wintering grounds . While migration has been well studied, information is missing on the migratory patterns of many species, notably the relationship between breeding locations and arrival timing at stopover sites, areas where migrants rest and refuel while en route . Spatiotemporal patterns of migration among populations have repercussions for resource competition at stopover sites with increasing conspecific densities [3, 4]. If populations overlap en route, then migrants will experience increased intraspecific competition in addition to other challenges, such as predation , unfamiliar habitat , interspecific competition , and weather .
Migratory movements involve both spatial and temporal components, which have been succinctly defined by Smith and colleagues . Spatially, a species can exhibit: 1) a chain migration pattern, where northern and southern populations show the same spatial pattern on the breeding and wintering grounds; or 2) a leap-frog migration pattern, where southern breeding populations winter further north than northern breeding populations. In terms of timing, species can show: 1) a type 1 migration pattern, where southern breeding populations initiate migration prior to northern breeding populations; or 2) a type 2 migration pattern, in which northern breeding populations begin to migrate before southern breeding populations. It is possible for species to show any combination of the spatial and temporal patterns described (e.g., “type 1 chain migration”, “type 2 chain migration”, etc.; see figure 1 in ), due to a variety of possible mechanisms. For example, northern breeding populations may face inclement conditions earlier in the season and initiate migration before southern breeding populations. Another possible driver is seasonal resource availability, which may influence initiation of migration between regional populations. Pienkowski and colleagues  hypothesized two additional factors to influence migration patterns, these include intraspecific competition and the cost of migration. If regional populations exhibit differential body size, then migration patterns may have evolved to alleviate intraspecific competition between populations. Finally, the cost of migration (the amount of time and energy it takes to travel between breeding and wintering locations) can influence migration timing of regional populations.
Ruby-throated hummingbirds (Archilochus colubris), hereafter “ruby-throats”, are a species with an expansive breeding range, which allows opportunities to examine variation in timing of migration. Ruby-throats are latitudinal migrants traveling between tropical wintering locations (Mexico and Central America) and temperate breeding grounds (United States and central Canada) along the same degree of longitude . Preliminary analysis of spatial migration patterns in ruby-throats by Hutcheson and colleagues  suggests a chain migration pattern (sensu [12, 13]). Our objective is to expand on their work to examine the temporal aspect of migration (type 1 chain migration or type 2 chain migration) in ruby-throats passing through a stopover site along the northern coast of the Gulf of Mexico using stable hydrogen isotope ratios. We hypothesize that timing of rubythroats’ migration will exhibit a strong temporal relationship with distance between the breeding location and our stopover site (e.g., latitude) since their southbound migration seems tightly tied to resource availability en route (see ). Specifically, we predict that ruby-throats show a type 1 chain migration pattern given that peaks in migration phenology seem to co-occur at stopover sites in the northern and southern portions of the species range [14, 15]. We also analyze factors identified by Pienkowski and colleagues  that may drive temporal migration patterns. We expect individuals from northern latitudes to have higher fat-free body masses (i.e. larger body size) compared to individuals from southern latitudes. We also expect lower fuel stores from individuals originating from higher (northern) latitudes since they would have travelled a longer distance compared to more southerly individuals when arriving at our study site.
We captured ruby-throats in Fort Morgan, AL (30°13′49″ N, 88°0′13″ W; see  for a description of the study site) between August 25 and November 1 in 2010, 2011, and 2014 using nylon mist nets (see  for a complete description of capture effort). Netting effort was both active (baiting some mist-nets with artificial feeders) and passive. We banded ruby-throats with a USGS aluminum band as well as sexed and aged (hatching year or after-hatching year) by bill corrugation, plumage, and morphology . We estimated subcutaneous fat , measured natural wing chord, recorded mass, and collected two outer rectrix (R4) feathers, one from each side of the tail. We stored feather samples individually in sealed and labeled paper envelopes until analysis of stable hydrogen isotope ratios.
Feather sampling and stable isotope analysis
For each year, we randomly selected juvenile (hatching year) individuals with feather samples (n = 25 for each sex per year; n = 150 total) throughout the autumn migration season to relate natal origin with date of passage. The feathers of hatching year birds during this phase of the annual cycle should reflect the stable hydrogen isotope ratios of the natal locations in which they were grown. We sent raw feathers to the Colorado Plateau Stable Isotope Laboratory (CPSIL; Flagstaff, Arizona, USA) for preparation and analysis of stable hydrogen isotope ratios. To prepare feathers for isotope analysis, CPSIL first cleaned the feathers with a phosphate-free detergent as well as a 2:1 chloroform methanol solution and rinsed them with deionized water before drying feathers at 50 °C overnight. CPSIL placed clipped feather material (0.350 mg; range: 0.330–0.370 mg) into silver capsules (3.5 × 5 mm) for analysis. CPSIL conducted sample pyrolysis with a Thermo Scientific TC/EA and hydrogen analysis via a Thermo Scientific Delta Plus IRMS configured through a Thermo Scientific CONFLO IV for automated continuous-flow analysis. The three normalization standards analyzed with feather samples included powdered forms of Keratin (SC Lot SJ; mean ± SD: −120.1 ± 1.0 ‰; expected: −121.6 ‰; n = 25), CBS – caribou hoof (mean ± SD: −197.5 ± 1.0 ‰; expected: −197.0 ‰; n = 8), and KHS – kudo horn (mean ± SD: −56.1 ± 1.6 ‰; expected: −54.1 ‰; n = 3). All δ2H values are reported in parts per mil (‰).
North America isotope precipitation model
We developed a δ2H isoscape model based on North American precipitation values analyzed in IsoMAP ([18, 19]; see  for results). We considered an extensive temporal range (1980–2009) in relation to four summer months (May–August) and based our δ2H isoscape inferences on the following independent variables: temperature, elevation, latitude and longitude [21, 22]. For the variable “temperature”, we used min, max, and average temperatures provided by Climate Research Unit . We limited the geographical extent of the δ2H isoscape model to a spatial range compatible with the Ruby-throated Hummingbird breeding distribution . This set of variables and parameters did not show significant spatial autocorrelation necessary for developing a δ2H geostatistical model of reference. Therefore, we adopted the δ2H statistical model, which is based on multiple linear regressions between δ2H precipitation values and the independent variables mentioned above .
To determine the natal origins of migrants stopping over along coastal Alabama, we built the species feather δ2H isoscape using a known-origin calibration dataset of published ruby-throat feather isotope data (n = 186) obtained from ten breeding populations [11, 25]. Then, we used the previously generated precipitation isoscape to extract environmental δ2H values at the same sampling locations reported by Hutcheson and colleagues ([26, 27]; see Additional file 1: Table S1). We resampled feather and precipitation δ2H values 1000 times using the site isotopic mean and SD values of the calibration dataset to create the rescaling function from our known-origin data. We used the mean and standard deviation of these bootstrapped regressions to convert the precipitation δ2H isoscape to feather δ2H values prior to assignment (see below and “rescale function” and “raster conversion function” in Additional file 2).
We adopted a likelihood approach to obtain natal location assignment probability surfaces for each individual sampled (n = 150) at Fort Morgan, AL in 2010, 2011, and 2014 and to generate the corresponding assignment map for each migrant using their δ2H feather values. This allowed us to obtain the individual assignment probabilities, which virtually connected each migratory ruby-throat to its putative natal site. Our assignment computations included the tissue-specific δ2H values from the raster conversion step, as well as the SD values for the calibration raster and for the original raster (regression model) computed in IsoMAP (for the computation details, see  and Additional file 2). We normalized the assignment probabilities with the cellStats function in the R package “raster”  and constrained the probability surfaces to the Ruby-throat breeding range (see Additional file 3: Figure S1, Additional file 4: Figure S2, Additional file 5: Figure S3, Additional file 6: Figure S4, Additional file 7: Figure S5 and Additional file 8: Figure S6) using species distribution data .
Assignment probabilities and latitudinal correlations
We implemented a Generalized Additive Mixed Model to investigate the linear or non-linear relationship(s) among predicted natal latitude (response variable) and non-correlated (Spearman’s rho: < 0.06) predictor variables, which included arrival day, fat-free body mass, and fuel load; year was included as a random factor. We assumed birds arrived on the date of first capture and transformed calendar date to an ordinal day. We calculated the sex-specific fat-free body mass of each individual using the data and regression method presented by Zenzal and Moore  for ruby-throats at Fort Morgan, AL, which is based the methods of Ellegren , as well as Owen and Moore . The fat-free body mass of each individual allowed us to 1) estimate a proxy for each individual’s body size based on fat-free mass, and 2) estimate the fuel load of each individual by subtracting the fat-free mass from the mass at capture. We conducted our analysis using the “mgcv” package [25, 31, 32] in the R statistical language (version 3.3.3; ).
As expected based on peaks in phenology occurring simultaneously at northern and southern stopover sites [14, 15], our results support a type 1 temporal migration pattern in juvenile ruby-throats – individuals from southern latitudes initiate migration before individuals from northern latitudes. This temporal migration pattern has also been found in sharp-shinned hawks  and two species of wood warbler (Parulidae; ) in the western United States and Canada during autumn migration. In our system, individuals from the lower-southeast (~ 30–35° N) passed through our study site in September followed by most individuals originating from higher latitudes (~ 40–50° N) in October. If we assume the same temporal pattern is exhibited as animals arrive on the wintering grounds and take into consideration previous work on the connectivity between breeding and wintering locations , then it is likely ruby-throats use a type 1 chain migration pattern during autumn (see figure 1a in ).
There are several mechanisms that have been hypothesized to drive the evolution of spatiotemporal migration patterns and we tested two purported by Pienkowski and colleagues . First, we found no relationship between ruby-throat body size (fat-free body mass) and natal origin, which does not support the hypothesis that differences in body size across the range reduce intraspecific competition. This finding is similar to that of American redstarts (Setophaga ruticilla; ) and contrasts with that of sharp-shinned hawks, which show a similar temporal migration pattern . However, Zenzal and Moore  found two distinct peaks in the phenology of juvenile ruby-throats, one peak in September and a separate peak in October, that coincidentally match when individuals from northern and southern latitudes passed through our study site as revealed by isotope analysis in this study. It is possible these peaks represent a southern population and northern population using stopover habitats along the northern coast of the Gulf of Mexico at different times, a strategy that would alleviate intraspecific competition. However, more study is needed to determine the origin of birds arriving during these separate peaks.
The second factor we tested, fuel load, did not support our prediction that there would be differences in the cost of migration across the range, rather individuals from both northern and southern latitudes arrived with similar fuel loads. These results are similar to other studies [8, 35–37] that found no relationship between fuel load and migratory distance. Moreover, most individuals in this study, as well as ruby-throats from this location previously described by Zenzal and Moore , typically arrived with some fuel stores. This may be due to birds being in close proximity to an ecological feature (Gulf of Mexico), where migrants tend to put on substantial fuel stores prior to negotiating a crossing (e.g., ). It is also possible given their small size  and high metabolism , ruby-throats always carry additional fuel stores to safeguard against times of energy shortfalls. Additionally, northern populations show synchrony between migratory movements and resource availability allowing birds to arrive on the Gulf coast in good energetic condition .
While beyond the scope of our study, two other mechanisms may influence migratory patterns. One possibility is that individuals show constant breeding area resident time (sensu ), allowing southern birds to breed earlier and hence migrate sooner than birds from northern latitudes. Another possibility is the availability of resources en route, which may influence when species time their migratory movements (e.g., [14, 41–44]) and this is strongly supported in hummingbirds [14, 41–43]. The timing of individuals from northern latitudes in our study tightly corresponded to the availability of Impatiens capensis and Lobelia cardinalis described by Bertin . The evolution of migration patterns may be tied to the ability to obtain resources en route, which has important conservation and habitat management implications. While ruby-throats are not a species of conservation concern, information on timing can be used to inform the general public when to focus on providing supplemental feeding – an activity with positive economic impacts [45, 46] and likely positive impacts on survival during migration.
Our results reveal that ruby-throats from across the breeding range used our coastal Alabama stopover site, with individuals from northern and southern latitudes arriving at different periods over the migratory season. When coupled with the spatial pattern described on the wintering grounds , ruby-throats show a type 1 chain migration pattern wherein individuals from southern latitudes initiate migration before individuals from northern latitudes. However, we were unable to identify the mechanism that might drive this pattern. We found no relationship between body size or migration cost and natal latitude, indicating other factors such as resource availability, competition, or constant breeding area resident time may be responsible for this migration pattern. Future research should work to find the mechanism(s) that may have influenced the evolution of the type 1 migration pattern in ruby-throats.
We thank the 2010-2014 Fort Morgan banding station field technicians and members of the Migratory Bird Research Group for their help and support. We also thank the Bon Secour National Wildlife Refuge and Fort Morgan State Historic Site for assistance and permission to use their properties. We would also like to thank C. Hutcheson, L. Hendrix, J. Moran, and, L. Wassenaar for providing isotope data from their work as well as the banders that provided feathers in their studies: L. Beall, F. Bassett, S. Campbell, A. Chartier, T. Gordon, C. Hutcheson, M. Niedermeier, R. Leberman, R. Mulvihill, E. Womack, R. Yunick, and L. Walewski. We thank the University of Southern Mississippi, University of Illinois at Urbana-Champaign, and University of Oklahoma for support, computing time, and copy services. This manuscript was improved through comments and critiques from H. Vander Zanden, members of the Ward/Brawn lab – especially J.K. Kelly and T.J. Jones, as well as four anonymous reviewers.
This work was supported by the National Science Foundation (NSF) (IOS 1147096; GK-12 0947944); National Geographic Committee on Research and Exploration (Awards 8971–11 and 9155–12); Birmingham Audubon Society; Alabama Ornithological Society; Wilson Ornithological Society; Louisiana Ornithological Society; Inland Bird Banding Association, and Eastern Bird Banding Association. T.J.Z. was supported by a NSF GK-12 Program Award (0947944), the Mississippi-Alabama Sea Grant Consortium (NA14OAR4170089), and the NOAA RESTORE Act Science Program (NA17NOS4510092). Funding agencies had no role in the study design, collection of data, analysis, interpretation, or writing of the manuscript.
Availability of data and materials
The datasets used and analysed during the current study are available from the corresponding author on reasonable request.
TJZ, FRM, and JFK conceived the project, TJZ collected the data, TJZ and AJC analyzed the results and wrote the paper, FRM and JFK provided substantial edits on the paper. All authors read and approved the final manuscript.
Ethics approval and consent to participate
All research activities were approved by the University of Southern Mississippi Institutional Animal Care and Use Committee (Protocol #11092210), U.S. Geological Survey Bird Banding Laboratory (Permit #21221), and the Bon Secour National Wildlife Refuge.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Rappole JH. The ecology of migrant birds: a Neotropical perspective. Washington, D.C.: Smithsonian Institution Press; 1995.Google Scholar
- Kelly JF. Stable isotope evidence links breeding geography and migration timing in wood warblers (Parulidae). Auk. 2006;123:431–7.View ArticleGoogle Scholar
- Moore FR, Yong W. Evidence of food-based competition among passerine migrants during stopover. Behav Ecol Sociobiol. 1991;28:85–90.View ArticleGoogle Scholar
- Kelly JF, DeLay LS, Finch DM. Density-dependent mass gain by Wilson’s warblers during stopover. Auk. 2002;119:210–3.View ArticleGoogle Scholar
- Cimprich DA, Moore FR. Energetic constraints and predation pressure during stopover. Proc Int Ornithol Congr. 1999. p. 834–46.Google Scholar
- Németh Z, Moore FR. Unfamiliar stopover sites and the value of social information during migration. J Ornithol. 2007;148:369–76.View ArticleGoogle Scholar
- Newton I. Weather-related mass-mortality events in migrants. Ibis. 2007;149:453–67.View ArticleGoogle Scholar
- Smith RB, Meehan TD, Wolf BO. Assessing migration patterns of sharp-shinned hawks Accipiter striatus using stable-isotope and band encounter analysis. J Avian Biol. 2003;34:387–92.View ArticleGoogle Scholar
- Pienkowski MW, Evans PR, Townshend DJ. Leap-frog and other migration patterns of waders: a critique of the Alerstam and Högstedt hypothesis, and some alternatives. Ornis Scand. 1985;16:61–70.Google Scholar
- Weidensaul S, Robinson TR, Sargent RR, Sargent MB. Ruby-throated Hummingbird (Archilochus colubris). Rodewald PG, editor. Birds N. Am. Online [Internet]. 2013 [cited 2017 Feb 27]; Available from: https://doi.org/10.2173/bna.204.
- Zenzal TJ, Fish AC, Jones TM, Ospina EA, Moore FR. Observations of predation and anti-predator behavior of ruby-throated hummingbirds during migratory stopover. Southeast Nat. 2013;12:N21–5.View ArticleGoogle Scholar
- Salomonsen F. Evolution and bird migration. Acta XI Congr Int Ornithol. 1955;337:339.Google Scholar
- Boulet M, Norris DR. Introduction: the past and present of migratory connectivity. Ornithol Monogr. 2006;61:1–13.View ArticleGoogle Scholar
- Bertin RI. The ruby-throated hummingbird and its major food plants: ranges, flowering phenology, and migration. Can J Zool. 1982;60:210–9.View ArticleGoogle Scholar
- Zenzal TJ, Moore FR. Stopover biology of ruby-throated hummingbirds (Archilochus colubris) during autumn migration. Auk. 2016;133:237–50.View ArticleGoogle Scholar
- Pyle P. Identification guide to North American birds: Columbidae to Ploceidae. Bolinas: Slate Creek Press; 1997.Google Scholar
- Helms CW, Drury WH. Winter and migratory weight and fat field studies on some North American buntings. Bird-Band. 1960;31:1–40.View ArticleGoogle Scholar
- Welker JM. Isotopic (δ18O) characteristics of weekly precipitation collected across the USA: an initial analysis with application to water source studies. Hydrol Process. 2000;14:1449–64.View ArticleGoogle Scholar
- IAEA/WMO. Global network of isotopes in precipitation [Internet]. GNIP Database. 2011. Available from: http://www.iaea.org/water.
- Contina A. IsoMAP job 59969, Isoscapes modeling, analysis and rediction (version 1.0). [Internet]. IsoMAP Proj. 2017. Available from: http://isomap.org.
- U. S. National Geophysical Data Center. ETOPO-5 five minute gridded world elevation. Boulder: NGDC; 1998.Google Scholar
- Mitchell TD, Jones PD. An improved method of constructing a database of monthly climate observations and associated high-resolution grids. Int J Climatol. 2005;25:693–712.View ArticleGoogle Scholar
- BirdLife International, Handbook of the Birds of the World. Bird Species Distrib. Maps World. 2016. Available from: http://datazone.birdlife.org/species/requestdis.
- Bowen GJ, Liu Z, Vander Zanden HB, Zhao L, Takahashi G. Geographic assignment with stable isotopes in IsoMAP. Methods Ecol Evol. 2014;5:201–6.View ArticleGoogle Scholar
- Wood SN. Stable and efficient multiple smoothing parameter estimation for generalized additive models. J Am Stat Assoc. 2004;99:673–86.View ArticleGoogle Scholar
- Hutcheson CA, Hendrix L, Moran JA. An isotopic analysis of migratory connectivity in ruby-throated hummingbirds. North Am Bird Bander. 2010;35:5–11.Google Scholar
- Hutcheson CA, Wassenaar LI, Hendrix L. A preliminary examination of the use of hydrogen isotope ratios in estimating the natal latitudes of hatching-year Ruby-throated Hummingbirds. North Am Bird Bander. 2007;32:68-74.Google Scholar
- Hijmans RJ. raster: Geographic data analysis and modeling. R package version 2.5–8 [Internet]. R Found. Stat. Comput. Vienna 2016. Available from: https://CRAN.R-project.org/package=raster.
- Ellegren H. Estimated effects of age and sex on the fat-free body mass of autumn migrating bluethroats Luscinia s. svecica. Ardea. 1992;80:255–9.Google Scholar
- Owen JC, Moore FR. Seasonal differences in immunological condition of three species of thrushes. Condor. 2006;108:389–98.View ArticleGoogle Scholar
- Wood SN. Generalized additive models: an introduction with R. Boca Raton: CRC Press; 2006.Google Scholar
- Wood SN. Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. J R Stat Soc Ser B Stat Methodol. 2011;73:3–36.View ArticleGoogle Scholar
- R Core Team. R: A language and environment for statistical computing [Internet]. Vienna, Austria: R Foundation for Statistical Computing; 2017. Available from: https://www.R-project.org.
- Norris DR, Marra PP, Bowen GJ, Ratcliffe LM, Royle JA, Kyser TK. Migratory connectivity of a widely distributed songbird, the American redstart (Setophaga ruticilla). Ornithol Monogr. 2006;61:14–28.View ArticleGoogle Scholar
- Hobson KA, Van Wilgenburg S, Wassenaar LI, Moore F, Farrington J. Estimating origins of three species of neotropical migrant songbirds at a gulf coast stopover site: combining stable isotope and GIS tools. Condor. 2007;109:256–67.View ArticleGoogle Scholar
- Wilson S, Hobson KA, Collister DM, Wilson AG. Breeding destinations and spring migration patterns of Swainson’s Thrush (Catharus ustulatus) at a Costa Rican stopover site. Auk. 2008;125:95–104.View ArticleGoogle Scholar
- Shneor O, Perlman G, Balaban A, Yom-Tov Y, Huchon D, Geffen E. Origin of passerine migratory waves: evidence from the blackcap at a stopover site. Isr J Ecol Evol. 2010;56:135–51.View ArticleGoogle Scholar
- Caldwell LD, Odum EP, Marshall SG. Comparison of fat levels in migrating birds killed at a central Michigan and a Florida Gulf Coast television tower. Wilson Bull. 1963;75:428–34.Google Scholar
- Lasiewski RC. Oxygen consumption of torpid, resting, active, and flying hummingbirds. Physiol Zool. 1963;36:122–40.View ArticleGoogle Scholar
- Thorup K, Tøttrup AP, Rahbek C. Patterns of phenological changes in migratory birds. Oecologia. 2007;151:697–703.View ArticlePubMedGoogle Scholar
- Tyler W. Ruby-throated hummingbird habitats. In: Bent AC, editor. Life Hist. North Am. Cuckoos goatsuckers hummingbirds their allies. Washington D.C.: United States National Museum Bulletin; 1940. p. 332–50.Google Scholar
- Grant KA, Grant V. Hummingbirds and their flowers. New York: Columbia University Press; 1968.Google Scholar
- Austin DF. Bird flowers in the eastern United States. Fla Sci. 1975;38:1–12.Google Scholar
- Aborn DA. Correlation between raptor and songbird numbers at a migratory stopover site. Wilson Bull. 1994;106:150–4.Google Scholar
- Glowinski SL. Bird-watching, ecotourism, and economic development: a review of the evidence. Appl Res Econ Dev. 2008;5:65–77.Google Scholar
- U.S. Department of the Interior, U.S. Fish and Wildlife Service, U.S. Department of Commerce, U.S. Census Bureau. National Survey of fishing, hunting, and wildlife-associated recreation. 2011.Google Scholar