Published Date
For further details log on website :
http://www.nature.com/ncomms/2016/160728/ncomms12306/full/ncomms12306.html
- Nature Communications
- 7,
- Article number:
- 12306
- doi:10.1038/ncomms12306
- Received
- Accepted
- Published
Data
For each sampling location or site in the PREDICTS database (November 2014 version), we calculated within-sample species richness, total abundance of individuals, rarefied richness (based on the fewest individuals at any site within each study) and community weighted mean log10geographic range size—the inverse of which was then plotted to give our endemicity measure. Each species’ range size was derived from its global occurrence in the Global Biodiversity Information Facility database. We recognise biases in the Global Biodiversity Information Facility data, but these are mitigated to some extent by our hierarchical modelling approach and our estimates compare reasonably well with estimates based on other data sources, listed in full in the Supplementary Information. Land use was classified using the study authors’ description for each site; this method has been shown to be repeatable33. Sites were considered to be protected if their geographical coordinates fell inside protected areas from the World Database on Protected Areas34 (see Supplementary Methods). We then derived two datasets: the first included all studies with sites inside and outside protected areas (all-sites data; Fig. 1b); the second retained only those sites from each study for which land use could be matched across the protected area boundary (matched-sites data 2; Fig. 1d). All sources of biodiversity data are listed in the Supplementary Information.Analyses
We used generalized linear mixed-effects models to account for differences in response variables due to study-specific methodologies and the spatial structure of sites46. The PREDICTS data present a rare opportunity to compare sites inside and outside protected areas, but do not have the geographic coverage required for a stricter counterfactual approach14, 17 in which sites are individually matched. To reduce the risk that any differences observed between sites inside and outside protected areas were caused by biases in the location of protected areas27, we considered elevation47 and derived slope at c. 1 km2 resolution and agricultural suitability48 at 10 km2 resolution as covariates in all models (see Supplementary Information for further details). To ensure independence of all variables in the model, we intentionally included only these three confounding variables that we considered to be fully independent of the presence of a protected area. For example, distances to roads and markets are affected by the presence of protected areas so are not independent confounding factors (see Supplementary Information for details). We sequentially compared models with and without each fixed effect and at each step dropped the term with the highest P-value, until all terms had P<0.05 (ref. 49).Assessing protection effects
We tested for biodiversity differences between sites inside and outside protected areas using the all-sites data, treating protection status (inside vs outside a protected area) as a fixed effect. We then tested whether biodiversity measures differed between management category groups by re-coding IUCN category as a four-level factor: unprotected, IUCN category III–VI, IUCN category unknown, and IUCN category I and II.Assessing protection effects within and among land uses
We used two approaches to test whether biodiversity differences between protected and unprotected sites varied with land use. First, using the all-sites data, we modelled the response of each biodiversity measure to protection status, land use, and their interaction. We also tested for the three-way interaction between land use, protection and either use intensity, latitudinal zone or taxonomic group. Second, using the matched-sites data, we re-ran models with protection status, and then with management category group as a fixed effect. We also split the matched-sites data by latitudinal zone and taxonomic groups to assess whether these factors influenced the effect of protection. Finally, we tested whether the site-level biodiversity response to protection varied with the size/age class of the protected area [four-level factor with all combinations of young (<20 years), old (20–85 years), small (<400 km2) and large (400–12,000 km2); these thresholds between categories were selected to give a similar number of sites in each group].Estimating global protected area effectiveness
The global effectiveness of protected areas (e) was estimated from e=1−(1−i)/(1−o), where modelled site-level biodiversity inside (i) and outside (o) protected areas are expressed as a proportion of that under ‘pristine’ conditions. We calculated the ratio of i/o from the model estimates for biodiversity inside relative to outside protected areas in each land use (Fig. 3), where each land-use parameter was weighted by the proportion of global terrestrial area within that land-use type. This value of i/ocould then be used to solve an equation expressing the global state of site-level biodiversity: 1−r=ai+(1−a)o, where r is the estimated global average loss of site-level biodiversity relative to pristine46 and a is the fraction of the total land area that is protected50. Solving this equation for i and o allowed us to estimate e. Finally, by using estimates for the effect of protection in IUCN categories I and II (Fig. 2a,b) to give i/o, we estimated e under the more restrictive management scenario. By rearranging the equations we estimated the total protected area (a) needed to obtain the same average local biodiversity outcome (1−r) inferred under this more restrictive management scenario. See Supplementary Information for more details.Data availability
The biodiversity data that support the findings of this study are available in the Natural History Museum data portal (data.nhm.ac.uk) with the identifier dx.doi.org/10.5519/0095544. R scripts are available at http://github.com/claudialouisegray/PREDICTS_WDPA.Additional Information
- How to cite this article: Gray, C. L. et al. Local biodiversity is higher inside than outside terrestrial protected areas worldwide. Nat. Commun. 7:12306 doi: 10.1038/ncomms12306 (2016).
- References
- Effectiveness of terrestrial protected areas in reducing habitat loss and population declines. Biol. Conserv. 161, 230–238 (2013). et al.
- http://www.unep-wcmc.org/resources-and-data/protected-planet-report-2014 (UNEP-WCMC, 2014). et al. Protected Planet Report 2014. Available at
- CBD. Decision X/2, The strategic plan for biodiversity 2011–2020 and the Aichi Biodiversity Targets, Nagoya, Japan, 18 to 29 October 2010. Available at http://www.cbd.int/decision/cop/default.shtml?id=13164 (2010).
- A mid-term analysis of progress toward international biodiversity targets. Science 346, 241–244 (2014). et al.
- A global analysis of protected area management effectiveness. Env. Manag. 46, 685–698 (2010). , , , &
- Averting biodiversity collapse in tropical forest protected areas. Nature489, 290–294 (2012). et al.
- Mapping change in human pressure globally on land and within protected areas. Conserv. Biol. 28, 1604–1616 (2014). , &
- Large mammal population declines in Africa’s protected areas. Biol. Conserv. 143, 2221–2228 (2010). et al.
- WWF. Living Planet Report 2014: species and spaces, people and places. Available at http://wwf.panda.org/about_our_earth/all_publications/living_planet_report. (WWF, (2014).
- The performance and potential of protected areas. Nature 515, 67–73 (2014). , , &
- Creating a safe operating space for iconic ecosystems. Science 347, 1317–1319 (2015). et al.
- Financial costs of meeting global biodiversity conservation targets: current spending and unmet needs. Science 338, 946–949 (2012). et al.
- Global-scale mapping of economic benefits from agricultural lands: Implications for conservation priorities. Biol. Conserv. 140, 40–49 (2007). &
- Global protected area impacts. Proc. R. Soc. B 278, 1633–1638(2011). &
- Increasing isolation of protected areas in tropical forests over the past twenty years. Ecol. Appl. 15, 19–26 (2005). , , &
- On the protection of ‘protected areas’. Proc. Natl Acad. Sci. USA 105, 6673–6678 (2008). , &
- Governance regime and location influence avoided deforestation success of protected areas in the Brazilian Amazon. Proc. Natl Acad. Sci. USA 110, 4956–4961 (2013). , , &
- Securing tropical forest carbon: the contribution of protected areas to REDD. Oryx 44, 352–357 (2010). et al.
- Effectiveness of strict vs. multiple use protected areas in reducingtropical forestfires: A global analysis using matching methods. PLoS ONE 6, e22722(2011). &
- Dudley N. (ed). Guidelines for Applying Protected Area Management Categories IUCN (2008).
- Analysis of deforestation and protected area effectiveness in Indonesia: a comparison of Bayesian spatial models. Glob. Environ. Change 31, 285–295 (2015). et al.
- Local scale comparisons of biodiversity as a test for global protected area ecological performance: a meta-analysis. PLoS ONE 9, e105824(2014). , &
- More strictly protected areas are not necessarily more protective: evidence from Bolivia, Costa Rica, Indonesia, and Thailand. Environ. Res. Lett. 8, 025011 (2013). et al.
- Protected area types, strategies and impacts in Brazil’s Amazon: public protected area strategies do not yield a consistent ranking of protected area types by impact. Phil. Trans. R. Soc. B 370, 20140273 (2015). , , &
- Extent and ecological consequences of hunting in Central African rainforests in the twenty-first century. Phil. Trans. R. Soc. B 368, 20120303 (2013). , , , &
- Quantifying biodiversity: procedures and pitfalls in the measurement and comparison of species richness. Ecol. Lett. 4, 379–391 (2001). &
- High and far: biases in the location of protected areas. PLoS ONE 4, e8273 (2009). &
- Protection reduces loss of natural land-cover at sites of conservation importance across Africa. PLoS ONE 8, e65370 (2013). et al.
- Protected area effectiveness in reducing conversion in a rapidly vanishing ecosystem: The Brazilian Cerrado. Conserv. Lett. 7, 216–223(2014). , , &
- Measuring the effectiveness of protected area networks in reducing deforestation. Proc. Natl Acad. Sci. USA105, 16089–16094 (2008). , , , &
- Roads, land use, and deforestation: a spatial model applied to Belize. World Bank Econ. Rev. 10, 487–512 (1996). &
- Determinants of deforestation and the economics of protection: an application to Mexico. Am. J. Agric. Econ. 84, 943–960 (2002). &
- The PREDICTS database: a global database of how local terrestrial biodiversity responds to human impacts. Ecol. Evol. 4, 4701–4735 (2014). et al.
- IUCN and UNEP. The World Database on Protected Areas (WDPA), July 2014. (Available at www.protectedplanet.net (UNEP-WCMC, 2014).
- Change the IUCN protected area categories to reflect biodiversity outcomes. PLoS Biol. 6, e66 (2008). et al.
- Biodiversity hotspots for conservation priorities. Nature 403, 853–858 (2000). , , , &
- Used planet: a global history. Proc. Natl Acad. Sci. USA 110, 7978–7985(2013). et al.
- Aligning conservation priorities across taxa in Madagascar with high-resolution planning tools. Science 320, 222–226 (2008). et al.
- Aichi targets: protect biodiversity, not just area. Nature 526, 195 (2015).
- Making parks make a difference: poor alignment of policy, planning and management with protected-area impact, and ways forward. Phil. Trans. R. Soc. B 370, 20140280 (2015). , &
- Targeting global protected area expansion for imperiled biodiversity. PLoS Biol. 12, e1001891 (2014). et al.
- Synergies and trade-offs in achieving global biodiversity targets. Conserv. Biol. 30, 189–195 (2016). et al.
- Replacing underperforming protected areas achieves better conservation outcomes. Nature 466, 365–367 (2010). et al.
- Global biodiversity indicators reflect the modeled impacts of protected area policy change. Conserv. Lett. 9, 14–20 (2015). et al.
- Bolder thinking for conservation. Conserv. Biol. 26, 1–4 (2012). et al.
- Global effects of land use on local terrestrial biodiversity. Nature 520, 45–50(2015). et al.
- Global multi-resolution terrain elevation data 2010 (GMTED2010). US Geological Survey Open File Report 2011–1073 (2011). &
- Global Agro-Ecological Assessment for Agriculture in the 21st Century: Methodology and Results. Plate 46: Suitability for rain-fed crops (maximizing technology mix). International Institute for Applied Systems Analysis and Food and Agriculture Organization of the United Nations (2002). , , &
- http://lme4.r-forge.r-project.org/lMMwR/lrgprt.pdf (2010). lme4: Mixed-effects modeling with R. Available at
- Shortfalls and solutions for meeting national and global conservation area targets. Conserv. Lett. 8, 329–337 (2015). et al.
- http://thematicmapping.org/downloads/world_borders.php (2016). World Borders Dataset 0.3. Available at
- We thank the hundreds of data contributors; all PREDICTS project volunteers, masters and PhD students that collated records; and Adriana De Palma, Helen Phillips, Diego Juffe-Bignoli, Neil Burgess, Max Gray, Daniel Ingram, Valerie Kapos, Naomi Kingston, Sarah Luke and the protected areas team at UNEP-WCMC for comments and assistance. We thank the School of Life Sciences at the University of Sussex for support and the Natural History Museum for a GIA travel award. The PREDICTS project is funded by the UK Natural Environment Research Council (NERC, grant number: NE/J011193/2). PREDICTS is endorsed by the Group on Earth Observations Biodiversity Observation Network (GEO BON). This is a contribution from the Imperial College Grand Challenges in Ecosystem and the Environment Initiative, and the Sussex Sustainability Research Programme.
Competing financial interests
The authors declare no competing financial interests.
For further details log on website :
http://www.nature.com/ncomms/2016/160728/ncomms12306/full/ncomms12306.html
No comments:
Post a Comment