
==== Front
Nature
Nature
Nature
0028-0836
1476-4687
Nature Publishing Group UK London

36859548
5690
10.1038/s41586-022-05690-1
Article
Tropical deforestation causes large reductions in observed precipitation
http://orcid.org/0000-0002-2705-8398
Smith C. ee13c2s@leeds.ac.uk

http://orcid.org/0000-0002-3720-4758
Baker J. C. A.
http://orcid.org/0000-0002-7551-4597
Spracklen D. V.
grid.9909.9 0000 0004 1936 8403 School of Earth and Environment, University of Leeds, Leeds, UK
1 3 2023
1 3 2023
2023
615 7951 270275
14 4 2022
15 12 2022
© The Author(s) 2023
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
Tropical forests play a critical role in the hydrological cycle and can influence local and regional precipitation1. Previous work has assessed the impacts of tropical deforestation on precipitation, but these efforts have been largely limited to case studies2. A wider analysis of interactions between deforestation and precipitation—and especially how any such interactions might vary across spatial scales—is lacking. Here we show reduced precipitation over deforested regions across the tropics. Our results arise from a pan-tropical assessment of the impacts of 2003–2017 forest loss on precipitation using satellite, station-based and reanalysis datasets. The effect of deforestation on precipitation increased at larger scales, with satellite datasets showing that forest loss caused robust reductions in precipitation at scales greater than 50 km. The greatest declines in precipitation occurred at 200 km, the largest scale we explored, for which 1 percentage point of forest loss reduced precipitation by 0.25 ± 0.1 mm per month. Reanalysis and station-based products disagree on the direction of precipitation responses to forest loss, which we attribute to sparse in situ tropical measurements. We estimate that future deforestation in the Congo will reduce local precipitation by 8–10% in 2100. Our findings provide a compelling argument for tropical forest conservation to support regional climate resilience.

A pan-tropical analysis using satellite, station-based and reanalysis datasets shows that deforestation causes reduced precipitation, and demonstrates that the effect increases with spatial scale.

Subject terms

Climate change
Tropical ecology
issue-copyright-statement© The Author(s), under exclusive licence to Springer Nature Limited 2023
==== Body
pmcMain

Tropical forests play an important role in moderating local, regional and global climate through their impact on energy, water and carbon cycles3. Crucially, tropical forests control local-to-regional rainfall patterns1,2. Evapotranspiration from tropical forests is a strong driver of regional precipitation4,5 contributing up to 41% of basin mean rainfall over the Amazon and up to 50% over the Congo6. Evergreen tropical forests are dependent on high annual rainfall for their survival and productivity7, and forest–rainfall feedbacks have been highlighted as an important determinant of tropical forest stability4,5,8, amid concerns that the exacerbating impacts of droughts and deforestation could threaten their viability9.

Rapid loss of forests is occurring across the tropics10. Tropical deforestation warms the climate at local-to-global scales by changing the surface energy balance and through emissions of carbon dioxide3. The impact of tropical deforestation on precipitation is less certain with a range of processes operating at different scales. Small-scale deforestation over the southern Amazon has been shown to increase precipitation frequency11,12 owing to thermally13 and dynamically12 induced circulations. At larger scales, deforestation reduces precipitation recycling leading to a reduction in precipitation1,14. Over Indonesia, deforestation has been linked to declining precipitation15, and exacerbation of El Niño impacts16. Global and regional climate models predict annual precipitation declines of 8.1 ± 1.4% for large-scale Amazonian deforestation by 2050 (ref. 17), but an observational study of the impacts of tropical deforestation on precipitation across spatial scales is lacking.

Here we present a pan-tropical assessment of the impact of forest loss on precipitation based on measurements. We use a satellite dataset of forest cover change over the period 2003–2017 to identify areas of forest loss, with a focus on evergreen broadleaf forests of the Amazon, Congo and Southeast Asia (SEA; Fig. 1). To provide a robust assessment of the impacts of deforestation on precipitation, we analysed 18 different precipitation datasets, including satellite (n = 10), station-based (n = 4) and reanalysis (n = 4) products (Extended Data Table 1). We compared the precipitation change over pixels experiencing forest loss with neighbouring pixels that had experienced less forest loss (Methods). This comparison against neighbouring pixels that will have experienced similar climate change focuses our analysis on the changes due to forest loss. To explore the impact of forest loss across scales, we analysed the impacts of forest loss on coincident precipitation at a series of spatial resolutions ranging from roughly 5 km to 200 km (0.05°, 0.1°, 0.25°, 0.5°, 1.0° and 2.0°).Fig. 1 Tropical evergreen broadleaf forest cover loss from 2003 to 2017.

a–f, Forest cover change at 0.05° (a), 0.1° (b), 0.25° (c), 0.5° (d), 1.0° (e) and 2.0° (f) resolution. The Amazon Basin, Congo Basin and SEA regions used in this study are outlined in purple. Maps of the different regions generated using Cartopy and Natural Earth51. Forest loss data from ref. 10.

Source Data

Precipitation response to forest loss

Observed precipitation responses to tropical forest loss across multiple spatial scales and precipitation products are presented in Fig. 2. Satellite-based precipitation datasets suggest that tropical forest loss causes statistically significant (P < 0.05) declines in median annual mean precipitation at all scales analysed. At larger scales (>0.5°), reductions exceed 0.03 mm per month for each percentage point loss of forest cover (Fig. 2d–f). The largest changes are observed at the 2.0° scale (approximately 220 km at the Equator; Fig. 2f), for which each percentage point reduction in forest cover causes 0.25 ± 0.1 mm per month reduction in annual precipitation.Fig. 2 Reductions in precipitation over regions of tropical forest loss.

a–r, Bars indicate the median absolute change in annual precipitation (millimetres per month) per percentage point of forest loss over 2003 to 2017 in each region (tropics (a–f), Amazon (g–l), Congo (m–r), SEA (s–x)) for each precipitation dataset category (satellite, station and reanalysis). Results are shown for forest loss scales of 0.05° (a,g,m,s), 0.1° (b,h,n,t), 0.25° (c,i,o,u), 0.5° (d,j,p,v), 1.0° (e,k,q,w) and 2.0° (f,l,r,x). Statistically significant (*P < 0.05; **P < 0.01) and nonsignificant (NS) differences in changes in mean precipitation (calculated as a multi-annual mean over 2003–2007 compared with 2013–2017) over deforested regions compared with control regions are indicated. Error bars show ±1 standard error from the mean. Datasets used in this analysis are detailed in Extended Data Table 1. ΔP, precipitation change.

Source Data

Analysis of precipitation change as a function of forest loss confirms larger reductions in precipitation for larger reductions in forest cover (Extended Data Fig. 1), although with considerable variability, as seen in the modelled response18. Observed reductions in precipitation are consistent across satellite datasets, with all ten satellite precipitation products agreeing on the sign of the rainfall response at 2° over the tropics (Extended Data Fig. 2). At the 2° scale, significant (P < 0.05) reductions in annual mean precipitation with forest loss were observed across all tropical regions (Fig. 2). Reductions in precipitation at 2° based on satellite datasets ranged from 0.48 ± 0.36 mm per month in SEA to 0.23 ± 0.12 mm per month in the Amazon, and 0.21 ± 0.19 mm per month in the Congo for each percentage point loss in forest cover, with at least 8 out of 10 satellite datasets agreeing on the sign of the response within each region (Extended Data Fig. 2). In SEA, it has been suggested that proximity to the ocean and the replacement of tropical forest with plantations as opposed to pasture or cropland may cause reduced sensitivity of precipitation to deforestation1. Our analysis suggests that forest loss in SEA causes reductions in precipitation consistent with or greater than reductions in precipitation in the Amazon and Congo.

Station-based datasets and reanalysis products exhibit contrasting annual mean precipitation responses to deforestation at 2.0° (Fig. 2). Across the tropics, station-based and reanalysis datasets showed no statistically significant changes in annual mean precipitation due to forest loss (Fig. 2f), and there was little agreement with satellite datasets at the regional scale (Fig. 2l,r,x), with some non-satellite precipitation products showing small increases in annual mean precipitation due to forest loss. Sparse in situ measurements across the tropics, particularly over regions of forest loss, mean that station-based datasets provide a weak constraint on precipitation changes. A comparison of station-based precipitation datasets revealed higher levels of uncertainty in the tropics, including the Amazon19. In regions of sparse data such as tropical forests20, interpolation methods may mask precipitation changes driven by forest loss. Reanalysis products, which are numerical models constrained by empirical data, are also expected to be less reliable in regions where in situ data are limited21. Our results indicate that precipitation data based on satellite remote-sensing measurements may have an advantage over tropical forest regions where in situ measurements are sparse or unavailable. For these reasons, we focus our analysis on satellite-based datasets and identify where agreement between datasets exists.

Our results are robust (Extended Data Fig. 3) to a range of methodological assumptions including the length of analysis period, the choice of start and end period and the spatial extent of control pixels (Methods). Our analysis period includes the 2015–2016 El Niño that resulted in negative precipitation anomalies over many tropical land regions (Supplementary Fig. 1). We found that the precipitation response to forest loss was robustly negative during both El Niño and non-El Niño years (Extended Data Fig. 3). Over the Amazon and SEA, we see a stronger reduction in precipitation over regions of forest loss during El Niño years. The relative impact of El Niño on precipitation is smaller in the Congo22, and correspondingly we do not see a stronger reduction here. A stronger precipitation response to forest loss in regions and periods impacted by El Niño is probably due to higher transpiration rates observed in tropical forests during El Niño years23 and because rainfall is more sensitive to reductions in moisture recycling during drought years5,24. Climate change is expected to lead to increased droughts over many tropical regions25, which may be further exacerbated by ongoing deforestation.

Seasonal precipitation reductions

Changes in precipitation due to forest loss during the dry, wet and transition seasons are nearly consistently negative (Fig. 3). For the tropics, absolute changes in precipitation with forest loss are greatest in the wet season (Fig. 3a, up to −0.6 mm per month per percentage point forest loss) whereas relative changes of precipitation with forest loss are similar (−0.2% per percentage point) across dry, wet and transition seasons (Supplementary Fig. 2). In the Amazon, deforestation causes the largest reductions in precipitation during the transition season (Fig. 3b) as has been found previously18,26,27.Fig. 3 Changes in seasonal precipitation due to forest loss.

a–d, Bars indicate the median change in precipitation (millimetres per month) per percentage point forest cover loss for satellite datasets during 2003–2017 for tropics (a), Amazon (b), Congo (c) and SEA (d). Error bars indicate ±1 standard error from the mean. Statistically significant (*P < 0.05; **P < 0.01) and nonsignificant (NS) differences in changes in mean precipitation over deforested regions compared with controls are indicated. Results are shown for the wettest 3 months (wet), the driest 3 months (dry) and the transition months (remaining 6 months). Datasets used in this analysis detailed in Extended Data Table 1.

Source Data

Previous case studies have indicated that dry-season precipitation can increase over deforestation in the Amazon11,28,29. We observed a nonsignificant increase in dry-season precipitation due to forest loss in the Amazon at 2° as well as increases in the Congo at 1° and 2° (Fig. 3). In SEA, forest loss causes reductions in dry-season precipitation across all scales (Fig. 3d). The mechanism through which forest loss impacts precipitation is likely to change with both season and spatial scale. At the smallest scales (5 km), thermally driven impacts are likely to dominate, shifting to dynamically driven impacts through reductions to surface roughness, then to reductions in moisture fluxes and precipitation recycling at the largest scales12,30. Our observation of greater reductions in precipitation due to deforestation at larger spatial scales is consistent with a reduction in moisture recycling emerging as the dominant mechanism1.

Comparison with climate models

A meta-analysis of climate model studies (predominantly global models with >2° resolution) found that forest loss in the Amazon resulted in a mean reduction in annual mean precipitation of 0.16 ± 0.13% per percentage point17, overlapping with our value of 0.25% per percentage point (Supplementary Fig. 2). Fewer simulations have been conducted for the Congo, with models predicting a reduction in precipitation of 0.16 ± 0.17% per percentage point2, similar to our reduction of 0.15% per percentage point (Supplementary Fig. 2). The large range of model estimates highlights the substantial uncertainty in model predictions. Our observationally derived analysis provides support for models that predict reductions in precipitation under regional deforestation at global climate model scales.

Our observational analysis documents the impacts of deforestation on precipitation across the tropics. Applying linear scaling to the reductions in precipitation observed in our analysis would suggest that complete deforestation could result in reductions in annual precipitation of 10–20%. Previous estimates of the impact of complete deforestation on precipitation range from a 16% (ref. 17) to 55–70% (ref. 31) reduction in the Amazon and an 18% (ref. 2) to 50% (ref. 32) reduction in the Congo.

Impacts of future deforestation

To further explore how future deforestation might modify precipitation, we combined our observationally derived estimates of precipitation responses to forest cover loss with future projections of land cover change from a high-deforestation scenario (Methods). We estimate that forest loss from 2015 to 2100 (Fig. 4a) could lead to reductions of annual mean precipitation of up to 16.5 ± 6.2 mm per month in the Congo (Fig. 4b), equivalent to precipitation declines of 8–10%. Forest loss is projected to be greatest in the western and southern Congo (Fig. 4c), which will also experience the strongest reductions in precipitation (Fig. 4d).Fig. 4 Impact of projected future forest loss on annual mean precipitation.

a, Mean forest cover loss over 2015–2100 under Shared Socioeconomic Pathway 3–Representative Concentration Pathway 4.5 for the tropics, Amazon, Congo and SEA. b, Impact of projected forest cover loss on precipitation (P; ±1 standard error from the mean). c, Spatial pattern of forest cover loss. d, Predicted P change (∆P) in 2100 due to forest cover loss. Results are shown for 2.0° resolution. Maps of the different regions generated using Cartopy and Natural Earth51.

Source Data

The sensitivity of precipitation to the extent of forest loss is an uncertainty in our analysis, a result of the relatively short observational record, compounded by large spatial and temporal variability in precipitation. The response of precipitation to forest loss greater than 30%, a threshold beyond which large reductions in precipitation have been postulated1, is one such uncertainty. Restricting our analysis to the forest losses of 0–30% that are well sampled in our observational dataset (Supplementary Fig. 3), through capping the impacts of greater forest loss at that of 30%, results in projected annual mean precipitation reductions of 6.5 ± 2.6 mm per month in the Congo and 6.2 ± 2.5 mm per month in SEA (Supplementary Fig. 4). However, restricting our analysis in this way is likely to underestimate the precipitation impacts over regions projected to experience the most extensive deforestation, including the Congo, where mean forest cover is projected to decline by 40 percentage points between 2015 and 2100 (Fig. 4a).

Previous studies have identified both linear9,33 and nonlinear1,31 responses of precipitation to forest loss. Such nonlinear interactions and feedbacks have the potential to further amplify or moderate the responses predicted here14,34. Our analysis shows large reductions in precipitation for relatively small amounts of forest loss and evidence for reduced sensitivity of precipitation to additional amounts of forest loss (Extended Data Fig. 1). Assuming a nonlinear relationship between forest loss and precipitation (Methods) reduces our projected reductions in precipitation by around a factor of 2 (Supplementary Fig. 5). Our observationally based approach will miss tipping points in the climate system that might be reached as deforestation extent progresses further1. Such tipping points have been postulated for the Amazon under future global change25,35. Thus, the substantial declines in precipitation projected in our analysis should be viewed as a conservative estimate of the potential precipitation response to future deforestation. Nevertheless, our analysis suggests that deforestation can drive local and regional precipitation changes that may match or exceed those predicted due to climate change over the same period36,37.

Implications of precipitation reductions

Reductions in precipitation induced by forest loss have important implications for society and the sustainability of remaining tropical forest. Deforestation-induced reductions in precipitation affect agriculture1,14 and hydropower generation38. On average, crop yields decline by 0.5% for each percentage point reduction in precipitation39. Our results indicate that forest-loss-induced changes to annual precipitation (Supplementary Fig. 2) could cause crop yields to decline by 1.25% for each 10-percentage-point loss of forest cover, potentially exacerbating the impacts of climate change and future drought events. The maintenance of regional rainfall patterns due to forests in the Amazon has been valued at up to US$9 ha−1 yr−1 and US$1.84 ha−1 yr−1 through sustaining agricultural yields and hydropower generation, respectively40. Global cropland area increased by 9% in the past two decades, with even higher increases in South America and tropical Africa41 largely at the expense of natural ecosystems. Further agricultural expansion in tropical forest regions may lead to overall reductions in production if declines in yield due to deforestation-induced reductions in rainfall outweigh increased production from expanded agricultural area14.

Furthermore, reductions in rainfall over remaining areas of tropical forest are expected to lead to additional forest loss9 as well as impacting species composition22, carbon sequestration42 and fire frequency43. Reductions in dry-season precipitation pose a particular threat to forest viability by exacerbating seasonal droughts and potentially delaying the onset of the wet season and extending the length of the dry season. Increases in dry-season length over recent decades have previously been reported for the Amazon44 and the Congo45, possibly linked to land cover changes27.

Deforestation may also shift precipitation patterns, increasing dry-season rainfall immediately downwind of forest loss and decreasing rainfall in upwind areas12. Our approach is restricted to observing deforestation impacts up to scales of 200 km (Methods). At larger scales, insufficient pixels experienced forest loss during the relatively short period of satellite observations for a robust analysis. Deforestation is also likely to alter precipitation at these larger scales through reducing moisture recycling leading to reductions in rainfall downwind of forest loss4,5,9,35. The length scale of moisture recycling has been estimated at 500–2,000 km in the tropics46, with a median value of 600 km in the Amazon5. In regions downwind of extensive forests, such as the southwestern Amazon, up to 70% of precipitation could be sourced from upwind evapotranspiration47,48. Tropical forest loss could therefore have severe implications for precipitation in these regions that are hundreds to thousands of kilometres downwind of the forest loss5. Through missing the impacts at these larger scales, our analysis is likely to underestimate the full impacts of deforestation on rainfall.

Our results highlight the importance of remaining tropical forests for sustaining regional precipitation. Despite efforts to reduce deforestation, rates of tropical forest loss have accelerated over the past two decades49. Renewed efforts are needed to ensure recent commitments to reduce deforestation, including the New York Declaration on Forests and The Glasgow Leaders’ Declaration on Forests and Land Use made at the 26th UN Climate Change Conference of the Parties, are successful. Global efforts to restore large areas of degraded and deforested land could enhance precipitation50, reversing some of the reductions in precipitation due to forest loss observed here.

Methods

Datasets

We used 18 precipitation datasets, listed in Extended Data Table 1. All datasets were downloaded at the highest available spatial resolution, which for some datasets was 0.04°, or approximately 4 km at the Equator. Data were obtained as monthly means or converted to monthly mean using the Python package xarray52. We categorized precipitation datasets as satellite (n = 10), station (n = 4) and reanalysis (n = 4). Satellite datasets are those based primarily on data from satellite sensors and include datasets that have both satellite and station-based data (that is, merged datasets). Station datasets include only ground-based information from weather stations and rain gauges. Reanalysis products are models constrained by surface and satellite data. Precipitation datasets have been compared previously over the Amazon20 highlighting the limited station data over tropical forest regions. Time series of precipitation (Supplementary Fig. 1) reveal variability across the different datasets highlighting the need to analyse impacts of deforestation across multiple datasets.

To analyse the changes in forest canopy cover, we used data from the Global Forest Change (GFC) version 1.9 (ref. 10). GFC v1.9 provides forest canopy cover in the year 2000 and subsequent annual forest loss from 2001–2020 at 30-m resolution. We analysed forest cover and precipitation changes over the period 2003 to 2017, which was the period common to all datasets.

Analysis across multiple spatial scales

We analysed the impacts of forest loss across a range of scales (0.05°, 0.1°, 0.25°, 0.5°, 1.0° and 2.0°). Each precipitation dataset was analysed at its native resolution and at all lower resolutions across this range of scales. Spatial regridding was carried out using the Python package xESMF53 with a bilinear regridding scheme. Two alternative regridding methods (xESMF: conservative-normalized; and iris: area weighted) were tested and had little impact on our results. For GFC data, we calculated forest loss using the original 30-m data and converted the resulting values to each of the six spatial resolutions analysed by taking the sum of all 30-m pixels within each larger pixel. Change in canopy cover from 2003 to 2017 at each resolution is shown in Fig. 1.

Assessing impact of historical deforestation on precipitation

We used a moving-window nearest-neighbour approach54 to compare the forest loss and precipitation change of each pixel with that of its immediate neighbours. We tested the sensitivity of the analysis to the size of the moving window and found similar results for 3 × 3 and 5 × 5 (Extended Data Fig. 2) moving windows. Results from the 3 × 3 moving-window approach can been seen in the main paper. We calculated the forest loss of each deforested pixel relative to neighbouring control pixels as the forest loss of the deforested pixel minus forest loss of the control. We constrained our analysis to the tropical evergreen broadleaf biome using the Moderate Resolution Imaging Spectroradiometer land cover dataset55. To be included in the analysis, deforested pixels must have experienced 0.1% more forest loss over time than their neighbouring control pixels. The number of deforested pixels analysed varied between analysis resolutions as follows: 0.05°, n = 243,254; 0.1°, n = 58,660; 0.25°, n = 9,604; 0.5°, n = 2,303; 1.0°, n = 586; 2.0°, n = 123. We observed similar distributions of canopy change for all spatial resolutions analysed (Supplementary Fig. 6).

We calculated the precipitation change of the deforested pixel relative to the precipitation change of the control pixel (ΔP) as the precipitation change of the deforested pixel over the analysis period (for example, 2003–2017) minus the precipitation change over the control pixel. To reduce the impact of interannual variability in precipitation on our results, we calculated 5-yr means for periods at the start (2003–2007) and end (2013–2017; Extended Data Fig. 5) of the analysis period. We then calculated the change in precipitation as the difference between the start and end of these multi-year means. We report precipitation changes (ΔP) as a function of forest loss by dividing by the difference in forest loss between deforestation and control pixels (units of millimetres per month per percentage point). We also report precipitation change as the percentage change in precipitation (ΔP/P, in units of per cent) as a function of forest loss (in units of per cent per percentage point).

To ensure that control pixels and deforested pixels experience a similar background climate, we conducted a sensitivity test in which we restricted our analysis to pixels for which the pre-deforestation precipitation across the control and deforested pixels differed by less than 10%. Restricting our analysis in this way had little impact on our results (Supplementary Fig. 7) showing that our nearest-neighbour approach is effective even at the largest scales analysed here.

To explore the role of the analysis period on our results, we compared the results for 5-yr means to those for shorter 3-yr means (2003–2005 versus 2015–2017) and found consistent results (Extended Data Fig. 3). Our analysis period includes the strong 2015/2016 El Niño that resulted in reductions in precipitation over most tropical land regions, particularly in 2015 (Supplementary Fig. 1). To explore the potential impacts of the 2015/2016 El Niño on our analysis, we estimated the impact of forest loss on precipitation using 3-yr (2003–2005 versus 2018–2020) and 5-yr (2003–2007 versus 2016–2020) multi-annual means spanning an extended time period. The 3-yr analysis completely excludes the 2015/2016 ENSO, and the 5-yr analysis excludes 2015, which was the driest year (Extended Data Fig. 3). Two datasets (TRMM and UDEL) were not available after 2017, so they were removed from this sensitivity analysis.

Statistical analysis

For each category of precipitation data (satellite, station and reanalysis), precipitation change values were grouped together for all deforestation pixels and all control pixels. We found that precipitation changes for deforested pixels and control pixels, and the difference in precipitation change between deforested and control pixels (Extended Data Fig. 4), were normally distributed. Error bars (Figs. 2 and 3) show ±1 standard error from the mean calculated and displayed using the Python package Seaborn56. To test whether mean precipitation changes over regions of deforestation were statistically different from changes over the control areas, we used a Student’s t-test. We also used the Mann–Whitney test to test for significant differences in median precipitation change between control and deforested pixels and found similar results.

Seasonal analysis

For the satellite datasets alone, in addition to calculating precipitation changes at the annual timescale, we calculated changes for the dry season (driest 3 months of each year), wet season (wettest 3 months of each year) and transition season (remaining 6 months). The driest, wettest and transition months were identified for each pixel using each individual precipitation dataset. For each season and dataset, we calculated the median change in precipitation across all of the pixels within the region of interest (Supplementary Figs. 8–10).

Predicting future precipitation change due to forest loss

We used projections of forest cover change available at 0.05° from the Global Change Analysis Model (GCAM) for 2015–2100 based on the Shared Socioeconomic Pathway 3–Representative Concentration Pathway 4.5 scenario, which represents a high-deforestation future57. GCAM includes the impacts of climate and land use on future forest cover. We summed forest cover from all forest categories and calculated forest cover loss in each year compared to a 2015 baseline. Forest cover loss data were regridded to 2°. We estimated the impact of forest loss on future precipitation at the 2° scale through multiplying the projected percentage point forest loss for each pixel by the observed median change in precipitation (millimetres per month) per percentage point forest cover loss across the satellite datasets. To estimate the uncertainty in our predictions, we applied an upper and lower limit on the sensitivity of precipitation to forest loss based on the median value ±1 standard error from the mean (see error bars in Fig. 2) and rescaled by forest loss. This provides a range of estimated precipitation impacts of future forest loss. We also tested the impact on our results of capping future forest loss in each pixel at 30%, which is the upper range of forest loss that is well sampled in the observations (Supplementary Fig. 3). For each region, we applied the tropical satellite precipitation response to forest loss (Fig. 2f), meaning that our projected regional precipitation changes are a product of the regional canopy cover change and the median tropical precipitation response. Our approach assumes a linear precipitation response to forest loss, which recent work suggests could provide a conservative estimate of deforestation impacts31. We tested the sensitivity of assuming a linear response of precipitation to canopy cover loss. We fitted a nonlinear function to the data presented in Extended Data Fig. 1 through applying the median sensitivity of precipitation to forest cover loss (millimetres per month per percentage point) within each forest cover loss bin. We then scaled by the projected forest cover loss. This approach reduces the projected reduction in precipitation to 2.4 mm per month in SEA and 1.5 mm per month in the Congo (Supplementary Fig. 5).

Online content

Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at 10.1038/s41586-022-05690-1.

Supplementary information

Supplementary Figures This file Supplementary Figs. 1–10, which support the results presented in the main paper. Here we also present results from the individual precipitation datasets.

Peer Review File

Source data

Source Data Fig. 1

Source Data Fig. 2

Source Data Fig. 3

Source Data Fig. 4

Source Data Extended Data Fig. 1

Source Data Extended Data Fig. 2

Source Data Extended Data Fig. 3

Source Data Extended Data Fig. 4

Source Data Extended Data Fig. 5

Extended data figures and tables

Extended Data Fig. 1 Annual precipitation change as a function of forest loss.

Results are shown at 2° spatial resolution for all satellite precipitation (P) datasets calculated as the change in P over time for deforested data pixels minus change over time for control data pixels. Data is binned according to forest cover change (%) with an equal number of pixels in each bin. Points show the median and error bars show ± 1 standard error from the mean. Details of each data product are provided in Extended Data Table 1.

Source Data

Extended Data Fig. 2 Annual precipitation change due to forest loss for individual datasets.

Results are shown for 2003 – 2017 for 5 year averages and 3x3 moving window. Bars show the median absolute change in annual P (mm month−1) per percentage point tree cover loss in each region (Tropics (a-f), Amazon (g-l), Congo (m-r), SEA (s-x)). Each P dataset is shown separately and ordered and coloured by category: satellite (orange), station (yellow) and reanalysis (turquoise). The datasets are numbered; 1) CHIRPS, 2) CMORPH, 3) CPC, 4) CRU, 5) ERA5, 6) GPCC, 7) GPCP, 8) GPM, 9) JRA, 10) MERRA-2, 11) NOAA 12) PERSIANN-CCS, 13) PERSIANN-CCSCDR, 14) PERSIANN-CDR, 15) PERSIANN-NOW, 16) PERSIANN, 17) TRMM, 18) UDEL. Results are shown for forest loss scales of 0.05° (a,g,m,s), 0.1° (b,h,n,t), 0.25° (c,i,o,u), 0.5° (d,j,p,v), 1.0° (e,k,q,w), 2.0° (f,l,r,x). Details of each data product are provided in Extended Data Table 1.

Source Data

Extended Data Fig. 3 Changes in precipitation due to forest loss for different time periods and nearest neighbour comparisons.

Changes in annual mean precipitation at 2.0° resolution are shown for satellite (orange), station (yellow) and reanalysis (turquoise) datasets for the tropics (a-f), Amazon (g-l), Congo (m-r) and Southeast Asia (SEA, s-x). Columns show the sensitivity of our results to changes in the analysis period, number of years used to compute multi-annual means at start and end of the analysis period, and size of the moving window used for nearest neighbour comparisons: 2003-2017, 3-year averages and 3x3 nearest neighbour (Column 1, a,g,m,s); 2003-2017, 3-year, 5x5 (Column 2; b,g,n,t); 2003-2017, 5-year, 3x3 (Column 3; c,i,o,u); 2003-2017, 5-year, 5x5 (Column 4; d,j,p,v); 2003-2020, 3-year, 3x3 (Column 5; e,k,q,w); 2003-2020, 5-year, 3x3 (Column 6; f,l,r,x). Error bars show ± 1 standard error from the mean. Details of each data product are provided in Extended Data Table 1. Full results for all tested resolutions are available in an online repository [10.5281/zenodo.7373832].

Source Data

Extended Data Fig. 4 Change in precipitation over deforested, control and difference between deforested and control pixels.

Change in precipitation over 2003 to 2017 is shown for deforested (a, b), control (c, d) and difference between deforested and control pixels (e, f) for 0.05° (a, c, e) and 2.0° (b, d, f) resolution. Details of each data product are provided in Extended Data Table 1.

Source Data

Extended Data Fig. 5 Mean precipitation from satellite, station and reanalysis datasets.

For each class of dataset, satellite (a, d, g), station (b, e, h) and reanalysis (c, f, i), the median value for the 5-year multi-annual mean at the start (2003-2007; a, b, c) and end (2013-2017; d, e, f) of the analysis period as well as the change over the analysis period (end – start; g, h, i) is shown. Mean values across tropical evergreen broadleaf forests are shown in units of mm/month at the top of each panel. Maps of the different regions generated using Cartopy and Natural Earth51. Details of each data product are provided in Extended Data Table 1.

Source Data

Extended Data Table 1 Precipitation datasets used in this study refs 58–73

Precipitation datasets used in this study

Extended data

is available for this paper at 10.1038/s41586-022-05690-1.

Supplementary information

The online version contains supplementary material available at 10.1038/s41586-022-05690-1.

Acknowledgements

The research has been supported by funding from the European Research Council under the European Union’s Horizon 2020 research and innovation programme (DECAF project, grant agreement no. 771492), and the Newton Fund, through the Met Office Climate Science for Service Partnership Brazil.

Author contributions

All authors developed the concept of the study, and contributed to experimental design, interpretation of results and writing of the manuscript. C.S. and J.C.A.B. contributed to the analysis.

Peer review

Peer review information

Nature thanks Arie Staal and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Data availability

Full results for all tested resolutions used in this analysis are available through 10.5281/zenodo.7373832. The original datasets are freely available to download from the following repositories: CHIRPS from https://data.chc.ucsb.edu/products/?C=M;O=D, CMORPH from https://ftp.cpc.ncep.noaa.gov/precip/CMORPH_RT/GLOBE/data/, CPC from https://psl.noaa.gov/data/gridded/data.cpc.globalprecip.html, CRU from https://crudata.uea.ac.uk/cru/data/hrg/, ERA5 from https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=overview, GPCC from https://opendata.dwd.de/climate_environment/GPCC/html/download_gate.html, GPCP from https://disc.gsfc.nasa.gov/datasets/GPCPMON_3.1/summary?keywords=GPCPMON, GPM from https://gpm1.gesdisc.eosdis.nasa.gov/data/GPM_L3/, JRA from https://climatedataguide.ucar.edu/climate-data/jra-55 and https://jra.kishou.go.jp/JRA-55/index_en.html, MERRA-2 from https://disc.gsfc.nasa.gov/datasets?project=MERRA-2, NOAA (PREC/LAND) from https://psl.noaa.gov/data/gridded/data.precl.html, PERSIANN (CCS, CDR, CCS-CDR, PDIR-NOW) from https://chrsdata.eng.uci.edu/, TRMM from https://disc.gsfc.nasa.gov/datasets/TRMM_3B43_7/summary, UDEL from https://psl.noaa.gov/data/gridded/data.UDel_AirT_Precip.html. The GCAM model output used in this study is available from 10.25584/data.2020-07.1357/1644253. Source data are provided with this paper.

Code availability

The code used in this analysis is available through 10.5281/zenodo.7373832.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Lawrence D Vandecar K Effects of tropical deforestation on climate and agriculture Nat. Clim. Change 2015 5 27 36 10.1038/nclimate2430
2. Spracklen DV Baker JCA Garcia-Carreras L Marsham JH The effects of tropical vegetation on rainfall Annu. Rev. Environ. Resour. 2018 43 193 218 10.1146/annurev-environ-102017-030136
3. Bonan GB Forests and climate change: forcings, feedbacks, and the climate benefits of forests Science 2008 320 1444 1449 10.1126/science.1155121 18556546
4. Spracklen DV Arnold SR Taylor CM Observations of increased tropical rainfall preceded by air passage over forests Nature 2012 489 282 285 10.1038/nature11390 22951966
5. Staal A Forest-rainfall cascades buffer against drought across the Amazon Nat. Clim. Change 2018 8 539 543 10.1038/s41558-018-0177-y
6. Baker JCA Spracklen DV Divergent representation of precipitation recycling in the Amazon and the Congo in CMIP6 models Geophys. Res. Lett. 2022 49 e2021GL095136 10.1029/2021GL095136 35859721
7. Guan K Photosynthetic seasonality of global tropical forests constrained by hydroclimate Nat. Geosci. 2015 8 284 289 10.1038/ngeo2382
8. Staal A Hysteresis of tropical forests in the 21st century Nat. Commun. 2020 11 4978 10.1038/s41467-020-18728-7 33020475
9. Zemp DC Self-amplified Amazon forest loss due to vegetation-atmosphere feedbacks Nat. Commun. 2017 8 14681 10.1038/ncomms14681 28287104
10. Hansen MC High-resolution global maps of 21st-century forest cover change Science 2013 342 850 854 10.1126/science.1244693 24233722
11. Chagnon FJF Bras RL Contemporary climate change in the Amazon Geophys. Res. Lett. 2005 32 L13703 10.1029/2005GL022722
12. Khanna J Medvigy D Fueglistaler S Walko R Regional dry-season climate changes due to three decades of Amazonian deforestation Nat. Clim. Change 2017 7 200 204 10.1038/nclimate3226
13. Garcia-Carreras L Parker DJ How does local tropical deforestation affect rainfall? Geophys. Res. Lett. 2011 38 L19802 10.1029/2011GL049099
14. Leite-Filho AT Soares-Filho BS Davis JL Abrahão GM Börner J Deforestation reduces rainfall and agricultural revenues in the Brazilian Amazon Nat. Commun. 2021 12 2591 10.1038/s41467-021-22840-7 33972530
15. McAlpine, C. A. et al. Forest loss and Borneo’s climate. Environ. Res. Lett. 13, 044009 (2018).
16. Chapman, S. et al. Compounding impact of deforestation on Borneo’s climate during El Niño events. Environ. Res. Lett. 15, 084006 (2020).
17. Spracklen DV Garcia-Carreras L The impact of Amazonian deforestation on Amazon basin rainfall Geophys. Res. Lett. 2015 42 9546 9552 10.1002/2015GL066063
18. Jiang Y Modeled response of South American climate to three decades of deforestation J. Clim. 2021 34 2189 2203 10.1175/JCLI-D-20-0380.1
19. Harris I Osborn TJ Jones P Lister D Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset Sci. Data 2020 7 109 10.1038/s41597-020-0453-3 32246091
20. Fassoni-Andrade AC Amazon hydrology from space: scientific advances and future challenges Rev. Geophys. 2021 59 e2020RG000728 10.1029/2020RG000728
21. Haiden, T., Janousek, M., Vitart, F., Ferranti, L. & Prates, F. Evaluation of ECMWF Forecasts, Including the 2019 Upgrade. ECMWF Technical Memorandum No. 853 (ECMWF, 2019).
22. Esquivel-Muelbert A Compositional response of Amazon forests to climate change Glob. Change Biol. 2019 25 39 56 10.1111/gcb.14413
23. Brum M ENSO effects on the transpiration of eastern Amazon trees Philos. Trans. R. Soc. B 2018 373 20180085 10.1098/rstb.2018.0085
24. Bagley JE Desai AR Harding KJ Snyder PK Foley JA Drought and deforestation: has land cover change influenced recent precipitation extremes in the Amazon? J. Clim. 2014 27 345 361 10.1175/JCLI-D-12-00369.1
25. Wunderling N Recurrent droughts increase risk of cascading tipping events by outpacing adaptive capacities in the Amazon rainforest Proc. Natl Acad. Sci. USA 2022 119 e2120777119 10.1073/pnas.2120777119 35917341
26. Fu R Li W The influence of the land surface on the transition from dry to wet season in Amazonia Theor. Appl. Climatol. 2004 78 97 110 10.1007/s00704-004-0046-7
27. Leite-Filho AT de Sousa Pontes VY Costa MH Effects of deforestation on the onset of the rainy season and the duration of dry spells in southern Amazonia J. Geophys. Res. Atmos. 2019 124 5268 5281 10.1029/2018JD029537
28. Negri AJ Adler RF Xu L Surratt J The Impact of Amazonian deforestation on dry season rainfall J. Clim. 2004 17 1306 1319 10.1175/1520-0442(2004)017<1306:TIOADO>2.0.CO;2
29. Chagnon FJF Bras RL Wang J Climatic shift in patterns of shallow clouds over the Amazon Geophys. Res. Lett. 2004 31 L24212 10.1029/2004GL021188
30. Chambers JQ Artaxo P Biosphere–atmosphere interactions: deforestation size influences rainfall Nat. Clim. Change 2017 7 175 176 10.1038/nclimate3238
31. Baudena M Tuinenburg OA Ferdinand PA Staal A Effects of land-use change in the Amazon on precipitation are likely underestimated Glob. Change Biol. 2021 27 5580 5587 10.1111/gcb.15810
32. Duku, C. & Hein, L. The impact of deforestation on rainfall in Africa: a data-driven assessment. Environ. Res. Lett. 16, 064044 (2021).
33. Akkermans T Thiery W Van Lipzig NPM The regional climate impact of a realistic future deforestation scenario in the Congo basin J. Clim. 2014 27 2714 2734 10.1175/JCLI-D-13-00361.1
34. Staal, A. et al. Feedback between drought and deforestation in the Amazon. Environ. Res. Lett. 15, 044024 (2020).
35. Xu, X. et al. Deforestation triggering irreversible transition in Amazon hydrological cycle. Environ. Res. Lett. 17, 034037 (2022).
36. Kooperman GJ Forest response to rising CO2 drives zonally asymmetric rainfall change over tropical land Nat. Clim. Change 2018 8 434 440 10.1038/s41558-018-0144-7
37. Chen, Z. et al. Global land monsoon precipitation changes in CMIP6 projections. Geophys. Res. Lett. 47, e2019GL086902 (2020).
38. Stickler CM Dependence of hydropower energy generation on forests in the Amazon Basin at local and regional scales Proc. Natl Acad. Sci. USA 2013 110 9601 9606 10.1073/pnas.1215331110 23671098
39. Challinor AJ A meta-analysis of crop yield under climate change and adaptation Nat. Clim. Change 2014 4 287 291 10.1038/nclimate2153
40. Strand J Spatially explicit valuation of the Brazilian Amazon forest’s ecosystem services Nat. Sustain. 2018 1 657 664 10.1038/s41893-018-0175-0
41. Potapov P Global maps of cropland extent and change show accelerated cropland expansion in the twenty-first century Nat. Food 2022 3 19 28 10.1038/s43016-021-00429-z 37118483
42. Li Y Deforestation-induced climate change reduces carbon storage in remaining tropical forests Nat. Commun. 2022 13 1964 10.1038/s41467-022-29601-0 35413947
43. Aragão LEOC Interactions between rainfall, deforestation and fires during recent years in the Brazilian Amazonia Philos. Trans. R. Soc. B 2008 363 1779 1785 10.1098/rstb.2007.0026
44. Marengo, J. A. et al. Changes in climate and land use over the Amazon region: current and future variability and trends. Front. Earth Sci.10.3389/feart.2018.00228 (2018).
45. Jiang, Y. et al. Widespread increase of boreal summer dry season length over the Congo rainforest. Nat. Clim. Change10.1038/s41558-019-0512-y (2019).
46. Van Der Ent RJ Savenije HHG Length and time scales of atmospheric moisture recycling Atmos. Chem. Phys. 2011 11 1853 1863 10.5194/acp-11-1853-2011
47. Sorí R Nieto R Vicente-Serrano SM Drumond A Gimeno L A Lagrangian perspective of the hydrological cycle in the Congo River basin Earth Syst. Dyn. 2017 8 653 675 10.5194/esd-8-653-2017
48. van der Ent RJ Savenije HHG Schaefli B Steele-Dunne SC Origin and fate of atmospheric moisture over continents Water Resour. Res. 2010 46 W09525
49. Feng Y Doubling of annual forest carbon loss over the tropics during the early twenty-first century Nat. Sustain. 2022 4 441 451
50. Tuinenburg, O. A., Bosmans, J. H. C. & Staal, A. The global potential of forest restoration for drought mitigation. Environ. Res. Lett. 17, 034045 (2022).
51. Met Office. Cartopy: a cartographic python library with a Matplotlib interface 2010–2015. Met Office https://scitools.org.uk/cartopy (2022).
52. Hoyer, S. & Hamman, J. xarray: N-D labeled arrays and datasets in Python. J. Open Res. Softw.10.5334/jors.148 (2017).
53. Zhuang, J. xESMF. Zenodo10.5281/zenodo.1134365 (2022).
54. Baker, J. C. A. & Spracklen, D. V. Climate benefits of intact Amazon forests and the biophysical consequences of disturbance. Front. For. Glob. Change10.3389/ffgc.2019.00047 (2019).
55. Schaaf, C. & Wang, Z. MCD43A3 MODIS/Terra+Aqua BRDF/Albedo Daily L3 Global - 500m V006. NASA EOSDIS Land Processes DAAC10.5067/modis/mcd43a3.006 (2015).
56. Waskom M Seaborn: statistical data visualization J. Open Source Softw. 2021 6 3021 10.21105/joss.03021
57. Chen M Global land use for 2015–2100 at 0.05° resolution under diverse socioeconomic and climate scenarios Sci. Data 2020 7 320 10.1038/s41597-020-00669-x 33009403
58. Funk C The climate hazards infrared precipitation with stations—a new environmental record for monitoring extremes Sci. Data 2015 2 150066 10.1038/sdata.2015.66 26646728
59. Xie, P. et al. NOAA Climate Data Record (CDR) of CPC Morphing technique (CMORPH) high resolution global precipitation estimates, version 1. NOAA National Centers for Environmental Information10.25921/w9va-q159 (2019).
60. Xie P A gauge-based analysis of daily precipitation over East Asia J. Hydrometeorol. 2007 8 607 626 10.1175/JHM583.1
61. Hersbach H The ERA5 global reanalysis Q. J. R. Meteorol. Soc. 2020 146 1999 2049 10.1002/qj.3803
62. Elke, R., Hänsel, S., Finger, P., Schneider, U. & Ziese, M. GPCC Climatology Version 2022 at 0.25°: monthly land-surface precipitation climatology for every month and the total year from rain-gauges built on GTS-based and historical data. GPCC10.5676/DWD_GPCC/CLIM_M_V2022_025 (2022).
63. Huffman, G. J. A., Behrangi, R. F., Adler, D. T., Bolvin, E. J. & Nelkin, G. G. Introduction to the new version 3 GPCP monthly global precipitation analysis. GPCP https://docserver.gesdisc.eosdis.nasa.gov/public/project/MEaSUREs/GPCP/Release_Notes.GPCPV3.2.pdf (2022).
64. Hou AY The global precipitation measurement mission Bull. Am. Meteorol. Soc. 2014 95 701 722 10.1175/BAMS-D-13-00164.1
65. Kobayashi S The JRA-55 reanalysis: general specifications and basic characteristics J. Meteorol. Soc. Japan 2015 93 5 48 10.2151/jmsj.2015-001
66. Gelaro R The modern-era retrospective analysis for research and applications, version 2 (MERRA-2) J. Clim. 2017 30 5419 5454 10.1175/JCLI-D-16-0758.1
67. Chen M Xie P Janowiak JE Global land precipitation: a 50-yr monthly analysis based on gauge observations J. Hydrometeorol. 2002 3 249 266 10.1175/1525-7541(2002)003<0249:GLPAYM>2.0.CO;2
68. Nguyen P The CHRS data portal, an easily accessible public repository for PERSIANN global satellite precipitation data Sci. Data 2019 6 1180296 10.1038/sdata.2018.296
69. Ashouri H PERSIANN-CDR: daily precipitation climate data record from multisatellite observations for hydrological and climate studies Bull. Am. Meteorol. Soc. 2015 96 69 83 10.1175/BAMS-D-13-00068.1
70. Nguyen P Persiann dynamic infrared–rain rate (PDIR-now): a near-real-time, quasi-global satellite precipitation dataset J. Hydrometeorol. 2020 21 2893 2906 10.1175/JHM-D-20-0177.1 34158807
71. Sadeghi M PERSIANN-CCS-CDR, a 3-hourly 0.04° global precipitation climate data record for heavy precipitation studies Sci. Data 2021 8 157 10.1038/s41597-021-00940-9 34162874
72. Huffman GJ The TRMM Multisatellite Precipitation Analysis (TMPA): quasi-global, multiyear, combined-sensor precipitation estimates at fine scales J. Hydrometeorol. 2007 8 38 55 10.1175/JHM560.1
73. Matsuura, K. & Willmott, C. J. Terrestrial precipitation: 1900-2017 gridded monthly time series. Global Precipitation Archive http://climate.geog.udel.edu/~climate/html_pages/Global2017/README.GlobalTsP2017.html (2018).
