Shifts in regional water availability due to global tree restoration

of large-scale global tree restoration on water availability by accounting for both the local, direct Tree restoration is an effective way to store atmospheric carbon and mitigate climate change. However, large-scale tree-cover expansion has long been known to increase evaporation, leading to reduced local water availability and streamflow. More recent studies suggest that increased precipitation, through enhanced atmospheric moisture recycling, can offset this effect. Here we calculate how 900 million hectares of global tree restoration would impact evaporation and precipitation using an ensemble of data-driven Budyko models and the UTrack moisture recycling dataset. We show that the combined effects of directly enhanced evaporation and indirectly enhanced precipitation create complex patterns of shifting water availability. Large-scale tree-cover expansion can increase water availability by up to 6% in some regions, while decreasing it by up to 38% in others. There is a divergent impact on large river basins: some rivers could lose 6% of their streamflow due to enhanced evaporation, while for other rivers, the greater evaporation is counterbalanced by more moisture recycling. Several so-called hot spots for forest restoration could lose water, including regions that are already facing water scarcity today. Tree restoration significantly shifts terrestrial water fluxes, and we emphasize that future tree-restoration strategies should consider these hydrological effects.

I n June 2021, the United Nations declared the Decade on Ecosystem Restoration to prevent, halt and reverse the degradation of ecosystems worldwide. Large-scale tree restoration is key in climate change mitigation and for enhancing and protecting biodiversity and ecosystem services 1 . Under the current climate conditions, it is estimated that an additional 900 million hectares of tree cover could exist on Earth's surface 2 without encroaching on agriculture and urban areas. During the past decade, numerous global and regional initiatives were initiated to increase tree cover, and this will play an important role in shaping global land use over the next decades. Despite these ongoing initiatives and the claims that ecosystem restoration is beneficial to all of the Sustainable Development Goals 3 , the impact of tree planting on the water cycle and water availability is still poorly understood 4,5 . As a result, potential impacts of ecosystem restoration on ensuring water availability both downstream and downwind are often overlooked.
Tree-cover expansion impacts water availability locally through its effects on the radiation balance, infiltration and soil water storage, evaporation, streamflow and precipitation 6 . Traditionally, local impacts of forest cover on streamflow have been investigated mainly using a so-called paired catchment approach. These studies compare two nearby headwater catchments with similar characteristics over a prolonged period, during which one of the catchments underwent land-cover change while the other did not undergo change. These observational studies have, virtually without exception, concluded that tree planting increases annual evaporation and decreases streamflow [7][8][9][10][11][12] . This high evaporation is attributed to the deeper roots of trees (facilitating access to water during dry periods), higher leaf area index (increasing the precipitation interception and canopy conductance), lower snow-free albedo (increasing the energy available for evaporation) and higher aerodynamic roughness (facilitating turbulent exchange) compared with the other vegetation types 9 . Higher evaporation has been reported across different climate zones and tree species, but the magnitude of evaporation differs with climate, tree species and tree age 7,8 . From these studies, it was predicted that large-scale tree restoration will decrease annual mean water availability and streamflow locally 9,[13][14][15] .
In contrast to these small-scale river-basin studies, more recent, large-scale research suggests that the impacts of tree restoration on streamflow are more complex 4,6,[16][17][18][19] . Through atmospheric feedbacks and transport, the increased evaporation from restored trees will partly recycle back to the terrestrial surface (via so-called evaporation or moisture recycling) and thereby potentially increase downwind rainfall and water availability. Such effects of tree-cover change can reach far beyond the river basin or even continental level: tree-cover change in the Amazon forest could impact precipitation in Canada, Northern Europe and all the way into Eastern Asia 20 . A host of regional and global-scale research has integrated the effects of evaporation recycling in land-cover change studies 4,19,21,22 . These studies have shown that evaporation recycling has a major influence on the water availability and that evaporation recycling should be considered in future land-cover change studies.
No study has quantified the effects of large-scale global tree restoration on water availability by accounting for both the local, direct effect of increased evaporation and the large-scale, indirect effect of evaporation recycling. The recently published datasets of the 'global tree-restoration potential' 2 and evaporation recycling 23 open up an emerging opportunity for such analyses. In our idealized study, we calculate how large-scale tree restoration (defined as increasing tree cover in any region, independent of the land-use history) influences water availability (defined as precipitation water that is not lost through evaporation, the total water available for consumption, on a yearly basis) and streamflow (the amount of water flowing in a stream; in this study, the accumulation of the as-defined water availability on the river-basin scale). More precisely, we calculate how a recent estimate of the global tree-restoration potential 2 (Fig. 1a) would impact the fluxes of evaporation, precipitation and streamflow. The global tree-restoration potential dataset highlights where more trees could naturally grow without encroaching on agricultural and urban land. To determine how tree-restoration impacts the long-term partitioning of precipitation between evaporation and streamflow, we use an ensemble of six data-driven Budyko models available in the literature. These six models all include a vegetation parameter that was calibrated separately for forest and non-forest conditions at a 1 km spatial resolution 11,24-27 (Extended Data Fig. 1 and Extended Data Table 1 and 2). We validated the   . c-f, The absolute annual change in water fluxes after tree restoration: change in evaporation (c), precipitation (d), water availability without evaporation recycling (e) and water availability with evaporation recycling (f). Note that e is the inverse of c: without the feedback of evaporation recycling, the local increase in evaporation equals the local decrease in water availability. g, The histogram shows the distribution of the global changes in water availability without and with evaporation recycling; 89% (without recycling) and 91% (with recycling) of the data fall within the displayed range of -20 mm yr -1 to +10 mm yr -1 . All maps display the 0.1° mean values, except for c, which displays the 0.5° mean value.
multimodel streamflow under current climate and forest-cover conditions against large river-basin run-off (Fig. 1b). We found good agreement over orders of magnitude, except for the Murray-Darling and Colorado basins, where water extraction probably causes the models to overestimate observed run-off. In addition to the Budyko models, we use the recent UTrack dataset of global evaporation recycling 23,28 to calculate where, and to what extent, the increased evaporation could increase (downwind) precipitation. This dataset is created using a state-of-the-art Lagrangian moisture-tracking model 29 and presents the monthly climatological mean evaporation recycling (Extended Data Fig. 2). We assume that tree restoration would intensify the current evaporation-recycling patterns as presented in the UTrack dataset. This approach focuses on the regional distribution of evaporated water but does not consider the effect that land-cover change has on local precipitation or atmospheric circulation and recycling patterns. These assumptions are further addressed in the discussion.

Impact on water fluxes
Under the current precipitation and potential evaporation, large-scale tree restoration would lead to a direct local increase in terrestrial evaporation of on average 8.2 ± 5.5 mm yr -1 , which is an increase of 1.2% (Fig. 1c, Extended Data Fig. 3). The standard deviation characterizes the variation over the six Budyko models. The sign of the average evaporation effect is clear, but the underlying distribution is heavily skewed, and the direct increase in evaporation can exceed 250 mm yr -1 in the South American tropical forest. When we include the indirect effects of increased evaporation on the atmospheric moisture content and precipitation, we find that 68% of this extra evaporated water would rain out over land (Fig. 1d). Following tree restoration, terrestrial precipitation would increase, with 4.8 ± 3.1 mm yr -1 on average due to recycling of increased evaporation (an average increase of 0.7%). While the patterns in the direct evaporation effect are determined by the patterns in tree restoration and local climate, the indirect evaporation-recycling effect of tree restoration on precipitation happens at the continental scale and is determined by atmospheric circulation. Without the effect of evaporation recycling, global mean water availability would decrease by 8.2 ± 5.3 mm yr -1 (Fig. 1e). When we include the process of evaporation recycling on the water budget calculations (thus taking into account that the evaporated water will partly rain out over land), water availability would decrease by 5.3 ± 5.6 mm yr -1 (Fig. 1f). There is, however, a large spatial variability, with regions showing a decreasing (up to 38%) or increasing (up to 6%) water availability (Extended Data Fig. 3). The United Kingdom, for example, has a high tree-restoration potential and therefore a high increase in evaporation. There is, however, low evaporation recycling due to the dominant westerly moisture transport from the country towards Eurasia and due to a low tree-restoration potential upwind in the Atlantic Ocean. Therefore, streamflow in the United Kingdom will decrease more substantially compared with the other regions. The low latitudes and the Tibetan Plateau, however, show an increase in water availability. In these regions, local evaporation recycling is high, for example, due to strong convection above the tropical forest or strong orographic lift in the mountainous regions, and travel distances of the atmospheric moisture are relatively short 28,32 . The results suggest that for half of Earth's surface (47%), the indirect moisture recycling effects of large-scale tree restoration could offset the direct evaporation effects, thus resulting in slight increases in water availability rather than decreases (Fig. 1g).

Varying effects on streamflow
To evaluate the impacts of tree restoration at the river-basin scale, we aggregate the direct effects (via increased evaporation) and indirect effects (through increased precipitation) for 21 large river basins. For all the river basins, enhanced evaporation reduces streamflow (up to 9%), but streamflow reduction can be close to zero, when increasing precipitation outweighs the increasing evaporation ( Fig. 2 and Extended Data Table 3).  The increased precipitation is a combination of local recycling (when evaporated water rains out within the same river basin) and recycling of water from other, further upwind, regions. For several tropical basins with a high local evaporation recycling, the high recycling nearly completely compensates for the loss of water through evaporation. Most river basins with a low wetness index (water-limited basins) have a low tree-restoration potential simply because arid regions can only support tree growth if there is significant groundwater convergence 33 . These arid river basins have a small absolute change in evaporation and precipitation. Overall, following tree restoration, streamflow will decrease for most of the world's important river basins despite the indirect effect of evaporation recycling.

Implications for water availability
Our study shows that large-scale tree restoration will shift water availability regionally. The hot spots for forest restoration could face a strong reduction in the water availability (Fig. 3a) because water supply by evaporation recycling is not enough to compensate for water loss through evaporation. However, regions with no or limited tree-restoration potential could see an increase in water availability due to the lateral transport of evaporated water into the region. This shows the importance of proper strategic planning of forest restoration projects 35 . Forest restoration could locally enhance precipitation through land-atmosphere interactions 18 , but these effects are complex and not included in our calculations. This enhanced precipitation could partly offset the high water loss in regions with high tree restoration.
Another important finding is that this tree-restoration scenario could further decrease annual water availability in several regions that are already facing water scarcity (Fig. 3b). Rainfall supply for some of these regions is largely of oceanic origin 28,36 , and these regions, therefore, do not, or only to a small extent, see an increase in precipitation after restoration. However, this tree-restoration scenario has the potential to increase water availability in other currently water-scarce regions. Although several regions could see a decrease in annually aggregated water availability, local water availability in the dry season could increase. Trees are known to increase soil porosity and soil organic carbon 37 . Therefore, tree restoration, especially in treeless regions, promotes the infiltration capacity and the soil water storage capacity and reduces the overland flow 6,38-40 . After tree restoration, the soils could thus potentially store more water, which would slowly be released during the dry season 35 . The true impact of tree restoration on water availability in the water-scarce regions may thus not be necessarily negative, and the impact depends on a combination of various aspects with high seasonal variability.

Discussion and conclusions
The presented results are calculated using state-of-the-art data and methods, but several feedbacks and processes are not included in the analyses. Our approach does not consider the potential impacts of increased tree cover on atmospheric circulation and land-atmosphere interactions. Tree restoration could locally enhance convergence, cloud cover and precipitation and change the travelling direction and distance of atmospheric moisture 18,[41][42][43] . Research suggests that forests could even impact large-scale wind patterns and draw atmospheric moisture from the oceans to the continents 44,45 , although the importance of this effect is still debated. These different feedbacks are poorly understood and difficult to incorporate in the present study because most evaporation-recycling models rely on meteorological reanalysis data, which are valid only under current land-cover conditions. The land-atmosphere interactions and potential changes in atmospheric circulation could impact the exact location and magnitude of increased precipitation. These processes could partly, but not fully, compensate for the enhanced evaporation, and this is unlikely to affect the main results of this study. In addition to data-driven studies, coupled land surface models have been used to reconstruct the effects of land-cover change on water availability 46,47 . Because of the large number of model parameters and feedbacks, their estimates of local impacts may be less constrained by observations and therefore more uncertain.
Furthermore, this study presents the impact of tree restoration under current climate conditions. However, global warming, and the tree restoration itself, will shift temperature and precipitation patterns 48,49 , and these are not considered in our analyses. Higher temperatures could reduce the global tree-restoration potential by 25% towards 2050 2 . Furthermore, higher temperatures will increase the potential evaporation and annual precipitation 50 , which causes a vertical shift in the Budyko framework (Extended Data Fig. 1b). Climate warming will impact the atmosphere's dynamics, which will possibly result in a poleward shift of the Hadley cells and storm  Fig. 3 | Implications of the percentage restored area and the implications for water scarcity. a, Water availability could decrease in 'hot spots for tree restoration'. The lower, middle and upper boundaries of the boxes display the 25th, 50th and 75th percentiles of the data (the interquartile range), respectively. The lines extend to a maximum of 1.5 times the interquartile range. b, Water availability will decrease in several regions that currently face water scarcity: the colour indicates whether local water availability will decrease (orange) or increase (blue) following the studied tree-restoration scenario. The cross-hatched regions currently face freshwater scarcity for at least three months per year, and the bold hatched regions face freshwater scarcity for at least nine months per year 34 . Note that we quantified the change in water availability on an annual timescale. The sign of change could vary seasonally, and water availability could, for example, decrease on a yearly basis but remain equal, or even increase, in the dry season.
tracks and a shift in magnitude of the jet streams [51][52][53][54] . The shifting global circulation will impact the patterns and magnitude of evaporation recycling. This study estimated the effects of global large-scale tree restoration on water fluxes, including evaporation recycling. The results show that tree restoration shifts water availability and streamflow at both local and continental scales and that the patterns of increasing and decreasing water availability are complex throughout the globe. Following tree restoration, the increased precipitation could increase water availability locally; however, we find a global net loss of water availability over the continents. These conclusions are consistent with other data-driven or model-based studies of the effects of land-cover change on water fluxes 18,46 . We stress that future tree-restoration strategies should consider the hydrological effects. Smart tree-restoration strategies are necessary to secure the water supply locally and remotely or enhance the water supply in water-scarce regions.

online content
Any methods, additional references, Nature Research 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 https://doi.org/10.1038/ s41561-022-00935-0.

Overview.
We calculate how a realistic scenario of large-scale tree restoration (restoring 900 million hectares of tree cover 2 ) would impact the fluxes of evaporation (E), precipitation (P) and streamflow (Q). The methodology includes four steps (Extended Data Table 1). (1) We use Budyko models to calculate E and Q before restoration. (2) We use Budyko models to calculate E and Q after restoration without evaporation recycling. (3) We use the UTrack dataset to calculate P after restoration on the basis of the increase in E calculated from steps 1 and 2. (4) We use Budyko models to calculate E and Q after restoration with evaporation recycling. The calculation of E and Q without evaporation recycling (step 2) is a local approach that assumes that evaporation increases with tree restoration but does not consider the effects of evaporation recycling. The calculation of E and Q with evaporation recycling (step 4) takes into account that the evaporated water partly recycles back to Earth's land surface and increases P.
Input data to calculate E and Q. The values of E and Q before tree restoration (step 1) are calculated using climatological P (Multi-Source Weighted-Ensemble Precipitation (MSWEP) v.2.8 55 ), climatological potential evaporation (WorldClim v.2 56 ) and current tree cover (Hansen tree cover v.1.7 57 ). For E and Q after restoration (steps 2 and 4), we use tree cover after restoration (Hansen tree cover 57 + tree-restoration potential 2 , where we set the maximum tree cover after restoration to 1). The values of E and Q with evaporation recycling (step 4) are calculated using P after restoration (MSWEP P + P from evaporation recycling calculated in step 3). We assume that tree restoration does not affect potential evaporation. All input data were resampled to 0.00833° spatial resolution (~1 km at the Equator).
Budyko models to calculate E and Q. The values of E and Q before and after restoration are calculated using six different Budyko-type models 11,[24][25][26][27] (Extended  Data Table 2). These models reflect global annual patterns in P partitioning between E and Q (Extended Data Fig. 1). The general form of the equations is: where E P is the fraction of precipitation partitioned into evaporation, E0 P is the aridity index and ω is a model parameter. Previous studies have suggested that the parameter ω is closely related to vegetation type or vegetation coverage 11 . Six models were available in the literature that have a calibrated ω parameter for forest and grassland cover. These models were calibrated using river-basin or lysimeter data from different climate zones. The parameter ω trees was calibrated with data from both natural forests (the major part of the data) and plantations. The vegetation parameter for grassland, ω grass , was calibrated using grassland data only or using both grassland and cropland data (Extended Data Table 3). The semi-empirical Budyko models reflect the yearly, catchment-integrated, effect of differences in interception, the plant available water, evaporation, the soil water storage capacity and the energy balance between grassland and forest.
For each of the six Budyko models, E and Q before and after restoration (without and with evaporation recycling) are calculated as: where E trees and Q trees are E and Q calculated for trees, TC is the tree cover, E grass and Q grass are E and Q calculated for grassland and NTC is the no tree cover, where: NTC, thus, includes, among others, bare land, shrubland and built-up area. Evaporation from these land-cover types is different from evaporation from grassland and cropland. We expect that these differences are negligible because bare land is generally very arid, where E is close to zero (both the Budyko E and true E), and built-up area has a limited extent. The Budyko vegetation parameters for shrubland (when calculated) closely resemble the parameters for grassland and cropland 25,26 .
Mean yearly Q before restoration was validated against streamflow data for 19 river basins (Fig. 1b). These 19 river basins are of the studied basins for which streamflow measurements and the river-basin outline were available in the Global Runoff Data Centre (GRDC) dataset 30 .
The changes in E and Q were calculated from the differences in E and Q before and after restoration. The results were resampled to 0.5° spatial resolution for the calculation of evaporation recycling.
Evaporation recycling-calculation of precipitation. We used the UTtrack dataset 23,28 of global atmospheric moisture recycling (at 0.5° spatial resolution) to calculate where the additional evaporated water from restored trees would rain out (P after restoration in step 3). The UTrack dataset was created using the state-of-the-art Lagrangian moisture-tracking model UTrack 29 . The dataset presents the monthly climatological mean evaporation recycling for each pixel. We aggregated the dataset's monthly mean evaporation recycling to yearly mean evaporation recycling to be able to combine it with the yearly Budyko models. This yearly aggregation of the UTrack dataset ignores the seasonality in direction and magnitude of atmospheric moisture transport that exists in some regions 18,28 . This aggregation was preferred above a disaggregation of the Budyko results because the Budyko model assumes zero change in water storage and can therefore be applied only on a multi-year mean timescale. Restoring tree cover could increase soil storage capacity and increase dry-season water availability 38 . Thus the Budyko models cannot be used to calculate evaporation at a monthly timescale. As a comparison, evaporation recycling was also calculated using a monthly time step for the river basins of the Mississippi, Orange and Ganges-Brahmaputra (Extended Data Fig. 2). When using the monthly time step, the total yearly evaporation is disaggregated over 12 months relative to the magnitude of evaporation for each month in the European Centre for Medium-Range Weather Forecasts Reanalysis v.5 data 58 . The monthly and yearly calculations have similar patterns and values, but there are differences of up to 1 mm yr -1 .
In our approach, we assume that increased evaporation will amplify the current evaporation-recycling patterns. It is commonly assumed that land-cover change impacts evaporation recycling (patterns), but the effects of land-cover change on evaporation-recycling processes are poorly understood and difficult to determine 41 . There is currently no possibility to include these processes at high resolution at global scale, and therefore we use the UTrack dataset as the best available option.