1. Introduction
Protected Areas (PAs) are a key strategy for safeguarding global biodiversity and ecosystem services. As of 2018, there were more than 230,000 terrestrial PAs worldwide, protecting 14.9% of the earth’s surface and inland waters outside of Antarctica [
1], and 5.2% of global tree cover (
Figure 1). The extent of PAs has increased substantially since the 1990s, and under the Convention on Biological Diversity, nations committed to further increasing the land area in PAs to 17% by 2020. Despite this, PAs are under increasing threat from anthropogenic activities, including encroachment for settlements, agriculture, mining, logging, and poaching and bushmeat hunting [
2]. Worldwide almost one-third of PAs are under intense human pressure, determined as the combined influence of built environments, agriculture, human population, and transportation infrastructure [
3]. Moreover, less than half of PAs are free of any human pressure, and this pressure has increased since the 1990s [
4].
Globally, between 2001 and 2012, 3% of PA forests and 5% of all forests were converted to other land cover types [
6]. While the lower rate of forest loss in PAs relative to the global average may suggest that PAs are effective in preventing some, if not all, forest loss, several studies have shown that PAs are preferentially located in areas that have a lower risk of deforestation [
7,
8,
9,
10]. The locations of PAs are biased towards areas with lower potential agricultural revenues and limited access, in order to minimize conflict with extractive industries and thus reduce the cost of acquisition and establishment [
11]. Nevertheless, studies controlling for these confounding factors generally demonstrate that PAs do provide additional protection beyond what would have been expected in the absence of their designation [
12,
13].
Here, we examine global tree cover loss in PAs globally over 2001–2018, providing the most up-to-date report on forest conversion trends in PAs. Previous analyses of forest loss have been restricted to national or regional scales (e.g., References [
8,
14]), and/or have not been recently updated (e.g., Reference [
7]). We highlight regions and countries where PAs are succeeding and failing to prevent forest loss, a proxy for their ability to safeguard intact habitats and protect other provisioning and regulating ecosystem services. Next, we examine the land cover following tree cover loss in PAs, and variations in follow-up land cover over space and time. In particular, we measured the magnitude (extent) of forest to agricultural land conversion in PAs, as agriculture has been shown to be a dominant driver of deforestation in the tropics [
15,
16] and globally [
17]. Our analysis provides useful information about what may be causing forest loss in PAs and informs the development of PA management and enforcement strategies that are tailored to the agents of change on the ground. Our global assessment should be considered a screening tool to identify priority regions for further detailed investigation of threats to PAs.
2. Materials and Methods
To conduct our analysis, we took advantage of three recently published, or recently updated, spatially explicit datasets (
Table 1), (1) protected areas from the World Database on Protected Areas (WDPA) [
5], (2) 30 m resolution Global Forest Change (GFC) data representing tree cover loss annually, from 2001–2018 [
18], and (3) 300 m resolution land cover maps for the years 2005, 2010, and 2015 from the European Space Agency–Climate Change Initiative (ESA-CCI) [
19].
From the WDPA, we excluded marine PAs, PAs which have been proposed but not formally designated, and those without spatial information (e.g., only provided as point data). We examined forest loss trends in a given PA beginning the year after which it was formally designated, and in the ~10% of cases where the establishment year was not provided, we assumed that the PA had been established prior to 2001. We included all PA types in our analysis (
Supplementary Table S1), including those designated in the International Union for Conservation of Nature (IUCN) categorization system as “Not Applicable”, “Not Assigned”, and “Not Reported”. These categories include some important types of protection, including lands managed by indigenous communities in Brazil and United Nations Educational, Scientific and Cultural Organization (UNESCO) Biosphere Reserves in Guatemala. However, these may be left uncategorized according to the IUCN typology due to reporting errors. We isolate trends in PAs which are more strictly protected according to the IUCN categorization system (strict nature reserves, wilderness areas, and national parks), in recognition of the fact that less stringent PA categories may support forest management and other types of sustainable land use change that would result in forest loss. Within IUCN category IV, it is recognized that active management, or modifications to the ecosystem (e.g., halting natural succession, providing supplementary food, or artificially creating habitats) will take place. Specifically, IUCN Category IV sites allow sustainable management of natural resources to maintain culturally defined ecosystems with unique biodiversity, but are not designed for industrial harvest levels [
5], and Category VI areas allow the sustainable use of natural resources to promote ecosystem services. In the case where PA polygons overlapped, we assumed the most stringent level of protection. We converted the vector shapefile to a raster grid with spatial resolution of 30 m.
The GFC dataset [
18] mapped tree cover loss, defined as the conversion from forest to non-forest, during the 2000–2018 period. We refer to this as forest loss, under the assumption that loss in PAs is predominantly natural forest loss as opposed to loss of planted trees or plantations. We restrict our analysis to areas with greater than or equal to 50% canopy cover in the year 2000 [
18]. We tabulated forest loss through three time periods: 2001–2004, 2005–2009, and 2010–2014. We then categorized each forest loss pixel to a land cover class in the year immediately following each of these periods: 2005, 2010, and 2015, respectively. Previous research demonstrated that the land cover following forest loss does not change substantially in a 1–10-year period after the loss occurred. Our approach is based on a period of 1–4 years after forest loss [
20]. For a given PA, we excluded any tree cover loss that occurred prior to the year of PA establishment. We additionally report loss in PAs from 2015 to 2018, though we cannot assign a follow-up land cover to this loss, as the most recent land cover map is from the year 2015.
We use the land cover type following forest loss to categorize the cause of deforestation. We acknowledge that subsequent land cover is only a proxy for the complex and dynamic causes of deforestation, but more detailed investigation of these underlying causes is not possible at the scale of our analysis. We reclassified the 22 land cover categories presented in the European Space Agency Climate Change Initiative (ESA-CCI) land cover dataset to seven categories (
Supplementary Table S2) [
21]. Our reclassification schema consolidated forest categories (e.g., broadleaf tree cover, needleleaf tree cover), shrubland categories (e.g., shrubland, mosaic herbaceous cover), grassland, and ‘other’ land cover types (e.g., urban, bare land, water bodies, snow cover). We retained three separate agriculture categories: cropland (including both rainfed and irrigated cultivated crops), mosaic cropland (>50% cropland mixed with trees, shrubs, and herbaceous cover), and mosaic vegetation (<50% cropland and >50% mixed trees, shrubs, and herbaceous cover). The mapped cropland land cover categories have reported accuracies of 73%–89%, except for the mosaic vegetation category. This land cover type has a reported accuracy of just 59%, largely due to commissions of the other cropland categories. Notably, ESA-CCIs agriculture category includes areas used for crop cultivation but does not include areas used for livestock grazing (pasture land or managed grasslands). ESA-CCI does include a grassland category, but does not differentiate natural grasslands from managed grasslands, as this is difficult at global scales using mid-resolution satellite imagery [
22]. We resampled the landcover dataset to a 30 m pixel raster grid, matching the resolution of the forest loss map.
To gauge the robustness of our approach to classifying the land cover following forest loss using global-scale data, we compared our results to two previous studies which investigated drivers of deforestation using nationally and regionally specific datasets (
Supplementary Table S3). We aggregated several land cover categories in order to facilitate comparison according to
Supplementary Table S3. In Indonesia, the authors of Reference [
15] found that about 67% of deforestation nationally was followed by agriculture and 40% of deforestation events in PAs were followed by agriculture. We found a similar proportion of forest loss to agriculture nationally (61%) and in PAs (38%). In South America, the authors of Reference [
16] reported that 20% of deforestation was followed by agriculture and another 69% followed by pasture lands from 2001 to 2005. We estimated that 44% of forest loss was followed by agriculture, and another 44% was followed by grassland. The differences in Brazil may be partially explained by the challenge of differentiating natural grasslands from managed grassland or pastureland. The results of these robustness checks provide confidence that we can broadly track agriculture as a driver of forest loss using the global land cover dataset, though differentiating pastureland as a driver of forest loss remains difficult with currently available land cover maps (
Supplementary Tables S4 and S5). We address the implications of this challenge in more detail in the Discussion Section.
We examined case studies of forest loss in PAs in four countries: Brazil, Indonesia, Democratic Republic of Congo, and the US (
Supplementary Figure S1). We selected these countries because each has some of the highest forest loss globally and we wanted to include representation across continents and biomes (
Supplementary Table S4). In each case, we investigated forest loss trends in a PA with one of the highest rates of forest loss nationally. We visualized forest cover and loss for each selected PA, examined the Landsat time series from 2000 to 2016 in Google Earth, and available high-resolution satellite imagery in the PA over the study period. We do not aim to provide a systematic validation of our approach to tracking forest loss and following land cover in PAs, but rather use these case studies to explore areas of concern in more detail with higher resolution satellite imagery.
4. Discussion
Between 2001 and 2018, two trends took place simultaneously. The absolute area of protected tree cover increased due to countries designating additional land as PAs. Conversely, the annual rate of tree cover loss inside PAs nearly doubled during this time period. The highest loss in tree cover within PAs occurred in 2016, when 0.44% of protected forests experienced loss. Globally, it seems that forests in PAs face the same economic, natural, and social pressures as non-protected forests, as shown by the consistency in trends of forest loss between the two forests categories (
Figure 2). This is despite PAs being in areas which should experience fewer human pressures of deforestation [
9,
11].
In the tropics, the extent of forest loss in PAs increased notably over the study period, and occurred largely in Indonesia, Brazil, and the Democratic Republic of Congo. Global land cover maps demonstrate that shrublands and grasslands were the dominant land cover following forest loss in PAs in the tropics. Our case study analysis demonstrated that this loss corresponds to areas impacted by fires, for example in Indonesian peat lands, and may also correspond to pasture land, for example in Brazil. In many countries in the tropics, agriculture was also a dominant land cover following forest loss, particularly in Sub-Saharan Africa.
In the northern hemisphere, the United States, Canada, and Russia contribute large total amounts of protected tree cover loss. Unlike in the tropics, forest loss in PAs in these countries did not increase noticeably over time. But similarly, shrublands were also the dominant land cover following forest loss over the study period. It is likely that most of this loss corresponds to natural occurrences such as fire, pests, or storm blowdown, which have been shown to be dominant drivers of deforestation in these regions outside of PAs [
17]. Indeed, our US case study identified wildfire as a dominant driver of loss.
PAs that do not fit within IUCN categorization schema, which comprise roughly one-third of all PAs, have the highest rates of tree cover loss. These PAs include indigenous lands and UNESCO reserves. There is evidence that indigenous land tenure recognition is effective at preventing deforestation [
27]. On the other hand, the majority of uncategorized PAs do not have any active management authority. It is possible that a lack of clear authority over these PAs may be one reason for higher rates of tree cover loss in these uncategorized PAs as a whole.
An important limitation of our assessment of land cover following forest loss in PAs is the reliance on a global mid-resolution land cover dataset. We used a global approach to allow for direct comparisons across regions, and to identify specific regions (within and across countries) where further investigation is needed. However, global land cover datasets do not necessarily address land use and may struggle, for example, to differentiate grazing and pasture lands from shrublands or grasslands [
22]. This limits our ability to reliably track pasture expansion into PAs in some geographies where it is important, including in Brazil. Also, mosaic land cover classes such as shrubland/mosaic natural vegetation according to the global map could actually be agroforestry or mixed cropland/agroforestry. This limits our ability to track small-scale and mixed agriculture classes in geographies where those are dominant land cover transitions, such as Central Africa.
We used a case study approach to gauge availability and usefulness of the additional information available via high-resolution imagery from Google Earth to identify and track drivers of forest loss in PAs. This was not intended as a validation of the global approach, but rather an exploration of the potential utility of this emerging technology for more in-depth examination of drivers of forest loss in hotspots of deforestation or priority conservation areas. We found that the spatial and temporal resolution of the imagery available on Google Earth helped inform possible reasons for forest loss, including wildfire, small-scale agriculture, and pasturelands.
Finally, the ESA-CCI dataset represents land cover at 300 m resolution, so pixels with a small proportion of a given land cover category may not be represented, even if they are identified in the 30 m resolution tree cover loss dataset (also discussed in Reference [
20]). This will impact our results in areas with highly heterogenous land cover, or small and isolated deforestation events such as targeted logging operations. Recent studies report that logging is the most common driver of loss in intact, but not necessarily protected, forests globally [
28]. Logging is difficult to detect via satellite imagery because in many cases, sufficient canopy cover remains following logging that land cover is still classified as forest. High-spatial resolution and frequent satellite imagery may be able to detect the most evident indications of logging, including access roads, skid trails, and tree fall gaps. However, research suggests that these may comprise as little as 20% of the total area impacted by logging activities [
29]. Because we use a relatively coarse resolution land cover map, very small-scale or ephemeral forest disturbances—even isolated tree cover loss events in the GFC loss map—will be reported as followed by forest cover. Future research with higher resolution imagery could support investigation of the role of logging in PAs globally.
Despite limitations, by calculating tree cover loss at the PA level, we now have a comprehensive global dataset that can be used to compare outcomes across PAs, to identify the specific characteristics of PAs which limit the rate of tree cover loss over time, and to evaluate the impact of PAs on reducing tree cover loss. There is a growing literature aimed at measuring the impact of human pressure on PAs [
3]. This dataset can complement future studies which aim to assess the impacts of socio- and macro-economic factors on ecosystem degradation within PAs. Also, as the land use sector is increasingly recognized for its important role in stabilizing future climate, this research can inform assumptions of land available for agriculture. Given that agriculture is occurring in PAs in some regions, despite their designation, researchers and modelers may not want to assume that all protected land will remain in a natural state to more accurately represent land cover dynamics globally.