Regional patterns of wild animal hunting in African tropical forests

Main
Biodiversity loss is occurring at an unprecedented rate globally1. Hunting has been identified as a major pressure on wild animals, and unsustainable hunting (overexploitation) can drive species towards extirpation or extinction2,3,4,5. The loss of species and populations (defaunation6) can lead to seemingly ‘empty’ habitats7,8,9 and cascading effects on ecosystem services and functions10. Hunting, consumption and trade of wild animals are pervasive throughout much of the world, particularly in tropical regions11,12. Hunted animals contribute to village subsistence economies and offer revenues for those involved in trade. Declines in wild animal populations caused by overexploitation can thus trigger increased food and economic insecurity among people relying on meat and other body parts of wild animals13,14 (hereafter ‘wild meat’ refers to the meat of wild terrestrial and aquatic animals, excluding fish12) and can lead to increased land conversion to agriculture to provide new sources of food and income15. Furthermore, access to wild meat plays important socio-cultural roles in some human communities16,17,18, which declines in wild animals may negatively impact. Ensuring that all extraction of wild meat is sustainable is a key component of achieving several of the United Nations Sustainable Development Goals (for example, Goal 2 (Zero Hunger) and Goal 15 (Life on Land)12,19) and the Convention on Biological Diversity Global Biodiversity Framework targets (particularly Targets 5 and 9 (ref. 20)) and for slowing or even reversing the loss of biodiversity21.
Despite decades of research and increasing interest in global policy22, the governance and management of wild animal hunting and trade remain substantial global challenges, with little evidence of interventions that have achieved sustainable and equitable wild meat hunting in tropical regions12. These challenges are perceived to be severe in sub-Saharan Africa, where current impacts on wildlife from overexploitation are high and likely to grow given the expected doubling of the human population between 2019 and 205023. West and Central African countries, which contain some of the most biodiverse ecosystems on Earth, are at risk of wildlife overexploitation given the current high levels of hunting pressure for wild meat, ultimately due to the need for food and income22,24. Over the past few decades, hunting practices have evolved, with traditional methods (for example, spears and nets) largely being replaced or accompanied by technologies that improve the efficiency of individual hunting such as wire snares25,26, guns27, mobile phones22 and LED flashlights28. Given these rapidly changing pressures, there is a need to understand the spatial and temporal patterns of wild meat extraction to guide transitions to sustainable food systems at scale29.
For decades, many researchers have quantified the number of wild animals hunted for meat during a given period (that is, offtake) at individual sites across West and Central Africa, predominantly to understand local hunting patterns and drivers30,31. By collating and analysing these studies, we aim to better understand the scale, magnitude and patterns of wild meat offtake across the region. Past efforts to investigate wild meat offtake patterns by analysing site-level hunting studies14,32,33,34,35 have been limited by the number of sites that were included (14–33 sites, excluding hunting camps), survey effort (including very low effort, ≤10 days), data compatibility (including opportunistic sampling) and the analytical methods used (not accounting for study differences). Here we collate and analyse data from 83 studies, representing 115 settlements, to investigate the socio-cultural, economic and landscape variables associated with key components of wild meat offtake, using a Bayesian modelling approach that accounts for differences among studies. We investigate factors associated with variation in (1) wild meat offtake rate per hunter (mean kg per hunter per day), (2) the proportion of hunted animals that are sold (as a proxy for the use of wild meat for income versus food within hunter households), (3) the proportion of all hunted animals killed by gun hunting (a measure of hunting technology use) and (4) four taxonomic composition metrics of the offtake (as proxies for wild animal depletion). Throughout, we also investigate the influence of hunting technology and survey effort (measured as both the number of hunters surveyed and the number of monitoring days) on our results. By improving our understanding of the scale and patterns of wild meat hunting, our results can inform relevant sustainability policy in the region, especially by providing insights for monitoring and intervention design.
Results
Summary of studies
In total, we collated data from 83 hunter offtake studies conducted in 78 sites across seven countries in West and Central Africa between 1991 and 2020 (Fig. 1, Extended Data Fig. 1 and Supplementary List 1). The sites accounted for data from 115 settlements; the number of settlements per site, including single villages, grouped villages and towns, ranged from 1 to 14. The studies monitored a mean of 23 ± 20 hunters for a mean of 178 ± 173 days (mean ± s.d.). Across all studies, 85,214 individual animals from 143 unique species and 33 orders were reported to be captured, accounting for a total live wild animal mass of >837,900 kilograms hunted over 462,435 hunter-days (including elephants). Mammals were the most frequently hunted taxa, accounting for a mean of 93 ± 5% of the total number of animals hunted in each study, with the remainder comprising reptiles (4 ± 4%), birds (2 ± 2%) and amphibians (>0 ± 0%). Across all studies, hunters sold a mean of 57 ± 17% of the animals hunted (range, 17–94%). The studies were biased towards areas with greater forest cover: the mean forest cover within 20 km of sites was 90 ± 8% (Fig. 1), whereas the mean forest cover across the study region was 50 ± 33%. The study sites were located closer to protected area (PA) boundaries (35 ± 30 km away) compared with the mean across the study region of 48 km (excluding points inside PAs), and in areas with lower human population density (7 ± 13 individuals per km2) than the mean across the study region of 40 ± 86 individuals per km2 in 2020 (excluding urban areas).

Map of West and Central Africa, showing forest cover (green shading) and sites where hunting data were collected (white circles). Circle edge colours indicate the proportion of hunters monitored that were village hunters (with a gradient from orange to purple representing only forest hunter-gatherers to only village hunters, respectively). The data were collated from seven countries in West and Central Africa and represent hunters from 115 settlements at 78 sites. Forest cover data are from Hansen/UMD/Google/USGS/NASA119. The country outlines were obtained using the wrld_simpl data in the maptools R package105. Three-digit country codes denote countries with studies available: Cameroon (CMR), Central African Republic (CAF), Democratic Republic of the Congo (COD), Equatorial Guinea (GNQ), Gabon (GAB), Liberia (LBR) and the Republic of the Congo (COG).
Modelling offtake and composition
Our analyses showed that studies with more hunters surveyed were associated with lower estimates of offtake (kg per hunter per day; Fig. 2, Extended Data Fig. 2a and Supplementary Table 1). Within the range of data in the model, studies that surveyed 84 hunters produced mean per-hunter offtake estimates less than one sixth the size of those that surveyed 2 hunters, when all other variables were held at means or medians. We found that studies with larger mean daily hunter offtakes were associated with a greater proportion of animals sold (that is, hunting for income)—a change from 17% to 94% of animals sold was associated with a more-than-four-times increase in mean daily hunter offtake (Fig. 3a). Mean daily hunter offtake was greater in areas with better forest condition (that is, lower forest degradation from infrastructure, agriculture or deforestation), such that the average hunter removed over six times more wild meat per day at sites with healthy forests (10/10 forest intactness index) than at sites with degraded forests (3.7/10; Fig. 3b). Mean daily hunter offtakes from studies that monitored only village hunters (traditionally agricultural rural communities—for example, many Bantu peoples) were nearly double those of studies that monitored only traditionally mobile forest hunter-gatherers (for example, Aka and Baka; Fig. 3c). Furthermore, hunter offtakes more than doubled farther away from towns (population >10,000 people) across the range of our data (Fig. 3d). We also found marginal support (89% uncertainty interval (UI)) for mean daily hunter offtakes being greater closer to PAs. While not quite at the marginal 89% level of support, our model showed some evidence suggesting that lower levels of individual hunter offtake per day may be correlated with higher levels of subnational human development (estimate, −0.13; lower UI bound, −0.28; upper UI bound, 0.01).

Bayesian generalized linear multilevel model standardized coefficient estimates and posterior density curves across the following response variables: mean daily hunter offtake (OFFTAKE; kg per hunter per day), the proportion of total animals that were sold (SOLD), the proportion of animals hunted by gun (GUN), species richness (RICHNESS), the ratio of larger to smaller duikers (DUIKER RATIO), the ratio of ungulates to rodents (U:R RATIO) and the proportion of total animals that were primates (PRIMATES). Possible explanatory variables, on the basis of a priori hypotheses, are the number of hunters surveyed (Hunters), the Subnational Human Development Index (SHDI), the distance to a protected area (PA), the number of monitoring days (Days), the proportion of animals hunted by gun (Gun), the mean human population density within 20 km (HPD), the travel time to the nearest town >10,000 people in minutes (Town), the proportion of monitored hunters who were village hunters rather than traditional forest hunter-gatherers (Village), the mean forest condition within 20 km (Forest; index from poor condition (0) to good (10)), the proportion of animals that were sold (Sold), the year (Year) and an interaction between year and the proportion of village hunters (Year:Village). For each variable, the vertical line shows the mean estimate, the outer light purple line represents the 95% uncertainty interval (UI) and the darker purple shaded area represents the 89% UI. Possible variable associations are not well supported in the model when the UI overlaps with 0 (the vertical light purple line). The variables are ordered by mean coefficient estimate within each model, from negative to positive.

a–f, Bayesian generalized linear multilevel model predicted effects of an association between hunter offtake and the proportion of animals sold (a), forest condition (b), the proportion of village hunters (c) and the travel time from a site to the nearest town with a population >10,000 (d). The proportion of offtake that is sold is associated with year (e) and the proportion of village hunters (f). The purple ribbons show the 89% (dark purple) and 95% UIs (light purple), while the black line shows the global average marginal effect. The points show individual studies (scaled by the number of hunters surveyed) and are coloured according to the proportion of hunters monitored that were village hunters (with a gradient from orange to purple representing only forest hunter-gatherers to only village hunters, respectively). The points are semi-transparent to show point density.
Our model showed that the mean proportion of animals sold within studies increased from 34% to 72% of the total offtake between 1991 and 2020 (Fig. 3e). The proportion sold was also shown to be greater among village hunters than among forest hunter-gatherers—a mean of 61% of the total catch was sold in studies that monitored only village hunters compared with 42% in studies that monitored only forest hunter-gatherers (Fig. 3f). The model showed marginal support for a positive association between the number of monitoring days and the proportion sold.
We found that the proportion of animals killed by guns was twice as high in studies of village hunters than in those of forest hunter-gatherers. Our model provides marginal statistical support for an increase in the proportion of animals hunted by guns over time, as well as marginal support for an interaction between time and the proportion of monitored hunters that were village hunters (Fig. 2). Specifically, the increasing trend in the proportion of animals killed by guns over time was most evident in studies where all the hunters monitored were village hunters, in which it increased from 25% to 76% between 1991 and 2020 (Fig. 4). Our model also marginally supported a negative association with the number of hunters surveyed and a positive association between gun-hunting and the proportion of animals that were sold.

Top, Bayesian generalized linear multilevel model predicted changes in the proportion of the total number of individual animals caught that were hunted by gun over time, plotted as the interaction term when the proportion of hunters monitored that were village hunters is 1 (that is, purple points only). Bottom, modelled relationships of the proportion of animals hunted by gun with the ungulate:rodent ratio (left) and with the proportion of animals hunted that were primates (right). The purple ribbons show the 89% (dark purple) and 95% UIs (light purple), while the black line shows the global average marginal effect. The points show individual studies (scaled by the number of hunters surveyed) and are coloured according to the proportion of hunters monitored that were village hunters (with a gradient from orange to purple representing only forest hunter-gatherers to only village hunters). The points are semi-transparent to show point density. Credit: animal silhouettes, PhyloPic.org under a Creative Commons license CC0 1.0.
The taxonomic composition of the hunter offtake was associated with several variables (Fig. 2), depending on the metric under consideration. We found that hunted species richness more than doubled with the number of monitoring days, from 20 to 50 species over the range of the data (Extended Data Fig. 2b), holding other variables constant. The ratio of larger- to smaller-bodied duikers was positively associated with travel time to towns (with a fourfold increase in this duiker ratio as travel time increased from 0 to 15 hours; Extended Data Fig. 3) and was marginally negatively associated with the proportion of animals hunted by gun. The ratio of ungulates to rodents increased from 2.5 to 17.1 with an increasing proportion of animals hunted by guns (from 0% to 99%; Fig. 4) and from 0.33 to 17.1 with increasing forest condition index. Finally, the proportion of primates in the offtake nearly tripled to 20% as the proportion of animals hunted by gun increased (Fig. 4).
Discussion
We analysed factors associated with regional hunting patterns across African tropical forests, which previously had been explored rigorously at only the local site level36,37. Using a large compilation of African hunting data, our analyses provide insights into regional patterns of the scale of wild meat hunting, possible wildlife depletion and effects of survey effort, which together have implications for wild meat monitoring, management and policy over large spatial scales.
Socio-economic variables associated with hunter offtake
Our results show that the proportion of offtake sold is positively correlated with daily hunter offtake levels. Mean daily hunter offtakes were also positively correlated with the proportion of hunters who were village hunters, as opposed to traditional forest hunter-gatherers. Village hunters typically sold a greater proportion of their offtake than traditional hunter-gatherers. These results could be because village hunters (probably including hunters coming from elsewhere specifically to hunt) tend to be more integrated into the market economy than traditional forest peoples, and have better access to modern sales technologies such as mobile phones and mobile money25,38,39. Furthermore, though we acknowledge that our dataset is not longitudinal across the same hunters/sites, our model suggests that the proportion of animals sold has increased over time. This result could be because the demand for wild meat, particularly in urban areas, also drives the proliferation of commercial hunters and encourages some subsistence hunters to hunt for income22. These results are probably comparable to those in other parts of the tropics such as in African savannahs where wild meat is sold for income8, and in East, South and Southeast Asia40 and Latin America22 where wild meat is openly sold in markets. As urban areas are a major destination for traded wild meat12, initiatives to reduce consumer demand and the price of wild meat in urban areas could disincentivize commercial hunting and in turn ease pressure on wildlife41. However, given that one outcome of reducing urban demand would be to reduce the immediate economic viability of hunting in rural areas, campaigns targeting urban areas will probably need to be conducted together with initiatives to develop other, more sustainable sources of rural income42. Where rural communities have legal rights to hunt wildlife, a substantial challenge to local hunting management and sustainability could be that communities have no legal means of excluding external commercial hunters, who may target areas with intact forest and/or close to PAs.
While more data are needed to understand whether there is a concrete relationship between hunting and development, development assistance focusing on building resilient local economies that are not dependent on wild meat hunting (for example, by raising health, standard of living and education levels), coupled with increased access to and diversification in employment and alternative proteins, may reduce hunting pressure over time. This will depend on trends in urbanization, consumer preferences, wealth and human population density. However, future development could lead to unintended impacts such as increased ecosystem degradation due to increasing agriculture, species declines43, industrial expansion, logging roads and potentially new markets for wild meat if it becomes a luxury commodity44.
Landscape variables associated with hunter offtake
At the regional scale, we found greater mean daily hunter offtake in areas with better forest condition. This result could be caused by hunters in areas with better forest condition hunting more because they have fewer alternative sources of meat and income and/or because areas with better forest condition are less depleted of wildlife, so hunters can hunt larger-bodied species. We found that mean daily hunter offtakes were greater in more remotely located sites and closer to PAs (although the latter relationship was marginally significant), mirroring the local-scale patterns for greater levels of hunting45 and wild meat consumption recorded in East Africa46. This suggests that remote areas or PAs may have a role as refuges and source areas for wild animals47, and/or that PAs are established in more remote intact landscapes where hunters may rely more on wild meat for food and income, and/or that these areas have yet to reach the depleted state of more accessible sites. Our results underscore the need to sustainably manage functional forest ecosystems and PAs, which provide local people with potentially renewable natural resources.
Implications of gun hunting
Our results suggest that village hunters are increasingly using guns to hunt animals rather than other methods, albeit at marginal statistical significance (see ‘Limitations’). This could result from multiple factors, including increased access to guns by village hunters26,48,49, increased access to aerial and arboreal animals that guns provide to hunters, and the increased selectivity of guns relative to traps, making guns the preferred weapon for those who hunt for income26. Furthermore, our results suggest that where hunters use guns proportionally more than other methods, they also sell a greater proportion of the offtake. This association could be because guns are the weapon of choice for commercial hunters, so that they can target commercially viable species, and because guns require financial investment, which is more likely to be possible for commercial hunters.
Increasing gun use over time is likely to impact fauna in different ways. Our results show that when hunters catch a greater proportion of animals by gun, the proportion of primates in the catch and the ungulate:rodent ratio increase. Consumers prefer ungulates and primates in some countries, and these taxa are commonly offered for sale at urban wild meat markets50,51,52. While gun hunting is usually more selective than trapping and could therefore be considered better for wildlife generally (for example, fewer carnivores are killed with guns than by traps53, and fewer animals are wasted54), increased gun hunting has implications for the conservation and management of ungulates and arboreal and aerial animals55. Primates, for example, may be particularly vulnerable to overhunting using guns and can be the target of specialized primate hunters given their high value56. In many cases, gun hunting without the correct permits (gun permit and game quotas) is already illegal, but quotas are extremely difficult to enforce, and hunting law is in general inadequately enforced. The scale of hunting for income, as well as gun proliferation, self-manufacture and use, means that purely enforcement-based management is likely to be challenging, especially in post-conflict areas55 and in areas where the rule of law is weak. Efforts to ensure that wildlife use is sustainable should simultaneously focus on enforcing laws that regulate the hunting of threatened species, including ensuring their survival in PAs, while also investigating and co-developing ways for local communities to manage wildlife sustainably, where possible57. Shifting the focus to spatial hunting management (for example, the rotation of hunting zones over time) may be more locally appropriate and resilient to changes within the hunting system (for example, changes in hunting technology).
Proxy-based evidence of depletion
Our results suggest that hunted fauna may be more depleted in areas more accessible from towns and in areas of worse forest condition. For example, the duiker ratio was lower in areas that are more accessible from towns, possibly indicating that the larger Cephalophus duikers have become more depleted in these areas58,59. This is corroborated by our findings that daily hunter offtake and the ungulate:rodent ratio were greater in areas with better forest condition. In depleted areas, we cannot discern whether the offtake has reached a post-depletion sustainability state in which fast-reproducing species are being hunted sustainably60,61. Together, and assuming that offtakes reflect the local abundance of wild animals to some extent55,62,63,64,65, these results underscore the role of remote areas with good forest condition as refuges for wildlife. Furthermore, Forest Stewardship Council (FSC)-certified logging concessions in Central Africa have been shown to hold much higher numbers of large-bodied animals than non-FSC certified concessions66 and thus are of great importance for sustainable hunting management across much of forested Africa67. Given that poor forest condition is characterized by greater levels of infrastructure, agriculture and deforestation, the effects of hunting pressure and forest degradation may be conflated in our models. However, this would not change the premise that remote areas serve as refuges for wildlife and possibly source populations for hunted areas.
Limitations
Our results have the following limitations. First, due to the disparate nature of our dataset, our results largely focus on past drivers and hunting patterns. In the future, hunting patterns may be influenced by other factors, which would make it challenging to predict them from this dataset. However, we suspect that many factors will remain important in the future, such as (1) accessibility and distance to PAs, given that they will probably remain an important refuge for wildlife; (2) forest integrity, due to the importance of source–sink dynamics and the cumulative effect of multiple stressors on wildlife; and (3) the proportion of wild meat sold, given increasing access to markets and the cash-based economy. Furthermore, it was not possible to include important contextual information that was not reported in our data sources, such as hunter typology (for example, commercial and/or migrant hunters travelling long distances), use of hunting camps, indicators of wealth or prosperity, governance, regulation, conflicts and PA effectiveness. Another limitation is that, given the absence of time-series studies at multiple locations across the region68, our analyses use space-for-time substitution out of necessity and therefore assume that comparisons across space represent trends over time (we recognize that correlation does not imply causation). Finally, while our results cannot be used to determine the sustainability of hunting practices at individual sites directly (see the discussions in ref. 69), they do provide important context for developing local sustainable management and monitoring plans. Given the now substantial catalogue of evidence regarding the drivers and impacts of hunting achieved over 20 years of research12, we urge researchers and practitioners to focus their efforts on identifying and implementing ways to sustainably manage wildlife.
Research, policy and action
The analyses presented here are an important step towards better understanding some of the factors associated with offtake patterns and could be used to guide policy, monitoring and management strategies. For example, we identified that a substantial proportion of offtake is destined for trade, which suggests a role for policies and actions to reduce demand in urban areas and to regulate wild meat supply chains. However, fundamental moral and ethical implications should be considered when designing related national policies or local interventions. We suggest that where local communities rely on wild meat for nutrition, and where hunting builds resilience in the local economy and remains an important component of cultural identity, management efforts should focus on attaining sustainable offtake under local governance, protecting threatened species and developing an effective hunting monitoring framework. Furthermore, given that we highlight how gun hunting may disproportionately impact primates, research is needed to conduct sustainability and population assessments of hunted primate species, followed by appropriate management and policy. We present evidence of wild animal depletion in areas accessible from towns and in areas with poor forest condition, which emphasizes the importance of active conservation efforts that protect and restore functioning ecosystems and that provide people with important natural resources.
Our results reflect hunting patterns in forested areas of West and Central Africa. In our analyses, the number of studies in several large countries (including the Democratic Republic of the Congo and the Central African Republic), in non-forest areas and in large areas of West Africa was limited. This means we cannot be certain about hunting patterns in these areas, so targeting future research in these areas will be valuable to fill these knowledge gaps. We highlight that longer-term monitoring of hunting activities, including ensuring representative samples of hunters, is needed to ensure robust conclusions about hunting patterns, species richness, sustainability and intervention effectiveness. While our results and others64,65,70 suggest that the duiker ratio and the ungulate:rodent ratio may be useful proxies for wildlife availability and abundance, our results also suggest that the indicators are sensitive to changes in hunting technology (for example, switches from snaring to gun hunting), which needs to be considered when interpreting these indicators. Research across several biomes using standardized survey methods and indicators (for example, see www.wildmeat.org), where survey effort is incorporated into survey design, monitoring and evaluation, needs to be undertaken. The development of a transparent and rigorous regional wild meat hunting, consumption and trade monitoring network would provide standardized data over time that could be used to inform and monitor the impacts of national and regional strategies and policies for wild meat management and governance.
Together, our findings provide information at a scale that is relevant for decision-making across the region. Our study shows the strong need for robust monitoring and management frameworks for hunters and traders, which are currently lacking in most West and Central African countries and in other parts of the tropics. Without monitoring of hunted species populations and hunter offtakes, the sustainability of hunting systems remains unknown. Such frameworks could also be supported through organizations such as the Central African Forests Commission (COMIFAC), an intergovernmental organization tasked with managing Central African forest sustainably. Assessing sustainability is key to monitoring progress towards Kunming-Montreal Global Biodiversity Framework Targets 5 and 9 on ensuring the sustainable use and management of wildlife by 2030, adopted by 196 countries in December 2022 (ref. 20). In conjunction with the necessary resources and suitable wild meat policies and legislation12, strong political and public will is crucial if such goals towards a more sustainable and just future for people and wildlife are to be achieved globally.
Methods
Wild meat data collation and inclusion criteria
Data on the number of animals caught by hunters in a given location (hereafter ‘hunter offtakes’) across West and Central Africa were extracted from published journal articles, grey literature (for example, non-governmental organization reports and university degree theses) and unpublished data (deemed suitable when they had been collected using published methods). Using the following methods, we searched for potential data sources between 1 August 2019 and 17 December 2021. Suitable sources were first identified from lists provided in two studies30,31 that used systematic searching techniques. We then conducted regular additional literature searches using search engines and bibliographic databases (Google Scholar (including using alerts) and the ISI Web of Knowledge (All Databases)). We performed searches using Boolean logic with the following words: (bushmeat, game, gibier, hunt*, offtake, poach*, viande de brousse, wild meat, wildmeat) AND (Africa, Afrique centrale, Afrique de l’Ouest, Afrique occidentale, Bassin du Congo, Central Africa, Congo Basin, West Africa). While this method was effective at finding sources, we were able to find further sources mostly through subsequent searching of reference lists and contacting wild meat researchers (both akin to chain-referral sampling71). We did not exclude sources on the basis of the year or biome within which data were collected, or the ethnic group of the hunters monitored.
Sources may contain one study or more, in which hunter offtakes were monitored at sites (settlements including villages and towns), over a known period. For seven studies that monitored hunting at more than one settlement in a limited geographic area, data were unavailable for each settlement separately; thus, some sites represent more than one settlement. Hunter offtake data were deemed suitable if (1) the number of individual wild animals hunted per taxonomic group was reported (or estimated, as in one study72); (2) the hunters monitored were deemed representative of hunters at a given site (that is, the source or the data owner did not mention monitoring a specified minority group of hunters—however, we could not determine from studies whether the communities that were monitored were chosen randomly); (3) the survey method was reported; (4) the coordinates of the site(s), or a map from which the coordinates could be acquired by georeferencing against Google Earth, were available; (5) the source included information on the number of days that hunting was surveyed and the number of hunters whose offtakes were monitored; (6) the settlement at which the hunting was monitored was a permanent dwelling (that is, not a hunting camp); and (7) the hunting method and the use of the hunted animals (that is, sold, eaten or other) were reported. We either extracted hunter offtakes directly from the sources or contacted the authors for their raw data.
Wild meat data curation
Species names were extracted directly from the sources when reported. When animals were not identified to the species level, we allocated them to the most resolved taxonomic level possible—for example, if an animal was identified as a ‘genet’, it was allocated to Genetta. Taxonomic identity was spell-checked across all sources and harmonized against the IUCN Red List (www.iucnredlist.org), and each unique taxonomic identity was then given a unique numeric identifier. So that we could convert the number of animals hunted into the total biomass extracted (see below), we also allocated each species a body mass. Given that sources rarely contained data on the mass of each carcass, we instead matched the average adult body mass for each species from Myhrvold et al.73. For Dendroaspis, Naja and Python, we matched animals to adult body mass data reported in Abugiche74. When animals were identified to a lower taxonomic resolution than species (for example, genus), we calculated and allocated the mean body mass of those species within that taxonomic unit present in our dataset. On rare occasions when an animal was identified to a taxonomic unit for which no other animals had been identified to species level within that unit in our dataset, we allocated the mean body mass of all species in that unit from Myhrvold et al.73.
Individual hunter offtakes typically vary considerably within a settlement75. To compare offtakes across time, space and study, we calculated a standardized mean hunter offtake rate (per day) for each study–site combination. Offtake rates were standardized by first multiplying the number of animals caught per species by their respective body mass and then summing the species biomasses within each study–site. As many studies did not have offtake data available for each hunter separately, we calculated the offtake variable used in our analyses—that is, mean hunter offtake per day (kg per hunter per day)—by dividing the total study–site biomasses by both the number of hunters that were continuously monitored and the number of days they were surveyed in each study. Elephants (Elephantidae) rarely appeared in hunting records given that elephant hunting is usually a specialist, illegal activity typically requiring expensive high-velocity rifles and as such is not conducted by most wild meat hunters76. Additionally, their enormous size means that elephant records disproportionately affect mean offtakes. To avoid overestimating typical hunter offtakes, we removed all elephants before calculating mean hunter offtake per day64.
The studies varied in the types of hunters monitored. While most studies solely monitored village hunters from typically agricultural rural communities (63 studies; most commonly Bantu peoples), 17 studies included or solely monitored traditional forest hunter-gatherers—that is, those whose subsistence depends on a forager and hunter-gatherer lifestyle (for example, Aka and Baka). For each study, we calculated the proportion of hunters monitored that were village hunters: a study that monitored only forest hunter-gatherers would have a proportion of 0, while a study that monitored only village hunters would have a proportion of 1. This proportion was then used directly as a variable in the analyses. For three studies located in the same area, it was unclear whether any forest hunter-gatherers were monitored, so we assigned these the mean proportion across all other studies.
The composition of offtake has been shown to reflect the diversity and relative abundance of hunted wild animals in an area55,62,63,64,65,70. We calculated four taxonomic composition metrics of the offtake as proxies for wild animal depletion. First, for each study, we calculated the minimum species richness in the hunter offtake, measured as the count of unique species hunted. In several studies, not all animals could be identified at the species level, but they could be identified to other taxonomic resolutions. Animals not identified to the species level often belonged to groups that are hunted relatively infrequently and were less likely to be already represented in a count of species richness based solely on animals identified to the species level. Therefore, to reduce the underestimation of species richness, we manually checked through all studies and added one to the richness count each time it was clear that individuals identified to a lower resolution were not already represented in the count of a given study. For example, if five animals were identified only as Genetta in a study, and no animals were identified to a specific species of Genetta, we added one to the study richness count.
Second, as a measure of possible depletion of larger-bodied ungulates, which are preferentially hunted in West and Central Africa and are potentially more sensitive to hunting (at least in the case of some species), we adapted a regional indicator called the “B/M ratio”58. Using this indicator, we explored the suggestion that the ratio of small, fast-breeding blue duikers (Philantomba monticola; B) to slower-breeding, medium-sized red duikers (Cephalophus spp., not including the larger-bodied C. silvicultor; M) in the offtake is indicative of the status of the fauna in a given area58. The ratio is based on the observation that larger, slower-reproducing duiker species are depleted first as hunting pressure increases, while smaller, faster-reproducing duikers increase in density under the same regime until hunting pressure becomes too high and they are too depleted58,70,77. This was supported in another study, which found that the proportion of red duikers was greater in more intact forests in Central Africa59. We adapted the B/M ratio in the following ways: (1) as our studies encompassed both West and Central Africa and blue duikers are not present in most of West Africa, we expanded the ratio to the genus level to include the highly similar Philantomba species from West Africa (P. maxwellii); and (2) so that a high ratio may be indicative of better faunal status of large-bodied animals sensitive to hunting pressure, we expressed our results as the M/B ratio (including both medium- and large-bodied duikers (that is, C. silvicultor) in the Cephalophus genus), with a larger value thus indicating greater relative abundance of larger-bodied duikers.
Third, we calculated the ratio of ungulates (here even-toed ungulates including deer, antelopes and pigs; hereafter ‘ungulates’) to rodents among all animals hunted, which has often been considered a metric of wild animal depletion in the tropics32,70, as increases in rodent offtake can coincide with the depletion of preferred larger-bodied animals65,78. The ungulate:rodent ratio was calculated as the number of ungulates divided by the number of rodents caught. A larger value of this indicator represents a greater relative abundance of ungulates than rodents.
Finally, we calculated the relative proportion of primates in the total number of animals caught to investigate the impact of gun hunting on a preferred group of species.
Processing and extraction of spatial layers
Several spatial layers were acquired or calculated for use as variables in our analyses (Supplementary Table 2). For some layers, we were able to directly download appropriate spatial layers, including the SHDI v.5.0 (ref. 79) and the estimated travel time to the nearest town80. The SHDI is an index constructed to represent the dimensions of education, health and standard of living and is created from expected years of schooling, mean years of schooling, life expectancy and gross national income per capita79. For the SHDI, we extracted the identity of each subregion from the SHDI layer for each site and then matched the index value to the same year as the study. For travel time to the nearest town, we used the 2015 layer80 to calculate the travel time from each grid cell to the nearest settlement with 10,000+ people81. Travel time estimates consider movement by foot, roads, railways, rivers, water bodies, land cover class, slope, elevation and national borders80. We then extracted travel times for each of our sites by conducting value point extractions from the recalculated spatial layer. The cut-off point of 10,000 people was selected to represent settlements from small towns upwards, where wild meat markets may be present82. While we recognize that travel times will have varied across the period for which we have data (1991–2020), we did not have access to contextual data allowing us to adjust travel times to the year of data collection.
To calculate the distance to the nearest PA for each site, we downloaded PA polygons and point locations from the World Database of Protected Areas v.1.6 (https://www.protectedplanet.net/en). We excluded Ramsar and UNESCO Man and Biosphere sites because they do not necessarily have protective legislation (apart from any portions in overlapping nationally protected areas, which would be retained in our dataset). Then for each site, we calculated the Euclidean distance from the site to the nearest PA boundary. This underestimates the actual distance for a hunter or trader travelling on paths or roads83, but for analysis at a regional scale, we deemed this an adequate approximation given the lack of site-level information on the most direct paths/routes. Mean distance from a PA was estimated on the basis of the Euclidean distance of randomly selected points within a 1 km2 grid overlaid in the study area (N = 1 × 106). As all sites were in forested areas, we included a measure of forest condition in our models. We used the forest integrity index (‘forest condition’)84, which is constructed as an index (ranging from 0 to 10) based on observed human pressures (infrastructure, agriculture and tree cover loss), proximity to observed human pressures (inferred) and change in forest connectivity84. For each site, we calculated the mean forest condition within a 20 km radius; this radius was chosen to represent the furthest distance subsistence hunters are likely to travel24 and to include the area probably able to act as the source of most hunted species. To estimate the human population density around each site, we first downloaded all human population density layers from WorldPop v.2.0, available every five years between 2000 and 2020. For each site, we then allocated the human population density layer for the year nearest to the study and calculated the mean human population density within 20 km, to match the forest layers.
All spatial layers were reprojected into an Albers equal-area conic projection to conserve distance, and if not already at a 1 km grid-cell resolution, we resampled layers to a 1 km × 1 km grid-cell size so that all layers had the same projection, extent and cell size (see Supplementary Table 2 for the native resolution). For each study site, we extracted the point-level estimates of spatial variables. Sites that represented more than one settlement (where data could not be separated per settlement) were all located close to each other, so for spatial extractions for these sites, we used the coordinates of the settlement that had the highest offtake where possible or the centre point of the surveyed settlements. All spatial extractions and manipulations were conducted in R version 4.1.1 (ref. 85).
Modelling offtake and composition
To examine the variables associated with hunter offtake, use and hunting methods, and taxonomic composition, we constructed a series of Bayesian generalized linear multilevel models. We first checked for correlations between possible variables and did not include variables in the same models where Pearson’s r ≥ 0.7 (ref. 86). We used variance inflation factors to check the extent to which variables were influenced by their correlation with the other independent variables and found all to have limited influence (variance inflation factor < 2). All of the following models included two measures of survey effort as variables: the number of hunters surveyed and the number of monitoring days. For hunter offtake, we constructed a core Bayesian generalized linear multilevel structural equation model87, which included three primary response variables: offtake as mean hunter offtake per day (kg per hunter per day; see above for how this was derived), the proportion of the total number of animals hunted that were sold and the proportion of the total number of animals that were hunted by gun. Different combinations of variables were selected for different responses on the basis of their a priori hypothesized relationships (Supplementary Tables 3 and 4): the travel time to town, the forest condition within 20 km, the distance to the nearest PA, the year of study, the SHDI, the proportion of hunters monitored that were village hunters, the proportion of animals hunted by gun and the proportion of animals hunted that were sold. For the proportion-hunted-by-gun model, we included an interaction term between the year of study and the proportion of hunters monitored that were village hunters. Four further response variables were included as a multivariate response with offtake: species richness, the duiker ratio, the ungulate-to-rodent ratio and the proportion of hunted animals that were primates. All four taxonomic composition metric models included the same variables: the travel time to town, the forest condition within 20 km, the distance to the nearest PA and the proportion of animals hunted by gun. To standardize the comparison of the effects of variables, all values of continuous variables were centred by subtracting the mean and scaled by dividing by the standard deviation.
Offtake, duiker ratio and ungulate-to-rodent ratio were modelled with gamma distributions (log link), while species richness was modelled with a negative binomial distribution (log link). Prior to modelling, we added 1 to the number of rodents in one study where the original value was 0 so that we could calculate and subsequently model the ungulate-to-rodent ratio. The number of animals that were primates was modelled using a beta binomial distribution (logit link), accounting for overdispersion and the differences in study sample sizes through the response variable structure in the model (coded as successes|trials)88. The proportions of animals hunted by gun and the proportions that were sold were extracted from studies and modelled with beta distributions (logit link). Before modelling, we added 0.001 to the proportion-of-gun-hunting variable for one study where the original value was 0, which a beta regression could not model88. For all models, we included varying intercepts (but not slopes) for each unique reference ID (modelled as correlated across formulas), within which studies and sites are linked to account for variation due to differences in the methods employed by researchers. Given that multisite studies were often conducted in the same area, this varying intercept may partially account for possible spatial autocorrelation within studies; however, we conducted further distance-based autocorrelation checks (see below). We did not include varying intercepts for each country due to the low number of studies in some countries, which meant that country accounted for the same variation as the study reference in some cases. The weak explanatory power of countries in the data was further shown in preliminary analyses, in which including country caused convergence issues.
We fitted the models in the Stan probabilistic programming language89 through the R statistical environment (v.4.1.1)85 using the package brms90,91. The models ran four chains, each with 6,000 iterations, of which the first 1,500 were used to calibrate the sampler, leading to a total of 18,000 posterior samples. We used weakly informative normal (0,5) priors for all population-level parameters estimated. We assessed model convergence by visually assessing trace plots for sufficient mixing and by checking the Gelman–Rubin convergence diagnostic, (hat{R}) (ref. 92); in all our models, (hat{R}) values for the estimated parameters were 1.0, indicating adequate model convergence. Residual checks and distance-based autocorrelation using Moran’s I were assessed using the DHARMa and DHARMa.helpers packages (Supplementary Table 4)93,94,95. To calculate Moran’s I, the residuals were first aggregated at the level of each unique site. Model performance and fit were checked through graphical diagnostic checks of the predictive distribution in the R package bayesplot using standardized protocols96,97 (Supplementary Table 5), and we observed a close correspondence between the posterior predictive distributions and the observed data. The posterior means and UIs were used to interpret the strength and uncertainty of the estimated effects91. We considered variables to be meaningful when the 89% UI of their effect sizes did not overlap zero (marginal support) and strongly supported when their 95% UI did not overlap zero98,99. To further check the model effects, we also estimated the probability of direction using the bayestestR package, which represents “the certainty with which an effect goes in a particular direction”100. All graphs of model predictions were produced using the expected values, holding other variables at means or medians depending on the variable. Several R packages were used for spatial mapping, extraction and analysis (raster, rgdal, sf, gdalUtils and maptools101,102,103,104,105); correlations (GGally106); data manipulation (reshape2, stringr, broom.mixed, tidyr, dplyr, tidyverse and tidybayes107,108,109,110,111,112,113); and plotting (ggplot2, patchwork, cowplot, gridExtra and rphylopic114,115,116,117,118).
While the appropriateness of weighting models is subject to academic debate, we tested the influence of a version of the models that was weighted by one metric of survey effort, the number of hunters surveyed. The weighted model further accounted for differences in survey effort in our models by weighting the likelihood contributions of each observation by a scaled version of the number of hunters surveyed. To avoid over- or underweighting each observation, we rescaled the number of hunters surveyed so that all values fell between 1 and 5, given that weights in the brms package are taken literally—that is, an observation with weight 2 receives 2 times more weight than an observation with weight 1, essentially allowing one to put more importance on some data points. The weighting was not applied to the models where the response was the proportion of animals that were primates because these were already adjusted through their model structure (successes|trials). The scaled weighting variable was calculated as follows:
where a and b are the range of the scaled variable (here 1 and 5), respectively, and x is the survey effort (the number of hunters surveyed). The weighted model results and comparison are presented in Supplementary Table 6.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Responses