Published Date
Author
Abstract
For further details log on website :
http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0138918
- Published: October 14, 2015
- http://dx.doi.org/10.1371/journal.pone.0138918
Author
Abstract
In spite of the high importance of forests, global forest loss has remained alarmingly high during the last decades. Forest loss at a global scale has been unveiled with increasingly finer spatial resolution, but the forest extent and loss in protected areas (PAs) and in large intact forest landscapes (IFLs) have not so far been systematically assessed. Moreover, the impact of protection on preserving the IFLs is not well understood. In this study we conducted a consistent assessment of the global forest loss in PAs and IFLs over the period 2000–2012. We used recently published global remote sensing based spatial forest cover change data, being a uniform and consistent dataset over space and time, together with global datasets on PAs’ and IFLs’ locations. Our analyses revealed that on a global scale 3% of the protected forest, 2.5% of the intact forest, and 1.5% of the protected intact forest were lost during the study period. These forest loss rates are relatively high compared to global total forest loss of 5% for the same time period. The variation in forest losses and in protection effect was large among geographical regions and countries. In some regions the loss in protected forests exceeded 5% (e.g. in Australia and Oceania, and North America) and the relative forest loss was higher inside protected areas than outside those areas (e.g. in Mongolia and parts of Africa, Central Asia, and Europe). At the same time, protection was found to prevent forest loss in several countries (e.g. in South America and Southeast Asia). Globally, high area-weighted forest loss rates of protected and intact forests were associated with high gross domestic product and in the case of protected forests also with high proportions of agricultural land. Our findings reinforce the need for improved understanding of the reasons for the high forest losses in PAs and IFLs and strategies to prevent further losses.
Figures
Citation: Heino M, Kummu M, Makkonen M, Mulligan M, Verburg PH, Jalava M, et al. (2015) Forest Loss in Protected Areas and Intact Forest Landscapes: A Global Analysis. PLoS ONE 10(10): e0138918. doi:10.1371/journal.pone.0138918
Editor: Madhur Anand, University of Guelph, CANADA
Received: March 16, 2015; Accepted: September 4, 2015; Published: October 14, 2015
Copyright: © 2015 Heino et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: The work was funded by Maa- ja vesitekniikan tuki ry. In addition, MK was funded by the Academy of Finland project SCART (grant no. 267463) and MMa by the Academy of Finland project Green Economy and Policies (grant no. 257919). PHV thanks the ERC grant agreement nr. 311819 (GLOLAND) and the FP7 project LUC4C for support. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Forests play a crucial role in sustaining life on earth. They maintain ecological diversity, regulate climate, store carbon, protect soil and water and provide resources and livelihoods for the world’s population [1–4]. Despite the increasing awareness of the importance of these ecosystems, global deforestation rates have remained alarmingly high over the past decades [2].
In year 2010 forest covered around 40 million km2 or 31% of the global land area according to country reports [2]. Estimates of global forest loss rates within the past decade vary between 130,000 km2/yr [2] and 177,000 km2/yr [5]. At the same time, some reforestation and natural regrowth have occurred, leading to net forest loss reports ranging between 52,000 km2/yr and 115,000 km2/yr, respectively. Although forest loss is still remarkably high, it has shown some signs of decline: FAO [2] reports that during the 1990s, the forest loss rate was 160,000 km2/yr; whereas between 2000–2010, the rate was 130,000 km2/yr. According to FAO [2] the tropics were the only domain where the rate of forest loss increased in the first decade of the 21stcentury compared to 1990s: deforestation increased in the tropics by 2,101 km2/yr on average across the decade. Deforestation in the tropics accounted for 32% of global forest loss within the period of 2000–2012 [5]. However, there are also positive signs in the tropics: the rate of forest loss in the Brazilian Amazon has declined in recent years [6,7].
The main drivers of global deforestation are linked to expansion of agriculture, wood extraction, infrastructure extension, population growth, and expansion of agriculture [8–13]. The dominant drivers, however, vary among the regions [8–13]. In addition to agriculture and population growth, a meta-analysis of 117 deforestation studies by Ferretti-Gallon and Busch [9] suggests that deforestation is generally lower in high, steep and wet areas while it is higher in areas where forests are closer to roads and urban areas. Deforestation has also shifted from a dominantly state initiated to an enterprise driven process in the tropics between 1970 and 2000 [14].
Alarming deforestation rates combined with the increasing awareness of the importance of forests have resulted in exponential growth of the world’s protected areas (PA) over the past decades [15–17]. Schmitt et al. [18] estimate the global forest cover using earth observation satellite data from MODIS2005 and the extent of protected forest using the World Database of Protected Areas (WDPA) for the year 2008. They find that 7.7% of global forest cover fell within IUCN’s four strictest protection categories (I-IV) and 13.5% within IUCN’s all six protection categories (I-VI) (IUCN protection categories vary from strictly protected areas (I) to protected areas with sustainable resource use (VI); for more detailed category definition see [19]). Schmitt et al. [18] further conclude that forest protection varied greatly between different regions and forest types, and that forest protection in priority areas, such as biodiversity hotspots, was insufficient.
Considerable efforts of protection are targeted to the primary forests and large intact forest landscapes (IFLs, i.e. unbroken expanse of natural ecosystems). These forest areas play crucial roles in sustaining ecological diversity [2,20]. FAO [2] estimates that primary forests accounted 14 million km2 (36% of the global total forests) in year 2010, having decreased alarming 400,000 km2 over the period of 2000–2010 (annual rate 0.4%). Potapov et al. [20] report that the extent of IFL is 13.1 million km2. Vast majority of the IFLs are found either in dense Tropical and Sub-Tropical forests (45% of world total) or in Boreal Forests (44%). The lowest proportion of IFL was found in Temperate Broadleaf and Mixed Forests. Potapov et al. [20] further find that 18.9% of IFLs are under protection of IUCN protection categories I-VI and only 9.7% of IFLs is strictly protected under IUCN protection categories I-III.
The protection effect of PAs is, however, questioned. The PAs are considered in many cases to be biased in their location, meaning that PAs are located in areas that are unlikely to face land conversion pressures [21]. Joppa and Pfaff [22] reveal that a majority of the PA networks are located in high elevations, steep slopes and far from roads and cities. Joppa and Pfaff [21] argue that this bias has resulted in overestimations in the protection effect of PAs. Local case studies (e.g. [23]) support these global findings. Joppa and Pfaff [24] use ‘matching’ approach that attempts to avoid the overestimation of protection effect of PAs in 147 countries by comparing protected and non-protected areas with similar land characteristics. They find that matching reduced the protection effect in 80% of the countries compared to an assessment without matching. Altogether they find that the protection reduced conversion of natural land cover in 75% of the assessed countries. Protected areas are also reported to become increasingly isolated, especially in tropics [25]. This is alarming because smaller protected areas are often under great threat [26] and isolation of these areas restricts the habitat size, i.e. limits the survival of a great number of fauna and flora species (e.g. [27,28]).Although the understanding of forest loss at a global scale with an increasing spatial resolution is growing rapidly [5], the global forest loss in PAs and in IFLs are not assessed with detailed and uniform datasets that allow consistent forest extent comparisons over space and time. So far only regional analyses exist at this level, for example, for Indonesia [29]. Furthermore, the success of protection in preserving the IFLs is not yet well understood. Therefore, in this study, we aim to conduct a consistent and spatially explicit assessment of the global forest loss, particularly within protected and intact forests by focusing on the variation among the countries between 2000 and 2012. Additionally, we aim to study whether socio-economic indicators can offer potential explanation to these observed global forest losses in PAs and IFLs. We hypothesised that i) the extent of forest loss within PAs and IFLs varies strongly among the countries but is always less manifested within PAs than outside of them, and ii) country level indicators of population size, land use change and state economy can offer potential explanation to the observed global forest losses in PAs and IFLs. The key terminology used in this study is explained in Table 1.
Download:
Table 1. Key terminology used in the study.
Materials and Methods
To conduct the assessment, we combined four global datasets (Fig 1; Table 2): forest extent based on the Global Forest Change (GFC) data [5], forest loss based on the GFC data [5], the World Database of Protected Areas (WDPA) [30] and the global Large Intact Forest Landscapes (IFL) data [20]. From GFC we assessed the forest extent and forest loss while WDPA and IFL datasets were used to determine how much of that loss took place in protected areas and intact forest landscapes respectively. The resulting spatial data were aggregated to country-scale and analysed with Weighted Least-Squares (WLS) regression analyses to assess whether one or a series of socio-economic indicators are associated to forest loss patterns, possibly explaining their occurrence. The forest loss results are given at national and global scale and by geographical regions, while regression analyses were conducted at national scale to correspond with the socio-economic indicator data and the pursued forest governance level. Below the data and methods are described in more detail.
Download:
Fig 1. Aggregated (A-B) and converted (C-D) datasets used for the forest loss analyses.
A) Forest extent in percentage of grid area at 1 km resolution, based on Global Forest Change (GFC) data [5]; B) Forest loss in percentage of grid area at 1 km resolution, based on GCF data [5]; C) Protected areas based on World Database of Protected Areas (WDPA) [30], converted to 1 km resolution; and D) Large Intact Forest Landscapes (IFLs) based on dataset by Potapov et al [20], converted to 1 km resolution.
Download:
Table 2. A description of the datasets used in this study.
GFC stands for Global Forest Change dataset.
Data
We used the GFC dataset [5] to map the forest extent and forest loss (Fig 1A and 1B; Table 2). GFC data are based on Earth observation satellite data and have a resolution of ~30 meters at the equator (i.e. 1 arc-second). The forest extent in the GFC data are given as canopy cover percentage per output grid cell for all vegetation taller than 5 meters for the year 2000 [5]. The forest loss in the GFC dataset is defined as stand-replacing disturbance, or complete removal of tree canopy, which has occurred between years 2000 and 2012 [5]. In this study, we used the same definition for forest loss (Table 1). The GFC dataset provides gross and net forest loss data. We chose to focus on gross forest loss instead of net forest loss, as we believe the former reflects better the occurrence of forest disturbances in protected and intact forests than the latter. The use of net forest loss could have masked some of the disturbances through forest gain.
It should be noted that the GFC dataset does not specify the reasons behind forest loss and thus, the reported forest loss inside protected areas might be due to direct human actions (e.g. illegal logging, management strategies) as well as natural causes (e.g. forest fires, diseases, pests and storms) [2,5,31,32]. For example, natural forest fires are reported to be the most significant cause of forest loss in the boreal forests [33].
We used the WDPA dataset to map the protected areas of the globe (Fig 1C; Table 2). The dataset is a joint project of IUCN and UNEP and it is stated to be the most extensive dataset on protected areas containing both nationally and internationally designated protected areas [30]. The dataset is based on both authority and community sources. We included only the protected areas reported by authority sources into our assessment.
For the intact forest areas we used the global extent of IFLs (Fig 1D; Table 2) [20]. The dataset is created using various remotely sensed datasets and maps, including data from MODIS and Landsat. They define IFL as an unbroken area of natural ecosystems within forests, without signs of significant human activity, and having an area of at least 500 km2. The reference year for this dataset is 2000, i.e. the starting year to our analysis.
Methods: forest extent and forest loss
We conducted the analyses on 1 km resolution (at the equator; i.e. on 30 arc-second resolution). All datasets were either aggregated (forest extent and forest loss) or converted (PA and IFL) to this resolution. The conversion preserved the original resolution of the input data (Table 2) and the aggregation to coarser resolution was carried out in a way that original information was not distorted (see below for more details).
Before any data aggregation, we transformed the original GFC forest cover data into boolean type at 30 m resolution. For the transformation we used 20% tree canopy cover threshold, similarly to Potapov et al. [33], to define whether the 30 m resolution grid cells are classified as forest or no-forest. These boolean forest cover data were multiplied by a surface area of 30 m cell in question. The forest extent area data (in km2) were then aggregated into 1 km resolution by summing the surface area of forest from 30 m grid cells within each 1 km grid cell (note: the area of 30 m grid cells vary over latitudes in the original WGS84 projection and this was taken into account in area calculations). The GFC forest loss data are boolean data in nature (loss or no loss). These forest loss boolean data were multiplied by a surface area of 30 m cell in question and aggregated to 1 km resolution similarly as forest extent (see above). Possible forest loss was only accounted in those 30 m grid cells that had tree canopy cover over the 20% threshold in the forest extent dataset. As a result we got the forest extent and forest loss in km2 for each 1 km grid cell without losing any information in the aggregation process. To illustrate the aggregated datasets, the forest extent and forest loss in percentages of grid cell area are presented in Fig 1A and 1B.
IFL and WDPA datasets (Fig 1C and 1D; Table 2) were used to calculate the extent and loss of protected forest, intact forest and protected intact forest. We first converted the vector data of both datasets to raster data with 1 km resolution, being the original resolution of IFL data while resolution of WDPA data varies between the entries. After this we used the aggregated GFC data to calculate the area of forest extent and forest loss inside each type of spatial criteria within an analysis unit (country, region or global): i) total forest extent and loss, ii) forest extent and loss within protected areas, iii) forest extent and loss within intact forest landscapes, and iv) forest extent and loss within protected intact forest landscapes. We report both absolute and relative values for forest extent and loss.
By using these results, we analysed the forest protection effect of protected areas. The effect was estimated by calculating anomaly ratios for relative loss of protected and non-protected forest, intact and non-intact forest, and for protected intact and non-protected intact forest. The anomaly ratios were calculated at national scale by dividing the relative forest loss in the protected areas in question by the relative forest loss in the corresponding non-protected areas, e.g. [(relative loss of protected intact forest / relative loss of non-protected intact forest)– 1]. Thus the anomaly ratios provide information on the effect of protected areas on forest loss as they indicate whether the relative forest loss is lower (<0) or higher (>0) in the protected forest than in the non-protected forest.
We further identified similarities in forest loss patterns between countries by clustering the country scale results with k-mean clustering method. In this method the clustering of the observations (i.e. country-values) is based on minimizing the sum of the squared Euclidean distances between each observation and the cluster-centroid-value. We assessed goodness of fit in terms of variation within clusters. To define the optimal number of clusters we used the SSE curve. We performed two different cluster analysis, each considering four different parameters. In the first cluster analysis, we identified clusters based on absolute and relative losses in total forest and protected forest areas while in the second cluster analysis, we considered absolute and relative forest losses in intact forest landscapes and protected intact forest landscapes.
It should be noted that while for forest extent and intact forest datasets the reference year is 2000, for protected areas the coverage in year 2010 was taken. The reference year 2010 for protected areas was chosen because harmonization of the reference year was not feasible. The WDPA dataset does not report the establishment year for all protected areas. This however did not compromise the analyses as we found that the WDPA data from the year 2010 adequately reflect the protection status of the forest over the study period: only 8% of the worlds protected areas (where an establishment year is provided) were established after the year 2000.
Methods: the relationship between forest loss and socio-economic indicators
To study the relationship between observed forest loss and socio-economic attributes, we employed WLS (i.e. Weighted Least Squares) regression analysis, which allowed the total extent of forest type in question to be used as weights. The degree to which the socio-economic attributes can account for the relative forest loss was analysed for each of the four forest loss measures considered in this study: total forest loss, loss in protected forest, loss in intact forest, and loss in protected intact forest. Altogether 11 commonly used socio-economic indicators were selected and they are listed in Table 3. The socio-economic indicators were, moreover, selected to represent the widely acknowledged causes for global deforestation: population size, land-use change, state governance and state economy. The selection was carried out with a requirement for available homogeneous data at the global scale. We briefly describe the WLS regression analysis below while step-by-step description is given in S3 Appendix.
Download:
Table 3. The socio-economic indicators used in the regression analysis.
The indicators were used as independent variables in the assessment of the forest loss drivers with Weighted Least-Squares (WLS) regression analysis.
The dependent variables were log10 –transformed prior to analyses to ensure the normality of their distributions. Similarly, the independent variables were transformed with a link function in cases it improved the linearity of the relation with the analysed forest loss measures. The link function of each independent variable was defined separately for each four dependent variables. Despite the transformations, the data included some statistical outliers. Exclusion of these outliers did not affect the significance of the results and thus the reported results were gained without omitting the statistical outliers. The multiple regression models were built by starting with the full model comprising of all 11 independent variables and ending with a model where no multicollinearity was detected. Multicollinearity was determined by the Variation Inflation Factor (VIF, a set limit to VIF < 4) and the collinear independent variable with the lowest explanatory power was always excluded from the model.
Results
Our analyses based on GFC data [5] revealed that forest covered one third of the total land area (approximately 43 million km2) in the year 2000, and 19% of this forest (over 8 million km2) was under some form of protection (Table 4). According to our analyses 25% (11 million km2) of global forests were intact, of which almost 35% was protected. Areas that are both intact and protected covered thus a total of almost 3.7 million km2 (or 8.7% of the total forest extent) (Table 4).
Download:
Table 4. Forest extent in 2000 and its loss over the period of 2000–2012.
The forest extent and loss are calculated for geographical regions in absolute [103 km2] and relative values [%].
The global forest loss based on GFC data was 2.14 million km2 between the years 2000 and 2012, being over 5% of the total forest extent in year 2000 (Table 4). This global forest loss rate corresponds to calculations by Hansen et al. [5]. Further, according to our analyses, over 10% of the total forest loss occurred in protected areas: 219,000 km2 of protected forest (3% of the total protected forest) was lost within that time period. We further found that 269,000 km2(2.5%) of intact and 55,000 km2 (1.5%) of protected intact forest were lost during the years 2000–2012 (Table 4). In the following sections the results are presented in more detail. The numeric results are also given for each country in S2 Table of the supplementary information. It is worth to highlight that our estimates on total forest extent and total forest loss on a global scale presented in this section are given only to validate our analyses (i.e. comparison to Hansen et al. [5]) and to put our results on PAs and IFLs into a context.
Forest extents and their protection status
In year 2000 the absolute forest extent was largest in Russia, China, Indonesia, Democratic Republic of the Congo, Brazil, United States, and Canada, exceeding 1 million km2 (Fig 2A). Together these countries contained over 60% of the global forest extent. When examining relative forest extent (Fig 2B), particularly Scandinavia and Finland (Nordic countries from here on), and countries in central Africa and Southeast Asia had a remarkable portion of their territory under forest. The countries with the largest protected forest areas were mainly the same countries with the greatest extent of forest (Fig 2C). However, when assessing the ratio of protected forest extent and the total forest extent, the largest protected forest areas were located in Central Europe, Australia, and in some countries of South America, Sub-Saharan Africa and the Middle East (Fig 2D), where over 30% of the forests were protected (though these often represent small areas of forest).
Download:
Fig 2. Absolute and relative forest extent by country in year 2000.
Absolute forest extents (left column) are presented in km2 and relative forest extents (right column) are in %. A and B) Total forest extent; C and D) Protected forest extent; E and F) Intact forest extent; G and H) Protected intact forest extent.
The countries with the largest total forest extent had also the largest areas of intact forest (over 400,000 km2), except Indonesia, United States and China (Fig 2E). On the other hand, a large part of European, Middle Eastern, North African and Southern African countries did not have any intact forest landscapes (Fig 1D). When comparing the area of intact forest to the area of total forest, we found that in many countries of North and South America, Southeast Asia, Middle Africa and Australia shares of intact forest of total forest were very high (Fig 2F). In contrast to these countries, Nordic countries had a large proportion of their land under forest (Fig 2B) but only less than 2% of that forest was classified as intact forest (Fig 2F).
Only Brazil had over 1 million km2 of protected intact forest (Fig 2G). Russia, Canada, United States and northern countries of South America had over 100,000 km2 of their intact forests under protection. When examining the relative extent of protected intact forest (extent of protected intact forest versus extent of total intact forest), we found that majority of intact forest, over 70%, were protected in countries such as Thailand, New Zealand, Japan, Madagascar, and Ethiopia as well as Nordic countries (Fig 2H). On the contrary, in some countries less than 20% of the intact forest was under protection, including Russia, Canada and Democratic Republic of the Congo (Fig 2H).
Loss of total and protected forest
The largest absolute forest loss values were encountered in Brazil, Canada, United States, Russia and Indonesia. Each of them experienced forest losses of over 100,000 km2 during the period of 2000–2012 (Fig 3A). The largest relative losses, over 10% of their forest extent, were found in some countries in north-western and southern Africa, and Southeast Asia (Fig 3B). In these countries the relative losses were high in their protected forests too (Fig 3D). Relatively high (> 5%) protected forest losses occurred in Australia, United Stated and some European countries. South Asia, Middle East and Central Asia showed, on the other hand, small numbers of both absolute and relative protected forest losses (Fig 3C and 3D).
Download:
Fig 3. Absolute and relative forest losses by country between 2000 and 2012.
Absolute forest losses (left column) are presented in km2 and relative forest losses (right column) are in %. A and B) Total forest loss; C and D) Loss in protected forests; E and F) Loss in intact forest; G and H) Loss in protected intact forests.
We were able to identify altogether seven distinct groups of countries differing by the nature of the total forest loss and loss in protected forest within the national territories (see Fig 4A). United States, Brazil and Russia formed a cluster with large areas of absolute loss in both total and protected forests while their relative losses were moderate (cluster #5 in Fig 4A; see also Fig 3A–3D). For the largest cluster with 100 countries, both absolute and relative forest losses were very small (cluster #1 in Fig 4A; see also Fig 3A–3D). The characteristics of another large cluster, consisting of 61 countries (cluster #7 in Fig 4A), were rather similar to the previous one, but the absolute and relative losses in both total and protected forests were slightly larger. One distinct cluster was formed also by Canada and Indonesia, together with some smaller islands in Southeast Asia, which all had moderate absolute and relative forest losses (cluster #2 in Fig 4A; see also Fig 3A–3D). The cluster with small absolute forest losses but moderately high relative losses includes countries such as Mongolia, South Africa and Paraguay (cluster #4 in Fig 4A; see also Fig 3A–3D).
Download:
Fig 4. Clustering of countries by their forest loss patterns.
A) country clusters based on absolute and relative losses in total and protected forest; and B) country clusters based on intact and protected intact forest. The cluster characteristics are illustrated by radar charts, of which radii represent the variables used in the clustering. The absolute range in each variable is given next to the variables. The cluster centroids on each variable are plotted on the radii, and these centroids are connected with lines. The countries within each cluster are mapped above the corresponding radar chart and the colours both, in the radar charts and their corresponding maps, indicate the cluster in question.
Loss of intact forest and protected intact forest
The largest losses of intact forest occurred in Canada with losses covering over 100,000 km2(over 5% of the total intact forest), while Russia lost over 50,000 km2 of its intact forests and United States, Brazil and Australia over 5,000 km2 (Fig 3E). In relative terms, Northern America, Australia, and Mongolia lost the largest areas of their intact forest, as over 5% of their intact forests were lost compared to the intact forest extent in year 2000 (Fig 3F).
The largest absolute losses of protected intact forest occurred in Canada, United States and Russia, as all of them lost over 10,000 km2 of protected intact forests (Fig 3G). Large relative losses of protected intact forest (over 5% of total) occurred in Australia, Mongolia, Canada and United States, while smallest relative losses occurred in Asia, Nordic countries and South America (Fig 3H).
When considering similarities and differences between the losses in total intact and protected intact forest, we found that the countries with large areas (Brazil, United States, Canada and Russia) formed each a cluster of their own (see cluster results in Fig 4B). Each of them have very specific pattern in the loss results. For example in the case of Brazil, its absolute intact forest losses were relatively high, but the relative forest losses were rather small when compared to the other countries with the largest forest extents (Fig 3E–3H). Nordic countries, many countries in northern South America, central parts of Africa and South and Southeast Asia fell in the same cluster, having all very small absolute and relative losses in intact forests (Fig 3E–3H). Paraguay and Australia on the other hand, formed a cluster with the rather large relative losses of intact and protected intact forests while their absolute losses were rather small (Fig 3E–3H).
The effect of protected areas
Our findings revealed that the anomaly ratio for loss within protected forest and loss within non-protected forest was largely positive in some African countries, Netherlands, Afghanistan and Mongolia (Fig 5A), indicating that relatively more forest was lost in protected areas than in non-protected areas. The anomaly ratios were around zero in Russia, some African countries, Australia, and many countries in Europe. On the contrary, in a large number of countries in Latin America, Sub-Saharan Africa, Central and Southeast Asia, and Nordic countries, relative forest loss in protected areas was much smaller than in non-protected areas (Fig 5A).
Download:
Fig 5. Comparison of forest loss between different forest types: A) Anomaly ratio for loss in protected forest and loss in non-protected forest; B) Anomaly ratio for loss in intact forest and loss in non-intact forest; C) Anomaly ratio for loss in protected intact forest and loss of non-protected intact forest (anomaly ratio <0 (>0) indicates that forest loss is lower (higher) in the protected forest type than in the non-protected forest type).
Note: all the anomalies were calculated with relative losses.
When examining the anomaly ratios between intact and non-intact areas it can be seen that in Mongolia, Nepal and the Dominican Republic, it was more likely that the lost forest was intact than non-intact (Fig 5B). Australia, Canada, and United States experienced as much relative loss in intact forests as in non-intact forests. In Russia the loss in intact forest was close to the loss in non-intact forest (anomaly ratio is –0.2) while in the rest of the world’s countries with intact forest, the anomaly ratio was below –0.5. In these countries much less intact forest than non-intact forest (relative to total area of forest in both classes) was lost.
In a few countries, including India, Japan, and Chile, the relative losses in intact forest were much larger in protected areas than non-protected areas while in Canada, United States, Colombia, Mongolia, Ethiopia and Russia this anomaly ratio was around zero (Fig 5C) (i.e. relative loss in protected intact forest was equal to relative loss in non-protected intact forest). The smallest anomaly ratios (i.e. relative loss in intact forest is much smaller in protected than non-protected areas) were found in Nordic countries, large parts of Asia, Middle Africa, and Latin America (Fig 5C).
The relationship between forest loss and socio-economic indicators
Each of the forest loss measures (total forest loss, loss in protected forest, loss in intact forest, and loss in protected intact forest), weighted by the forest area in question, had highly significant relationship with several socio-economic indicators (Table 5). All WLS regressions were highly significant and in the case of intact forest, the regressions explained over 75% of the weighted variations among countries in overall forest loss (Table 5).
Download:
Table 5. The associations of common socio-economic parameters to forest loss between 2000–2012.
The four WLS regression models have each a different dependent variable: (a) relative total forest loss, (b) relative loss in protected forest, (c) relative loss in intact forest, and (d) relative loss in protected intact forest. The weight used in the regression for the forest loss measures was the total forest extent in question. Standardised beta coefficient (B) is given for each independent variable included in each of the four models. The independent variables are explained in Table 3. The country specific data are given in S1 Table and S2 Table of the supplementary.
Each of the four forest loss measures represents specific dimensions of forest loss during 2000–2012 and was best explained by a somewhat different set of socio-economic indicators. This variation could not be explained within the scope of the analyses, but it is expected that part of it originates from the different configuration of WLS regression models for different forest loss measures: the identity and number of countries varied between WLS regression models (Table 5). In addition, some of the used dependent variables were cross-correlated leading to a different choice of variables upon small differences between the dependent variables. The lack of socio-economic indicator data for all the countries also reduced the final number of countries included in the WLS regression models (see the NaN values for the socio-economic indicators in the S1 Table of the supplementary information).
Population density growth had, surprisingly, a negative relationship with weighted loss in intact and protected intact forests (Table 5). This indicates that an increase in population density is associated with a smaller forest loss. Population density itself had, similarly to population density growth, a negative relationship with weighted forest loss, although this relationship was only found with the loss in protected and protected intact forests. This indicates a stronger protection effect to prevent forest losses in densely populated countries. We found highly significant positive relationship between Gross Domestic Product (GDP) per capita and weighted loss in protected, intact, and protected intact forests: the higher the GDP, the higher the weighted forest loss. Moreover, polity (i.e. level of democracy) was also found to have a positive relationship with weighted forest loss (i.e. higher developed governance level was linked with higher forest loss) although, the statistical significance was observed in all forest loss types except in protected forest. Human Development Index (HDI) and corruption correlated highly with GDP (Tables A-D of S1 Appendix), and although they can as well offer explanatory power over forest loss, they had to be omitted from the final models.
More straightforwardly, agricultural land area had a significant positive relationship with total weighted forest loss and weighted loss in protected and protected intact forests. Similarly the agricultural land area growth had a positive relationship but only with weighted loss in protected intact forest (Table 5). In S2 Table of the supplementary information, we provide the country level data of our assessment that can be used for further analysis.
Discussion
In our assessment of global forest loss between 2000 and 2012, the loss in protected and protected intact forests was found to be 3% and 1.5%, respectively, on global scale with a high variation among countries and geographical regions. The largest relative losses in protected forests occurred in North America (6%) as well as in Australia and Oceania (5%) and Central America (5%). When assessing the relative losses in intact forests, the largest regional losses were found in North America (6%) and in Eastern Asia (5%), where also protection did not prevent the high losses. Protection of forests and intact forests did not reduce forest loss in approximately one quarter of the analysed countries when the occurring loss within the protected areas was compared to outside of those areas. Looking at the relation between socio-economic indicators and global forest loss in PAs and IFLs, we found metrics of population density, GDP and agricultural land extent to provide probable causes for forest loss, although sometimes in a surprising direction.
Forest extent and its loss
Our estimate for global forest extent based on data from Hansen et al. [5] show relatively consistent values when compared with previous studies (see a detailed comparison in S2 Appendix): our forest extent estimates are very close to those by FAO [2] and Schmitt et al. [18]. Our estimates, nevertheless, for global protected forest extent somewhat exceed the previous findings, as we included protected areas that are not assigned with IUCN protection category, while Schmitt et al. [18] focus on protected areas with IUCN protection category. When the same IUCN protection categories were used in the comparison, our findings would be very close to those by Schmitt et al. [18]. When focusing on the estimates of intact forests, the findings in this study indicate a smaller extent compared to the previous estimates by Potapov et al. [20] and FAO [2], who use primary forest extent as a somewhat equivalent measure for intact forests. When in turn focusing on the protected IFL extent estimates, our results indicate almost double the area compared to the findings by Potapov et al. [20]. But as above, the different definitions for protected areas by the two studies can cause the large difference (for more information, see S2 Appendix).
Looking at the forest loss, the annual average global forest loss estimates of FAO [2] are 34% smaller compared to this study and to the study by Hansen et al. [5], whose data we used. On contrary, while we found that 22,400 km2 intact forest was lost per year (Table 4), FAO [2] reports annual lost of 36,400 km2 of primary forest. Similarly to the comparison of forest extent, a direct comparison of forest loss between the present study and the study by FAO [2] is constrained by the different definitions used in the studies. Despite the problematic comparison between studies, the use of different definitions for assessed parameters related to global forest loss, expand the understanding of the issue.
The findings presented above thus reinforce the current knowledge of the global forest extent and loss concerning especially protected and intact forests. When focusing on these areas, we found that 19% of total forest extent and 34% of intact forest extent were within protected areas in year 2000. By year 2012 the extent of these protected forests had decreased by 3% and 1.5%, respectively (Table 4). These forest losses, furthermore, varied considerably between geographical regions supporting partly our first hypothesis. For example, loss in protected forest exceeded 5% in Australia and Oceania, and North and Central America, while the largest intact forest and protected intact forest losses were found in Eastern Asia and North America between 2000 and 2012 (Table 4; Fig 3B–3D).
Effectiveness of protection
We used our analysis in protected areas to estimate the effectiveness of protection to prevent forest loss in those areas. Supporting further our first hypothesis, we found that in 146 of 191 assessed countries (76%), relative forest loss was smaller in protected areas than outside them (anomaly ratio < –0.125: Fig 5A). The protection, moreover, seemed to be particularly effective for reducing forest loss in South America, Southeast Asia and Sub-Saharan Africa (Fig 5A). Nevertheless, the effect of protection on forest loss was the opposite in some countries in Middle Africa, Europe, Middle East and Mongolia in Asia and thus rejects our first hypothesis of a unified positive protection effect. On country scale, our findings suggest positive news for Brazilian Amazon where protection was found to be very effective in preventing forest loss. This is in agreement with the earlier research findings on the Brazilian Amazon [6,39]. In the case of intact forest, the global results were similar, i.e. protection was effective in 43 of 60 of assessed countries (72%) (Fig 5B). Largest deviations to the global majority were shown by Mongolia, Nepal and Australia. It is worth to acknowledge that the PAs are not necessarily always aimed to protect the forest stands per se, instead protection can have other goals. Prior to our study, the effectiveness of protection is globally analysed by Joppa and Pfaff [24]. Our study further strengthens those findings, yet our findings are retrieved data-wise in more consistent manner by using temporally comparable forest extent datasets (Table 2), which only became available after the study by Joppa and Pfaff [24].
Drivers of forest loss
Importantly, we also explored the drivers of forest loss by looking at the statistical linkage between socio-economic indicators and weighted forest loss (Table 5). Weighing the country specific forest losses by their total forest extent in question allowed our analysis to indicate the drivers that are expressive for the global forest extent instead of focusing on drivers explaining the inter-country variation. This globally uniform assessment showed forest loss having a strong connection with agricultural land extent. Moreover agricultural land expansion is recognised as one of the most important processes causing forest loss (see e.g. [8–13]) and our analysis agrees well with this, particularly concerning the protected intact forests (Table 5). However, it is worth to recognise that agriculture related forest loss involves much more complex dynamics than our regression analyses are able to capture (see e.g. [40]). We further found strong connection between losses in protected and/or intact forests and population density and GDP to fully support our second hypothesis.
While several global scale studies have found that higher population density increases forest loss (e.g. [41,42]), these all report total forest loss. We did not find significant relationship between total forest loss and population density, whereas our findings for losses in intact and protected intact forests (Table 5) show opposite direction compared to these existing studies on total forest loss. Our findings also seem to disagree with a meta-analysis by Porter-Bolland et al. [43]. However, a direct comparison between our results and those of Porter-Bolland et al. [43] is discouraged by the different scopes and approaches in the studies. Unlike our study, the meta-analysis assessed the effectiveness of only 40 protected areas in the tropics. Additionally, population densities are measured by different manners: in our study population densities were country averages while the meta-analysis based its data on population density in the immediate neighbourhood of protected areas. Our desire for single statistical analysis including multiple explanatory variables determined the selection of country averages on the population metric to match with the available level of the other variables. All in all, our findings together with those by Porter-Bolland et al. [43] point to a need for further research on the role of population dynamics in forest loss.
In terms of GDP, we did not find a statistically significant relationship with total forest loss (Table 5). This finding is in line with some previous studies [44–46], whereas some other studies suggest that lower GDP could indicate higher forest loss (e.g. [47,48]). While the effect of GDP on forest loss in protected areas has not been assessed previously on a global scale, it is evaluated by Nagendra [49] on 56 protected areas. Her findings indicate that the rates of land cover clearing do not differ significantly between countries with different levels of GDP per capita. We found, in turn, high GDP per capita to indicate high weighted forest loss rate in protected areas (Table 5). This relation between GDP and weighted forest loss on protected areas can possibly indicate different management strategies in countries with varying GDP, i.e. countries with higher GDP may have more versatile management strategies such as, those that aim to strengthen biodiversity or provide habitats for specific species, compared to countries with lower GDP.
Our assessment, i.e. looking at the forest loss drivers globally, naturally used rather limited set of socio-economic indicators on a coarse scale. This constriction on indicators was needed for the sake of homogeneous data. Another story would be revealed by more local and regional studies that could offer fruitful insights to this global study, as is for example shown by the global meta-analysis looking at the drivers behind the effectiveness of protection management [50]. Our study provides a global empirical assessment, weighted by the total forest extents in question, that can be used as a reference point for more detailed studies at other spatial scales.
Taking further–subnational results and comparison among datasets
We acknowledge that the country level results presented in this study do not reveal finer level heterogeneity. Therefore, as an example, we explored how the forest loss in PAs varies within a country. Taking Colombia as an example, in our results we showed that the mean rate of forest loss within all protected areas was 0.08% per year. In a finer scale assessment using the same dataset (i.e. GFC) we found that the forest loss varied among the individual Colombian PAs from zero to 1.48%, the observed rates being low over large PAs in the Amazon and generally higher for protected areas in the Andes. The mean rate of annual forest loss per PA was 0.12% with standard deviation of 0.19%. This highlights the importance to study further the impact of protection on forest loss in lower spatial scales.
Furthermore, we carried out the same analysis by using an alternative dataset, Terra-I [51], to explore possible differences in available datasets. Terra-i dataset deviates from GFC in its use of MODIS (250 m resolution) instead of Landsat (30 m resolution) and 2004–2014 instead of 2000–2012. The range in forest loss in the individual protected areas (0% to 1.72%) was similar to the finding made by GFC but the mean annual forest loss rate within the PAs was one third (mean rate of 0.043% with standard deviation 0.14%) of the found losses compared to GFC dataset findings. This comparison reflects differences in data type, methodology and period but also indicates the uncertainty associated with these datasets and, thus calls for further comparisons between the widely used datasets.
Methodological limitations and ways forward
It is important to expand the discussion towards the limitations created by the used datasets that also our study is naturally subjected to. For example, for the forest cover and loss we used the GFC dataset [5], which can give deviating results compared to other datasets mainly due to the differences in methods such as: definition of forest, separation of anthropogenic from naturally occurring forest loss, and handling of atmospheric distortions (e.g. cloud cover) in mapping the forest cover. Although, diving deeper into these differences and the effects they have on the findings is beyond the scope of this analysis, it would offer a fruitful topic for a comparative study. In addition, Lee [52] raises a concern on the IFL dataset [20] that its definition of intactness of forest may not capture the socio-ecological variation between different geographical regions and may thus introduce inaccuracies in country specific analyses. Inter-comparison of global scale datasets as well as more detailed analyses at the country/regional level would be needed to unravel the state of global classification of IFLs.
It should be noted that our study focused on forest loss, or more precisely on stand-replacing disturbance, ignoring the possible forest gain. This focus gave us more space to explore the loss of original habitats within protected areas and intact forest landscapes and the drivers behind these losses. A comparison of our findings to findings that would be based on the net forest loss, i.e. a loss assessment including the planting of forests and natural forest regeneration (see [2,5] and references therein for further information), would allow tracking of the forest canopy cover per se.
Our estimations on the effectiveness of protection to reduce forest loss within PAs compared to the forest loss outside those areas provide useful information the protection effect at the country level. At the country level, we were nevertheless not able to assess the leakage effect of individual protected areas, as done by e.g. Oliveira et al [53]. As highlighted by Ewers and Rodrigues [54] a comparison of forest loss in the vicinity of and inside protected areas is essential in understanding the effects of forest protection on a local scale. It is also worthwhile to recognise that forested areas might have been protected for other reasons than protect the forest stands per se. The global datasets used in this study could be further employed to this kind of assessment, comparing the leakage effect within and among regions and furthermore extending it all the way to global level. Most importantly, global data on forest management goals per each protection area would be highly needed to support more detailed assessment on the impact of protection on forest loss.
Conclusions
Our global scale analysis indicates that substantial areas of protected (3%) and intact forests (2.5%) were lost over the past decade, and large part of that loss occurred in countries with well-established protection (e.g. in Australia and Oceania, and North America). In terms of comparison between forest loss inside and outside of PAs and IFLs we could not confirm our first hypothesis, which stated that forest loss would always be less manifested within PAs than outside of them. Although the relative forest loss was smaller inside the PAs and protected IFLs in global averages and in majority of countries, there were numerous countries that showed an opposite effect. Also our second hypothesis was not entirely confirmed: while population size and other indicators were able to explain differences in weighted global forest loss within PAs and IFLs, a large part of the variation remained unexplained. Moreover, the statistical associations did not always correspond with our causal understanding of forest loss processes. Our findings thus highlight that protection of forested areas does not always guarantee a lower rate of forest loss. There is, indeed, a high geographical variation in the effectiveness of protection against forest loss, a variation that is probably at least to some extent explained by countries’ different means and intensities to combat against forest loss.
Supporting Information
S1 Appendix. Cross-correlation tables.
Cross-correlation tables of socio-economic indices in cases of different forest loss categories.
doi:10.1371/journal.pone.0138918.s001
(PDF)
S2 Appendix. Comparison of our findings to previous studies: forest extent and loss.
doi:10.1371/journal.pone.0138918.s002
(PDF)
S3 Appendix. Weighted Least Squares (WLS) regression analysis steps.
doi:10.1371/journal.pone.0138918.s003
(PDF)
S1 Table. Socio-economic data at country level.
Socio-economic data used in the WLS regression analysis.
doi:10.1371/journal.pone.0138918.s004
(XLSX)
S2 Table. Country level results.
Forest extent and forest loss results presented at country level.
doi:10.1371/journal.pone.0138918.s005
(XLSX)
Acknowledgments
We are grateful for CSC—IT Center for Science for providing the computational facilities for the analyses. We also thank the reviewers for their constructive comments and Prof. Olli Varis, Dr. Joseph Guillaume, and other members of our research team at Aalto University for their encouragements and fruitful discussions.
Author Contributions
Conceived and designed the experiments: MH MK M. Mulligan TAR. Performed the experiments: MH MK M. Mulligan MJ TAR. Analyzed the data: MH MK M. Makkonen M. Mulligan MJ TAR. Contributed reagents/materials/analysis tools: MH MK M. Makkonen M. Mulligan MJ TAR. Wrote the paper: MH MK M. Makkonen M. Mulligan PHV TAR.
References
- 1.Bonan GB (2008) Forests and Climate Change: Forcings, Feedbacks, and the Climate Benefits of Forests. Science 320: 1444–1449 doi: 10.1126/science.1155121. pmid:18556546
- 2.FAO (2010) Global forest resource assessment 2010. Rome: Food and Agrculture Organisation of the United Nations.
- 3.Gibson L, Lee TM, Koh LP (2011) Primary forests are irreplaceable for sustaining tropical biodiversity. Nature 478: 378–381. doi: 10.1038/nature10425. pmid:21918513
- 4.Pan Y, Birdsey RA, Fang J, Houghton R, Kauppi PE, Kurz WA, et al. (2011) A Large and Persistent Carbon Sink in the World’s Forests. Science 333: 988–993 doi: 10.1126/science.1201609. pmid:21764754
- 5.Hansen MC, Potapov PV, Moore R, Hancher M, Turubanova SA, Tyukavina A, et al. (2013) High-resolution global maps of the 21st-century forest cover change. Science 342: 850–853. doi: 10.1126/science.1244693. pmid:24233722
- 6.Nepstad D, Soares-Filho BS, Merry F, Lima A, Moutinho P, Carter J, et al. (2009) The End of Deforestation in the Brazilian Amazon. Science 326: 1350–1351. doi: 10.1126/science.1182108. pmid:19965742
- 7.Godar J, Gardner TA, Tizado EJ, Pacheco P (2014) Actor-specific contributions to the deforestation slowdown in the Brazilian Amazon. Proceedings of the National Academy of Sciences 111: 15591–15596.
- 8.DeFries RS, Rudel T, Uriarte M, Hansen M (2010) Deforestation driven by urban population growth and agricultural trade in the twenty-first century. Nature Geoscience 3: 178–181.
- 9.Ferretti-Gallon K, Busch J (2014) What drives deforestation and what stops it? A meta-analysis of spatially explicit econometric studies. Washington, DC: Center for Global Development.
- 10.Geist HJ, Lambin EF (2002) Proximate causes and underlying driving forces of tropical deforestation. BioScience 52: 143–150.
- 11.Hosonuma N, Herold M, De Sy V, De Fries RS, Brockhaus M, Verchot L, et al. (2012) An assessment of deforestation and forest degradation drivers in developing countries. Environmental Research Letters 7.
- 12.Rudel TK (2005) Tropical Forests. Regional Paths of Destruction and Regeneration in the Late Twentieth Century. New York: Columbia University Press.
- 13.Rudel TK, DeFries R, Asner GP, Laurance WF (2009) Changing drivers of deforestation and new opportunities for conservation. Conservation Biology 23: 1396–1405. doi: 10.1111/j.1523-1739.2009.01332.x. pmid:20078640
- 14.Rudel TK (2007) Changing agents of deforestation: From state-initiated to enterprise driven processes, 1970–2000. Land Use Policy 24: 35–41.
- 15.Juffe-Bignoli D, Burgess ND, Bingham H, Belle EMS, de Lima MG, Deguignet M, et al. (2014) Protected Planet Report 2014. Cambridge, UK: UNEP-WCMC.
- 16.Naughton-Treves L, Holland MB, Brandon K (2005) The role of protected areas in conserving biodiversity and sustaining local livelihoods. Annu Rev Environ Resour 30: 219–252.
- 17.Radeloff VC, Beaudry F, Brooks TM, Butsic V, Dubinin M, Kuemmerle T, et al. (2013) Hot moments for biodiversity conservation. Conservation Letters 6: 58–65.
- 18.Schmitt CB, Burgess ND, Coad L, Belokurov A, Besançon C, Boisrobert L, et al. (2009) Global analysis of the protection status of the world’s forests. Biological Conservation 142: 2122–2130.
- 19.Dudley N (2008) Guidelines for Applying Protected Area Management Categories. Gland, Switzerland: IUCN.
- 20.Potapov P, Yaroshenko A, Turubanova S, Dubinin M, Laestadius L, Thies C, et al. (2008) Mapping the world's intact forest landscapes by remote sensing. Ecology and Society 13.
- 21.Joppa LN, Pfaff A (2010) Reassessing the forest impacts of protection: The challenge of nonrandom location and a corrective method. Annals of the New York Academy of Sciences 1185: 135–149. doi: 10.1111/j.1749-6632.2009.05162.x. pmid:20146766
- 22.Joppa LN, Pfaff A (2009) High and far: Biases in the location of protected areas. PloS ONE 4: e8253–e8273.
- 23.Verburg PH, Overmars KP, Huigen MGA, de Groot WT, Veldkamp A (2006) Analysis of the effects of land use change on protected areas in the Philippines. Applied Geography 26: 153–173.
- 24.Joppa LN, Pfaff A (2011) Global protected area impacts. Proceedings of the Royal Society B 278: 1633–1638. doi: 10.1098/rspb.2010.1713. pmid:21084351
- 25.DeFries R, Hansen A, Newton AC, Hansen MC (2005) Increasing isolation of protected areas in tropical forests over the past twenty years. Ecological Applications 15: 19–26.
- 26.Maiorano L, Falcucci A, Boitani L (2008) Size-dependent resistance of protected areas to land-use change. Proc R Soc B June 7, 2008 275 1640 1297–1304 275: 1297–1304.
- 27.Brook BW, Sodhi NS, Bradshaw CJA (2008) Synergies among extinction drivers under global change. Trends in Ecology & Evolution 23: 453–460.
- 28.Turner IM (1996) Species loss in fragments of tropical rain forest: A review of the evidence. Journal of Applied Ecology 33: 200–209.
- 29.Margono BA, Potapov PV, Turubanova S, Stollle F, Hansen MC (2014) Primary forest cover loss in Indonesia over 2000–2012.
- 30.IUCN, UNEP-WCMC (2014) The World Database on Protected Areas (WDPA). In: UNEP-WCMC, editor. Cambridge, UK. Available: www.protectedplanet.net. Accessed 7 July 2014.
- 31.Bonda WJ, Keeleyb JE (2005) Fire as a global ‘herbivore’: the ecology and evolution of flammable ecosystems. Trends in Ecology and Evolution 20: 387–394. pmid:16701401
- 32.Logan JA, Régnière J, Powell JA (2003) Assessing the impacts of global warming on forest pest dynamics. Frontiers in Ecology and the Environment 1: 130–137.
- 33.Potapov P, Hansen MC, Stehman SV, Loveland TR, Pittman K (2008) Combining MODIS and Landsat imagery to estimate and map boreal forest cover loss. Remote Sensing of Environment 112: 3708–3719.
- 34.UN (2014) World Population Prospects: The 2012 Revision. UN Population Division. Available: http://esa.un.org/unpd/wpp/Excel-Data/population.htm. Accessed 7 August 2014.
- 35.World Bank (2014) World Bank Development Indicators. World Bank. Available: http://data.worldbank.org/indicator. Accessed 7 August 2014.
- 36.Transparency International (2010) Corruption Perception Index. Transparency International. Available: http://www.transparency.org/cpi2010/. Accessed 14 June 2014.
- 37.Marshall MG, Gurr T (2014) Polity IV Project, Political Regime Characteristics and Transitions, 1800–2013. Center for Systematic Peace. Available: http://www.systemicpeace.org/inscrdata.html. Accessed 14 June 2014.
- 38.UNDP (2014) Human Development Index 2012. United Nations Development Programme. Available: https://data.undp.org/dataset/Table-1-Human-Development-Index-and-its-components/wxub-qc5k. Accessed 7 August 2014.
- 39.Nolte C, Agrawal A, Silvius KM, Soares-Filho BS (2013) Governance regime and location influence avoided deforestation success of protected areas in the Brazilian Amazon. PNAS 110: 4956–4961. doi: 10.1073/pnas.1214786110. pmid:23479648
- 40.Pagnutti C, Bauch CT, Anand M (2013) Outlook on a Worldwide Forest Transition. PLoS ONE 8: e75890. doi: 10.1371/journal.pone.0075890. pmid:24130750
- 41.Mather AS, Needle CL (2000) The relationships of population and forest trends. Geographical Journal 166: 2–13.
- 42.Mather AS, Needle CL, Fairbairn J (1998) The human drivers of global land cover change: the case of forests. Hydrological Processes 12: 1983–1994.
- 43.Porter-Bolland L, Ellis EA, Guariguata MR, Ruiz-Mallen I, Negrete-Yankelevich S, Reyes-Garcia (2012) Community managed forests and forest protected areas: An assessment of their conservation effectiveness across tropics. Forest Ecology and Management 268: 6–17.
- 44.Damette O, Delacote P (2012) On the economic factors of deforestation: What can we learn from quantile analysis? Economic Modelling 29: 2427–2434.
- 45.Koop G, Tole L (1999) Is there an environmental Kuznets curve for deforestation? Journal of Development Economics 58: 231–244.
- 46.Tole L (1998) Sources of Deforestation in Tropical Developing Countries. Environmental Management 22: 19–33. pmid:9419282
- 47.Ewers RM (2006) Interaction effects between economic development and forest cover determine deforestation rates. Global Environmental Change 16: 161–169.
- 48.Smith RJ, Muir RDJ, Walpole MJ, Balmford A, Leader-Williams N (2003) Governance and the loss of biodiversity. Nature 426: 67–70. pmid:14603318
- 49.Nagendra H (2008) Do parks work? Impact of protected areas on land cover clearing. Ambio 37: 330–337. pmid:18828277
- 50.Leverington F, Lemos Costa KL, Pavese P, Lisle A, Hockings M (2010) A Global Analysis of Protected Area Management Effectiveness. Environmental Management 46: 685–698. doi: 10.1007/s00267-010-9564-5. pmid:20859627
- 51.Reymondin L, Jarvis A, Perez-Uribe A, Touval J, Argote K, Rebetez J, et al. (2012) A methodology for near real-time monitoring of habitat change at continental scales using MODIS-NDVI and TRMM. Available: www.terra-i.org. Accessed 16 December 2014
- 52.Lee PG (2009) Caution against using intact forest-landscapes data at regional scales. Ecology and Society 14: r1. Available: http://www.ecologyandsociety.org/vol14/iss11/resp11/.
- 53.Oliveira PJC, Asner GP, Knapp DE, Almeyda A, Galván-Gildemeister R, Keene S, et al. (2007) Land-Use Allocation Protects the Peruvian Amazon. Science 317: 1233–1236. pmid:17690260
- 54.Ewers RM, Rodrigues ASL (2008) Estimates of reserve effectiveness are confounded by leakage. Trends in Ecology and Evolution 23: 113–116. doi: 10.1016/j.tree.2007.11.008. pmid:18280005
For further details log on website :
http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0138918
No comments:
Post a Comment