Abstract
Subsistence hunting of neotropical birds is common and widespread in the tropical forests of Latin America. Although its sustainability under different scenarios is subject to debate, hunting has already contributed to the decline and local extirpation of several taxa and is considered to be a significant threat to a range of large-bodied species. Gaining a better understanding of the variability of hunting patterns, as well as the factors that can potentially be used to predict them, is important if we are to develop conservation strategies that target the species most likely to be experiencing declines. In this article, we examine the avian hunting profiles of 65 communities in the neotropics. We describe their variability and look at the relationship between a hunting profile and (a) its geographical location, (b) the community’s age, (c) the community’s population size, and (d) the year in which the survey was carried out. We find that there is a significant but weak relationship between a community’s geographic location and the composition of its bird hunting profile, and that prey profiles can be considerably different even among close neighbors. We found no relationship between a community’s age or size and the mean biomass of bird prey hunted. Our results challenge the assumption that older and larger settlements have a predictable impact upon avian prey communities and suggest that cultural preferences or the starting availability of prey species can change rapidly over short distances.
Introduction
Subsistence hunting of birds, mammals, and reptiles for their meat and other body parts is common and widespread among both indigenous and settler communities in the tropical forests of Latin America. Some species have historically provided an important source of protein (Begazo & Bodmer, 1998; Lenselink, 1972; Smith, 1976; Vickers, 1980) and fats (Sirén & Machoa, 2008), whereas others are sources of tools, traditional medicine, and adornments that play a prominent role in the culture of many indigenous groups (Bezerra, de Araujo, Alves, & Alves, 2013; Fernandes-Ferreira, Mendonça, Albano, Ferreira, & Alves, 2012; Mena, Stallings, Regalado, & Cueva, 2000; Santos-Fita, Naranjo, & Rangel-Salazar, 2012). Hunting for the above purposes, however, is considered to be one of the major drivers of global biodiversity loss and has led to population declines in several species whose life history traits or perceived value makes them vulnerable to overexploitation (Brooks & Strahl, 2000; Cullen, Bodmer, & Valladares Pádua, 2000; Nuñez-Iturri & Howe, 2007; Peres, 1990, 2000; Rosin & Swamy, 2013). A recent synthesis of 176 studies found that bird and mammal abundances are an average of 58% and 83% lower in hunted sites than nonhunted sites respectively (Benítez-López et al., 2017), with higher declines associated with areas that are close to access points (such as roads) and have greater access to urban markets (Benítez-López et al., 2017). Rapid increases in the rate of infrastructure construction in Latin America (Fraser, 2014), as well as the potential development of commercial wild meat markets similar to those monitored in Africa (see Suárez et al., 2009), are raising important questions about the sustainability of future wild meat harvests, as well as the implications that species declines will have for the ecology of tropical forests.
Much of the hunting literature in the neotropics has focused on mammals, and the susceptibility of species such as tapirs, peccaries, and large ateline primates to declines. Birds, however, have received less attention, though hunting is recognized to be a major threat to several neotropical species (BirdLife International, 2017a; Galetti, Martuscelli, Olmos, & Aleixo, 1997), and some communities prefer their meat over that of mammals (Mena et al., 2000; Santos-Fita et al., 2012). Most ethnozoological studies on the use of birds have additionally focused on their use as pets or on their cultural importance (Haenn, Schmook, Reyes, & Calmé, 2014; Nascimento Rodrigues, Czaban, & Nóbrega Alves, 2015; Nóbrega Alves, De Farias Lima, & Araujo, 2012; Roldán-Clarà, López-Medellín, Espejel, & Arellano, 2014), and not their role as game. Hunting is considered to be a particularly serious threat to several species belonging to the Cracidae (including guans, curassows, and chachalacas), as a result of their large body size (which increases their appeal to hunters) and their slow rate of population recovery (Brooks & Fuller, 2006; Thiollay, 2005). Of the 56 neotropical species of Cracid, 41% are in a threatened IUCN category, whereas an additional 9% are considered near threatened (IUCN Galliform Specialist Group, 2017). Hunting is believed to be a significant threat to all six Cracid species classed as critically endangered and is considered to have exacerbated the decline of the Alagoas currasow (
The decline or extirpation of these large, frugivorous species is likely to have knock-on effects on the ecology of tropical forests (Peres, Thaise, Schietti, Desmoulieres, & Levi, 2016; Wright, Hernandéz, & Condit, 2007). Toucans, for example, are the most important seed dispersers of some species of
The diversity of avian species that are targeted by hunters in each particular community will be the product of two influences: (a) how easy each species is to hunt (itself influenced by the abundance and behavior of that species in a particular area) and (b) its perceived value, in terms of its taste or cultural prestige. Although the outcome of this balance has been recorded in detail for a number of communities, we lack an understanding of how variable avian hunting profiles are on a cross-continental scale. The first aim of this article is therefore to assess the variability of bird hunting profiles across Latin America and to investigate whether an easy-to-obtain attribute of communities—their geographical location—is a good predictor of hunter preferences. Understanding the extent to which profiles are universal and the extent to which they vary on a case-by-case basis is important if we are to develop management plans that address the taxa most likely to be impacted; understanding the predictive power of geography alone will also allow us to assess whether we are able to accurately predict the likely preferences of communities that have not received research attention.
The second part of this article aims to assess whether avian hunting profiles change in a predictable way in response to high or prolonged of hunting pressure. Hunters in each community are principally assumed to act in a way that tries to maximize the amount of meat harvested for a given amount of hunting effort (although the cultural value of different prey species is expected to introduce small levels of skew; Jerozolimski & Peres, 2003). They may therefore be initially expected to focus on species that are both large and relatively abundant, which has been shown to be the case in multiple studies (e.g., Franzen, 2006; Mena et al., 2000; Ohl-Schacherer et al., 2007; Sirén, 2004). However, if larger species become scarce through overexploitation and the effort involved to hunt them becomes uneconomic, hunters may be expected to shift to a wider range of smaller prey species. We might therefore expect the prey items of settlements that are larger (and thus would be expected to exert a higher hunting pressure on their surrounding forest) or have been hunting for a longer amount of time to have a lower mean biomass. Although this phenomenon has been widely observed in fisheries data (Friedlander & DeMartini, 2002; Pauly, Christensen, Dalsgaard, Froese, & Torres Jr, 1998; Pauly & Palomares, 2005), evidence of it occurring in terrestrial systems remains more tentative. In the neotropics, Jerozolimski and Peres (2003) found that the mean biomass of mammalian prey hunted declined as the age of settlements increased but not as their populations increased. Similarly, a study by Ingram et al. (2015) on hunting profiles from Central and West Africa found that the mean body mass of mammalian prey hunted had decreased significantly between 1966 and 2010. Contrary to expectations, however, the mean body mass of birds harvested increased between 1975 and 2010; possibly as a result of changes in demand for large birds and their bills such as the black-casqued hornbill. The authors suggest that an indicator based on the mean body mass of prey hunted (the mean body mass index) could be used as part of a suite of measures to monitor hunting pressure in the absence of comprehensive monitoring systems. Whether patterns such as these can be detected in neotropical bird hunting data remains untested. In this article, we therefore explore whether the age and size of communities can be used to predict a decline in the mean biomass of birds hunted; we also investigate the relationship between the mean body mass of prey and the year in which the study monitoring hunting activity took place in order to allow comparison of our data to the Ingram et al.’s (2015) study. We expect the mean biomass of prey items to decrease as settlements get older, larger, or had their hunting study conducted in later years.
Finally, it is important to consider the interactions between harvests of avian prey and the harvest of other taxa including mammals and reptiles, as shortages in preferred species belonging to one taxon might well be compensated for by increasingly targeting large species in a different group of animals. We briefly assess whether communities that hunt a large total biomass of birds hunt a lower biomass of mammals or reptiles. Through our analyses, we hope to (a) better characterize the variability of neotropical hunting profiles across space and time, (b) assess whether high or prolonged hunting pressure has a predictable impact on the mean biomass of birds hunted, and (c) identify data gaps that could better inform these assessments.
Methods
Literature Search and Data Compilation
We performed a literature search for studies listing hunting records of neotropical birds, using papers returned from the standard search terms:
Geographic Variability
For each community, we calculated the percentage of total kills and percentage of total biomass accounted for by each avian order. In cases where the biomass of a particular species was not recorded, we averaged their recorded biomass across all studies that had actively measured weights (i.e., did not rely on literature values), and multiplied by the number of individuals harvested. If information on the biomass of species was not available (because it had not been hunted elsewhere or because it was only included as part of a large species grouping such as
By using such a high taxonomic level to measure similarity, we will inevitably be losing resolution and the ability to detect preferences for different families. However, because not all families occur across all study locations (Phasianidae, for example are only represented by a single species [
To test whether communities that were closer together geographically had similar hunting profiles, we used our data set to generate two distance matrices: one based on the similarity of hunting profiles as defined by the percentage of individuals hunted belonging to each avian order, and one based on geographical proximity using the latitude and longitude coordinates of each settlement. The hunting similarity matrix was calculated using a Bray-Curtis dissimilarity index in the
Correlates of Average Biomass
We used Spearman’s rank correlations, as in, Spearman's rank correlations to investigate the relationship between the mean biomass of birds hunted and (a) the age of communities, (b) the population size of communities, and (c) the year in which each hunting study was carried out. Average biomass was calculated by dividing the total biomass of birds hunted by the total number of individuals recorded in the profile. Information on the population of communities at the time of study was available for
Interactions Between Bird, Mammal, and Reptile Hunting
We investigated whether the total biomass of birds in each profile was correlated to the total biomass of mammals and reptiles recorded using a linear model with
Results
Overview
Number of Species Belonging to Each Order Recorded Across

Number of mammal and bird genera reported in
Galliformes was the dominant order hunted by both the number of animals harvested and by biomass (Figure 2) and appeared in the highest proportion of lists (97%). They accounted for an average of 57% of the total number of individuals harvested in the profiles of communities where at least one member of the order was hunted and were the only order of birds targeted in six communities (Figure 2). Within the Galliformes, in communities where it is included in prey lists, (a) Percentage of kills accounted for by the top six avian orders reported in offtake lists using only communities where a member of that order is listed as being hunted and (b) percentage of kills accounted for by all Galliform genera, using only those where a genus is included in an offtake list. The number of prey profiles featuring each genus is given below the genus name, alongside its average biomass (± standard error), calculated from all studies which weighed carcasses and did not use literature values. Red diamonds show the average percentage kills for each order or genus. Whiskers extend to the largest value no further than 1.5*IQR (inter-quartile range) from each hinge.
Geographic Variability
The geographic proximity of settlements and the shape of their hunting profile (in terms of the percentage of kills accounted for by different avian orders) were significantly but weakly correlated (Mantel test, Proportion of (a) kills and (b) biomass belonging to each order, from Comparison of trees generated using the Ward method of hierarchical clustering, showing the relationship of settlements according to (a) their geographical proximity and (b) the similarity of their hunting profile according to the percentage of kills belonging to each avian order.

Several clusters of settlements that are geographically close have similar hunting profiles. Settlements in Yucatán and Quintana Roo, for example, almost solely hunt Galliformes, which is reflected in their tight clustering in the hunting profile similarity tree. Settlements in the Ecuadorian Oriente tend to have more diverse profiles that are principally composed of Galliformes and Piciformes, and most of them grouped together in our profile similarity tree. However, there are also several groups of communities that are geographically close but have very different profiles. Notably the Nicaraguan communities in our sample, though clustered tightly geographically, are scattered throughout the hunting offtake similarity tree; indicating that in certain cases, preferences can change rapidly over short distances.
Correlates of Avian Prey Biomass
Figure 5 shows the average biomass of bird kills versus (a) the age of communities, (b) the size of communities, and (c) the starting year of each hunting study. We found no significant relationship between any of the three variables and the average biomass of birds hunted (Table 2). There was a significant positive relationship between the biomass of mammals and the biomass of birds harvested across communities but not between the biomass of birds and reptiles (Table 3, Figure 6).
Average biomass of bird kills versus (a) settlement age ( Spearman’s Rank Correlation Coefficients (Probability, Sample Size) for the Age of Communities, Size of Communities, and Year of Study Versus the Mean Undressed Biomass of Bird Prey Included in Their Offtake List. Linear Model of Bird Biomass Hunted Versus Mammal and Reptile Biomass Hunted (Number of Observations = 54, Biomass of birds hunted versus the biomass of mammals hunted in 

Discussion
This article aimed to give an overview of the variability of bird hunting profiles on a pan-continental scale and to investigate whether different aspects of a community’s hunting profile could be predicted by four variables: its location, its age, its population size, and the year in which the hunting survey was conducted. We found that there was a significant but weak correlation between a community’s location and the structure of its hunting profile in terms of the proportion of kills belonging to each avian order. In addition, there was no correlation between the age of a community, its population size, or the year in which the hunting study was conducted and the average biomass of individual birds hunted. Our results raise several interesting questions for further research but also highlight important gaps in the literature pertaining to bird hunting which, if filled, could allow more sophisticated analysis of its causes and impacts.
Although the data sets included in this study covered a wide geographical area, they were tightly clustered and revealed sizeable gaps in the availability of avian hunting offtake data for large regions of Central America, the Amazon, and the Guianan shield. We did not find any published data for studies carried out in the state of Amazonas, Brazil, or Loreto, Peru, despite them covering a significant proportion of their respective country’s Amazonian territory. Likewise, we were unable to find hunting offtake lists from any communities in Costa Rica, Guatemala, Honduras, Belize, or Guyana. Nevertheless, our data showed that neotropical bird offtake lists tend to be less diverse (in terms of genera) than mammal lists, and that they are dominated by five orders which were recorded in over 50% of profiles (Galliformes, Tinamiformes, Psittaciformes, Gruiformes, and Piciformes). This narrower scope could potentially reflect the smaller selection of bird genera that reach a size where hunting them is considered worthwhile. Carcasses of relatively large genera such as
Geographic Variation
Several studies have shown that what a community hunts will not only be influenced by the efficiency of harvesting different species in that particular community’s catchment but also the cultural and personal considerations of the
A third confounding factor to consider is the wide time span over which the data included in this study were gathered. The 42 years over which the data included in this study were collected have seen many communities undergo considerable cultural and technological changes (Roosevelt, 2013). Many communities have become better connected to market towns and cities (Apaza et al., 2002; Stafford, Alarcon-Valenzuela, Patiño, Preziosi, & Sellers, 2016; Suárez et al., 2009), new technologies such as guns, outboard motors, and lamps have become more widely available (Hames & Vickers, 1982; Yost & Kelley, 1983), new laws have been passed that modify the legality and cost of different practices (e.g., Sirén & Wilkie, 2014), and many communities are experiencing population growth (Holt, 2005). This raises important questions about the validity of comparing community preferences that may well have changed in the intervening years. Controlling for year of study in our analysis testing whether communities that were closer together had similar hunting profiles made little difference to the outcome of the analysis. Nevertheless, more longitudinal data sets are needed to search for trends and to assess the comparability of data sets covering a wide time span.
Correlates of Prey Profiles
We found no evidence that the mean biomass of individual birds hunted decreased as settlements get older or larger, or of decreased mean biomass in studies that were conducted in later years. Taken at face value, our findings indicate either (a) that larger, preferred bird species are, in general, being hunted at sustainable rates or (b) that it is not worth switching to smaller species of birds once the larger ones become scarce, and a better economic choice may be to switch to other sources of wild meat such as mammals or fish. The former explanation is not supported by data from studies which have found a lower abundance of Cracids, macaws, and toucans in hunted forests (Begazo & Bodmer, 1998; Peres & Nascimento, 2006; Thiollay, 2005), although some studies suggest that the availability of nearby unhunted areas to act as a source of migration to replenish stocks may be enough to protect some bird species from depletion (Ohl-Schacherer et al., 2007), and not all species hunted in these studies were harvested unsustainably. Investigating the second possibility is difficult because hunting profiles are rarely accompanied by a full inventory of species that were theoretically available to hunt, and we could not find any studies that explicitly tested how preferences may switch between groups. Longitudinal data sets that track preferences of all groups (including fish) as well as the abundance of different species through time may present a better avenue under which to investigate these patterns.
A third possibility is that the three proxies that we selected for our analysis—the age of settlements, size of settlements, and year of the hunting study—are not sufficiently informative of actual hunting pressure. The true hunting pressure an area experiences will be influenced by several factors, including the proportion of a settlement’s population that actively hunts, the availability of alternative protein sources such as domestic animals and fish, the size of hunting catchments, and their governance. In addition, the ability of each particular species to persist in a catchment will be influenced by the productivity of forests (which is known to vary across the continent; Emmons, 1984; Haugaasen & Peres, 2005) and the availability of nearby unhunted areas that can replenish hunted stocks (Ohl-Schacherer et al., 2007). Although previous studies found that decreases in biomass could be detected without adjusting for these different sources of variation (Jerozolimski & Peres, 2003), the discord between our failure to find a pattern and the declines of large species documented in the literature suggests a need to investigate whether accounting for these factors can give a better picture of the impacts of hunting in cross-sectional (as opposed to longitudinal) data. However, even simple measures such as the age and population size of communities were only available for 41% and 72% of settlements, respectively. Restricting our analysis further to include only studies where these additional variables were available would have significantly reduced the available pool of hunting profiles.
Finally, our data suggest that communities that hunt large quantities of mammal biomass also hunt high quantities of bird biomass, suggesting that they do not tend to specialize in either one or the other. It would be interesting to investigate whether a similar relationship exists with the diversity of species hunted (i.e., whether communities that hunt a wide diversity of mammal species hunt a wide diversity of bird species), but the lack of identification to species or genus level that was commonly observed in our bird offtake lists meant we were not able to investigate this possibility.
Implications for Conservation
Our finding that geography is a significant but weak predictor of the relative proportions in which avian orders are hunted suggests that profiles in communities that have not been monitored will not necessarily follow the same patterns as their neighbors. This stresses the need for management plans or conservation interventions (if they are deemed necessary) to carefully consider the particulars of each individual community without assuming it will follow a similar pattern to communities nearby. In addition, our analysis suggests that the mean biomass of bird prey hunted does not decrease as settlements get older or larger. Assuming that the age and size of communities are good proxies for hunting pressure (which may not be the case), better data are needed to understand how mammal, bird, reptile, and fish hunting interact with each other, in order to assess whether our observed pattern is indicative of sustainable practices or likely to be caused by diminishing returns of one prey type (birds) being compensated for by an increased focus on other taxa such as reptiles, mammals, and fish.
Footnotes
Acknowledgments
The authors are thankful to the authors of the studies included in this paper and to two anonymous reviewers for their valuable comments on the manuscript.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Stafford is funded by a Natural Environment Research Council (NERC) studentship (grant numbers NE/J500057/1 & NE/L501591/1).
Supplementary Material
Supplementary material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
