Introduction {#Sec1}
============

Climate change and land-use change are among the major drivers of biodiversity loss^[@CR1],[@CR2]^. Impacts of climate change are predicted to be particularly important in already fragile areas with low ecosystem resilience^[@CR3]^. Anthropogenic land conversion and climate change may often act in synergy, triggering a negative feedback whereby habitats may rapidly deteriorate, resulting in wildlife population declines and possible extinctions^[@CR4]^. Such synergistic impacts may be particularly important in most regions of the world where large swaths of land are under some form of production regimes (such as farmland) and are at high vulnerability from climate change^[@CR3],[@CR5]^.

Among those regions under rapid anthropogenic conversion and of high vulnerability to climate change are the arid savannahs of Africa^[@CR5]^. This biome was recently found as highly endangered, as most of it is currently unprotected^[@CR6]^. The arid savannahs of Africa support important taxonomic and phylogenetic diversity^[@CR7]^, provide key ecosystem services, such as food and fiber, and generate profit from ecotourism that contribute to improve local livelihoods^[@CR8]^. These areas, once largely natural, are now being widely converted to farmland, either for livestock or for crop production^[@CR9]^. In these man-managed arid savannah landscapes, relatively small shifts towards harsh conditions triggered by climate change may strongly impact on the life-history of some species^[@CR10]--[@CR13]^.

Land-use and environmental change may have particularly disproportionate impacts on specialist species at the top of the food chain, such as obligate scavengers which rely on multiple trophic levels for their survival and reproduction^[@CR14]^. As the sole obligate scavengers, vultures play an important role in the savannah ecosystem, and across many other terrestrial systems of the world^[@CR15]^. By comprising an exclusive functional guild, vultures contribute to clean the environment from decomposing organic matter and keep the ecosystem in balance^[@CR15],[@CR16]^. Carrion, the primary food source for vultures, constitutes a largely unpredictable and localized resource that is also used by competitor facultative scavengers^[@CR14],[@CR15]^.

While anthropogenic activities are causing widespread vulture declines^[@CR15],[@CR16]^, in the absence of anthropogenic-driven mortality, vulture populations are generally limited by resource availability and their slow reproductive rates^[@CR14]^ which are typical of large-bodied long-lived species often representing K life history strategy^[@CR17]^. The ranging behavior and distribution of vultures are strongly related to the availability of carrion^[@CR18]--[@CR21]^. Ultimately, food availability for vultures will be determined by the abundance of wild ungulates (or domestic livestock) and their mortality rate, which may vary in space and time in response to environmental conditions and management practices^[@CR22]^. Carrion availability will also depend on the presence and abundance of mammalian competitors, as well as carcass visibility and accessibility^[@CR23],[@CR24]^.

Environmental conditions may directly affect the body condition of nestlings (e.g., through thermoregulation and water stress^[@CR11]^) or indirectly through their influence on foraging success and movement patterns of their parents^[@CR25]^. Nestling body condition may ultimately determine breeding success or survival during the post-fledging period and consequent recruitment to the breeding population^[@CR26]^. For example, environmental conditions experienced early in life were found to affect future survival prospects of long-lived vulture species^[@CR27]^. There is a large body of scientific work on the movement and distribution of vultures in relation to the coverage of protected areas as well as spatio-temporal variation in environmental conditions^[@CR19],[@CR22],[@CR28],[@CR29]^. However, studies on the effects of environmental and climate-driven drivers on the breeding biology of obligate scavengers are still relatively scarce.  More generally, studies addressing the interactive impacts of land-use and environmental change on species life-histories are still very rare, hindering the development and implementation of appropriate adaptation measures under global and local changes^[@CR30]^.

In this study we explore the independent and interactive impact of extrinsic factors on nestling body condition of two obligate scavenger species, the lappet-faced (*Torgos tracheliotus*) and African white-backed (*Gyps africanus*; hereafter white-backed vulture) vulture, in an arid ecosystem of Africa.  For the white-backed vulture, only independent impacts were explored due to collinearity (see methods).  We hypothesize that forage availability for ungulates (here quantified using a measure of vegetation greenness) could affect ungulate mortality (e.g. during prolonged periods of poor forage availability) and thus carcass availability for vultures, in turn affecting nestling condition.  We predict that harsh, dry conditions assumed to increase mortality rates of ungulates will benefit carrion eating vultures.  We used a unique dataset of ringing observations collected on the populations of these two species in Namibia (Fig.  [1](#Fig1){ref-type="fig"}).  The area experiences rather extreme environmental conditions, e.g. prolonged and intense droughts.  The large dataset, accumulated over 14 years and across protected and unprotected areas provides a unique spatio-temporal coverage for testing our hypothesis.  Specifically, we aim to quantify how nestling body condition is affected by (i) the cover of protected areas around a nest and (ii) environmental conditions (i.e. climate-driven forage availability for ungulates).  The spatio-temporal extent of the dataset also allows us to test (iii) the sensitivity of the results to different temporal scales of changes in environmental conditions (i.e. short or prolonged periods of below-average forage availability assumed to affect ungulate mortality, and in turn nestling body condition).  Protected and unprotected areas represent contrasting environments which may be associated with different ungulate and livestock mortality rates under stressful environmental conditions.  For example, dry periods may reduce water and forage availability for ungulates in protected areas, resulting in higher mortality there as compared to unprotected farmland areas where constant water and forage are provided to livestock.  Therefore, we also tested (iv) for interactive effects of climate-driven environmental conditions and protected areas on nestling body condition. Figure 1Map of Namibia with the location of the three study regions (delimited with thick black line).  Left panel shows also the cover of protected areas (grey shades) and the locations of nest where nestling measurements for the lappet-faced (empty circles) and white-backed (empty triangles) vultures included in this study.  The right panel shows the normalized difference vegetation index (NDVI) across the study landscape, in this case representing the situation as an yearly average NDVI value for the year 2007.  Figure created in ArcGIS 10.1.0 software (<http://desktop.arcgis.com/en/>).  The polygons of administrative borders and protected areas of Namibia were retrieved from http://www.nnf.org.na/eis/30.html a free online information resource.  The NDVI dataset used was acquired from NASA LP DAAC.  MODIS/Terra + Aqua Vegetation Indices 16-Day L3 Global 250 m (MOD13Q1).  USGS Earth Resources Observation and Science (EROS) Center, Sioux Falls, South Dakota (<https://lpdaac.usgs.gov>), accessed January 31, 2017. 

Results {#Sec2}
=======

Body condition of lappet-faced vulture nestlings was negatively related to NDVI (Table [1](#Tab1){ref-type="table"} and Fig.  [2](#Fig2){ref-type="fig"}).  This implies that body condition of lappet-faced vulture nestlings was higher when environmental conditions prior to nestling measurement were drier (lower vegetation greenness, or productivity).  This pattern was not sensitive to the period considered for calculating the NDVI, as it was similar for all five temporal scales (from 1 to 36 months; see Table [1](#Tab1){ref-type="table"} and Fig.  [2](#Fig2){ref-type="fig"}).  Moreover, we found an interactive effect of the protected area coverage and NDVI on body condition (Table [1](#Tab1){ref-type="table"}).  This effect was strongest, i.e. the upper and lower credible intervals of the posterior mean did not overlap zero, for the NDVI referring to the shortest periods (one and three months) preceding the body condition measurements (Fig.  [3](#Fig3){ref-type="fig"}, Table [1](#Tab1){ref-type="table"} and Fig.  [S2](#MOESM1){ref-type="media"}).  Specifically, the effect of protected area on nestling body condition was most positive after dry periods, i.e. low NDVI during preceding months, but became negative after wet periods (Fig.  [3](#Fig3){ref-type="fig"}).  Body condition of white-backed vulture nestlings decreased with the cover of protected areas around the nest, but was not affected by NDVI at any of the five temporal scales considered (Table [1](#Tab1){ref-type="table"}, Figs [2](#Fig2){ref-type="fig"} and [S3](#MOESM1){ref-type="media"}).  Moreover, white-backed vulture nestling body condition declined markedly over the 14 years of study (Table [1](#Tab1){ref-type="table"} and Fig.  [S3](#MOESM1){ref-type="media"}).  A similar but weaker decline in condition seems also apparent for lappet-faced vulture nestlings. Table 1Summary of the posterior marginal distributions for Bayesian hierarchical models exploring the relationship between the body condition of lappet-faced and white-backed vulture nestlings in relation to the proportion of area protected (PA cover), and the normalized difference vegetation index (NDVI).VariableLappet-faced vultureWhite-backed vultureMeanLower CIUpper CIMeanLower CIUpper CIPA cover0.07−0.010.15**−0.39−0.70−0.07**NDVI 1 m**−0.12−0.20−0.03**−0.08−0.220.07NDVI 3 m**−0.10−0.18−0.01**0.13−0.030.29NDVI 12 m**−0.14−0.23−0.05**−0.01−0.140.13NDVI 24 m**−0.14−0.23−0.05**0.06−0.120.23NDVI 36 m**−0.15−0.25−0.06**0.14−0.070.34Time−0.07−0.140.00**−0.16−0.26−0.07**NDVI 1 m \* PA cover**−0.06−0.11−0.01**NDVI 3 m \* PA cover**−0.06−0.10−0.01**NDVI 12 m \* PA cover−0.04−0.090.01NDVI 24 m \* PA cover−0.04−0.090.01NDVI 36 m \* PA cover−0.04−0.090.01We calculated NDVI at five temporal scales representing the situation during the one month prior to the measurement (NDVI 1 m) and up until 3 years prior to the measurement (NDVI 36 m) and modeled these separately.  PA cover and the NDVI refer to the area within a 50 km radius around the nest (see methods).  The variable named time aimed to depict continuous inter-annual trends in body condition.  Each variable and related statistics are derived from separate models including only that variable and time as continuous in the fixed part of a model (see methods and Appendix S1).  The statistics for the main variables effects were derived from models without interactions.  Values represent the mean of the posterior distribution with upper and lower 95% credibility intervals.  Variables for which the credibility intervals do not overlap zero are shown in bold font. Figure 2The relationship between the body condition of lappet-faced (left panels) and white-backed (right panels) vulture nestlings and the normalized difference vegetation index (NDVI).  We explored NDVI effects across five temporal scales representing the situation during the one month prior to nestling measurement (uppermost row) and up until 3 years prior to the measurement (lowermost row).  Each panel shows the raw data (circles) and the fitted linear relationship obtained from the Bayesian hierarchical models (straight black line, with 95% credibility intervals shown as grey shade; see methods and Table [1](#Tab1){ref-type="table"}).Figure 3Predicted body condition (here in the form of the scaled body mass index) of lappet-faced vulture nestlings in relation to the normalized difference vegetation index (NDVI) and the cover of protected areas around the nest.  The different panels show results obtained by using NDVI calculated at five temporal scales - representing the situation during the one month prior to the measurement and up until 3 years prior to the measurement.  These predictions are derived from the model including the interaction term between protected area cover and each of the five NDVI measures (see methods for further details).  Warm colors indicate better condition compared to cold colors. 

Discussion {#Sec3}
==========

We show that the effects of land-use, i.e. coverage of protected areas, and climate-driven conditions, i.e. NDVI, on vulture nestling body condition are species-specific and can be interactive.  Partly in accordance with our initial hypothesis, we found that, for lappet-faced vulture nestlings only, body condition is increased after periods of low availability of forage for ungulates (below average NDVI, presumably linked to increased ungulate mortality), and that this effect is mediated by land-use, i.e. relative coverage of protected areas.  Specifically, at nests largely surrounded by protected areas, body condition increased after periods of low availability of ungulate forage, but body condition worsened following periods of high forage availability.  White-backed vulture nestlings, however, were not influenced by variation in forage availability, their body condition was lower at largely protected nests and also declined through the study period. 

Body condition of altricial nestlings is typically mediated by the amount and quality of food that the parents can provide^[@CR17]^. The finding that white-backed vulture nestlings are in better condition in unprotected as compared to protected areas is interesting.  Unprotected areas in Namibia are largely dominated by livestock and game farms.  As obligate scavengers, vultures may rely on livestock mortality as well as ungulates killed by predators and for trophy hunting at game farms, a growing industry in Namibia's unprotected land^[@CR31]^. As white-backed vultures are social foragers that roam large areas in search of food^[@CR14],[@CR20]^, they may more readily take advantage of the carrion availability offered by livestock and game farmland as compared to the lappet-faced vultures that typically forage in pairs^[@CR20]^. These contrasting foraging strategies may likely contribute to explain the result reported above.  This assertion is also in line with previous findings from South Africa suggesting that white-backed vultures can readily take advantage of carrion provided by expanding game and cattle farms^[@CR32]^.

Climatic conditions were found to be important determinants of nestling body condition and breeding success in a number of species^[@CR11],[@CR33],[@CR34]^, with reported effects of climate being direct, e.g. causing physiological stress^[@CR11]^, or indirect, causing changes in food availability^[@CR34],[@CR35]^. Ungulate populations of the African savannah are strongly regulated by primary productivity, which is mediated by climatic factors, such as rainfall, and in turn impacts on the food base for vultures, with repercussions on their movement, reproduction and survival (see e.g.^[@CR36]^).  In arid savannah landscapes, such as those of most of Namibia, below average rainfall likely triggers additional ungulate mortality, due to reduced forage availability^[@CR22],[@CR37],[@CR38]^. In this study, NDVI was used as a proxy for ungulate forage availability and is presumably inversely correlated with their mortality.  The increased ungulate mortality resulting from below average forage availability (i.e. low NDVI) likely leads to increased carrion availability for the lappet-faced vultures and improved nestling condition in those years. 

Putting our findings into a broader context of longer time periods, we note that although low ungulate forage availability within protected areas may improve lappet-faced vulture nestling condition in the short-term through the boosted food availability from ungulate mortality, this benefit may be ephemeral.  Arid savannah landscapes, such as those of Namibia, are expected to experience heavy droughts under climate change^[@CR3],[@CR5],[@CR39]^. In the long run, this may cause desertification, with consequent reduction in the carrying capacity of the environment and overall declines of ungulate numbers^[@CR5],[@CR40]^, ultimately reducing the food base for vultures. 

The fact that NDVI (i.e. ungulate forage availability) did not seem to affect white-backed vulture nestling condition is surprising.  It is possible that adult white-backed vultures are able to better buffer against climate-driven environmental variability as compared to the lappet-faced vultures.  For example, white-backed vultures typically roam across larger foraging areas, which may support adequate food resources during favorable and adverse years, as compared to the lappet-faced vultures^[@CR20]^. White-backed vultures also rely on larger animals as a food source compared to lappet-faced vultures^[@CR14]^. Larger ungulates may be more resistant to harsh environmental conditions (e.g. decreased food quality) as compared to smaller ungulates which may be also impacted by increased predation in dry years^[@CR40]^. Consequently, years of poor forage availability for ungulates may selectively increase small ungulate mortality, and favor the lappet-faced vulture more than the white-backed vulture.  Moreover, and perhaps most importantly, the white-backed vultures are social foragers, which enables them to scan larger areas more efficiently as compared to the lappet-faced vultures that typically forage in pairs^[@CR14],[@CR20]^. This feature of the white-backed vultures may allow them to efficiently gather around and take advantage of already discovered carcasses, irrespective of climate-driven yearly variation in carrion abundance.  The finding that white-backed vulture nestlings seem to buffer their body condition against varying environmental conditions requires further ad-hoc investigations beyond the scope of this study. 

Interestingly, body condition of lappet-faced vulture nestlings was strongly affected by the interaction between protected area cover and climate-driven ungulate forage availability.  Nestling body condition was most improved after years of poor forage availability for ungulates but only in largely protected sites.  Indeed, the effect of variation in forage availability on nestling body condition within protected areas is rather large, amounting to a gain of approximately one kg of body mass when environmental conditions change from the highest to the lowest forage availability (corresponding to highest to lowest NDVI values in Fig.  [3](#Fig3){ref-type="fig"}).  The gain of one kg described above appears remarkable considering the average weight of just above six kg of the lappet-faced vulture nestlings considered here.  We propose that environmental stressors for ungulates may be most pronounced within ecosystems found in protected areas as compared to man-managed unprotected landscapes.  In the latter, water and pasture scarcity may be buffered by the food and water supplementation or by artificial reduction in stock density.  Thus, while in protected areas lower forage availability may lead to additional available carrion for local lappet-faced vultures (as suggested by the improved nestling condition under low NDVI), this is not the case in farmland dominated unprotected landscapes. 

Outside of the protected areas of the world, the human footprint is increasing, causing changes across species, ecosystems and habitats^[@CR4]^, including space use patterns of various species^[@CR41]^. Under these changing conditions, our finding that lappet-faced vultures in unprotected land did not benefit from increased resource availability like their counterparts in protected areas highlights an important prediction for future dynamics of this species.  Predicted changes in land-use (i.e., farmland expansion and/or changes in farm management practices) accompanied by climate change will further restrict species ability to benefit from harsh conditions and to cope with environmental change.  This suggests that while protected areas are very important under current conditions, in the near future they will play an even more crucial role in preserving healthy ecosystems where vultures can benefit from periods of favorable conditions (and unfavorable for ungulates). 

In the long term, however, climate change may trigger vulture distribution range shifts, pushing them progressively outside of the current protected area network, as shown by a recent study on the cape vulture (*Gyps coprotheres*)^[@CR42]^. Therefore, an immediate challenge is to ensure the preservation of healthy ecosystems within the current protected areas and address the threats to vultures in unprotected landscapes dominated by livestock and game farms.  Here, the use of poisons to kill carnivores^[@CR43],[@CR44]^ as well as environmental contaminants^[@CR45]^, such as lead present in hunted game, represent impelling threats that need to be urgently addressed.  This would allow for the short as well as long term conservation of vulture populations across arid savannah landscapes under rapid environmental change. 

This work represents one of the very few cases (but see^[@CR33]^) where climate effects, acting in synergy with land-use, are shown to impact nestling body condition, with potential, yet currently unknown, demographic consequences.  Juvenile vultures are typically characterized by low survival as compared to adults^[@CR46],[@CR47]^. Nestling body condition was found to have carryover effects on survival of several species^[@CR48],[@CR49]^. Therefore, the improved nestling condition of lappet-faced vultures in protected areas under dry conditions, and the importance of livestock and game farmland for white-backed vulture nestlings, carry important conservation implications.  They highlight the need for implementing conservation actions within and outside protected areas to maintain viable vulture populations in the long term.  A broad landscape level approach is crucial for species like obligate scavengers that roam very large areas, often well beyond the protected area networks^[@CR19]^.

The worsening trend in body condition of white-backed vulture nestlings over the 14 years study period may be explained by deteriorating environmental conditions, particularly outside of protected areas.  Here, while the recent expansion of the game farm industry may provide an increased food source for vultures, it may also negatively affect their condition through increased exposure to lead from hunting ammunitions^[@CR45]^. Another factor that may explain the negative trend in body condition is the progressive bush encroachment (i.e. the increase in woody vegetation cover), that may reduce the detectability and accessibility of carcasses to vultures.  This growing phenomenon is occurring particularly across Namibia's livestock farmland^[@CR50]^ and was indicated as one of the causes of the decline of the national cape vulture population^[@CR51]^. The fact that the decline in body condition is strongest for the white-backed vulture as compared to lappet-faced vulture nestlings may be due to their different foraging ecology and recent population trends.  Populations of these two vulture species have declined over the past decades in the study region^[@CR43]^, as well as across Africa, largely due to poisoning^[@CR52]^. These declines in number may have affected the foraging efficiency of white-backed vultures, which are social foragers and rely on social information to detect carrion to a much larger extent than the lappet-faced vultures^[@CR20]^. The resulting decline in foraging efficiency of the white-backed vultures when its numbers become scarce likely translates in less food provided to the nestling.  This may contribute to the decline in body condition observed here.  The direct causal link portrayed above would however warrant further ad-hoc investigations. 

While we show here that vulture nestlings' body condition may vary according to the coverage of protected areas and environmental conditions, our correlative study does not allow to establish causation in the patterns found.  Nevertheless, these results underscore that the variation in environmental conditions, in combination with land-use, may trigger measurable changes in nestling body condition, as shown for the lappet-faced vulture (Fig.  [3](#Fig3){ref-type="fig"}).  These changes may carry over to the successive stages of the birds' life, but the consequences are largely unknown.  However, we can speculate that an excess of one kg in body mass may make a disproportionate difference for a nestling of a species that largely relies on occasional and unpredictable food, such as carrion.  In such extreme conditions, extra body resources may allow a nestling to survive through periods of food scarcity.  Thus, our findings can pave the way to future work which could explore the carry-over effects of nestling body condition on survival, recruitment and ultimately fitness, ideally using a semi-natural experimental design.  Such investigations would require very long-term and individual-based data given the long lived and slow reproducing features characterizing vultures.  They would however be particularly relevant in light of the observed worsening I body condition of the white-backed vulture, an IUCN (International Union for Conservation of Nature; [www.iucn.org](http://www.iucn.org)) Critically Endangered species.  Understanding potential demographic consequences of the observed trend in nestling condition will allow informing the conservation management for this threatened species. 

Information on nestling sex was also not available in the dataset used here.  Sexing vultures would require DNA extraction, because the sexes are phenotypically alike in both studied species^[@CR14]^. Biased sex ratio at hatching has been found to be related to environmental conditions in broods of species with a clear sexual-size dimorphism, such as the Tengmalm's owl (*Aegolius funereus*)^[@CR53]^. Given the lack of sexual-size dimorphism in the two vulture species considered here, we deem unlikely that nestling sex would represent a confounding factor in our analyses. 

Climate-driven environmental change and land-use change may act independently but often interact in their impact on biodiversity^[@CR30]^. Our findings carry timely implications for preserving healthy savannah ecosystems through the adoption of adaptive management strategies that would ensure the preservation of viable vulture populations within and beyond protected areas under global change.  Ultimately, this can only be achieved if the drivers of change, and their interactive impacts, are well understood. 

Understanding the drivers of change and the mechanisms associated to them is crucial for predicting and managing impacts on biodiversity under global change.  Here we provide evidence of species-specific interactive impacts of different drivers of change on the life-history of a critical functional guild^[@CR15]^. The finding that two important drivers (such as land-use and climate-mediated vegetation productivity) act in synergy in their influence on lappet-faced vulture nestling condition further highlights the importance of considering interactive effects of multiple drivers of change. 

Methods {#Sec4}
=======

Ethics Statement {#Sec5}
----------------

The work was conducted in accordance with all relevant national and international guidelines.  The handling and measuring of nestlings were carried out by experienced bird ringers holding a valid SAFRING (South African Bird Ringing Unit) ringing license approved by the Namibian Ministry of Environment and Tourism, and following the guidelines for ringing provided by SAFRING (http://safring.adu.org.za/). 

Study species and study landscape {#Sec6}
---------------------------------

The study was based in Namibia, a country dominated by arid savannah landscapes, with low yearly rainfall largely concentrated during the period October to March^[@CR54]^. While most of Namibia is covered by sparse savannah woodland, the coastal areas form a barren and typically sandy desert, such as that protected within the Namib-Naukluft National Park.  In this study we focused on three discrete regions where most of the ringing data of vulture nestlings are available (see Fig.  [1](#Fig1){ref-type="fig"}).  A northern region (hereafter named Etosha), is partly included within the Etosha National Park.  A central region (hereafter Central) which lies on the highlands of central Namibia, and is largely represented by commercial livestock farmland.  A third region is partly included within the Namib-Naukluft National Park (hereafter named Namib; Fig.  [1](#Fig1){ref-type="fig"}).  Unprotected areas within the study region are typically managed for livestock farming, game farming or a mixture of the two, whereas crop farming is scarce in these regions due to the arid conditions^[@CR54]^. Part of these unprotected areas belong to commercial farmers, and part are used by communal subsistence farmers which typically belong to conservation conservancies^[@CR8]^.

We focused on the two most common breeding species of vultures in Namibia, the lappet-faced and the white-backed vulture^[@CR55]^. The sexes of both species are phenotypically alike and can only be separated by means of DNA analyses^[@CR14]^. These two species typically nest on trees in the dry savannah, mostly represented by Acacia species^[@CR55]^. The lappet-faced vulture and the white-backed vulture are classified as Endangered and Critically Endangered by IUCN owing to their rapid population declines that are likely to continue into the future^[@CR52]^. The main threats to both species are intentional and unintentional poisoning, habitat degradation and collision mortality with infrastructures^[@CR43],[@CR44],[@CR52]^. In Southern Africa, breeding occurs during the Austral winter, with hatching taking place during the middle of the dry season, typically between July and August.  Both species lay one egg and the nestling period lasts about four months^[@CR14]^. This period coincides with progressively drying conditions in Namibia, which typically exacerbate towards the second half of the nestling period, between September and December^[@CR54],[@CR55]^. In this way, vultures synchronize the most energy demanding phase of their life cycle with a pulse in food availability resulting from increased ungulate mortality induced by low availability of nutritious vegetation during dry periods^[@CR22],[@CR36]^.

Nestling measures and calculation of body condition {#Sec7}
---------------------------------------------------

We use nestling ringing data for the period 2003--2016 when ringing effort was high and constant throughout.  Nests have been typically located either through systematic surveys from the ground or from the air (as in the case of Etosha and the Namib regions).  Alternatively, nests were located by opportunistic sightings of the nest with adults by the ringer, farmer or farm-workers (as in most of the observations for the Central study area).  Nestlings were ringed at the age between two months and before fledging.  Ringers typically reached the nest using a ladder.  Handling involved placing a metal ring on the leg and measuring the wing length and the weight of the individual chick.  The chick was then placed back to the nest, with the whole procedure lasting half an hour on average.  Overall, the dataset includes 608 white-backed and 899 lappet-faced vulture complete data points (i.e. including body weight, wing length and coordinates of each observation).  Of these, 292 and 316 white-backed vulture data points were available for the Central and Etosha regions, respectively.  No ringing data for this species are available for the third region, the Namib, as white-backed vultures did not breed there.  For the lappet-faced vulture, 730, 26 and 143 ringing data points were available from the Namib, Central and Etosha regions, respectively. 

Body condition was estimated for each nestling from the mass and wing length measures and using the scaled mass index method^[@CR56]^. This index scales each individual body mass to the expected value if all nestlings had the same body size.  The scaled mass index is deemed highly robust to measurement errors and allows reliable comparability of body condition from measurements taken on subjects with different age^[@CR56]^. It also provides a more reliable indication of the true condition of an animal compared to other approaches, such as the common one based on residuals calculation from an ordinary least square regression of body mass against length^[@CR56]^. The formula for calculating the scaled mass index used here is as follows: body mass × (464.16/wing length)^0.69^ for the lappet-faced vulture nestlings, and body mass × (410.94/wing length)^0.65^ for the white-backed vulture nestlings.  In these formulas, the values within brackets represent the mean wing length of the population of each species (in mm) to which all body mass values are adjusted to, respectively, and the upper case coefficients are the slopes of standardized major axis regressions coefficients (see^[@CR56]^ for more details).  In practice, because the index scales the body mass of each individual to the average body size of the population under study, the scaled mass index can be compared among individuals for which the measurements were taken at different ages.  The scaled mass index has been successfully used to estimate body condition of adults and nestlings of several species (see e.g.^[@CR57]^).  For simplicity and clarity, we will use the term body condition instead of scaled mass index throughout. 

Predictors of body condition {#Sec8}
----------------------------

In order to test the above mentioned hypotheses we extracted the cover of protected area (hereafter PA; see below) and normalized difference vegetation index (NDVI, here used as a proxy for climate-driven ungulates forage availability; see below) from an area with a circular buffer of 50 km radius around the nest location.  The rationale for choosing the 50 km radius was based on known movement ecology of the two study species in Namibia indicating that most of the daily foraging movements occur within this distance^[@CR20]^. Although this is a somewhat arbitrary spatial threshold, we deem this distance to represent an optimal compromise between type I and II errors given the available data quality.  On one hand, it is a large enough radius to likely incorporate most of the foraging area of the two species during the breeding season (presumably minimizing omission of utilized area, i.e. Type I error).  On the other hand, it is restricted enough so that it includes only areas that are directly relevant for foraging, avoiding areas not utilized by the feeding parents (thereby confounding and possibly masking some genuine patterns, i.e. Type II error).  We found very strong spatial autocorrelation in the NDVI data (see Fig.  [1](#Fig1){ref-type="fig"} right panel) and thus we argue that radii of different sizes would most likely lead to very similar results as those obtained with the current selection. 

Protected areas were extracted as polygon files for GIS software from the World Database on Protected Areas (WDPA; <https://protectedplanet.net/>).  We selected only areas belonging to the IUCN protected area categories from I to IV that are currently designated (largely represented by National Parks, see Fig.  [1](#Fig1){ref-type="fig"} left panel).  PA cover was calculated as the percentage area protected within the 50 km buffer around each nest.  Due to restrictions on resource extraction (including hunting) and human land-use (e.g. farming), protected areas typically support ecosystems at a more "natural" state as compared to unprotected sites.  In PAs wild ungulate populations are at the mercy of the environmental variation (e.g. prolonged dry periods; although localized water points may exist in some of the PAs) and carcasses are left *in situ*, thus providing food for vultures, as compared to unprotected land (largely dominated by livestock farmland where water provision is regular) where dead animals may be removed^[@CR32],[@CR54]^. Thus, protected areas may positively affect the body condition of vulture nestlings. 

In addition, we also calculated the mean NDVI value for the same 50 km buffer area described above.  NDVI offers a remotely sensed measure of the greenness of the vegetation, which is in turn related to primary productivity and forage availability for wild and domestic ungulates^[@CR22],[@CR42],[@CR58]^. NDVI has also been related to ungulate population demography, including increased mortality under prolonged droughts^[@CR59],[@CR60]^. This may also apply to domestic livestock in Namibia^[@CR44]^. Therefore, we calculated NDVI for five different time periods, spanning backwards from nestling measurement date for one month, three months, one, two and three years.  This temporal gradient in the lag of NDVI effects allowed to explore whether prolonged periods of below average NDVI values, which means prolonged drought, positively correlate with vulture nestling condition through increased wild and domestic ungulate mortality.  The original NDVI values used in this study were available at 16-day composite values with a resolution of 250 m X 250 m and obtained from^[@CR61]^.

Statistical analyses {#Sec9}
--------------------

We used Bayesian hierarchical models to examine the relationship between nestling body condition and PA coverage and the five NDVI timescales, all as continuous variables.  Year was included as a continuous covariate in all models (in this case it was named *time* for clarity, see Table [1](#Tab1){ref-type="table"}) in order to account for unexplained temporal trends in the data.  Region identity was included as random factor (Etosha, Central and Namib).  This random structure was defined a priori to account for pseudo-replication due to potential measurements of nestlings from the same pair within the same region.  To deal with the potential confounding effects of spatial autocorrelation, we also included in the models a spatially correlated random intercept term.  The random intercept is assumed to be a Gaussian Markov Random Field (GMRF) with mean 0 and covariance matrix ∑. Because of the difficulty in estimating a covariance matrix for large datasets, an approach based on Continuous Domain Stochastic Partial Differential Equations (SPDE) is used to calculate the covariance matrix ∑ for the GMRF^[@CR62]^. This is a computationally efficient approach for large datasets, as it avoids the computation of a full covariance matrix for the GMRF.  The models were fitted using integrated nested Laplace approximations (INLA).  INLA was chosen as it allows efficient estimation regression parameters within a Bayesian framework, without the need to employ computationally intensive Markov Chain Monte Carlo algorithms.  Default priors were assigned for all fixed-effect parameters as recommended by^[@CR63]^, which are approximations of non-informative priors designed to have little influence on the posterior distribution.  Specifically, the default prior for fixed effects is specified by a Gaussian distribution with mean of 0 and a precision parameter τ = 0.001.  The prior for region hyperparameter is a log-gamma distribution with shape parameter of a = 1 and scale parameter b = 0.00005.  The spatial random effect component is defined in terms of two hyperparameters, κ and τ which are related with the range and scale of the spatial effect.  Their priors are centred values such that the range is about 20% of the diameter of the region and the variance is equal to 1, see^[@CR62]^.

Prior to analyses, we scaled all predictor variables (i.e. the five NDVI variables, PA and year), to their overall mean.  We also scaled and centered the response variable to aid convergence, because the original body condition values were in the order of thousands.  We reiterate here that while body measurements were taken from nestlings of different age, there is no need to account for age in the models because age was already explicitly taken into account in the response when calculating the scaled mass index (see above). 

We ran separate models for each of the two vulture species.  There was a relatively high degree of collinearity (r \> 0.7) among the NDVI variables and, for the white-backed vulture only, between NDVI variables and PA coverage (Fig.  [S1](#MOESM1){ref-type="media"} in Supplementary Methods).  Thus, each predictor was tested in turn in a separate model where year was always included as continuous (variable *time*, see above), and with a random structure as detailed above.  We also ran a separate model to test for possible interactive effects between each of the five NDVI variables and PA coverage for the lappet-faced vulture data where collinearity between NDVI and PA coverage was low (i.e. r \< 0.7; Fig.  [S1](#MOESM1){ref-type="media"} in Supplementary Material).  This was done because the magnitude of impact of changes in NDVI on nestling body condition may vary according to the PA coverage.  This interaction was not tested for the white-backed vulture data due to the strong collinearity (r \> 0.7) between PA coverage and NDVI in these data (Fig.  [S1](#MOESM1){ref-type="media"} in Supplementary Material).  Residual spatial autocorrelation was assessed by means of investigating correlograms, but no signs of spatial autocorrelation were found (see Fig.  [S4](#MOESM1){ref-type="media"} in Supplementary Material).  Analyses were performed in R software (version 3.4.0^[@CR64]^).  The cover of PA around each nest was extracted using ArcGIS 10.1.0 software (ESRI). 

Electronic supplementary material
=================================

 {#Sec10}

Supplementary Material

**Publisher\'s note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. 

Electronic supplementary material
=================================

**Supplementary information** accompanies this paper at 10.1038/s41598-018-31344-2. 

We thank Wilferd Versfeld, Holger Kolberg and Mark Boorman, along with the numerous vulture ringers for gathering and collating the ringing data.  We also thank the Ministry of Environment and Tourism of Namibia, as well as several sponsors, for providing financial support and logistics for the ringing activities.  For financial support of this study, we also thank the Kone foundation and Academy of Finland (Grant No.  307909 to A.S.).  O.S. was supported the Israeli Ministry of Science and Technology Shamir Postdoctoral Fellowship.  Two referees are also acknowledged for their insightful comments that helped improve this work. 

A.S. conceived the study, A.S. and P.B. collated the data, A.S., M.G. and O.S. designed the study, M.G. analysed the data.  A.S. and M.G. prepared the figures and tables.  A.S. wrote the manuscript, with contributions from O.S., M.G. and P.B.

The dataset generated and analysed during the current study is available in the figshare repository, 10.6084/m9.figshare.6876752. 

Competing Interests {#FPar1}
===================

The authors declare no competing interests.