Hostname: page-component-78c5997874-v9fdk Total loading time: 0 Render date: 2024-11-10T11:27:38.259Z Has data issue: false hasContentIssue false

Livelihood vulnerability increases human–wildlife interactions

Published online by Cambridge University Press:  27 September 2021

Joana Pereira*
Affiliation:
cE3c – Centre for Ecology, Evolution and Environmental Changes, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, 1749-016, Lisboa, Portugal
Luís Miguel Rosalino
Affiliation:
cE3c – Centre for Ecology, Evolution and Environmental Changes, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, 1749-016, Lisboa, Portugal CESAM – Centro de Estudos do Ambiente e do Mar, Departamento de Biologia, Universidade de Aveiro, Campus Universitário de Santiago, 3810-193Aveiro, Portugal
Serafino Mucova
Affiliation:
CESAM – Centro de Estudos do Ambiente e do Mar, Departamento de Biologia, Universidade de Aveiro, Campus Universitário de Santiago, 3810-193Aveiro, Portugal Faculty of Natural Sciences, Lúrio University, Pemba P.O. Box 958, Mozambique
Yasalde Massangue
Affiliation:
Faculty of Natural Sciences, Lúrio University, Pemba P.O. Box 958, Mozambique
Murchide Abdulrazak
Affiliation:
Faculty of Natural Sciences, Lúrio University, Pemba P.O. Box 958, Mozambique
Somar Vahossa
Affiliation:
Faculty of Natural Sciences, Lúrio University, Pemba P.O. Box 958, Mozambique
Mouzinho Selemane
Affiliation:
Faculty of Natural Sciences, Lúrio University, Pemba P.O. Box 958, Mozambique
Carlos Fonseca
Affiliation:
CESAM – Centro de Estudos do Ambiente e do Mar, Departamento de Biologia, Universidade de Aveiro, Campus Universitário de Santiago, 3810-193Aveiro, Portugal ForestWISE – Collaborative Laboratory for Integrated Forest & Fire Management, Quinta de Prados, 5001-801 Vila Real, Portugal
Maria João Santos
Affiliation:
Department of Geography, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
*
Author for correspondence: Joana Pereira, Email: jgopereira@fc.ul.pt
Rights & Permissions [Opens in a new window]

Summary

Human–wildlife interactions (HWIs) occur in many rural African communities, with potential impacts on livelihood vulnerability. High livelihood vulnerability may force communities to employ strategies that increase the risk of negative HWIs, yet the extent to which HWIs drive or are driven by vulnerability is unclear. We hypothesized that more vulnerable households are more likely to be exposed to wildlife and experience negative interactions. To test this hypothesis, we calculated the Livelihood Vulnerability Index (LVI) of rural households in and around Quirimbas National Park (north-eastern Mozambique) and assessed whether there is a link between livelihood vulnerability and HWIs. We found a two-way association between LVI and HWIs, with more vulnerable households indeed taking greater risks and encountering wildlife when fetching water from rivers, whereas less vulnerable households tended not to employ strategies likely to increase wildlife encounters. We also observed that HWIs exert a strong effect on livelihood vulnerability, suggesting that HWIs should be included as an exposure factor in vulnerability assessments for rural households. We recommend that livelihood strategies and community vulnerability should be considered when designing HWI mitigation schemes and implementing conservation measures.

Type
Research Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2021. Published by Cambridge University Press on behalf of Foundation for Environmental Conservation

Introduction

Human population growth can increase the vulnerability of local communities to food and water scarcity, poverty and climate change (Mondal Reference Mondal2019). Although strategies such as building a highly connected social network are very successful in reducing vulnerability (Chambers & Conway Reference Chambers and Conway1992), other strategies, such as collection of non-timber forest resources for commercial purposes, may alleviate vulnerability in one way but exacerbate it in others (Duffy et al. Reference Duffy, St John, Büscher and Brockington2016), such as by increasing the frequency of human–wildlife interactions (HWIs) (Khumalo & Yung Reference Khumalo and Yung2015). Globally, millions of people are at risk of negative HWIs, such as human fatalities or crop and livestock losses (Barua et al. Reference Barua, Bhagwat and Jadha2013). However, most vulnerability assessments have primarily focused on climate change risks (e.g., Hanh et al. Reference Hahn, Riederer and Foster2009) rather than the more immediate vulnerability risks associated with negative HWIs. Some recent studies have taken a more holistic approach to vulnerability assessments (e.g., Yadava & Sinha Reference Yadava and Sinha2020, Notelid & Ekblom Reference Notelid and Ekblom2021), including a variety of stress scenarios such as variability in crop yield and price (Jezeer et al. Reference Jezeer, Verweij, Boot, Junginger and Santos2019, Junquera & Grêt-Regamey Reference Junquera and Grêt-Regamey2020). Projected climate-driven changes in human migration and wildlife range shifts suggest that HWIs may become more frequent, so a better understanding of the relationship between HWIs and livelihood vulnerability is fundamental (Khumalo & Yung Reference Khumalo and Yung2015).

Some studies of HWIs have already explored the role of livelihood vulnerability in wildlife interactions for particular regions and social groups (Ogra Reference Ogra2008, Khumalo & Yung Reference Khumalo and Yung2015, Seoraj-Pillai & Pillay Reference Seoraj-Pillai and Pillay2017). It is widely recognized that living alongside wildlife brings many benefits and services to communities, such as through the provision of food and other resources (Cox & Gaston Reference Cox and Gaston2018), pest control (Morales-Reyes et al. Reference Morales-Reyes, Pe, Margalida, Dona and Lazcano2015), tourism revenue (Naidoo et al. Reference Naidoo, Fisher, Manica and Balmford2016) and cultural and recreational value (Bateman & Glew Reference Bateman and Glew2010). However, proximity to wildlife may also engender costs and disservices (Ceauşu et al. Reference Ceauşu, Graves, Killion, Svenning and Carter2019) through property damage (Lamichhane et al. Reference Lamichhane, Persoon, Leirs, Poudel, Subedi and Pokheral2018), competition for food and land (Treves Reference Treves2009), disease transmission (Blair & Meredith Reference Blair and Meredith2018) and human injury or loss of life (Ratnayeke et al. Reference Ratnayeke, Van Manen, Pieris and Pragash2014). These costs may differ with levels of wildlife exposure and with community vulnerability (Ogra Reference Ogra2008). Vulnerable communities may adopt livelihood strategies that escalate negative HWIs, as vulnerability can force their adoption of risky behaviours to obtain essential resources and/or to increase income (e.g., bushmeat hunting, poaching), in turn increasing wildlife exposure (Bevan Reference Bevan2000) and retaliation (Inskip et al. Reference Inskip, Ridout, Fahad, Tully, Barlow and Barlow2013). Moreover, socio-economic characteristics (e.g., age, gender, wealth, occupation, education, religion), individuals’ emotions (e.g., fear) and perceptions of risk, social norms and political inequalities (Bond Reference Bond2014, Bhatia et al. Reference Bhatia, Redpath, Suryawanshi and Mishra2019) can all contribute to escalating or alleviating HWIs. Communities and households with alternative sources of income and/or with a highly connected social network (e.g., who cooperate in crop damage prevention and mitigation (Stone et al. Reference Stone, Phalke, Warren, Av and Krishnan2019), share resources and participate in group protection) are less vulnerable to negative HWIs (Butt et al. Reference Butt, Shortridge and WinklerPrins2009).

In Africa, 60% of the population lives in rural areas, and of this, 43% live in extreme poverty (World Bank 2018) and are highly dependent on natural resources for subsistence. African wildlife populations are among the most diverse and dense in the world, and the overlap between wildlife habitat and human-occupied areas is substantial. Africa accounts for 66% of the negative HWIs reported globally (Seoraj-Pillai & Pillay Reference Seoraj-Pillai and Pillay2017). Damage by and encounters with wildlife have forced rural communities to adopt non-agricultural livelihoods to reduce their vulnerability (Gupta Reference Gupta2013). Since wildlife presence is an important driver of change for rural livelihoods (Gupta Reference Gupta2013), there is a need to better understand how livelihoods and vulnerability relate to HWIs.

In this study, we aim to understand whether and how HWIs and livelihood vulnerability are linked in a case study in Mozambique (southern Africa), an impoverished country with recent protected areas and where rural communities are living alongside wildlife (Merz et al. Reference Merz, Pienaar, Filk and Muyengwa2021). Mozambique’s plan for protection requires reduced use of protected areas as well as maintaining and improving wildlife populations and minimizing HWIs (Anderson & Pariela 2015). The country has experienced a growing incidence of HWIs as people scout the landscape for the natural resources necessary for their livelihoods (Le Bel et al. Reference Le, Chavernac, Mapuvire and Cornu2014). This likely increases potential encounters with freely roaming wildlife (Dunham et al. Reference Dunham, Ghiurghi, Cumbi and Urbano2010). We hypothesized that more vulnerable households are more likely to be exposed to wildlife and to experience negative interactions. More specifically, individuals and communities displaying higher sensitivity (the degree to which individuals and communities may be affected by climate change and HWIs) and lower adaptive capacity (a lower ability to take action to reduce exposure and sensitivity) will be more vulnerable to negative HWIs (Majale Reference Majale2002). We explore how livelihood vulnerability varies across communities and evaluate which drivers have the greatest influence on such vulnerability, its components (i.e., exposure, sensitivity and adaptive capacity) and, ultimately, HWIs. Finally, we assess how HWIs and livelihood vulnerability influence each other by examining whether HWIs are either external drivers or internal components of that vulnerability. Our approach expands upon previous vulnerability analyses by considering the interaction between HWIs and livelihood vulnerability. This approach is crucial for understanding the links between livelihood vulnerability and biodiversity conservation (García-Frapolli et al. Reference García-Frapolli, Ayala-Orozco, Oliva and Smith2018), as well as for guiding biodiversity management and sustainable development.

Methods

Study area

We studied communities living within and close to Quirimbas National Park (QNP; Province of Cabo Delgado, north-eastern Mozambique; –12°30’0 S, 39°24’0 E) (Fig. 1a). There are 153 villages in QNP, harbouring a total of c. 200 000 people, 57% of whom live within the Park, and the remainder inhabit its buffer zone (Appendix A; MITADER 2012). We also included three villages located outside QNP encompassing an estimated combined population of 8786 residents.

Fig. 1. (a) Location of Quirimbas National Park (QNP) in the Republic of Mozambique on the African continent. (b) Sampled villages within and close to QNP.

Data collection

Sampling design

We sampled villages inside the Park (n = 9), in the buffer zone (n = 2) and bordering or outside the Park (n = 3) (Fig. 1b). We used structured questionnaires (n = 224) that included both closed and open-ended questions. We selected villages according to their geographical location so that they would be representative of socio-demographic characteristics (e.g., religion and ethnicity), as well as for accessibility to QNP and its surroundings (see Appendix A). Upon arrival in a village, we first consulted with the community leaders to explain the purpose of the project and to obtain permission to visit households. We used two questionnaires (Appendix B), one for community leaders (n = 14) and another for households (n = 210), enabling us to capture the respective leaders’ specialized knowledge about village dynamics, characteristics and problems, as well as the perceptions of householders. We randomly selected 15 households and interviewed the heads of the households (typically men) to obtain information on livelihoods and vulnerability. To ensure a balanced sex ratio (Inskip et al. Reference Inskip, Ridout, Fahad, Tully, Barlow and Barlow2013), we also interviewed the partners of the heads of the households, resulting in a sample set of 118 men (56%) and 94 women (44%). Participants were informed about the purpose of the study, its anonymity and that participation was voluntary, before requesting signed or fingerprinted consent. We limited identifying information to village name and questionnaire number. The study was approved by the Ethics Committee for the Collection and Protection of Scientific Data (‘Comissão de Ética para Recolha e Protecção de Dados de Ciências’ – CERPDC) of the University of Lisbon, Portugal. Each interview lasted 35 minutes on average (range = 15–45 minutes; SD = 0.006) and had a response rate of 100%. Interviews were conducted in Portuguese by the lead author or in Makua by local native speakers hired for the project, or they were translated into the local dialect by a person from the community hired to work on the survey. All interviewers were trained in the sampling design, survey technique, confidentiality and participation consent protocol.

Questionnaire structure

For community leaders, we collected information on population size and village infrastructure (e.g., schools, healthcare centres, hospitals, markets, religious places, water fountains, electricity and telephone networks). For households, we obtained information to calculate the LVI (see the ‘Livelihood Vulnerability Index’ section; Fig. 2 & Appendix C). Our household questionnaire was based on that of Hanh et al. (Reference Hahn, Riederer and Foster2009), which was adapted for our context and field conditions. The questionnaire included seven sections: household demographics, livelihood strategies, social networks, health and health services, food security, access to water and community problems. For instance, we removed questions from the Hanh et al. (Reference Hahn, Riederer and Foster2009) questionnaire that referred to sensitive topics (e.g., percentage of households with orphans) or that we considered irrelevant to our study (e.g., average malaria exposure). We included questions about HWIs within the sections on health, food security and access to water by enquiring about household risk perception of zoonotic diseases (Decker et al. Reference Decker, Evensen, Siemer, Leong, Riley and Wild2010), crop damage by wildlife, attacks on domestic animals by wildlife and wildlife encounters. We also collected socio-demographic data from individuals with respect to gender, ethnicity, religion, household size and number of family members with any level of education. To better describe the local context, we calculated the village development index (Sahn & Stifel Reference Sahn and Stifel2003) and the intensity of healthcare requirements (Chambers & Conway Reference Chambers and Conway1992).

Fig. 2. The Livelihood Vulnerability Index (LVI) components and specific indicators, organized according to contributory factors to the Intergovernmental Panel on Climate Change’s (IPCC) Vulnerability Assessment Framework (i.e., exposure, sensitivity and adaptive capacity). The framework considers human–wildlife interactions as a possible component of exposure in the LVI, as examined in our analysis. The figure has been generated according to the IPCC’s Fourth Assessment Report (AR4) on climate change vulnerability.

Additional data

Village characteristics

We obtained the most recent geographical information for QNP (2012) from the National Administration for Conservation Areas (ANAC) that manages QNP to measure village distance to roads, distance to the nearest strict protection area and location (within or outside the park or in the buffer zone). Distance to the nearest strict protection area was measured as the Euclidean distance from the centroid of the village to the centroid of the nearest protected area (Madsen & Broekhuis Reference Madsen and Broekhuis2018). We used the centroid instead of the edge because animal population densities are likely be higher at the core of protected areas (Kiffner et al. Reference Kiffner, Stoner and Caro2013). We also calculated a village accessibility index as a combination of road class (primary (connecting provincial capitals), secondary (connecting primary roads and economic centres) and tertiary (connecting secondary roads and residential areas); INE 2017), road type (asphalt or unpaved), road condition (good, reasonable) and distance to the nearest primary road (<5 km, 5–10 km, >10 km). For road class, we also considered a subcategory of road types (national (N; major intercity roads), regional (R; connecting towns/localities)). All distance metrics were calculated in QGIS 3.6.1 (QGIS Development Team 2018).

Climate data

Climate data at 1-km2 resolution were obtained from the global database ‘WorldClim version 2’ (http://worldclim.org/). We used variables that reflected monthly variation in temperature and rainfall during 1970–2000, namely: temperature seasonality (BIO4), mean temperature of warmest quarter (BIO10), mean temperature of coldest quarter (BIO11), precipitation seasonality (BIO15), precipitation of warmest quarter (BIO16) and precipitation of coldest quarter (BIO17). We calculated the average of each bioclimatic variable for the area of the village and used it as an indicator of ‘climate variability’.

Data analysis

Livelihood Vulnerability Index

We used the LVI developed by Hahn et al. (Reference Hahn, Riederer and Foster2009) to measure livelihood vulnerability for the selected communities. The LVI is an additive indicator combining seven components deemed to influence livelihoods, namely: ‘socio-demographic profile’, ‘livelihood strategies’, ‘social networks’, ‘health’, ‘food’, ‘water’ and ‘climate variability’ (Fig. 2 & Appendix C). The LVI builds on the Intergovernmental Panel on Climate Change (IPCC) vulnerability assessment framework, which considers exposure, sensitivity and adaptive capacity as factors contributing to vulnerability (IPCC 2001). We calculated the LVI and its components as standardized scores following the approach of Hahn et al. (Reference Hahn, Riederer and Foster2009), which considers exposure as ‘climate variability’ (e.g., bioclimatic variables); sensitivity as a combination of access to ‘health’ (e.g., distance to health facility, percentage of family members with chronic disease), ‘food’ (e.g., crop diversity, number of months without food) and ‘water’ (e.g., distance to water source, percentage of households without daily water availability); and adaptive capacity as a combination of ‘livelihood strategies’ (e.g., livelihood diversification, percentage of households solely dependent on agriculture), ‘social network’ (e.g., percentage of households lending and/or borrowing money, percentage of households asking for help from community leaders) and ‘socio-demographic profile’ (e.g., percentage of female heads of households, dependency ratio) (Hanh et al. Reference Hahn, Riederer and Foster2009). We further explored the possibility of adding HWIs as a factor of exposure (see the ‘Relationship between the LVI and HWIs’ section).

Human–wildlife interactions

Four parameters relating to HWIs (percentage of households whose agricultural fields were damaged by wildlife, percentage of households who became ill due to wildlife-transmitted diseases, percentage of households whose livestock production had been affected by wildlife, percentage of households who encountered wildlife at water sources) were used to calculate standardized HWI scores and a respective average per household.

Drivers of the LVI and HWIs

We used a non-parametric Kruskal–Wallis analysis of variance (Zar Reference Zar2010) to test whether the variance in the LVIs differed between villages and districts. We performed a post hoc comparison of the pair-wise means using Tukey’s honest significant difference test (Tukey Reference Tukey1949).

First, we used a principal component analysis (PCA; Jolliffe Reference Jolliffe2002) to examine whether the LVI components and HWIs are related in a multidimensional space. Then, we used a set of linear mixed-effects models (LMMs; Zuur et al. Reference Zuur, Ieno, Walker, Saveliev and Smith2009) to quantify the effects of potential drivers on the LVI and its components. We explored three sets of variables (socio-demographics, village characteristics and HWIs) as potential drivers of the LVI. Similarly, we explored the same sets of variables as drivers of the LVI components, except that we added additional LVI components to HWIs. In this ‘HWI + other LVI components model’, we tested the effects of the LVI components on each other. For example, if ‘food’ was the response variable, then the components ‘water’, ‘health’, ‘socio-demographic profile’, ‘social network’, ‘livelihood strategies’, ‘climate’ and ‘HWI’ were added as predictors. Lastly, for the model with ‘HWI’ as the response, we considered four sets of variables (i.e., socio-demographics, village characteristics, LVI components and ‘LVI’). The variables ‘village’ and ‘district’ were included in all of the models as nested random effects (i.e., ‘village’ nested within ‘district’).

We used a multiple-stage modelling approach, whereby initially we built independent sets of models and then built a combined model encompassing the variables included in the best-performing models from the previous stage for which coefficients had been reliably estimated (Morin et al. Reference Morin, Yackulic, Diffendorfer, Lesmeister, Nielsen and Reid2020). We tested for collinearity between variables using a variance inflation factor (Zuur et al. Reference Zuur, Ieno and Smith2007), prompting us to remove one variable (‘local zone’) from the analysis. We selected the best-performing models from each set of models using Akaike’s information criterion corrected for small sample sizes (AICc; Akaike Reference Akaike1974). We chose ΔAICc < 5 to identify the best models, and final parameter and error estimates were calculated by model averaging of the best model(s) (Burnham & Anderson Reference Burnham, Anderson, Burnham and Anderson2002). We only report variables for which it was possible to reliably estimate an effect (i.e., the 95% confidence interval (95% CI) around the respective coefficient (β) did not encompass zero; see Appendix D for additional results). For the variables included in the best models, we estimated relative importance (RI) as the cumulative model weights of all models that included those variables (Arnold Reference Arnold2010).

Relationship between the LVI and HWIs

For the LVI components having an effect on the LVI (i.e., those contributing to the best models), we examined whether their indicators were related to HWIs. We developed an additional LMM in which ‘HWI’ acted as the response variable and individual indicators were the predictors, and we followed the same procedure as detailed above.

We also examined whether HWIs exerted an indirect (as an external driver) or direct (as a component within the LVI) effect on the LVI. We assessed an indirect HWI effect by generating a LMM to test the impact of HWIs on exposure, sensitivity and adaptive capacity individually. Next, we tested whether HWIs are a direct component of the LVI by adding HWIs as an indicator of exposure. We built three LMMs to test the effect of this expanded parameter of exposure on sensitivity. We used three variable sets for exposure – (1) HWI, (2) Climate and (3) HWI + Climate – and used the same modelling approach as detailed above.

Statistical analysis was performed in R software (R Core Team 2019) using the packages Hmisc (Harrell & Frank Reference Harrel and Frank2015), lme4 (Bates et al. Reference Bates, Maechler, Bolker and Walker2015), MuMIn (Bartón 2019), FactoMiner (Le et al. Reference Le, Josse and Husson2008), Factoextra (Kassambara & Mundt Reference Kassambara and Mundt2016) and missMDA (Husson & Josse Reference Josse and Husson2016).

Results

Almost all (96%) of the interviewed households reported that their crops had been damaged by wildlife, a majority of which (59%) indicated that it was a daily occurrence. The interviewees identified six crop-raiding species, including baboons (Papio cynocephalus; 47%), bush pigs (Potamochoerus porcus) and/or warthogs (Phacochoerus aethiopicus; 32%), vervet monkey (Chlorocebus pygerythrus; 13%), elephants (Loxodonta africana; 7%) and a single case of hippopotamus (Hippopotamus amphibius; 1%). Apart from crop-raiding events, 35% of the households mentioned encountering wildlife when collecting water from nearby rivers. Only 3% of respondents associated zoonotic diseases with family illness, and another 4% reported attacks by wildlife on domestic animals.

Livelihood vulnerability across QNP

Livelihood vulnerability was similar among districts (LVI range = 0.33–0.39; Table C.1 & Appendix E) and villages (LVI range = 0.29–0.43), but with high intra-village variability (LVI range = 0.21–0.53; Fig. 3a, Table C.2 & Appendix E). We observed considerable variability in LVI components across villages (Fig. 3b). In particular, the parameters ‘socio-demographics’ (range: 0.008–0.8), ‘livelihood strategies’ (range: 0.04–1.0) and ‘water’ (range: 0–0.8) presented the widest ranges.

Table 1. Averaged model coefficients for variables included in the best combined models (ΔAICc < 5). The variables with a 95% confidence interval not including zero are highlighted in bold.

AICc = Akaike’s information criterion corrected for small sample sizes; CI = confidence interval; HWI = human–wildlife interaction; LVI = Livelihood Vulnerability Index; RI = relative importance.

Fig. 3. (a) The Livelihood Vulnerability Index (LVI) for sampled villages within and close to Quirimbas National Park (QNP). Bubble size illustrates index values, ranging from 0.29 (less vulnerable) to 0.43 (more vulnerable). (b) The radar graph displays the scores for the LVI components in different colours. Cumulative scores per village are represented by individual bars. The villages have been grouped by district. Significantly different mean LVI scores after the Tukey’s honest significant difference test for each district are represented by asterisks and are grouped by letter case.

We identified a link between ‘HWI’, ‘water’ and ‘social network’, as these three variables were all positioned on the positive side of PCA axis 1 (Fig. 4). In contrast, ‘livelihood strategies’, ‘health’ and ‘food’ appeared on the negative side of the PCA axis 1. PCA axis 2 further separated ‘livelihood strategies’ from ‘social network’. Together, the two PCA axes accounted for 36.7% of the total variance in the data (PCA1 = 20.4% and PCA2 = 16.3%). The variables most contributing to the two PCA axes were ‘water’ (21.1%), ‘HWI’ (20.7%), ‘livelihood strategies’ (18.6%) and ‘food’ (16.7%). Village segregation across the two PCA axes implies a contribution of different livelihood vulnerability components to data variability (Figs 3b & 4).

Fig. 4. Principal component analysis plot in which primary contributing components to the Livelihood Vulnerability Index are displayed as arrows, questionnaires per district are represented by different symbols and ellipses represent the 95% confidence interval. The blue dots are the centroids of the villages.

Drivers of livelihood vulnerability across QNP

Our LMM results showed that livelihood vulnerability was best predicted by ‘HWI’, ‘gender’, ‘accessibility’ and ‘population size’. More vulnerable households tended to be more exposed to ‘HWI’ and men were less vulnerable than women. Villages with better access to main roads were also more vulnerable, whereas villages with larger population sizes were less vulnerable. Of all the variables we considered, ‘HWI’ and ‘gender’ had the strongest positive effect on the LVI (Table 1, LVI model; Appendix F Fig. F.1).

Drivers of components of livelihood vulnerability

Exposure

None of our variables constituted good predictors of ‘climate variability’. The null model was the only model selected from a total of 190 models.

Sensitivity

‘Food’ was negatively associated with ‘age’ of the household head, implying that older people experienced lower food scarcity. ‘Food’ was also positively associated with ‘livelihood strategies’, indicating that households with greater food scarcity were less resilient (Appendix F Figs F.1a & F.1b & Table F.1). We also found that access to water was positively associated with the interaction between ‘household size’ and ‘HWI’. In other words, larger households with limited access to water had greater exposure to HWIs (Appendix F Fig. F.1c & Table F.2). None of the variables were good predictors of ‘health’ sensitivity (only the null model was selected).

Adaptive capacity

The ‘socio-demographic profile’ of households was negatively associated with ‘gender’ of the head of household and positively associated with village ‘accessibility’ (Appendix F Figs F.1d & F.1e & Table F.3), with women prone to having a more vulnerable socio-demographic profile. ‘Food’ was also positively associated with ‘livelihood strategies’, as food scarcity enhanced the likelihood of vulnerable ‘livelihood strategies’ (Appendix F Fig. F.1f & Table F.4). Finally, the ‘social network’ of households was negatively affected by ‘socio-demographic profile’ and ‘livelihood strategies’ (Appendix E Figs F.1f & F.1h & Table F.5).

Drivers of HWIs

More vulnerable households, in particular those with more restricted access to water, were also more exposed to HWIs. We found that ‘HWI’ was best explained by ‘LVI’ and ‘water’, with both of those variables exerting a positive effect on ‘HWI’ (Table 1, HWI model; Appendix G Fig. G.1 & Table G.1).

The type of ‘water source’ (river/lake/lagoon or village well) was the only variable with an effect on the average best model (Appendix G Fig. G.2 & Table G.2), revealing that HWI risk was greater when households fetched water from a river, lake or lagoon than when they used wells in the village.

HWIs as a component of livelihood vulnerability

We identified a positive association between ‘HWI’ and livelihood ‘sensitivity’, with households more exposed to ‘HWI’ being the ones displaying higher sensitivity. In fact, the variable ‘sensitivity’ was present in the best four models, out of a total of eight models, and it was reliably estimated (Table 1, see HWI vulnerability categories). ‘Exposure’ and ‘adaptive capacity’ were present in just two of the eight best models, but the respective coefficients were not reliably estimated (Appendix H Table H.1). Furthermore, the models of ‘sensitivity’ as a function of ‘HWI’ and of ‘sensitivity’ as a function of ‘climate’ were both weaker than models of sensitivity as a function of ‘climate’ plus ’HWI’, with this latter reflecting exposure (Fig. 5).

Fig. 5. Sensitivity as a function of Exposure to ‘Climate’, ‘Climate + HWI’ and ‘HWI’. Model performance metrics include Akaike’s information criterion corrected for small sample sizes (AICc) and ΔAICc. The null model had an AICc = –385.43. HWI = human–wildlife interaction.

Discussion

We show the extent to which HWIs and livelihood vulnerability are linked in QNP, Mozambique. We found that high sensitivity and low adaptive capacity render individuals and communities more vulnerable to risks from wildlife. Furthermore, we show that HWIs, such as crop raiding or wildlife encounters when fetching water from rivers, were prevalent and almost daily occurrences in the study area. Although in many cases these encounters were not serious, some involved potentially dangerous species such as elephants (33%), crocodiles (Crocodylus niloticus; 13%), carnivores (lion (Panthera leo), leopard (Panthera pardus), hyenas (Crocuta crocuta); 4%) and hippopotamus (only one report). Moreover, when we included HWIs as a factor of exposure, we found that it was strongly related to sensitivity and overall livelihood vulnerability.

On average, communities living in QNP (LVI = 0.37) seem to have higher livelihood vulnerability than those inhabiting some other districts of Mozambique (e.g., Moma, LVI = 0.316; Mabote, LVI = 0.326; Hahn et al. Reference Hahn, Riederer and Foster2009), but they are likely less vulnerable than communities in Limpopo National Park (also in Mozambique). In Limpopo, communities were qualitatively classified as exceedingly vulnerable to climate change and very exposed to wildlife damage (Notelid & Ekblom Reference Notelid and Ekblom2021). These results suggest that a larger number of factors can increase vulnerability, and in some cases the contribution of protected areas and HWIs to the overall livelihood vulnerability might be minor.

Fundamentally, the livelihood vulnerability of communities in QNP is enhanced by water scarcity, as communities take risks while fetching water. Rivers are attractive areas for wildlife and they are used as dispersal corridors. Accordingly, it is not surprising that we uncovered a higher probability of HWIs when communities used these water sources (Madsen & Broekhuis Reference Madsen and Broekhuis2018). To counteract this problem, we advocate increasing the number of communal closed wells near villages, thereby reducing the need to retrieve water from nearby rivers and, consequently, decreasing the risk of potentially negative HWIs. In addition, the association between the variables ‘HWI’, ‘water’ and ‘social network’ (Fig. 4) indicates that households with less social support struggled to manage their water resources, so they displayed an increased need to retrieve water from rivers, exposing them to HWIs. In other words, households with a strong social network were more likely to fetch water in groups, reducing their exposure to wildlife encounters and boosting security. Indeed, a critical role for social ties in reducing vulnerability to water scarcity has been demonstrated previously, such as by enabling water resources to be shared among households in ways that are beneficial to the social network and the wider community (Faurès & Santini Reference Faurès and Santini2009) through increased communication (e.g., on weather) and partnerships (Dickson et al. Reference Dickson, Schuster-Wallace and Newton2016). Our modelling results have revealed that the relationship between water scarcity and HWIs is stronger for larger households, most likely because such households need more water (Angoua et al. Reference Angoua, Dongo, Templeton, Zinsstag and Bonfoh2018). Moreover, since access to wells in QNP villages is limited (i.e., a limited number of community wells and a large number of users), a greater number of households must fetch water more frequently, engendering a greater risk of HWIs.

Surprisingly, and despite 96% of householders reporting crop-raiding events, we did not find evidence for a relationship between food security and HWIs. This outcome contrasts with the findings of other similar studies (e.g., Barua et al. Reference Barua, Bhagwat and Jadha2013). We believe that this result could be due to the indicator of crop damage that we used (percentage of households whose agricultural fields were damaged by wildlife) not being sufficiently sensitive to the intensity or extent of such events. Since every respondent household suffered crop damage at least once a month, with a consequent reduction in food supply, all households displayed maximum vulnerability for this indicator.

We found that women had higher livelihood vulnerability than men. More specifically, young women who acted as the head of a household were the most vulnerable. Young women in rural communities are less empowered, are less likely to have undertaken formal education and are more restricted in their activities due to rigidly defined gender roles (Khumalo & Yung Reference Khumalo and Yung2015). Thus, most women are unlikely to have a direct source of income, hindering their potential to invest in alternative activities to diversify their livelihoods (Khumalo & Yung Reference Khumalo and Yung2015). Accordingly, they are less able to recover from or adopt strategies to overcome crop losses (Naughton-Treves Reference Naughton-Treves1997). Moreover, these vulnerable women are also more exposed to HWIs, especially as they are frequently responsible for fetching water and for the maintenance of agricultural fields very often visited by wildlife (Mwangi et al. Reference Mwangi, Akinyi, Maloba, Ngotho, Kagira, Ndeereh and Kivai2016), and they are less effective at guarding against and preventing crop damage (Naughton-Treves Reference Naughton-Treves1997).

Our analyses have also revealed that households in villages with more access to nearby cities suffered greater livelihood vulnerability. This outcome is contrary to the common expectation that remote villages are more highly vulnerable because they lack assets and business opportunities (Salerno Reference Salerno2016). However, villages closer to primary roads may suffer more rapid resource depletion, with negative consequences for the livelihoods of their communities (Mwangi et al. Reference Mwangi, Akinyi, Maloba, Ngotho, Kagira, Ndeereh and Kivai2016). We found that villages with greater road access presented more vulnerable socio-demographic profiles and lower development indices, reinforcing the evidence for a negative influence of road proximity (Appendix I.1). Furthermore, the presence of a road system does not necessarily mean greater mobility between villages (Asafo-Adjei & Iyer-Raniga Reference Asafo-Adjei and Iyer-Raniga2017). For the communities in QNP, mobility remains very limited, since almost no villager owns a car or motorcycle and public transportation is non-existent. Our data show that villages with larger population sizes tended to display a lower LVI than smaller villages, with this finding potentially being related to the greater capacity for development and infrastructure of the former (International Monetary Fund 2008) and their greater access to health facilities and water sources (Appendix I.2).

Although the LVI incorporates different components of livelihood vulnerability, it still has some limitations. The LVI oversimplifies the complex reality of livelihoods into indices and indicators that are difficult to validate (Hanh et al. Reference Hahn, Riederer and Foster2009). Furthermore, in the case of HWIs, our results suggest caution regarding the type of indicator of HWIs used in LVI assessments, such as crop damage intensity. We propose that spatial modelling of the probability of HWIs and their effects on the LVI could improve our understanding of the HWI–LVI bidirectional relationship. However, currently available data on wildlife ranges and habitat use (MITADER 2012) are insufficient for proper analyses of species distributions. Other relevant indicators for the study area may also be lacking (e.g., human migration) because the underlying socio-economic, cultural and ecological contexts are still poorly understood.

Conclusions

Understanding the relationship between livelihood vulnerability and HWIs is fundamental to attaining improved management and conservation of wildlife populations and the sustainable development of rural populations throughout Africa. We found evidence of a strong relationship between livelihood vulnerability and HWIs, implying that this dependency needs to be better addressed if both of these goals are to be met. Based on our results, simple measures such as enhancing cooperation in water and food activities within communities and adding closed wells near villages could quickly diminish exposure to HWIs. In the long term, households could be encouraged to invest in crops that are more resistant to wildlife damage and to switch from agriculture-based activities, although this latter would have land-use impacts elsewhere. We also recommend that community livelihoods and sensitivities must be considered when designing HWI mitigation schemes and implementing conservation actions. Furthermore, our study highlights that the problem of HWIs cannot be addressed without addressing the root causes of livelihood vulnerability. Accordingly, we propose that HWIs should be included as a component of exposure by extending the LVI framework proposed by Hanh et al. (Reference Hahn, Riederer and Foster2009). Finally, in order to address the multiple dimensions of livelihood vulnerability, specifically where HWIs have a profound impact on poor rural communities and given that HWIs affect millions of people globally, communities should be provided with the capacity to respond to the risks of HWIs.

Supplementary material

To view supplementary material for this article, please visit https://doi.org/10.1017/S037689292100028X.

Acknowledgements

The authors thank all participants in the study from the villages of Bilbiza, Chimpembe, Linde, Meluco, Miagane, Mirate, Muaja, Muaguide, Nairoto, Namange, Quissanga, Ravia, Sitate and Soco – the interviewees, community chiefs, guides and translators. The study would not have been possible without the logistical support of Lúrio University in Pemba, especially from Marcelino Caravela, Cristovão Nanvonamuquitxo and Isildo Nganhane, nor without the authorization from the National Administration for Conservation Areas (ANAC) to work in QNP.

Financial support

This work was supported by a National Geographic Society early career grant (EC-51226C-19) and a Fundação para a Ciência e Tecnologia (FCT) PhD fellowship (SFRH/BD/137682/2018). LMR is grateful to FCT/MCTES for the financial support to CESAM (UID/AMB/50017/2019) and cE3c (UIDB/00329/2020), through national funds, and for the co-funding by the FEDER, within the PT2020 Partnership Agreement and Compete 2020. MJS was supported by the University of Zurich University Research Priority Program in Global Change and Biodiversity. CF was funded by the Programme ‘Apoio à Contratação de Recursos Humanos Altamente Qualificados’ (NORTE-06-3559-FSE-000045), supported by Norte Portugal Regional Operational Programme (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, and also is grateful to FCT/MCTES for the financial support to CESAM (UIDP/50017/2020, UIDB/50017/2020), through national funds.

Conflict of interest

None.

Ethical standards

The authors assert that all procedures contributing to this work comply with applicable ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.

References

Akaike, HAI (1974) A new look at the statistical model identification. IEEE Transactions on Automatic Control 19: 716723.CrossRefGoogle Scholar
Anderson, JL, Pariela, F (2005) Strategies to Mitigate Human–Wildlife Conflicts Mozambique (Report). Rome, Italy: FAO.Google Scholar
Angoua, ELE, Dongo, K, Templeton, MR, Zinsstag, J, Bonfoh, B (2018) Barriers to access improved water and sanitation in poor peri-urban settlements of Abidjan, Côte d’Ivoire. PLoS ONE 13: e0202928.CrossRefGoogle ScholarPubMed
Arnold, TW (2010) Uninformative parameters and model selection using Akaike’s information criterion. Journal of Wildlife Management 74: 11751178.CrossRefGoogle Scholar
Asafo-Adjei, CU, Iyer-Raniga, GAM (2017) Mobility and access to transportation for the rural poor: examining the role of intermediate means of transport (IMTS) in Ghana. Presented at: 13th International Postgraduate Research Conference. University of Salford, Manchester, UK, 14–15 September.Google Scholar
Barua, M, Bhagwat, SA, Jadha, S (2013) The hidden dimensions of human–wildlife conflict: health impacts, opportunity and transaction costs. Biological Conservation 157: 309316.10.1016/j.biocon.2012.07.014CrossRefGoogle Scholar
Bateman, IJ, Glew, D (2010) Tigers, markets and palm oil: market potential for conservation. Oryx 44: 230234.CrossRefGoogle Scholar
Bates, D, Maechler, M, Bolker, B, Walker, S (2015) Fitting linear mixed-effects models using lme4 . Journal of Statistical Software 67: 14.CrossRefGoogle Scholar
Bevan, P (2000) Poverty in Ethiopia. Presented at: DFID Ethiopia Seminar. London, UK, 6 November.Google Scholar
Bhatia, S, Redpath, SM, Suryawanshi, K, Mishra, C (2019) Beyond conflict: exploring the spectrum of human–wildlife interactions and their underlying mechanisms. Oryx 54: 104108.Google Scholar
Blair, AG, Meredith, TC (2018) Community perception of the real impacts of human–wildlife conflict in Laikipia, Kenya: capturing the relative significance of high-frequency, low-severity events. Oryx 52: 497507.CrossRefGoogle Scholar
Bond, J (2014) A holistic approach to natural resource conflict: the case of Laikipia County, Kenya. Journal of Rural Studies 34: 117127.10.1016/j.jrurstud.2014.01.008CrossRefGoogle Scholar
Burnham, KP, Anderson, DR (2002) Formal inference from more than one model: multimodel inference (MMI). In: Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach, eds Burnham, KP, Anderson, DR (p. 149). New York, NY, USA: Springer-Verlag.Google Scholar
Butt, B, Shortridge, A, WinklerPrins, A (2009) Pastoral herd management, drought coping strategies, and cattle mobility in southern Kenya. Annals of the Association of American Geographers 99: 309334.CrossRefGoogle Scholar
Ceauşu, S, Graves, RA, Killion, AK, Svenning, JC, Carter, NH (2019) Governing trade-offs in ecosystem services and disservices to achieve human–wildlife coexistence. Conservation Biology 33: 543553.CrossRefGoogle ScholarPubMed
Chambers, R, Conway, G (1992) Sustainable Rural Livelihoods: Practical Concepts for the 21st Century. Brighton, UK: Institute of Development Studies.Google Scholar
Cox, DTC, Gaston, KJ (2018) Human–nature interactions and the consequences and drivers of provisioning wildlife. Philosophical Transactions of the Royal Society B: Biological Sciences 373: 20170092.10.1098/rstb.2017.0092CrossRefGoogle ScholarPubMed
Decker, DJ, Evensen, DTN, Siemer, WF, Leong, KM, Riley, SJ, Wild, MA et al. (2010) Understanding risk perceptions to enhance communication about human–wildlife interactions and the impacts of zoonotic disease. ILAR Journal 51: 255261.CrossRefGoogle ScholarPubMed
Dickson, SE, Schuster-Wallace, CJ, Newton, JJ (2016) Water security assessment indicators: the rural context. Water Resources Management 30: 15671604.CrossRefGoogle Scholar
Duffy, R, St John, FAV, Büscher, B, Brockington, D (2016) Toward a new understanding of the links between poverty and illegal wildlife hunting. Conservation Biology 30: 1422.CrossRefGoogle Scholar
Dunham, KM, Ghiurghi, A, Cumbi, R, Urbano, F (2010). Human–wildlife conflict in Mozambique: a national perspective, with emphasis on wildlife attacks on humans. Oryx 44: 185193.CrossRefGoogle Scholar
Faurès, J-M, Santini, G (2009). Water and the Rural Poor Interventions for Improving Livelihoods in Sub-Saharan Africa (Report). Rome, Italy: FAO Water and Land Division.Google Scholar
García-Frapolli, E, Ayala-Orozco, B, Oliva, M, Smith, RJ (2018) Different approaches towards the understanding of socio-environmental conflicts in protected areas. Sustainability 10: 117.CrossRefGoogle Scholar
Gupta, AC (2013) Human–wildlife conflict and rural livelihoods around Chobe elephants, safety nets and agrarian culture: understanding National Park, Botswana. Journal of Political Ecology 20: 14.CrossRefGoogle Scholar
Hahn, MB, Riederer, AM, Foster, SO (2009) The Livelihood Vulnerability Index: a pragmatic approach to assessing risks from climate variability and change – a case study in Mozambique. Global Environmental Change 19: 7488.CrossRefGoogle Scholar
Harrel, JR, Frank, E (2015) Hmisc: Harrel miscellaneous [www document]. URL https://cran.r-project.org/web/packages/Hmisc/index.html.Google Scholar
INE (2017) 4° Censo Geral da População e Habitação 2017. Instituto Nacional de Estatística de Moçambique [www document]. URL www.ine.co.mz Google Scholar
Inskip, C, Ridout, M, Fahad, Z, Tully, R, Barlow, A, Barlow, CG et al. (2013) Human–tiger conflict in context: risks to lives and livelihoods in the Bangladesh Sundarbans. Human Ecology 41: 169186.CrossRefGoogle Scholar
International Monetary Fund (2008) Lao People’s Democratic Republic: Poverty Reduction Strategy Paper International Monetary Fund (Country Report No. 03/310). Washington, DC, USA: International Monetary Fund.CrossRefGoogle Scholar
IPCC (2001) Climate Change 2001: Impacts, Adaptation, and Vulnerability. Contribution of Working Group II to the Third Assessment Report. Cambridge, UK: Cambridge University Press.Google Scholar
Jezeer, R, Verweij, PA, Boot, RGA, Junginger, M, Santos, MJ (2019) Influence of livelihood assets, experienced shocks and perceived risks on smallholder coffee farming practices in Peru. Journal of Environmental Management 242: 496506.CrossRefGoogle ScholarPubMed
Jolliffe, IT (2002) Principal Component Analysis, 2nd edition. New York, NY, USA: Springer.Google Scholar
Josse, J, Husson, F (2016) missMDA: a package for handling missing values in multivariate data analysis. Journal of Statistical Software 70: 10.18637/jss.v070.i01.CrossRefGoogle Scholar
Junquera, V, Grêt-Regamey, A (2020) Assessing livelihood vulnerability using a Bayesian network: a case study in northern Laos. Ecology and Society 25: 131.CrossRefGoogle Scholar
Kassambara, A, Mundt, F (2016) Factoextra: extract and visualize the results of multivariate data analyses [www document]. URL https://cran.r-project.org/package=factoextra Google Scholar
Khumalo, KE, Yung, LA (2015) Women, human–wildlife conflict, and CBNRM: hidden impacts and vulnerabilities in Kwandu Conservancy, Namibia. Conservation and Society 13: 232243.CrossRefGoogle Scholar
Kiffner, C, Stoner, C, Caro, T (2013) Edge effects and large mammal distributions in a national park. Animal Conservation 16: 97107.CrossRefGoogle Scholar
Lamichhane, BR, Persoon, GA, Leirs, H, Poudel, S, Subedi, N, Pokheral, CP et al. (2018) Spatio-temporal patterns of attacks on human and economic losses from wildlife in Chitwan National Park, Nepal. PLoS ONE 13: e0195373.CrossRefGoogle ScholarPubMed
Le, BS, Chavernac, D, Mapuvire, G, Cornu, G (2014) FrontlineSMS as an early warning network for human–wildlife mitigation: lessons learned from tests conducted in Mozambique and Zimbabwe. Electronic Journal of Information Systems in Developing Countries 60: 113.Google Scholar
Le, S, Josse, J, Husson, F (2008) FactoMineR: an R package for multivariate analysis. Journal of Statistical Software 25: 118.CrossRefGoogle Scholar
Madsen, EK, Broekhuis, F (2018) Determining multi-species site use outside the protected areas of the Maasai Mara, Kenya, using false positive site-occupancy modelling. Oryx 3: 395404.Google Scholar
Majale, M (2002) Towards pro-poor regulatory guidelines for urban upgrading. Presented at: Workshop on Regulatory Guidelines for Urban Upgrading. Bourton-on-Dunsmore, UK, 17–18 May.Google Scholar
Merz, L, Pienaar, EF, Filk, T, Muyengwa, S (2021) Tolerance of wildlife in protected areas borderlands [www document]. URL https://www.biorxiv.org/10.1101/2021.03.27.436188 Google Scholar
MITADER (2012) Plano de Maneio 2011–2021 Parque Nacional das Quirimbas. Pemba, Mozambique: Ministry of Tourism.Google Scholar
Mondal, MSH (2019) The implications of population growth and climate change on sustainable development in Bangladesh. Jàmbá Journal of Disaster Risk Studies 11: 535.CrossRefGoogle ScholarPubMed
Morales-Reyes, Z, Pe, JM, Margalida, A, Dona, A, Lazcano, C (2015) Supplanting ecosystem services provided by scavengers raises greenhouse gas emissions. Scientific Reports 5: 7811.CrossRefGoogle ScholarPubMed
Morin, DJ, Yackulic, CB, Diffendorfer, JE, Lesmeister, DB, Nielsen, CK, Reid, J (2020) Is your ad hoc model selection strategy affecting your multimodel inference? Ecosphere 11: 116.CrossRefGoogle Scholar
Mwangi, DK, Akinyi, M, Maloba, F, Ngotho, M, Kagira, J, Ndeereh, D, Kivai, S (2016) Socioeconomic and health implications of human–wildlife interactions in Nthongoni, eastern Kenya. African Journal of Wildlife Research 46: 87102.CrossRefGoogle Scholar
Naidoo, R, Fisher, B, Manica, A, Balmford, A (2016) Estimating economic losses to tourism in Africa from the illegal killing of elephants. Nature Communications 7: 13379.CrossRefGoogle ScholarPubMed
Naughton-Treves, L (1997) Farming the forest edge: vulnerable places and people around Kibale National Park, Uganda. Geographical Review 87: 27.CrossRefGoogle Scholar
Notelid, M, Ekblom, A (2021) Household vulnerability and transformability in Limpopo National Park. Sustainability 13: 2597.CrossRefGoogle Scholar
Ogra, MV (2008) Human–wildlife conflict and gender in protected area borderlands: a case study of costs, perceptions, and vulnerabilities from Uttarakhand (Uttaranchal), India. Geoforum 39: 14081422.CrossRefGoogle Scholar
QGIS Development Team (2018) QGIS Geographic Information System. Open Source Geospatial Foundation Project [www document]. URL http://qgis.org Google Scholar
Ratnayeke, S, Van Manen, FT, Pieris, R, Pragash, VSJ (2014) Challenges of large carnivore conservation: sloth bear attacks in Sri Lanka. Human Ecology 42: 467479.Google Scholar
R Core Team (2019) R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing.Google Scholar
Sahn, DE, Stifel, D (2003) Exploring alternative measures of welfare in the absence of expenditure data. Review of Income and Wealth 49: 463489.CrossRefGoogle Scholar
Salerno, J (2016) Migrant decision-making in a frontier landscape. Environmental Research Letters 11: 044019.CrossRefGoogle Scholar
Seoraj-Pillai, N, Pillay, N (2017) A meta-analysis of human–wildlife conflict: South African and global perspectives. Sustainability 9: 34.CrossRefGoogle Scholar
Stone, M, Phalke, S, Warren, N, Av, DK, Krishnan, A (2019) Gauging farmers’ acceptance of ‘social barrier’ mechanisms for preventing elephant crop raids. Gajah 50: 2328.Google Scholar
Treves, A (2009) Hunting for large carnivore conservation. Journal of Applied Ecology 46: 13501356.CrossRefGoogle Scholar
Tukey, J (1949) Comparing individual means in the analysis of variance. Biometrics 5: 99114.CrossRefGoogle ScholarPubMed
World Bank (2018) World Development Indicators, Rural population [www document]. URL https://data.worldbank.org/indicator/SP.RUR.TOTL.ZS?locations=ZG Google Scholar
Yadava, RN, Sinha, B (2020) Vulnerability assessment of forest fringe villages of Madhya Pradesh, India for planning adaptation strategies. Sustainability 12: 1253.CrossRefGoogle Scholar
Zar, JH (2010) Biostatistical Analysis. Upper Saddle River, NJ, USA: Pearson Prentice Hall.Google Scholar
Zuur, AF Ieno, EN Smith, GM (2007) Analysing Ecological Data. New York, NY, USA: Springer.CrossRefGoogle Scholar
Zuur, AF, Ieno, EN, Walker, NJ, Saveliev, AA, Smith, GM (2009) Mixed Effects Models and Extensions in Ecology with R. New York, NY, USA: Springer.CrossRefGoogle Scholar
Figure 0

Fig. 1. (a) Location of Quirimbas National Park (QNP) in the Republic of Mozambique on the African continent. (b) Sampled villages within and close to QNP.

Figure 1

Fig. 2. The Livelihood Vulnerability Index (LVI) components and specific indicators, organized according to contributory factors to the Intergovernmental Panel on Climate Change’s (IPCC) Vulnerability Assessment Framework (i.e., exposure, sensitivity and adaptive capacity). The framework considers human–wildlife interactions as a possible component of exposure in the LVI, as examined in our analysis. The figure has been generated according to the IPCC’s Fourth Assessment Report (AR4) on climate change vulnerability.

Figure 2

Table 1. Averaged model coefficients for variables included in the best combined models (ΔAICc < 5). The variables with a 95% confidence interval not including zero are highlighted in bold.

Figure 3

Fig. 3. (a) The Livelihood Vulnerability Index (LVI) for sampled villages within and close to Quirimbas National Park (QNP). Bubble size illustrates index values, ranging from 0.29 (less vulnerable) to 0.43 (more vulnerable). (b) The radar graph displays the scores for the LVI components in different colours. Cumulative scores per village are represented by individual bars. The villages have been grouped by district. Significantly different mean LVI scores after the Tukey’s honest significant difference test for each district are represented by asterisks and are grouped by letter case.

Figure 4

Fig. 4. Principal component analysis plot in which primary contributing components to the Livelihood Vulnerability Index are displayed as arrows, questionnaires per district are represented by different symbols and ellipses represent the 95% confidence interval. The blue dots are the centroids of the villages.

Figure 5

Fig. 5. Sensitivity as a function of Exposure to ‘Climate’, ‘Climate + HWI’ and ‘HWI’. Model performance metrics include Akaike’s information criterion corrected for small sample sizes (AICc) and ΔAICc. The null model had an AICc = –385.43. HWI = human–wildlife interaction.

Supplementary material: File

Pereira et al. supplementary material

Pereira et al. supplementary material

Download Pereira et al. supplementary material(File)
File 1.4 MB