Abstract
Conflicts between forest conservation and socio-economic development in the Brazilian Legal Amazon (BLA) have persisted for years but the effects of Indigenous territory (ITs) and protected area (PAs) status on deforestation there remain unclear. To address this issue, we analysed time-series satellite images and qualified annual forest area in the BLA under different governance and management regimes. Between 2000 and 2021, areas classified as ITs or PAs increased to cover 52% of forested areas in the BLA while accounting for only 5% of net forest loss and 12% of gross forest loss. In the years (2003â2021) after establishment, gross forest loss fell 48% in PAs subject to âstrict protectionâ and 11% in PAs subject to âsustainable useâ. However, from 2018 to 2021 the percentage rate of annual gross forest loss in ITs/PAs was twice that of non-designated areas. Our findings reveal the vital role of, and substantial progress achieved by, ITs and PAs in Amazonian forest conservation as well as the dangers of recent weakening of Brazilâs forest policies.
Similar content being viewed by others
Main
The Brazilian Legal Amazon (BLA) has the largest tropical rainforest area, the highest biodiversity and the largest amount of aboveground biomass in the world1,2,3. Since 2000, the Indigenous territories (ITs) and the protected areas (PAs) in the BLA have increased substantially and by 2013 they accounted for 43% of the total land area and covered about half of the total forest area in the region1. ITs and PAs have important roles in forest and biodiversity conservation and climate mitigation in the BLA4,5. Many studies have investigated the effects of the ITs/PAs in reducing deforestation and strengthening forest conservation in the BLA1,6,7,8,9,10, but their results differ substantially as different approaches and datasets were used. A comparison of deforestation rates between the PAs and the non-protected areas (non-PAs) showed that the deforestation rates in the PAs were 1.6â2.2 times6 or ten times1,7,8 lower than those in the non-PAs, respectively. PAs are usually located in remote areas with small deforestation pressure (Extended Data Fig. 1), whereas the non-PAs often have high deforestation pressure. The strict-protected PAs, located in more remote areas, have more forest cover and less deforestation pressure in comparison to sustainable-use PAs and ITs (Extended Data Fig. 2). The comparison of deforestation rate between PAs and non-PAs may not fully show the effects of the PAs in reducing deforestation. Spatial matching methods were applied to assess the effects of PAs on deforestation in the BLA9,10; however, due to extensive deforestation and forest degradation11, it is challenging to find locations in the non-PAs with conditions similar to those in each PA to serve as the controls for the spatial matching methods. Various studies have assessed the effects of ITs/PAs using the official Brazilian deforestation dataset (PRODES)6,7,8, which reports deforestation of primary forests based on the reference forest map in the 1980s and does not include secondary forests that are important for the carbon cycle and biodiversity conservation12.
The conflicts between forest conservation and socio-economic development in the BLA persist but vary over the years13,14,15,16. Forest conservation in ITs/PAs has encountered increasing threats from loosened environmental laws and regulations, changing governmental policies and massive economic development17,18, especially after 201219. About 100âÃâ106âha, including at least 20% of the area in strict protection PAs and ITs, have applications pending for mineral prospecting or mining operations20,21. The satellite data show that illegal mining in the BLA substantially increased from 2000 to 2012 (an average area of 2,500âhaâyrâ1) and hit a record high area in 2020 (10,000âhaâyrâ1) amid widespread protests from Indigenous people22. Primary forest loss reached a decade high in 202023. Between March and September 2020, Brazil passed 27 legislative acts that weakened environmental protection24. Fines for violation of environmental and conservation laws dropped by 72% from March to August 2020, despite an increase in deforestation24. Besides, COVID-19 started to spread widely in the Amazon in early 2020, causing high mortality rates among Indigenous groups25, creating favourable conditions for encroachment of Indigenous lands by illegal loggers and miners.
PAs in the BLA have two types of governance and institution (national versus state PAs) and two types of management objective (strict protection versus sustainable use; Extended Data Figs. 1 and 2). To what degree were deforestation dynamics in the BLA and ITs/PAs under different governance and management affected by the COVID-19 pandemic and changing policies in recent years? Here, we investigate and report on (1) interannual changes in forest areas and deforestation in the BLA, ITs, PAs and non-PAs (that is, non-designated areas) and (2) the effects of ITs, national PAs (nPAs) and state PAs (sPAs) in reducing deforestation from 2000 to 2021 under the changing environmental laws and governance, natural disturbances and COVID-19. We generated annual evergreen forest maps by using daily images acquired by moderate resolution imaging spectroradiometer (MODIS) onboard the Terra satellite from 2000 to 2021 and the forest mapping algorithm1,2.
Results
Increased areas of ITs/PAs for forest conservation
The number and area of ITs/PAs in the BLA increased from 1980 to 2018 with noticeable characteristics (Fig. 1a). The ITs area rose slowly from 0.2âÃâ106âha in 1980 to 4âÃâ106âha in 1988 but started to increase rapidly after 1988 and reached 74âÃâ106âha in 2000 and 115âÃâ106âha in 2016. This large and rapid increase in numbers and areas of ITs was driven by the 1988 Constitutionâs requirement that the government demarcate all Indigenous lands within 5âyears. Although the Constitutionâs requirement has yet to be fulfilled, it stimulated substantial expansion of ITs, including 40âÃâ106âha demarcated with financial support from the G7 Pilot Program to Conserve the Brazilian Rainforest. In 2000, the Brazilâs National Protected Areas System was officially established to expand PAs and better manage the PAs for forest conservation. The nPAs area increased from 9âÃâ106âha in 1980 to 27âÃâ106âha in 2000, 62âÃâ106âha in 2008 and 66âÃâ106âha in 2018. The sPAs area increased from less than 1âÃâ106âha in 1980â1989 to 25âÃâ106âha in 2000, 56âÃâ106âha in 2007 and 60âÃâ106âha in 2018. When lumped together, the total area of ITs/PAs substantially expanded from 10âÃâ106âha in 1980 to 126âÃâ106âha in 2000 and 241âÃâ106âha in 2018. The numbers of ITs, nPAs and sPAs also had similar changes over time and substantially increased from 2, 11 and 2 in 1980 to 227, 70 and 95 in 2000 and 387, 146 and 191 in 2018, respectively.
As most ITs/PAs are covered by forests, the large increase of ITs/PAs strengthens forest conservation substantially in the BLA. Figure 1b shows the interannual change in forest area in ITs, nPAs and sPAs during 2000â2018. Forest area in ITs increased from 67âÃâ106âha in 2000 to 105âÃâ106âha in 2016, an increase of 55%. Forest area in nPAs increased from 25âÃâ106âha in 2000 to 62âÃâ106âha in 2018, an increase of 146%. Forest area in sPAs increased from 18âÃâ106âha in 2000 to 50âÃâ106âha in 2018, an increase of 169%. We overlaid the 2000 annual forest map (a total of 394âÃâ106âha of forest) with the 2018 boundary maps of ITs/PAs and these ITs/PAs in 2018 covered 206âÃâ106âha forest, accounting for 52% of the total forest area in 2000, which clearly shows the importance of the ITs/PAs for forest conservation.
Interannual change of forest area during 2000â2021
We used our annual forest maps to quantify the interannual changes in forest area during 2000â2021 (Fig. 2a). The total forest area in the BLA decreased substantially from 394âÃâ106âha in 2000 to 366âÃâ106âha in 2021, a loss of 28âÃâ106âha (~7% of the forest area in 2000 or an area loss larger than Brazilâs state of Rondônia). The 2018 ITs/PAs boundary maps were overlaid on the annual forest maps to calculate the interannual change of forest area in ITs/PAs during 2000â2021. Total forest area in ITs, nPAs and sPAs decreased from 105.7âÃâ106âha, 62.6âÃâ106âha and 50.5âÃâ106âha in 2000 to 104.9âÃâ106âha, 62.2âÃâ106âha and 50.0âÃâ106âha in 2021, respectively, with the average loss rates of 0.04%âyrâ1, 0.03%âyrâ1 and 0.10%âyrâ1 over their forest areas in 2000 (Fig. 2bâd). As some of the ITs/PAs overlap each other, we combined the ITs/PAs together; the total forest area in the combined ITs/PAs decreased slightly from 206.4âÃâ106âha in 2000 to 204.9âÃâ106âha in 2021, an average annual loss rate of 0.08âÃâ106âhaâyrâ1 (0.04%âyrâ1). To our surprise, ITs/PAs together had a slightly larger forest area in 2018â2021 (204.6âÃâ106âha) than in 2014â2017 (204.1âÃâ106âha), which may be related to the recovery of forests after severe damage in 2015/2016 and, to a much lesser extent, from tree-planting projects. In comparison, the forest area in the non-PAs (Fig. 2e) decreased from 187.8âÃâ106âha in 2000 to 161.3âÃâ106âha in 2021, with an average annual loss rate of 1.3âÃâ106âhaâyrâ1 (0.7%âyrâ1), about 14 times more than the ITs/PAs. The loss of 27âÃâ106âha forest area during 2000â2021 in the non-PAs accounted for ~95% of the total loss of forest area in the BLA over the same period.
Geographically, the interannual change of forest area in ITs, nPAs and sPAs from 2000 to 2021 had noticeable spatiotemporal patterns. The trend analysis of forest areas (Fig. 2fâi) showed divergent dynamics among ITs/PAs. ITs/PAs in the southern and eastern portions of the BLA (the âarc of deforestationâ) had the largest losses of forest area. The ITs/PAs in the northern portion of the BLA either had no significant change in forest area or had increased forest area. In total, 39.8% of nPAs (59 out of 146 nPAs), 40.4% of ITs (154 out of 387 ITs) and 44.0% of sPAs (84 out of 191 sPAs) had significant forest area loss from 2000 to 2021 (Pâ<â0.1).
Spatiotemporal dynamics of primary forest area loss
We used the 2001 forest map as the reference map and identified the first year the forest pixels in 2001 were classified as non-forest pixels during 2002â2021 (Fig. 3a,b) and we counted the number of pixels with a change from forest to non-forest (forest loss) in a year as annual gross forest area loss over the BLA, ITs, PAs and non-PAs from 2002 to 2021 (Fig. 3câg). The cumulative gross forest area losses during 2002â2021 were 49âÃâ106âha for the BLA, including 2.1âÃâ106âha for ITs, 1.2âÃâ106âha for nPAs, 2.8âÃâ106âha for sPAs and 43.1âÃâ106âha for non-PAs. The combined ITs/PAs had a 5.9âÃâ106âha gross forest area loss from 2002 to 2021, accounting for ~12% of total gross forest area loss in the BLA, which clearly indicates the critical role of the ITs/PAs in forest conservation.
The interannual change of gross forest area loss in the BLA (Fig. 3c) during 2002â2021 reveals three interesting results. First, annual gross forest area losses in 2005, 2007, 2010 and 2015 were substantially larger than those in previous and subsequent years. Years 2005, 2007, 2010 and 2015 were characterized by strong El Niño events (Extended Data Fig. 3), high air temperature or severe drought26. Second, annual gross forest area loss in 2013 was least during 2000â2021. Year 2013 was a year with high air temperature without an El Niño or Atlantic Multidecadal Oscillation event. Third, when those years (2005, 2007, 2010, 2013 and 2015) are taken out, the temporal dynamics of annual gross forest area loss show five noticeable phases and each of these phases lasted a few years: (1) increased forest area loss in 2002â2004, (2) reduced forest loss in 2006â2009, (3) increased forest area loss in 2011â2014, (4) decreased forest area loss in 2016â2018 and (5) increased forest area loss in 2019â2021. In 2018â2021, the annual gross forest area loss rates increased 3.6 times in ITs/PAs, larger than the increase in non-PAs (1.6 times), indicating increasing deforestation pressure and an alarming signal in ITs/PAs.
Varying effects of ITs/PAs on annual forest area loss
To investigate the effect of ITs/PAs on reducing forest area loss after they were established, we selected those ITs/PAs that were established after 2002 and analysed the average annual gross forest area loss rates before and after their establishment years (Fig. 4). The results show that 72 ITs, 32 nPAs and 38 sPAs had substantially reduced annual forest area loss rates but 49 ITs, 21 nPAs and 27 sPAs still had small to moderate increases in annual forest area loss rates; the other 24 ITs, 2 nPAs and 7 sPAs had no change in forest area loss rates as they had little or no forest area loss (Fig. 4aâc). There was no clear geographical cluster among those ITs/PAs with reduced or increased forest loss rates (Fig. 4dâf). Most of the sPAs in the northern and western BLA had small increases in annual forest area loss rates. In terms of institution and governance, the average annual gross forest area loss rates were reduced substantially for the nPAs (36%, 13.7âÃâ103âha) and the ITs (30%, 10.7âÃâ103âha) but only slightly for the sPAs (5%, 2.7âÃâ103âha; Fig. 4g and Extended Data Fig. 4a). In terms of management objectives, the average annual gross forest area loss rates were reduced substantially for the strict-protection PAs (48%, 8.6âÃâ103âha) and the ITs (30%, 10.7âÃâ103âha) but only moderately for the sustainable-use PAs (11%, 7.7âÃâ103âha; Fig. 4h and Extended Data Fig. 4b).
Different deforestation dynamics inside different PA types
Forest area in nPAs with the strict-protection objective (Fig. 5a) varied slightly during 2000â2013 but decreased moderately in 2014â2016, with a net loss of 0.14âÃâ106âha (0.45%) from 2000 to 2021. Forest area in the nPAs with the sustainable-use objective (Fig. 5a) decreased continuously from 2000 (31.7âÃâ106âha) to 2021 (31.5âÃâ106âha), a net loss of 0.22âÃâ106âha (0.7%). Forest area in sPAs with the strict-protection objective (Fig. 5b) decreased from 5.91âÃâ106âha in 2000 to 5.88âÃâ106âha in 2021, a net loss of 0.03âÃâ106âha (0.51%). Forest area in sPAs with the sustainable-use objective (Fig. 5b) decreased from 44.65âÃâ106âha in 2000 to 43.53âÃâ106âha in 2016, a net loss of 1.12âÃâ106âha (2.5%) but it had a modest recovery by 2021 (44.19âÃâ106âha), thus a net loss of 0.46âÃâ106âha (1.0%) from 2000 to 2021.
Differences in annual gross forest area loss rates were small and not significant between the PAs with different governance and management (Fig. 5c,d). The sPAs with the sustainable-use objective had the highest average gross forest area loss rate (0.31â±â0.14%âyrâ1 of the forest area in 2001), followed by nPAs with the strict-protection objective (0.11â±â0.05%âyrâ1) and nPAs with the sustainable-use objective (0.09â±â0.04%âyrâ1). The sPAs with the strict-protection objective had the lowest average gross forest area loss rate (0.04â±â0.04%âyrâ1) due to these PAs being far away from the âarc of deforestationâ and therefore were under little deforestation pressure (Fig. 4f). The gross forest area loss rates increased from 2018 to 2021; these losses were probably related to the loosened forest conservation policies during the Bolsonaro presidential administration that began in January 201917,18,27. Gross forest area loss from 2018 to 2021 increased 1.5 times in nPAs with the strict-protection objective and 5 times in nPAs with the sustainable-use objective. Gross forest area loss from 2018 to 2021 increased 12.4 times in the sPAs with the strict-protection objective and 4.3 times in sPAs with the sustainable-use objective.
Discussion
Interannual change of forest area and deforestation
Annual forest maps from our mapping tools, PRODES28 and Global Forest Watch (GFW)29 were evaluated for South America30 and the BLA1, showing higher accuracy in our annual forest maps. The interannual changes of forest area and deforestation during 2000â2017 from these datasets were reported1. Here, we extended the data record from 2017 to 2021 (Fig. 6). We compared our forest area with the newly developed MapBiomas31 forest area and they had similar interannual change trends in forest area in the BLA (Extended Data Figs. 5 and 6). During 2000â2020, the average annual forest area from our dataset is ~4% lower than the MapBiomas and our forest area declined by 8% (32.7âÃâ106âha), close to the MapBiomas (7%, 27.9âÃâ106âha). As the MapBiomas dataset has not provided annual gross forest area loss data to the public yet, it was not used for detailed analysis in this study. Annual gross forest area loss from our dataset was 1.0âÃâ106âha in 2018 and increased to 2.9âÃâ106âha in 2021 (Fig. 6a). Annual gross forest area loss from the PRODES deforestation dataset was 0.8âÃâ106âha in 2018 and increased to 1.3âÃâ106âha in 2021 (Fig. 6b). Compared to the PRODES dataset, annual gross forest area loss from our dataset was higher but had a similar temporal trend. Our previous study explained the differences in forest area loss estimates between the PRODES dataset and our forest map, using annual forest maps in the BLA from 2002 to 20161. Annual gross forest area losses from the GFW dataset in 2018â2020 were slightly higher than those from our dataset (Fig. 6c). Annual forest area loss from the GFW dataset has been modified in numerous ways and these forest area loss data may not be comparable in those years32. Despite their differences in forest definition, in the satellite image data used and, in the forest mapping algorithms, all three datasets report that annual forest area loss in the BLA increased from 2018 to 2021. Together these results provide strong satellite-based evidence of increased forest area loss in 2020/2021.
Effects of ITs/PAs on forest conservation
Several studies analysed the PRODES deforestation datasets and assessed the effects of the ITs/PAs on forest conservation under changing laws, policies and climate change9,33,34. One study found that during 2000â2008 ITs/PAs in the states outside of the âarc of deforestationâ had little impact on deforestation within their boundaries but ITs and nPAs in the âarc of deforestationâ were more effective in reducing deforestation than were sPAs9. Another study reported that 91 sPAs established between 2005 and 2016 reduced deforestation both within their boundaries and in their adjacent surroundings during 2005â201733. The third study also reported that ITs/PAs reduced deforestation during 2000â2010 in the BLA34. The PAs with the strict-protection objective reduced deforestation more than PAs with the sustainable-use objective and the ITs were particularly effective at reducing deforestation in locations with high deforestation pressure34.
Our study uses a longer (2000â2021) and updated dataset to assess the effects of the ITs/PAs for forest conservation in the BLA. Our results over the 2000â2013 period agree with the findings from these previous publications and showed that ITs/PAs reduced deforestation within their boundaries. PAs with the strict-protection objective had small forest area losses in 2000â2013 but large forest area losses in 2013â2021 (Fig. 5). ITs/PAs in the states outside of the âarc of deforestationâ reduced deforestation within their borders (Fig. 4), which differs from that of ref. 9, possibly due to the different study periods or to limited data availability in PRODES caused by cloud cover1. Deforestation fronts also reached or encroached into some ITs/PAs (Fig. 3a and Extended Data Fig. 7). Those PAs with large deforestation areas were more likely to experience downgrading, downsizing and degazettement35. To avoid deforestation and reconcile the conflicts between forest conservation and socio-economic development, both intensive agriculture production in the deforested areas and forest restoration projects need to be fully explored. While these measures are important for environmental quality and for generating employment, they cannot be expected to have a âland sparingâ effect in reducing gross deforestation rates19. PAs with the sustainable-use objective, where small resident populations have low-impact uses of natural resources36, are designed to promote conservation. They have been demonstrated to offer a winâwin solution for biodiversity conservation and socio-economic development37; however, substantial deviations from the intended low-impact uses have sometimes occurred38.
Investigation of various driving factors and development of models that predict the effects of ITs/PAs on reducing deforestation is important for stake-holders and decision-makers. Several studies identified a few driving factors (for example, Indigenous land rights, institutional context and property rights of ITs) and different methods (spatial matching and regression discontinuity model)39,40,41,42. Our exploratory data analyses (Supplementary Information) reveal the importance and limitation of the governance (ITs, nPAs and sPAs), management (strict protection and sustainable use), sizes and locations of the ITs/PAs and forest areas within the ITs/PAs for reducing deforestation. Our preliminary study also highlights the need for casual analysis and more efforts by the research community, stake-holders and decision-makers, more data collection (for example, politics43, social-economic conditions and management practices in individual ITs/PAs), integrated analytics and models across multiple spatial and temporal scales.
Challenges for forest conservation in ITs/PAs
The changes in laws, policies, agriculture (soybean and beef cattle) and climate strongly affect forest conservation44,45 and recent changes in those factors suggest that forest conservation in the ITs/PAs and non-PAs could face increasing challenges in the coming years. As shown in Extended Data Fig. 8, the Action Plan for the Prevention and Control of Deforestation in the Legal Amazon (PPCDAm)23,46, which reduced deforestation, had been implemented since 2004 but was interrupted in 2019 by the Bolsonaroâs administration. The Forest Code and the credit restrictions to owners of deforested lands implemented by the Brazilâs Central Bank as of 200819, which also reduced deforestation, had also been weaken since 201247. The Bolsonaroâs administration is well known for its proximity to agribusiness18,48 and since 2019 has implemented various measures that weaken forest policies and laws and impede their enforcement. In 2021, there was an even stronger pushback against the Brazilian legal framework governing the PAs: five draft bills (PL 490/2007, PL 191/2020, PL 2633/2020, PL 2159/2021 and PLS 510/2021) would further loosen constraints to the economic activities in the ITs, reduce government authority over PAs and provide incentives for agricultural expansion in these areas. The Brazilian Supreme Court is currently trying a case (RE 1017365) that could restrict the constitutional provisions that favoured the demarcation of Indigenous lands and could result in only areas effectively possessed by Indigenous peoples in 1988 being granted demarcation as Indigenous lands. In addition, as the prices of soybeans and beef increased substantially in 2021 (Extended Data Fig. 9), more deforestation for agriculture expansion is likely to continue.
Our study also reveals that annual gross forest area loss rates were high in the El Niño and tropical Atlantic dipole drought years, such as 2005, 2007, 2010 and 2015 (Fig. 3), which may be related to the extensive tree mortality caused by the droughtâfire interactions45. Over the period 2002â2021, a total of 10.4âÃâ106âha forest was burned in ITs/PAs and 55.1âÃâ106âha forest in non-PAs and out of these burned forest areas deforestation occurred in 31% for ITs/PAs and 48% for non-PAs (Extended Data Fig. 10). As the Amazon is projected to have more frequent and severe droughts in the future49, how to manage and protect the remaining intact forests in the ITs/PAs continues to be a major concern.
Summary
The results of this study confirm the critical role of ITs/PAs in forest conservation and raise serious concerns about increased deforestation in the BLA in 2018â2021, especially in the ITs/PAs. The large infrastructure projects50 and the severe COVID-19 pandemic could further increase deforestation and forest degradation in the ITs/PAs. Both deforestation and forest degradation had severe potential impacts on biodiversity and carbon stock in the BLA2,51,52. More attention is urgently needed to strengthen the environmental policies and laws, uphold the existing legal protections and resist the changes being staged in Brazilâs National Congress. More investments from the Brazilian government, private companies and international organizations for the expansion and management of ITs/PAs are also critically needed53. Further reducing deforestation and forest degradation and supporting forest conservation and Indigenous people could prevent passing the tipping point for the Amazon forest ecosystems to flip into savanna ecosystems54.
Methods
Annual forest cover data from MODIS forest mapping tool
We published annual forest maps in the BLA during 2000â2017, which were generated by using the time-series MOD09A1 data product (8âd temporal resolution, 500âm spatial resolution) and the forest mapping algorithm1,2,55. The MOD09A1 8âd composite selects the best-quality observation within each 8âd period. We used cloud-free MOD09A1 observations based on the quality layer. The forest mapping algorithm is mainly based on the unique features of evergreen forest in terms of vegetation greenness, land surface water content and phenology1,2,55, specifically, all good-quality observations for an evergreen forest pixel in a year have the enhanced vegetation index (EVI)ââ¥â0.2 and the land surface water index (LSWI)ââ¥â0. Our annual MODIS forest maps in the BLA had high overall accuracy (>97%) when they were validated by three independent ground reference datasets1: (1) 18 5âÃâ5âkm2 sample blocks from the Global Land Cover Validation Reference Dataset, which was produced from analyses of very high spatial resolution images at 2âm spatial resolution (a total of 1,268 pixels at the 500âm spatial resolution); (2) 416 10âÃâ10âkm2 sample blocks from the TREES-3 forest and non-forest dataset at 30âm spatial resolution (a total of 262,514 pixels at the 500âm spatial resolution); and (3) 1,991 stratified random sample pixels at 500âm spatial resolution for forest changes generated by visual interpretation of time-series Landsat-5 TM, Landsat 7 ETM+ and Landsat 8 OLI images at 30âm spatial resolution. Our annual forest maps have area and spatial distribution similar to the forest maps derived from microwave images1, which are less affected by frequent clouds in the Amazon.
In this study, we made a minor improvement to the forest mapping algorithm, as a few pixels were contaminated by clouds or aerosols but may not be detected by the quality layer, which resulted in the EVI values dropping substantially although LSWI and the normalized difference vegetation index (NDVI) did not change much. Thus, we added another criterion to identify evergreen forest, that is, we classified a pixel as evergreen forest when it had >90% cloud-free observations that met the criteria of EVIââ¥â0.2 and LSWIââ¥â0 and had annual minimum LSWIââ¥â0. We analysed MOD091 data from February 2000 to December 2021 and we generated the annual evergreen forest maps from 2000 to 2021. We applied a 3âyr consistency check procedure to reduce the potential error in the annual maps of evergreen forest. We used the annual evergreen forest maps to generate two reports. One report is on forest area by year. In this report, we counted primary forest and secondary forest together as forest (no separation into these two categories). The second report is on âprimary forestâ (used 2001 as the reference year) and primary forest area loss. Here, we tracked which year deforestation first occurred for individual pixels. It is possible that some of forest pixels in 2001 were âreforested or recovered forestâ. As we have no data before 2000, we do not know how many pixels there were for this case, thus we kept this caveat during our data analysis and result interpretation.
Annual forest area data from the MapBiomas project
The MapBiomas project was launched in 2015 and generates annual land-cover and land-use maps in Brazil31. The algorithm theoretical basis document (ATBD) of MapBiomas presents the cross-reference of the MapBiomas land-cover and land-use classes with classes from other classification systems, including Food and Agriculture Organization of the United Nations (FAO), Brazilian Institute of Geography and Statistics (IBGE) and National Greenhouse Gas Emissions Inventory. The MapBiomas project (Collection 6) uses six steps to generate the annual land-cover and land-use maps by analyses of surface reflectance data from Landsat Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+) and the Operational Land Imager and Thermal Infrared Sensor (OLI-TIRS). The first step is to generate annual Landsat mosaics for specific temporal windows. The second step is to derive the spectral and temporal attributes from Landsat spectral bands to train the Random Forest classifier. The third step is to generate annual land-cover and land-use maps in each biome and cross-cutting theme using the random Forest algorithm (the classification of aquaculture, mining, irrigation, rice and citrus is based on the U-Net convolutional neural network classifier) and training samples. The fourth step is to apply spatiotemporal filters to reduce noise, including gap fill, spatial filter, a 3â5âyr temporal filter, frequency filter and incident filter. The fifth step is to merge the filtered land-cover and land-use maps of each biome and cross-cutting themes and apply the spatiotemporal filters again. The sixth step is the accuracy assessment based on the 75,000 independent samples per year from 1985 to 2018. At the level-1 (forest, non-forest natural formation, farming, non vegetation area and water) classes, the land-cover and land-use maps have 91% global accuracy. We used the level-1 class of forest from the MapBiomas land-cover and land-use maps in this study. For more information about the MapBiomas project see https://mapbiomas-br-site.s3.amazonaws.com/Metodologia/ATBD_Collection_6_v1_January_2022.pdf.
Annual forest cover loss data from the GFW
The 30âm GFW (v.1.8) annual forest cover loss in 2001â2020 was generated by using the decision-tree algorithms and time-series Landsat images acquired during the growing season29. In terms of the year in which a pixel experienced forest loss, the data producers reported that they are 75% confident that the forest loss occurred within the stated year and 97% confident that it occurred within the period from 1âyr before to 1âyr after the stated year.
Annual deforestation data from PRODES
The PRODES dataset28 is generated by Brazilâs National Institute for Space Research (INPE) and is the official Brazilian deforestation dataset. PRODES selected high spatial resolution images (tens of metres) from Landsat-5/7/8, ChinaâBrazil Earth Resources Satellite (CBERS-2/2B), Indian Remote Sensing Satellites (IRS-1) and United Kingdom-Disaster Monitoring Constellation-2 (UK-DMC-2) with two criteria: (1) acquisition dates close to the reference date (1 August) and (2) cloud cover as low as possible. Both visual interpretation and digital image classification were used to identify forest and non-forest pixels.
Active fire data during 2002â2021
We used the MODIS active fire data product: MODIS/Terra Thermal Anomalies/Fire 8-Day L3 Global 1âkm SIN Grid. We only used nominal and high-confidence observations.
The ITs and PAs data
We downloaded the data for ITs and PAs from the Amazonian Network for Socio-environmental Information (RAISG)56. The dataset covers the 1980â2018 period and its attributes include the boundary maps (shapefiles) of individual ITs/PAs, names, years of establishment, governance (by national or state government agency) and management objectives (strict protection or sustainable use; Extended Data Figs. 1 and 2). There are 387 ITs, 146 nPAs and 191 sPAs within the BLA in this study, including the establishment of 145 ITs, 55 nPAs and 72 sPAs after 2002.
The boundary map of BLA
BLA includes nine states: Amazonas, Pará, Mato Grosso, Amapá, Roraima, Acre, Rondônia, Tocantins and Maranhão.
Geospatial data analyses
We used the resultant annual evergreen forest maps, the boundary maps for the BLA, individual ITs/PAs in ArcGIS software to calculate annual forest area by the BLA, ITs and PAs. We also used the evergreen forest map for 2001 as the reference map and calculated the annual gross forest area loss from 2002 to 2021. We grouped individual PAs first by governanceâ(1) nPAs and (2) sPAsâand second by management objectivesâ(1) strict protection and (2) sustainable use. We then compared the forest area and deforestation dynamics among these four types of PA.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
MODIS evergreen forest maps are available at Figshare: https://figshare.com/articles/dataset/Annual_evergreen_forest_cover_maps_in_the_Brazilian_Amazon_from_2000_to_2021/21298497. MOD09A1 and MOD14A2 products are available at the US Geological Survey Land Processes Distributed Active Archive Center (LP DAAC): https://lpdaac.usgs.gov/. The PRODES project data are available from INPE: http://www.obt.inpe.br/OBT/assuntos/programas/amazonia/prodes. The GFW product is available from the University of Maryland: https://glad.earthengine.app/view/global-forest-change#dl=1;old=off;bl=off;lon=20;lat=10;zoom=3. The MapBiomas data are available from https://mapbiomas.org/en. The ITs/PAs boundary maps are available from RAISG: https://www.amazoniasocioambiental.org/en/.
Code availability
The Google Earth Engine code for evergreen forest mapping is available at Figshare: https://figshare.com/articles/software/Code_for_evergreen_forest_cover_mapping_in_the_Brazilian_Amazon_version_2_/21298725. Other code for data analysis in R x64 3.4.2, ENVI/IDL 5.2 and Matlab R2017a should be addressed to X.X. (xiangming.xiao@ou.edu).
References
Qin, Y. et al. Improved estimates of forest cover and loss in the Brazilian Amazon in 2000â2017. Nat. Sustain. 2, 764â772 (2019).
Qin, Y. et al. Carbon loss from forest degradation exceeds that from deforestation in the Brazilian Amazon. Nat. Clim. Change 11, 442â448 (2021).
Jenkins, C., Pimm, S. & Joppa, L. Global patterns of terrestrial vertebrate diversity and conservation. Proc. Natl Acad. Sci. USA 110, E2602âE2610 (2013).
Nogueira, E., Yanai, A., de Vasconcelos, S., de Alencastro, G. & Fearnside, P. Brazilâs Amazonian protected areas as a bulwark against regional climate change. Reg. Environ. Change 18, 573â579 (2018).
Ochoa-Quintero, J., Gardner, T., Rosa, I., Ferraz, S. & Sutherland, W. Thresholds of species loss in Amazonian deforestation frontier landscapes. Conserv. Biol. 29, 440â451 (2015).
Cabral, A., Saito, C., Pereira, H. & Laques, A. Deforestation pattern dynamics in protected areas of the Brazilian Legal Amazon using remote sensing data. Appl. Geogr. 100, 101â115 (2018).
Nepstad, D. et al. Inhibition of Amazon deforestation and fire by parks and Indigenous lands. Conserv. Biol. 20, 65â73 (2006).
Ricketts, T. et al. Indigenous lands, protected areas, and slowing climate change. PLoS Biol. 8, e1000331 (2010).
Herrera, D., Pfaff, A. & Robalino, J. Impacts of protected areas vary with the level of government: comparing avoided deforestation across agencies in the Brazilian Amazon. Proc. Natl Acad. Sci. USA 116, 14916â14925 (2019).
Jusys, T. Changing patterns in deforestation avoidance by different protection types in the Brazilian Amazon. PLoS ONE 13, e0195900 (2018).
Matricardi, E. et al. Long-term forest degradation surpasses deforestation in the Brazilian Amazon. Science 369, 1378â1382 (2020).
Silva, C. et al. Benchmark maps of 33 years of secondary forest age for Brazil. Sci. Data 7, 269 (2020).
Laurance, W. et al. The future of the Brazilian Amazon. Science 291, 438â439 (2001).
Laurance, W. et al. Development of the Brazilian Amazon. Response. Science 292, 1652â1654 (2001).
Silveira, J. Development of the Brazilian Amazon. Science 292, 1651â1654 (2001).
Kauano, Ã., Silva, J., Diniz, J. & Michalski, F. Do protected areas hamper economic development of the Amazon region? An analysis of the relationship between protected areas and the economic growth of Brazilian Amazon municipalities. Land Use Policy 92, 104473 (2020).
Silveira, F., Ferreira, M., Perillo, L., Carmo, F. & Neves, F. Brazilâs protected areas under threat. Science 361, 459â459 (2018).
Begotti, R. & Peres, C. Brazilâs indigenous lands under threat. Science 363, 592â592 (2019).
Fearnside, P. Deforestation of the Brazilian Amazon. Oxford Research Encyclopedias: Environmental Science (Oxford Univ. Press, 2017); https://doi.org/10.1093/acrefore/9780199389414.013.102
Ferreira, J. et al. Brazilâs environmental leadership at risk. Science 346, 706â707 (2014).
Villén-Pérez, S., Anaya-Valenzuela, L., Conrado da Cruz, D. & Fearnside, P. Mining threatens isolated indigenous peoples in the Brazilian Amazon. Glob. Environ. Change 72, 102398 (2022).
Tollefson, J. Illegal mining in the Amazon hits record high amid Indigenous protests. Nature 598, 15â16 (2021).
Silva, C. et al. The Brazilian Amazon deforestation rate in 2020 is the greatest of the decade. Nat. Ecol. Evol. 5, 144â145 (2021).
Vale, M. et al. The COVID-19 pandemic as an opportunity to weaken environmental protection in Brazil. Biol. Conserv. 255, 108994 (2021).
Charlier, P. & Varison, L. Is COVID-19 being used as a weapon against Indigenous Peoples in Brazil? Lancet 396, 1069â1070 (2020).
Davidson, E. et al. The Amazon basin in transition. Nature 481, 321â328 (2012).
Ferrante, L. & Fearnside, P. Brazilâs new president and âruralistsâ threaten Amazoniaâs environment, traditional peoples and the global climate. Environ. Conserv. 46, 261â263 (2019).
PRODES Legal Amazon Deforestation Monitoring System (INPE, 2020); http://www.obt.inpe.br/OBT/assuntos/programas/amazonia/prodes
Hansen, M. et al. High-resolution global maps of 21st-century forest cover change. Science 342, 850â853 (2013).
Qin, Y. et al. Annual dynamics of forest areas in South America during 2007â2010 at 50âm spatial resolution. Remote Sens. Environ. 201, 73â87 (2017).
Collection 6 of the Annual Land Use Land Cover Maps of Brazil (MapBiomas Project, accessed 10 July 2022); https://mapbiomas.org/en
Tree Cover Loss (Global Forest Watch, 2021); https://www.globalforestwatch.org/map/?modalMeta=tree_cover_loss
Fuller, C., Ondei, S., Brook, B. & Buettel, J. Protected-area planning in the Brazilian Amazon should prioritize additionality and permanence, not leakage mitigation. Biol. Conserv. 248, 108673 (2020).
Nolte, C., Agrawal, A., Silvius, K. & Soares, B. Governance regime and location influence avoided deforestation success of protected areas in the Brazilian Amazon. Proc. Natl Acad. Sci. USA 110, 4956â4961 (2013).
Tesfaw, A. et al. Land-use and land-cover change shape the sustainability and impacts of protected areas. Proc. Natl Acad. Sci. USA 115, 2084â2089 (2018).
OECD Environmental Performance Reviews: Brazil (OECD, 2015).
Campos-Silva, J. et al. Sustainable-use protected areas catalyze enhanced livelihoods in rural Amazonia. Proc. Natl Acad. Sci. USA 118, e2105480118 (2021).
Fearnside, P., Nogueira, E. & Yanai, A. Maintaining carbon stocks in extractive reserves in Brazilian Amazonia. Desenvolv. Meio. Ambie. 48, 446â476 (2018).
Nelson, A. & Chomitz, K. Effectiveness of strict vs. multiple use protected areas in reducing tropical forest fires: a global analysis using matching methods. PLoS ONE 6, e22722 (2011).
BenYishay, A., Heuser, S., Runfola, D. & Trichler, R. Indigenous land rights and deforestation: evidence from the Brazilian Amazon. J. Environ. Econ. Manag. 86, 29â47 (2017).
Bonilla-MejÃa, L. & Higuera-Mendieta, I. Protected areas under weak institutions: evidence from Colombia. World Dev. 122, 585â596 (2019).
Baragwanath, K. & Bayi, E. Collective property rights reduce deforestation in the Brazilian Amazon. Proc. Natl Acad. Sci. USA 117, 20495â20502 (2020).
Mangonnet, J., Kopas, J. & Urpelainen, J. Playing politics with environmental protection: the political economy of designating protected areas. J. Politics 84, 1453â1468 (2022).
Nepstad, D. et al. Slowing Amazon deforestation through public policy and interventions in beef and soy supply chains. Science 344, 1118â1123 (2014).
Brando, P. M. et al. Abrupt increases in Amazonian tree mortality due to droughtâfire interactions. Proc. Natl Acad. Sci. USA 111, 6347â6352 (2014).
West, T. & Fearnside, P. Brazilâs conservation reform and the reduction of deforestation in Amazonia. Land Use Policy 100, 105072 (2021).
Soares-Filho, B. et al. Cracking Brazilâs forest code. Science 344, 363â364 (2014).
Ferrante, L. & Fearnside, P. Military forces and COVID-19 as smokescreens for Amazon destruction and violation of indigenous rights. J. Geogr. Soc. 151, 258â263 (2020).
Jiménez-Muñoz, J. et al. Record-breaking warming and extreme drought in the Amazon rainforest during the course of El Niño 2015â2016. Sci. Rep. 6, 33130 (2016).
Ferrante, L. & Fearnside, P. The Amazonâs road to deforestation. Science 369, 634â634 (2020).
Feng, X. et al. How deregulation, drought and increasing fire impact Amazonian biodiversity. Nature 597, 516â521 (2021).
Aragão, L. et al. 21st century drought-related fires counteract the decline of Amazon deforestation carbon emissions. Nat. Commun. 9, 536 (2018).
Silva, J., Barbosa, L., Topf, J., Vieira, I. & Scarano, F. Minimum costs to conserve 80% of the Brazilian Amazon. Perspect. Ecol. Conserv. 20, 216â222 (2022).
Lovejoy, T. & Nobre, C. Amazon tipping point. Sci. Adv. 4, eaat2340 (2018).
Xiao, X., Biradar, C., Czarnecki, C., Alabi, T. & Keller, M. A simple algorithm for large-scale mapping of evergreen forests in tropical America, Africa and Asia. Remote Sens. 1, 355â374 (2009).
Natural Protected Areas and Indigenous Territories Maps in Brazil (RAISG, 2018); https://www.amazoniasocioambiental.org/en/
Acknowledgements
Y.Q. and X.X. were supported by NSF EPSCoR track 1 project (OIA-1946093), NASA âGeostationary Carbon Cycle Observatory (GeoCarb) Missionâ (GeoCarb Contract no. 80LARC17C0001) and NASA Science Team for the OCO Missions (80NSSC21K1077). P.F. is supported by Conselho Nacional de Desenvolvimento CientÃfico e Tecnológico (CNPq) (312450/2021-4), Financiadora de Estudos e Projetos/Rede Brasileira de Pesquisas sobre Mudanças Climáticas Globais (FINEP/Rede Clima) (01.13.0353-00) and Fundação de Amparo à Pesquisa do Estado do Amazonas (FAPEAM) (01.02.016301.000289/2021-33). We thank R. Giovanelli for his assistance on law-relevant work in Brazil.
Author information
Authors and Affiliations
Contributions
X.X. and Y.Q. designed the overall study plan. Y.Q. and X.X. prepared the annual evergreen forest maps. Y.Q., X.X. and F.L. carried out data processing and analysis. X.X., Y.Q. and F.L. interpreted the results. Y.Q. and X.X. drafted the manuscript and F.L., F.S.S., Y.S., E.A. and P.M.F. contributed to the writing and revision of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Sustainability thanks Ana Cabral, Kathryn Baragnawath and José da Silva for their contribution to the peer review of this work.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 Spatial distribution of created years for ITs, national PAs and state PAs.
(a) ITs. (b) National protected areas (nPAs). (c) State protected areas (sPAs). States in the Brazilian Amazon: PA (Pará), AM (Amazonas), RR (Roraima), MT (Mato Grosso), MA (Maranhão), RO (Rondônia), TO (Tocantins), AC (Acre) and AP (Amapá).
Extended Data Fig. 2 Spatial distribution of national PAs (strict protection, sustainable use) and state PAs (strict protection, sustainable use).
(a) National PAs (strict protection, sustainable use). (b) State PAs (strict protection, sustainable use).
Extended Data Fig. 3 Multivariate El NiñoâSouthern Oscillation (ENSO) Index Version 2 (MEI).
ENSO index values are provided by the National Oceanic and Atmospheric Administration (NOAA) Physical Sciences Laboratory (PSL).
Extended Data Fig. 4 Changes in average annual forest area loss rates before and after the ITs and PAs establishment.
(a) The types of ITs and PAs. (b) The management categories of ITs and PAs.
Extended Data Fig. 5 Interannual changes of forest area estimates in the BLA from 2000 to 2021 by this study and the MapBiomas dataset.
The forest area data in 2020 from the MapBiomas dataset is the newest dataset.
Extended Data Fig. 6 Interannual changes of forest area estimates in the nine states in the BLA from 2000 to 2021 by this study and the MapBiomas dataset.
(a) Acre. (b) Amazonas. (c) Roraima. (d) Rondonia. (e) Para. (f) Amapa. (g) Mato Grosso. (h) Tocantins. (i) Maranhao.
Extended Data Fig. 7 Two hotspot regions of annual gross forest area losses in the ITs and PAs from 2002 to 2021.
Each region covers an area of ~400âkmâÃâ400âkm.
Extended Data Fig. 8 Major events in laws, policies, governance and COVID-19 from 2000 to 2021, which are related to forest conservation in the BLA.
PPCDAm means the Action Plan for the Prevention and Control of Deforestation in the Legal Amazon.
Extended Data Fig. 9 Interannual changes of soybean price and live cattle futures.
(a) Soybean price. (b) Live cattle futures. The data are provided from https://www.macrotrends.net/.
Extended Data Fig. 10 Forest fire area and gross forest area loss in the ITs/PAs and BLA from 2002 to 2021.
(a) Interannual changes of forest fire areas and gross forest area loss with fire in the ITs and PAs. (b) Forest fire area, gross forest area loss with fire and gross forest area loss without fire in the ITs and PAs. (c) Forest fire area, gross forest area loss to fire and gross forest area loss without fire in the BLA.
Supplementary information
Supplementary Information
Supplementary data analysis.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Qin, Y., Xiao, X., Liu, F. et al. Forest conservation in Indigenous territories and protected areas in the Brazilian Amazon. Nat Sustain 6, 295â305 (2023). https://doi.org/10.1038/s41893-022-01018-z
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41893-022-01018-z
This article is cited by
-
Climate justice, forests, and Indigenous Peoples: toward an alternative to REDDâ+âfor the Amazon
Climatic Change (2024)
-
Protected areas reduce deforestation and degradation and enhance woody growth across African woodlands
Communications Earth & Environment (2023)