Introduction
Estimating the true, total number of birds per unit area is one of the central problems of avian population biology and addressing this question is of utmost importance, since population estimates are decision-support tools of critical value for conservation and management. Within the context of biological and ecological inquiry, this issue has been approached from small-to-large scales, with literally dozens of techniques employed to estimate population numbers in areas that range in size roughly from 1 × 102 to 1 × 106 m2 (e.g. Ralph and Scott Reference Ralph and Scott1981, Bibby et al. Reference Bibby, Burgess, Hill and Mustoe2000).
The precision of those estimates decreases as researchers increase the geographic extent of the study (e.g. Flather and Sauer Reference Flather and Sauer1996). Habitat heterogeneity and its resulting variation in density of occupation, the scale of sampling, number of individuals that are undetected but present, and many other sources of error, affect their precision.
At the other end of the spectrum, large-scale estimates ranging from state to continental (roughly from 1 × 108 – 1 × 1010 m2), have received considerable attention in recent years as benchmarks to evaluate the success of conservation measures (Fitzpatrick Reference Fitzpatrick2002).
The European bird conservation community initiated this approach and organized a group of ornithologists to calculate country estimates for all species. These first results were published 15 years ago (Tucker et al. Reference Tucker, Heath, Tomialojc and Grimmett1994), with a subsequent assessment published a decade later (Burfield and van Bommel Reference Burfield and van Bommel2004).
Global population estimates, such as those made for waterbirds (Delany and Scott Reference Delany and Scott2006), involve the contribution of data from multiple field projects that are then summarized into a single figure. Since Delany and Scott's (Reference Delany and Scott2006) publication is the fourth iteration of this work, they also provide calculations of population change.
In North America, the Partners in Flight coalition (PIF) released a North American Landbird Conservation Plan (Rich et al. Reference Rich, Beardmore, Berlanga, Blancher, Bradstreet, Butcher, Demarest, Dunn, Hunter, Iñigo-Elías, Kennedy, Martell, Panjabi, Pashley, Rosenberg, Rustay, Wendt and Will2004) that contains an appendix with continental-scale population estimates for 448 species derived predominantly from Breeding Bird Survey data (BBS). The confidence of these values is scored in six categories of “accuracy rating” and seven categories of “precision.” In the geographic scope of Rich et al. (Reference Rich, Beardmore, Berlanga, Blancher, Bradstreet, Butcher, Demarest, Dunn, Hunter, Iñigo-Elías, Kennedy, Martell, Panjabi, Pashley, Rosenberg, Rustay, Wendt and Will2004: 6), “North American” excludes México, and we adopt this convention for the purpose of this paper.
These first North American estimates are a useful reference point of central importance. The authors of them (Rich et al. Reference Rich, Beardmore, Berlanga, Blancher, Bradstreet, Butcher, Demarest, Dunn, Hunter, Iñigo-Elías, Kennedy, Martell, Panjabi, Pashley, Rosenberg, Rustay, Wendt and Will2004, Rosenberg and Blancher Reference Rosenberg, Blancher, Ralph and Rich2005) have acknowledged their conservative nature and have invited the contribution of data from other sources in refining them. Most of the reviews and discussion over the conservation plan, however, have addressed programmatic aspects (e.g. Faaborg Reference Faaborg2005, Vickery and Shriver Reference Vickery and Shriver2005), and only Thogmartin et al. (Reference Thogmartin, Howe, James, Johnson, Reed, Sauer and Thompson2006) have covered the topic of the population estimates, with suggestions to improve the analyses for deriving them. Although we strongly support the idea of improving precision through enhancements in analysis methods, we also think that an approach that incorporates other empirical data can help us reach the goal of better figures.
In this paper, we present a case study of new North American population estimates for seven species of waterbirds, vultures, and diurnal raptors based on our long-term research in Veracruz, México. We recognize that our migration-count-based suggested revisions are more compelling for species whose breeding range lies entirely north of Mexico and winter range south of Veracruz. In three cases (waterbirds), we use our field data to indicate that existing figures are underestimates and suggest a mechanism to produce more appropriate calculations. In the remaining four cases (vultures and hawks), we use our migration count data to produce new estimates. In the closing section of this paper, we discuss the importance of incorporating data from multiple types of surveys to improve global and continental estimates, particularly (but not exclusively) those of the PIF.
Methods
Since 1991, we have made standardized vulture, hawk, and waterbird migration counts during (boreal) spring and autumn seasons in Veracruz, México. This region is a geographic bottleneck for migrants. It is located at the intersection of two major mountain systems, the Sierra Madre Oriental and the Central Volcanic Belt, which constrain the width of the Gulf Coastal Plain at about 19° N, c. 30 km north of Veracruz City (Ruelas et al. Reference Ruelas, Goodrich, Hoffman, Tingay, Chancellor and Meyburg2000, Reference Ruelas, Hoffman, Goodrich, Ralph and Rich2005).
These migration counts are made annually from two fixed locations 11 km apart (Cardel, 19°22′N-96°22′W, elevation 29 m, and Chichicaxtle 19°21′N-96°28′W, elevation 120 m), along an east-to-west axis (perpendicular to the front of migration). Both sites are independent samples of the flow of birds migrating through the area, as the distance between sites does not allow double counts, and the combined totals from these localities are used for long-term population monitoring. Because the goal of this project is monitoring rather than an attempt to quantify the total number of migrants in transit through the region, we know we are missing a proportion of the population that migrates east and west of our field stations. Most of the waterbirds regularly recorded from our migration sites rarely fly over the Gulf, whereas hawks and vultures have never been recorded flying over this large water body, making the shoreline of the Gulf of Mexico an effective barrier delimiting the eastern extent of this flyway.
Daily autumn counts have been collected from 20 August to 20 November of 1992–2004, and spring counts (roughly) from 20 February through 15 May of 1991, 1994, 1995, and 2004, from 08h00–18h00. Field teams collected hourly weather, flight-recording conditions, and migrant count totals, following the protocol of the Hawk Migration Association of North America (Kerlinger Reference Kerlinger1989, HMANA 2009). For a detailed description of the methods in use, please see Ruelas (Reference Ruelas2005), and Ruelas et al. (Reference Ruelas, Hoffman, Goodrich, Ralph and Rich2005).
Precision of Veracruz migration counts
Accepting our suggested revisions also implies accepting a level of precision in Veracruz migration counts. Any person who has witnessed the autumn migration in Veracruz (e.g. the c. 500 ornithologists attending the IV North American Ornithological Conference that visited our two monitoring sites in October 2006) can easily reach the conclusion that accurate counts are impossible due to the large volume of migration (Figure 1).
The method of collecting this information in the field involves a team of three field biologists per site that are randomly assigned to teams and localities. Most birds occur in large flocks that use the thermal convection of the atmosphere for cross-country flights (Pennycuick Reference Pennycuick1998). These counts were made once the thermal-climbing period had finished and birds were in their straight-line gliding flights. Our field teams directly counted the number of birds when quantities were low (e.g., a few tens) and made estimations when numbers were large. Teams estimated the number of birds with the help of a hand tally-counter by counting a small group of birds, e.g. 10 individuals, and extrapolating the approximate area occupied by these 10 birds to the rest of the flock. Most commonly-used multiples in these estimations were 3, 5, 10, and 50 individuals.
This method has been applied consistently over the years, and the relatively low interannual variation in a 13-yr series of data (Table 1), suggests that the level of error of these estimates remains relatively constant (Smith et al. Reference Smith, Farmer, Hoffman, Lott, Goodrich, Simon, Riley, Ruelas, Bildstein, Smith, Ruelas I. and Veit2008). Variation in migration counts is lower than the variation reported among BBS routes for most raptor species (Farmer et al. Reference Farmer, Hussell and Mizrahi2007). We have performed some experiments using slides, bird-counting software (Wildlife Counts http://www.wildlifecounts.com/), and video, and concur with other sources of information (Erwin Reference Erwin1982, Hodges Reference Hodges1993) that our methods produce conservative estimates, particularly for large-flocking migrants. Although the magnitude of such underestimations has not been assessed quantitatively, calculating this underestimation would be a good avenue to pursue for future estimates. Estimations of lower quantities of birds, however, are more precise than those of higher magnitude.
1 This estimate includes populations from northeastern México
We compared our mean, minimum, and maximum, annual autumn count of each species with existing North American population estimates from different published sources. Most water bird species estimates originate from colony counts and aerial surveys, while those of vultures and raptors are primarily from the BBS. Population estimates made by the authors listed under each species section have recently been revised by Kushlan et al. (Reference Kushlan, Steinkamp, Parsons, Capp, Acosta, Coulter, Davidson, Dickson, Edelson, Elliot, Erwin, Hatch, Kress, Milko, Miller, Mills, Paul, Phillips, Saliva, Snydeman, Trapp, Wheeler and Wohl2002), Rich et al. (Reference Rich, Beardmore, Berlanga, Blancher, Bradstreet, Butcher, Demarest, Dunn, Hunter, Iñigo-Elías, Kennedy, Martell, Panjabi, Pashley, Rosenberg, Rustay, Wendt and Will2004), Delany and Scott (Reference Delany and Scott2006) and Partners in Flight (2007). It is important to notice that our migration count data includes both juvenile and adult individuals, as the two age classes are impossible to track separately in the field, while the estimates of all the remaining published sources reportedly include only reproductive individuals.
Results and Discussion
American White Pelican Pelecanus erythrorhynchos. The global estimated population of this species is variable: 106,690 individuals (53,345 nests × 2 fide Koonz Reference Koonz1987); 120,000 individuals (Kushlan et al. Reference Kushlan, Steinkamp, Parsons, Capp, Acosta, Coulter, Davidson, Dickson, Edelson, Elliot, Erwin, Hatch, Kress, Milko, Miller, Mills, Paul, Phillips, Saliva, Snydeman, Trapp, Wheeler and Wohl2002); 134,000 individuals “east of the Rocky Mountains” (King and Anderson Reference King and Anderson2005), and over 180,000 (Delany and Scott Reference Delany and Scott2006: 31). The mean Veracruz count (85,679 individuals, Table 1) is lower than any of the North American estimates. However, since the majority (over 80%) of the winter range for this species lies north of our study sites in both the Gulf and the Pacific coasts (Knopf Reference Knopf2004), and, since these populations do not migrate through Veracruz, it is reasonable to state that Delany and Scott's (Reference Delany and Scott2006) highest estimate is substantially low, possibly as much as 60%, assuming a homogenous winter density, and therefore this estimate needs to be revised. An appropriate method for a more accurate estimate should include colony counts and winter density estimations across the non-breeding range of the species, possibly aided with Christmas Bird Count (CBC) data. Migration counts from localities north and south of our field sites could also help refine these estimates.
Anhinga Anhinga anhinga. Current North American population estimates range from 20,000–34,000 individuals (10,000–17,000 pairs × 2 fide Frederick and Siegel-Causey Reference Frederick and Siegel-Causey2000 and Kushlan et al. Reference Kushlan, Steinkamp, Parsons, Capp, Acosta, Coulter, Davidson, Dickson, Edelson, Elliot, Erwin, Hatch, Kress, Milko, Miller, Mills, Paul, Phillips, Saliva, Snydeman, Trapp, Wheeler and Wohl2002) to 51,000 individuals for the south eastern United States alone (Delany and Scott Reference Delany and Scott2006: 39). A very large proportion of the population of this species overwinters in its breeding range (Frederick and Siegel-Causey Reference Frederick and Siegel-Causey2000), and therefore is not recorded during migration in Veracruz. Our mean autumn count of 31,633 (Table 1), may record less than 40% of the individuals that belong to the migrating population (this estimate is based on the range map in Frederick and Siegel-Causey Reference Frederick and Siegel-Causey2000), and therefore indicates that all previous population estimates are low and should also be revised. We are unable to determine how low this estimate might be, but suggest that migration counts, non-breeding season colony counts, CBC, and winter density estimations paired with non-breeding season geographic distribution data, could help obtain a more precise number.
Wood Stork Mycteria americana. This species is estimated to have 25,000 individuals in the United States breeding range (Coulter et al. Reference Coulter, Rodgers, Ogden and Depkin1999). Delany and Scott (Reference Delany and Scott2006: 56) estimate 28,000 individuals for the southeastern United States alone. A more recent figure of 46,000 individuals includes México and Central America (Kushlan et al. Reference Kushlan, Steinkamp, Parsons, Capp, Acosta, Coulter, Davidson, Dickson, Edelson, Elliot, Erwin, Hatch, Kress, Milko, Miller, Mills, Paul, Phillips, Saliva, Snydeman, Trapp, Wheeler and Wohl2002).
The mean number recorded in Veracruz is higher than any of these figures (56,975, Table 1) and nearly 20% larger than the highest published estimate. These results suggest that previous population estimates for North America are conservative, but in the case of those including Mexico and Central America (Kushlan et al. Reference Kushlan, Steinkamp, Parsons, Capp, Acosta, Coulter, Davidson, Dickson, Edelson, Elliot, Erwin, Hatch, Kress, Milko, Miller, Mills, Paul, Phillips, Saliva, Snydeman, Trapp, Wheeler and Wohl2002) these numbers are too conservative, up to 65% below, based on winter range distribution relative to our field sites, and therefore these should be revised.
Wood Storks recorded on southbound autumn migration in Veracruz include an unknown proportion of juveniles that migrate north along the Gulf coast from colonies in the wetlands of the Grijalva and Usumacinta rivers (Tabasco and Chiapas) and the region of Centla (Tabasco), and possibly from the states of Campeche and Yucatán, during their post-fledgling period from late May through July. These individuals may join southbound migrants of northern origin later in the autumn and the number of individuals recorded in autumn migration counts may not belong to populations originating exclusively in the southeastern United States, but include birds from southern México. An adequate knowledge of the post-breeding dispersal and migration ecology of this species, colony counts in the non-breeding season, and migration counts will also be helpful in producing more appropriate estimates.
Turkey Vulture Cathartes aura. Population estimates cover the broad and expanding range of this species throughout North America (Kirk and Mossman Reference Kirk and Mossman1998). The United States and Canada population estimate in Rich et al. (Reference Rich, Beardmore, Berlanga, Blancher, Bradstreet, Butcher, Demarest, Dunn, Hunter, Iñigo-Elías, Kennedy, Martell, Panjabi, Pashley, Rosenberg, Rustay, Wendt and Will2004) is 1,305,000 individuals (or 29% of 4,500,000 globally). Our mean annual count is higher (1,895,679, Table 1) exceeding the highest published figure by over 30%, yet, Veracruz counts detect only a portion of the North American breeding population. Much (very likely more than 50%) of the eastern population of Turkey Vultures does not migrate as far as Mexico during the winter period (Kirk and Mossman Reference Kirk and Mossman1998).
Our field sites are not sampling a proportion of the population migrating west of our study sites (10%), since thermal-soaring obligates such as this species, are restricted to the lowlands of the coastal plain where they find abundant thermals for energy-saving migration. Our rounded estimate of the North American and northern México population migrating through Veracruz is 2,945,000 individuals. An accurate estimation of the North American population requires more information on the geography of migration of this species.
Mississippi Kite Ictinia mississippiensis. The North American global population estimate by Rich et al. (Reference Rich, Beardmore, Berlanga, Blancher, Bradstreet, Butcher, Demarest, Dunn, Hunter, Iñigo-Elías, Kennedy, Martell, Panjabi, Pashley, Rosenberg, Rustay, Wendt and Will2004) is 190,000 individuals, very similar to the estimate of Parker (Reference Parker1999), who obtained this information from our count data. This is a total migrant (breeding and wintering ranges are disjunct), whose most important migration route is the Gulf coastal plain of México, and a very high percentage of the global population of this species is recorded in our count sites. The mean number recorded in Veracruz is 157,199 (Table 1).
The number of Mississippi Kites tracked in our field sites is a sample of the total population migrating through this area. This is a facultative thermal-soaring migrant in Veracruz, and a proportion of the flight that we estimate to be 15% (and this percentage is slightly larger than in species of thermal-soaring obligates, Ruelas Reference Ruelas2005), is missed west (inland) of our study sites; therefore, our rounded population estimate for this species is 181,000 individuals. This species is undergoing an explosive population increase, at a rate of > 15% per year (Smith et al. Reference Smith, Farmer, Hoffman, Lott, Goodrich, Simon, Riley, Ruelas, Bildstein, Smith, Ruelas I. and Veit2008).
Broad-winged Hawk Buteo platypterus. Global population estimates of 1.7 million individuals are based on research from our study site (Goodrich et al. Reference Goodrich, Crocoll and Senner1996). This estimation concurs with the one of Rich et al. (Reference Rich, Beardmore, Berlanga, Blancher, Bradstreet, Butcher, Demarest, Dunn, Hunter, Iñigo-Elías, Kennedy, Martell, Panjabi, Pashley, Rosenberg, Rustay, Wendt and Will2004) of 1,728,000 individuals (96% of 1,800,000). The mean number recorded in Veracruz (1,931,255, Table 1) is slightly larger than earlier North American estimates. This species is a total migrant and, with the exception of populations overwintering in southern Florida and individuals migrating into Cuba, a very large proportion of the global population is recorded during migration in Veracruz. The Broad-winged Hawk is a thermal-soaring obligate restricted to the coastal plain and we estimate that 10% of the transient population is not sampled west of our study sites. Our rounded estimate for North America is 2,125,000 individuals.
Swainson's Hawk Buteo swainsoni. This is another species of total migrant whose global population estimates are based on unpublished data from our research sites (Ruelas et al. Reference Ruelas, Goodrich and Hoffman1996). During the period 1995–1996, our most conservative estimate was 448,000 (England et al. Reference England, Bechard and Houston1997), which coincides with the figure provided by Rich et al. (Reference Rich, Beardmore, Berlanga, Blancher, Bradstreet, Butcher, Demarest, Dunn, Hunter, Iñigo-Elías, Kennedy, Martell, Panjabi, Pashley, Rosenberg, Rustay, Wendt and Will2004) of around 460,600 (94% of 490,000). The mean number recorded in Veracruz (1,197,850 birds, Table 1) is 56.6% larger than the previous highest population estimate.
The migration of this ‘western’ species follows a path across northern México into the Gulf coastal plain and most Swainson's Hawks overwinter in Argentina, with the exception of populations in California, western México (Binford Reference Binford1989), and Florida (England et al. Reference England, Bechard and Houston1997). The disjunct distribution of this species in the breeding and wintering area, the limited options for migration routes, and the outcome of migration counts and satellite telemetry studies suggest a very large proportion of the global population migrates through this bottleneck (e.g. Fuller et al. Reference Fuller, Seegar and Schueck1998). Swainson's Hawk is also a thermal-soaring obligate restricted to the coastal plain and we estimate our sites are missing 10% of the total population transient through the area. Our rounded population estimate is 894,000 individuals.
Population estimates as a product of a metanalysis
The main difference between the Global Waterbird (Delany and Scott Reference Delany and Scott2006), European (Burfield and van Bommel Reference Burfield and van Bommel2004), and North American (Rich et al. Reference Rich, Beardmore, Berlanga, Blancher, Bradstreet, Butcher, Demarest, Dunn, Hunter, Iñigo-Elías, Kennedy, Martell, Panjabi, Pashley, Rosenberg, Rustay, Wendt and Will2004) population assessments is that the first two use an approach that combines data from multiple types of surveys that cover the breeding and the non-breeding season, whereas the latter is based almost completely upon BBS data that most likely excludes all non-breeding adult and juvenile individuals (Rich et al. Reference Rich, Beardmore, Berlanga, Blancher, Bradstreet, Butcher, Demarest, Dunn, Hunter, Iñigo-Elías, Kennedy, Martell, Panjabi, Pashley, Rosenberg, Rustay, Wendt and Will2004: 79, Rosenberg and Blancher Reference Rosenberg, Blancher, Ralph and Rich2005).
The authors of the North American landbird estimates (Rich et al. Reference Rich, Beardmore, Berlanga, Blancher, Bradstreet, Butcher, Demarest, Dunn, Hunter, Iñigo-Elías, Kennedy, Martell, Panjabi, Pashley, Rosenberg, Rustay, Wendt and Will2004, Blancher et al. Reference Blancher, Rosenberg, Panjabi, Altman, Bart, Beardmore, Butcher, Demarest, Dettmers, Dunn, Easton, Hunter, Iñigo-Elías, Pashley, Ralph, Rich, Rustay, Ruth and Will2007) acknowledge the limitations of this single-survey-based approach. Including data from different surveys could enable us to derive more appropriate population estimates, since no survey is capable of generating abundance indices for all species.
This is particularly important in the case of raptors. BBS data have been shown to have limitations estimating trends and indices of abundance in species detected at rates of < 1.0 birds per route (Peterjohn et al. Reference Peterjohn, Sauer, Robbins, Martin and Finch1995), which is the case for all the species treated in this paper. Similarly, Farmer et al. (Reference Farmer, Hussell and Mizrahi2007) have shown that, for most species, BBS population change estimates are of lower precision than those of migration counts.
Towards a precise conservation tool
The word of caution offered by Rosenberg and Blancher (Reference Rosenberg, Blancher, Ralph and Rich2005: 65), that “using inaccurate population estimates to set conservation targets may lead to misdirected conservation actions” is of critical importance. One example of this is the recently recovered Peregrine Falcon (Falco peregrinus). The estimate of Rich et al. (Reference Rich, Beardmore, Berlanga, Blancher, Bradstreet, Butcher, Demarest, Dunn, Hunter, Iñigo-Elías, Kennedy, Martell, Panjabi, Pashley, Rosenberg, Rustay, Wendt and Will2004) for the North America population is 276,000 (23% of 1,200,000 worldwide). However, an estimate based on nest surveys by White et al. (Reference White, Clum, Cade and Hunt2002), suggests the North American population ranges between 52,000 and 62,000 individuals, less than one-quarter than reported in the PIF estimate. As we further monitor the recovery of this species and the impact of management activities on their regional populations (Millsap and Allen Reference Millsap and Allen2006, USFWS 2007), it becomes essential to work with the best population estimates possible.
The opportunity to have a revised database of population estimates as a conservation tool of high precision may be in front of us. The PIF population estimates are now available at the state level (Partners in Flight 2007). The PIF Science Committee use the BBS-based approach (Blancher et al. Reference Blancher, Rosenberg, Panjabi, Altman, Bart, Beardmore, Butcher, Demarest, Dettmers, Dunn, Easton, Hunter, Iñigo-Elías, Pashley, Ralph, Rich, Rustay, Ruth and Will2007) that will soon incorporate the suggestions of Thogmartin et al. (Reference Thogmartin, Howe, James, Johnson, Reed, Sauer and Thompson2006) at this level of detail (P. Blancher, pers. comm.). Now that these state-level, baseline data exist, the fine-grain resolution will facilitate reducing the error of these estimates and will be open to the contribution of more data from other surveys and researchers (Confer et al. Reference Confer, Serrell, Hager and Lahr2008, Lanctot et al. Reference Lanctot, Hartman, Oring and Morrison2008). An adjustment to these estimates, like the one in our case study using migration count data at a large scale, is more feasible at finer scales.
Generating these first baseline data should be considered a work in progress, rather than a final word, in the creation of this collaborative conservation tool (Lanctot et al. Reference Lanctot, Hartman, Oring and Morrison2008, Farmer Reference Farmer2008). We believe that (1) The contribution of data derived from diverse surveys, and (2) Sampling data from individual projects that involve large proportions of the population of interest, enhances the precision and value of North American estimates.
Acknowledgements
We thank Peter Blancher, Tony Diamond, and Stuart Marsden for reviewing this paper and for their suggestions for improvement. Wesley Hochachka, Stefan Hames, André Dhondt's Lab group (Bird Population Studies, Cornell Lab of Ornithology), Peter Vickery, and John Faaborg also provided criticism of this manuscript. We thank long-term sponsors of the Veracruz River of Raptors Project, National Fish and Wildlife Foundation, U.S. Fish and Wildlife Service's Neotropical Migratory Bird Conservation Act, Lannan Foundation, American Bird Conservancy, Fondo Mexicano para la Conservación de la Naturaleza, John D. and Catherine T. MacArthur Foundation, Lincoln Park Zoo, U.S. Agency for International Development mission in Mexico City, Hawk Mountain Sanctuary, HawkWatch International, Pronatura Veracruz, and many private donors. We extend our gratitude to our dedicated field crews. ERI is indebted to John Faaborg and the University of Missouri for their support for this project; this paper was written while supported by a Howard Hughes Medical Institute postdoctoral fellowship through Dartmouth College. This is Hawk Mountain Sanctuary contribution to conservation science No. 182.