Remote land use impacts on river flows through atmospheric teleconnections

. The effects of land-use change on river ﬂows have usually been explained by changes within a river basin. How-ever, land–atmosphere feedback such as moisture recycling can link local land-use change to modiﬁcations of remote precipitation, with further knock-on effects on distant river ﬂows. Here, we look at river ﬂow changes caused by both land-use change and water use within the basin, as well as modiﬁcations of imported and exported atmospheric moisture. We show that in some of the world’s largest basins, precipitation was inﬂuenced more strongly by land-use change occurring outside than inside the basin. Moreover, river ﬂows in several non-transboundary


Introduction
River flows (Q) are fundamental for ecosystems, nutrient transport, hydropower, navigation, and human well-being (Oki and Kanae, 2006). Land-use change (LUC) has been suggested to be the most important driver of both past (Piao et al., 2007;Sterling et al., 2012) and future (Betts et al., 2015;Milly et al., 2005) changes in river flows ( Q). Central to the analysis of Q is the river basin unit, and estimates of Q from LUC often assume that impacts occur exclusively within a basin (Gerten et al., 2008;Piao et al., 2007;Rost et al., 2008a, b;Sterling et al., 2012). Water governance is strongly focused on frameworks such as the Integrated River Basin Management (IWRM) and largely assumes that there is no land-atmosphere feedback, even in discussions of spatial misfit between institutions and hydrological realities (Hoekstra, 2010;Giordano et al., 2015). In fact, landatmosphere feedbacks are not incorporated into most recent literature on a wide range of topics of relevance for water management, such as virtual water (Dalin et al., 2017), the freshwater planetary boundary (Rockström et al., 2009;Steffen et al., 2015), water scarcity (Mekonnen and Hoekstra, 2016), the relative role of climate and LUC for water flows (Zheng et al., 2016), and land acquisition impacts on water (Johansson et al., 2016;Rulli et al., 2012).
However, studies on land-atmosphere interactions clearly show that changes in land surface properties can considerably influence precipitation (P ) and Q through landatmosphere feedback, sometimes well beyond the local scale (Badger and Dirmeyer, 2016;Garcia et al., 2016;Avissar and Werth, 2005). For example, general circulation model simu-Published by Copernicus Publications on behalf of the European Geosciences Union.
lations suggest that complete deforestation of central Africa may decrease February P by 35 % in the Great Lakes region (Avissar and Werth, 2005), and irrigation in India may support up to 40 % of the P in some arid regions in eastern Africa (de Vrese et al., 2016). Under a business-as-usual deforestation scenario, Q in the Xingu River basin in the Amazon was found to increase by 10-12 % without landatmosphere feedback, and decrease by 30-36 % when such feedback was taken into account (Stickler et al., 2013). Furthermore, statistical analyses of observed data suggest that irrigation in the US High Plains enhances downwind Q (Kustu et al., 2011), and coupled regional climate modelling shows that irrigation in the California Central Valley can be linked to an about 30 % increase in Colorado Q (Lo and Famiglietti, 2013). At the global scale, Q from future climate and LUC scenarios changed from decrease to increase by considering land-atmosphere feedback and by closing the water balance (Betts et al., 2015).
Land-atmosphere interactions can influence Q through thermal layer processes, terrestrial moisture recycling (TMR), and circulation perturbation (Goessling and Reick, 2011). First, thermal layer processes refer to the boundary layer and mesoscale circulation perturbation that may lead to a change in total terrestrial evaporation (E) and can locally lead to both positive and negative P responses (Guillod et al., 2015;Seneviratne et al., 2010;Koster et al., 2003). Local forest clearing has for example been shown to enhance P in downwind areas due to turbulence changes (Khanna et al., 2017;Saad et al., 2010). Second, TMR refers to the process of terrestrial E returning to land as P and is underpinned by the mass conservation of water (Brubaker et al., 1993). TMR is often the dominating land-atmosphere process at the regional to continental scale (D'Almeida et al., 2007;Spracklen et al., 2012;Lawrence and Vandecar, 2014;Tuinenburg, 2013). About 40 % of global terrestrial P (van der  originates from terrestrial E and the average distance travelled in the atmosphere is 500-5000 km (van der Ent and Savenije, 2011) -a distance likely to exceed the size of most river basins. Lastly, large-scale atmospheric circulation perturbation allows extreme LUC (e.g. complete tropical deforestation) to impact P in geographically remote regions and continents in unexpected ways (Avissar and Werth, 2005;Badger and Dirmeyer, 2016;Garcia et al., 2016;Lawrence and Vandecar, 2014). Monsoon regions are particularly sensitive to circulation perturbation, and irrigation may for example reduce P by weakening the monsoon onset (Tuinenburg, 2013).
The previous studies that illustrated the importance of remote LUC for basin P and Q did not examine the effect of taking moisture recycling into account for estimating LUC effects on Q and attributing them to influence from different nations, nor did they analyse the interplay between LUC within and outside the river basin. These effects are, however, important to disentangle since they can have profound water governance implications for, for example, riparian wa-ter rights and transboundary river basin treaties (Keys et al., 2017;Dirmeyer et al., 2009;Ellison et al., 2017). Thus, there is a missing interdisciplinary bridge between understanding the role of land-atmosphere feedback over large distances and its importance for water governance at the basin scale.
This study aims to (i) investigate the potential impacts of human LUC on Q worldwide accounting for TMR, (ii) disentangle the relative influence on Q from within-and extrabasin LUC, (iii) attribute potential human LUC impacts on Q to nation states, and (iv) discuss the potential implications for water governance. We focus on the TMR effect because it is transparent, closes the water balance, and explicitly links changes in land and water geographically. Given these advantages, similar TMR approaches have in recent years been used to analyse unexplored relations, e.g. LUC impacts of crop yields (Bagley et al., 2012), self-amplifying forest dieback from TMR changes (Zemp et al., 2017), and vulnerability to LUC-induced reductions in P Miralles et al., 2016). For a comparison of different methods for analysing LUC impacts on Q, see Table S1 in the Supplement.

Hydrological modelling
We used the process-based Simple Terrestrial Evaporation to Atmosphere Model (STEAM) hydrological model  to simulate water fluxes based on land cover and land use. STEAM partitions evaporation into five fluxes: vegetation interception, floor interception, transpiration, soil moisture evaporation, and open-water evaporation. STEAM uses the Penman-Monteith equation (Monteith, 1965) to estimate potential evaporation, the Jarvis-Stewart equation (Stewart, 1988) to compute stomatal resistance, and Jolly's growing season index (function of minimum temperature, soil moisture content, and daylight) to describe phenology (Jolly et al., 2005). STEAM operates at 1.5 • × 1.5 • and a 3 h resolution. Based on the long-term water balance, mean annual river flow (Q) is assumed to approximately equal the difference between mean annual P and E, i.e. Q = P − E. STEAM was validated in previous studies ) and compared well with recent observation-based analyses of evaporation partitioning by land-cover type (Wei et al., 2017). Modifications from the original version of STEAM  include (1) update of land-use classification, parameterization, and parameterization approach (Table S2 and Fig. S1 in the Supplement), (2) use of a temperature threshold of 0 • C for differentiating snowfall from rainfall, and (3) differences in input data (i.e. root zone storage capacity, land surface map, precipitation data source as Hydrol. Earth Syst. Sci., 22, 4311-4328, 2018 www.hydrol-earth-syst-sci.net/22/4311/2018/ described in Data). Evaluation against runoff data is shown in Fig. S2. Simulated land-use change effects on evaporation increase and decrease are compared with literature values in Table S3 and found to be in the conservative range. With the study period being 2000-2013, the years 1995-1999 were used as spin-up for STEAM.

Moisture tracking
Atmospheric moisture is tracked using the Eulerian Water Accounting Model-2 layers (WAM-2layers) moisture tracking scheme . WAM-2layers tracks atmospheric moisture from zero pressure to surface pressure in two layers. Within the layers, atmosphere is assumed to be well mixed. WAM-2layers tracks vapour flows by applying the water balance. For example, the following equation is used to track where evaporation from a given region falls as precipitation (i.e. forward tracking): where S tracked is the tracked atmospheric storage in an atmospheric column in one layer, t is time, u and v are wind components in the x zonal and y meridional directions, E tracked is tracked evaporation entering and P tracked is precipitation exiting an atmospheric column and layer, and F vertical,tracked is the tracked vertical moisture transport between the two layers. An analogous equation is used for tracking the source of precipitation to a given region (i.e. backward tracking). The spatial resolution of WAM-2layers is 1.5 • and input data are linearly interpolated to the 15 min time step to maintain numerical stability. WAM-2layers has been employed previously for analysing atmospheric moisture transport over terrestrial areas (Keys et al., 2012 and validated against other types of moisture tracking algorithms (van der Ent et al., 2013). We used the MATLAB version of WAM-2layers, but a Python version is also openly available on Github (van der Ent, 2016). With the study period being 2000-2013, the year 1999 is used as spin-up in forward tracking in WAM-2layers, and 2014 is used as spin-up for backward tracking in WAM-2layers.

Coupling of the moisture tracking scheme and the hydrological model
Hydrological flows in the current land-use scenarios are simply represented by current data and simulation. To obtain E and P under potential land cover, STEAM is coupled with WAM-2layers by (1) simulating present-day E in STEAM and forward tracking terrestrial E with WAM-2layers, meaning that the E tracked is equal to all evaporation from terrestrial surfaces, i.e. not belonging to the oceans, (2) simulating E in STEAM based on present-day P and potential land cover, and forward tracking the fate of terrestrial E with WAM-2layers, (3) calculating the change in P tracked , (4) updating the present-day P with the changes in P tracked , and (5) simulating E in STEAM based on updated P and potential land cover, and forward tracking the fate of terrestrial E with WAM-2layers; see Fig. 1. Steps 3-5 are iterated until the annual P change is below 1 % and the monthly P change is below 5 mm month −1 in every grid cell, which in our case ultimately resulted in four iterations in total. This procedure assumes that land-use induced changes in terrestrial E will result in proportional changes in P with terrestrial origin.

Land data
Land-use and land-cover data input to STEAM are based on the Ramankutty potential land-cover (Ramankutty and Foley, 1999) and current land-use scenarios (Ramankutty et al., 2008) for consistency. We further added permanent wetlands, permanent snow or ice, and urban or built-up areas from the Land Cover Type Climate Modeling Grid (CMG) MCD12C1 International Geosphere Biosphere Program (IGBP) land classification created from Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) data (Friedl et al., 2010) for the year 2005. Monthly irrigated rice and irrigation non-rice crops were obtained from the data set of Monthly Irrigated and Rainfed Crop Areas around the year 2000 (MIRCA2000) V1.1 (Portmann et al., 2010). The urban and irrigated areas were only added to the current land-cover map. In this merging procedure, MODIS is allowed to override the Ramankutty data sets, and MIRCA2000 is allowed to override the Ramankutty map as long as it does not extend over the cropland areas. The scenarios used are shown in Fig. S3 and the land-use change is illustrated in Fig. 2. The root zone storage capacity map is based on a climateobservation-based root zone storage capacity (S R ) (Wang-Erlandsson et al., 2016) derived from satellite and energy balance-based evaporation, gauge-based precipitation, and modelled irrigation. The best performing Gumbel normalized root zone storage capacity (S R,CRU-SM,merged ) was used. Root zone storage capacity for both current and potential land-cover and land-use scenarios was constructed from the mean of land-cover type and Köppen-Geiger climate class (Kottek et al., 2006). The mean root zone storage capacity of single land-cover types was used only in places where the combination of land-cover type and climate zone that exists in the potential land-cover scenario did not exist in the current land-use map.

Meteorological forcing and runoff data
Meteorological data used in WAM-2layers and STEAM, except for land precipitation, were taken from the Earth Retrospective Analysis Interim (ERA-I) from the European Cen- Figure 1. Model coupling schematic. Model coupling between STEAM and WAM-2layers based on current land use and potential vegetation scenarios. P stands for current precipitation; E stands for evaporation. Subscript t stands for terrestrial origin, pv denotes simulation with potential vegetation, cur denotes simulation with current land use, and n stands for the number of iterations.
tre for Medium-Range Weather Forecasts (ECMWF) (Dee et al., 2011). ERA-I meteorological forcings to STEAM are snowmelt, temperature at 2 m height, dew point temperature at 2 m height, wind speed (meridional and zonal vectors) at 10 m height, incoming shortwave radiation, and net longwave radiation. In addition, ERA-I evaporation data were used to downscale calculated daily potential evaporation in STEAM to the 3 h time step. ERA-I model level forcings used in WAM-2layers are specific humidity and wind speed at 6-hourly resolution, spanning from zero to surface pressure. Moreover, 3-hourly ocean evaporation is taken from ERA-I. The Modern-Era Retrospective analysis for Research and Applications (MERRA) reanalysis has in a previous study been used as input to WAM-2layers for comparison and generated similar persistent moisture recycling patterns, except in South America, where differences arise due to underestimation of precipitation in MERRA (Keys et al., 2014). Precipitation forcing for WAM-2layers and STEAM both come from the state-of-the-art Multi-Source Weighted-Ensemble Precipitation (MSWEP V1) product  that was specifically created for hydrological modelling. The use of MSWEP as forcing for STEAM resulted in runoff estimates that compare well to observed runoff data (Fig. S2). All meteorological forcing data cover temporally 1995-2014.
Runoff data used for benchmarking were taken from the composite (observed river discharge consistent with the water balance model) from the Global Runoff Data Centre (GRDC) (Fekete et al., 2002). The separate GRDC water bal-ance model runoff fields are included in the comparison for reference (Fig. S2).
The spatial coverage of all data used is 57 • S-79.5 • N latitudes at 1.5 • × 1.5 • resolution. MSWEP originally at 0.25 • and GRDC runoff at 0.5 • were aggregated to 1.5 • resolution by simple averaging.

Changes in hydrological flows
River flow change without TMR ( Q noTMR ) is where P cur is current-day precipitation data from MSWEP, E cur is current-day evaporation based on STEAM simulation, and E pv,1 results from STEAM simulation in the potential vegetation scenario and forced with current-day precipitation ( Fig. 1). River flow change after accounting for TMR ( Q) is where P pv,4 is the converged precipitation (i.e. meeting the convergence requirement of mean annual precipitation change < 1 % yr −1 and monthly precipitation change < 5 mm month −1 in every grid cell) achieved at the fourth iterative coupling between STEAM and WAM-2layers, and E pv,5 is the evaporation under the potential vegetation scenario simulated in STEAM with precipitation forcing P pv,4 . Change in tracked basin precipitation ( P tracked,basin ) occurring outside the river basin boundaries is referred to Changes in (a) land use (current land use is shown, with grid cells without major land-use change masked out) and (b) total area of difference between potential vegetation (y axis) and the current landuse map (colour legend).
as P import , whereas P tracked,basin originating from within the basin boundary is referred to as P basin-recycling . Internally recycled evaporation ( E basin-recycling ) corresponds to P basin-recycling and all other basin evaporation change is considered exported ( E export ).

Country influence on changes in river flows
The influence on river flow change in river basin b from country c without considering TMR (I b,c,noTMR ) is where E b,c is evaporation change in the part of river basin b located in country c. The influence on river flow change in basin b from country c with consideration of TMR (I b,c,TMR ) is where E b,c,export is the evaporation change exported from the part of basin b located in country c, and P b,c,import is the precipitation change imported to basin b from country c. Influences from countries below 5 % of total influences in a specific basin (I b,c,noTMR < 0.05 × I b,c,TMR ) were lumped into the category "Other".

LUC impacts on global water flows
Our results show that human LUC (from potential land cover to current land use) (Fig. 2) has led to reductions in E and P , and to increases in Q, in large parts of the world (Fig. 4bd). E has decreased primarily in Southwest China, Europe, western Africa, the south of Congo, and south-eastern South America, resulting from substantial pasture and agricultural expansion (Ramankutty et al., 2008). Following prevailing wind directions (Fig. 3c), subsequent P has decreased in all tropical regions, southern central China, the eastern US, and Europe.
Nevertheless, in some areas, E increased due to incremental irrigation -notably in India, the western US, Northeast China, and the Middle East (Fig. 4a, b). Due to the combination of heavy irrigation in India and orography, P has increased substantially along the Himalaya mountain ridge (Fig. 4b, c). Weak increases in P are observed in other downwind regions: the Sahel (i.e. downwind irrigation areas along the Nile) and in the western US. Continental precipitation recycling ratios are modified -with some exceptions -in a similar pattern to P (Fig. 4e, f). Large Q are seen in the La Plata basin in South America, the Zambezi in southern Africa, the Yangtze in China, and the Indus in northern India (Fig. 4g), and relative changes in Q are large in for example the Colorado basin in the US, the Odra basin in eastern Europe, and the Lake Chad river basin in Africa (Fig. 4g).

The role of TMR for Q
In aggregate (Fig. 5), when accounting for TMR, LUC changed global terrestrial E by −1251 km 3 yr −1 (−1.8 % from 69 211 km 3 yr −1 ), P by −586 km 3 yr −1 (−0.5 % from 107 800 km 3 yr −1 ), and Q by 664 km 3 yr −1 (+1.7 % from 38 589 km 3 yr −1 ). The estimated changes to Q tend to fall in the conservative end of previous estimates (Gerten et al., 2008;Piao et al., 2007;Rost et al., 2008a, b;Sterling et al., 2012) (Fig. 5). However, recent research  suggests that consumptive water use is severely underestimated in earlier studies (e.g. Döll et al., 2009;Sterling et al., 2012). Q with TMR corresponds to the difference between E and P change including TMR ( Fig. 5, solid bars), whereas Q without accounting for TMR simply corresponds to E without TMR (Fig. 5, hollow bars). Including TMR nearly halves the global Q estimate. This is because E returns as P over land and thus compensates for the initial water "loss" from the basin. This suggests that previous studies without TMR (e.g. Gerten et al., 2008;Piao et al., 2007;Sterling et al., 2012) may have substantially overestimated the net LUC impacts on Q. Our estimate of LUC impact on Q is slightly larger than some of the estimates of CO 2 fertilization (e.g. Alkama et al., 2010;Gerten et al., 2008), but substantially smaller than climate change and overall human impact (e.g. Alkama et al., 2010;Gerten et al., 2008) (Fig. 5).
Our river basin analysis shows that accounting for TMR considerably alters estimates of Q (Fig. 7a): in the Congo, Volga, and Ob basins, Q are reduced by more than half; in the Amazon, Q drops from 1630 to 270 m 3 s −1 ; and in the Yenisei, the sign of Q is reversed from an increase (150 m 3 s −1 ) to a decrease (−220 m 3 s −1 ).
At the basin level, the TMR effect on river flow change is estimated to be the largest in large and relatively wet basins such as the Amazon, Congo, and Yangtze River basins in terms of absolute volumes (Fig. 6a). Not accounting for TMR clearly generates the largest relative deviations in river flow change estimates in the Amazon (i.e. Q noTMR is approximately 5 times larger than Q), and large relative TMR effects are seen in many large basins worldwide, including e.g. the Congo ( Q noTMR is 150 % higher than Q), Yenisei ( Q noTMR is 165 % lower than Q), and Ob ( Q noTMR is 140 % higher than Q) river basins (Fig. 6b). The TMR effect relative Q cur (Fig. 6c) shows that TMR effect can be important also in more arid basins such as the Colorado, Niger, and the Yellow River.

The interplay between internal and external LUC
Furthermore, atmospheric moisture does not respect river basin boundaries (Fig. 7a, and spatial maps in Figs. 8, 9, S4, S5, S6, and S7). In fact, P over the basins has been modified more significantly by external than by internal LUC (change in imported precipitation P import > change in internally recycled precipitation P basin-recycling ) in some of the largest basins (Fig. 7a). Likewise, internally recycled evaporation changes (E basin-recycling ) (Fig. 7b, II) are substantially smaller than E, affecting P elsewhere ( E basin-recycling < change in exported evaporation E export ) for all selected river basins (Fig. 7a).
Internal moisture recycling (Fig. 7b, II) does not affect Q directly, but only indirectly if P basin-recycling affects subsequent E export under transient change ( Fig. 7b and Methods). Thus, provided steady state, Q simply corresponds to the difference between E export and P import (Fig. 7a). For example, Q in the Amazon is very small because the reduced P import is almost entirely offset by reduced E export . In Congo, about half of the within-basin LUC-induced Q increase is counteracted by extra-basin LUC (i.e. P import ≈ 0.5 E export ). The effect of TMR on Q ( Q noTMR − Q, where subscript noTMR denotes simulation without TMR) corresponds to total P (i.e. P import +  . Land-use change-induced changes in hydrological flows (current land-use-potential vegetation scenario): (a) absolute change in evaporation, (b) relative change in evaporation, (c) absolute changes in precipitation, (d) relative change in precipitation, (e) absolute change in continental precipitation recycling ratio (i.e. precipitation with terrestrial origin divided by total precipitation P tracked /P and converted to the unit of percent), (f) relative change in continental precipitation recycling ratio, (g) absolute change in river flows at outlet, and (h) relative change in river flows at outlet. P basin-recycling ) and any indirect E (i.e. E noTMR − E, not shown). In the Yangtze, the Q is mitigated mostly by P basin-recycling . The strong flow reduction in the heavily irrigated Indus, however, is only mildly compensated for by TMR (i.e. P import E export ). The pattern of overlapping precipitationsheds (i.e. P source regions)and evaporationsheds (i.e. E sink regions) illustrated in Fig. 7b and moderated by wind directions can also be clearly seen in the basin-specific precipitationshed and evaporationshed maps (Fig. 8). In the Amazon (Fig. 8a,   b), the moisture arrives from the east, is stopped up by the Andes, and changes direction towards the south-east. The hotspot of precipitation source and sink within the Amazon basin does not overlap, with major moisture providing spots located along the north-eastern border and the major moisture receiving spots located along the Andes in the west. In the Yangtze (Fig. 8c, d), the moisture comes from a large area in the south, and leaves in the direction of Japan in a relatively narrow band. In the Yenisei (Fig. 8e, f), the moisture follows the westerlies, coming in straight from the west, P07 S12 S12 R08b P07

R08a
R08a G08 G08 G08 G08 G08 A10 Figure (Table S4), where LUC implicitly accounts for consumptive water use. Note that while consumptive water use alone always reduces river flows, other human impacts have both positive and negative influences that are concealed by the global aggregate. and leaving straight towards the east. In the Niger (Fig. 8g, h), the moisture is mostly supplied from the east from terrestrial areas, and flows towards the west into the Atlantic. For precipitationsheds and evaporationsheds of other basins, see Figs. S4 and S6 respectively. While changes in precipitationsheds and evaporationsheds are conditioned by the original moisture flows, the resulting pattern ultimately depends on the distribution of LUCinduced hydrological change (compare Figs. 8 and 9). For example, although the Amazon precipitationshed is weak over Africa (Fig. 8a), the precipitationshed change is in fact relatively strong there due to strong LUC-induced hydrological change (Fig. 9a). In other cases, aggregated changes in Fig. 7 hide spatially heterogeneous increases and decreases in moisture flows. For example, agricultural activities and irrigation in India, the Sahel, and regions around the Nile increase moisture supply to the Yangtze, Yensisei, and Niger basins and offset deforestation-induced moisture supply decrease elsewhere (Fig. 9c, e, g). For changes in precipitationsheds and evaporationsheds of other basins, see Figs. S5 and S7 respectively.

Attributing influence on Q to nations
Typically, TMR attributes LUC influence on Q (methods described in Sect. 2.3.2) to a larger number of nations than when only basin boundaries are considered (Fig. 10). In several of the studied basins (such as the Amazon, Congo, Volga, Ob, Yenisei, and Niger basins; see Fig. 10a, b, e, f, i, l), the share of nations contributing less than 5 % to Q grows considerably when TMR is considered. In some cases, nations not considered key influencers of Q in fact influence Q by more than 5 % when TMR is accounted for: in the Mekong, India is only an important influencer (10 % influence) when TMR is considered (Fig. 10g); in the Yenisei, Mongolia falls below 5 % influence, while Kazakhstan (11 %) and China (6 %) climb considerably in influence (Fig. 10i); and in the Niger basin, Sudan/South Sudan (8 %) and Niger (5 %) replace Ivory Coast and Guinea as important influencers (Fig. 10l). Notably, basins geographically confined within one nation can be influenced by LUC taking place in foreign nations. This is for example the case in the Yangtze, Yellow, and Huai, where irrigation in India increases the basins' P (Fig. 10d, m, n). The TMR leads to a limited difference in nation influence only in the North American basins (Fig. 10h, o) and La Plata (Fig. 10c).

Interplay between TMR and LUC
At the global scale, Q as a response to LUC can be almost halved by taking TMR into account (Fig. 5)  pensate for the majority of the impact on Q from local LUC in other basins (e.g. Amazon, Fig. 7a) and even propose new transboundary relationships (e.g. Yangtze, Fig. 10d). From a TMR perspective, the impact on Q from within-basin LUC depends on the E exported from the basin as much as the P imported to the basin. Our analysis shows the importance of considering LUC on par with TMR to identify anthropogenic influence on water resources, beyond analyses of pure moisture exchanges (Dirmeyer et al., 2009;Keys et al., 2017). While Africa does not constitute a major moisture source of Amazonian P (7 % of all Amazon P , 13 % of Amazon P with continental origin; see also Fig. 8a), the spatial extent of E from LUC was sufficient to elevate the relative importance of African LUC for Amazonian Q (28 % of Amazon P ; see also Fig. 10a). Similarly, India is not identified as a major moisture source of the Yangtze (see Fig. 8c and Wei et al., 2012), but has about 10 % influence on Yangtze Q (Fig. 10d).

Potential governance relevance
Our results indicate that both precipitationsheds and evaporationsheds of river basins are relevant governance units. Previous studies of TMR for water management (Berger et al., 2014;Keys et al., 2017) have emphasized the importance of considering the P source region, i.e. the precipitationshed (Keys et al., 2012), which was introduced as a concept analogue to watershed for water resource management. This study finds that the evaporationshed (van der Ent and Savenije, 2013), i.e. the E sink region, is just as important when considering changes to Q.
LUC impacts Q through TMR in different ways depending on how precipitationshed, river basin, and evaporationshed are aligned. For example, where an evaporationshed has a limited overlap with river basin boundaries, reforesting a river basin may lead to unexpectedly large reductions in Q, if considerable deforestation simultaneously occurs in the precipitationshed outside the river basin.
The magnitude of TMR effects from remote LUC on Q can be comparable to managed water flows. For example, the Yangtze River provides 36 % of the country's surface water resources, and is subject to two of the world's most ambitious water engineering projects: the Three Gorges Dam and the South-to-North Water Diversion (CWRC, 2017). The overall TMR effect on mean annual LUC-induced Q is here estimated at 980 m 3 s −1 in the Yangtze basin, and the mean annual moisture change imported to the basin from foreign countries is estimated at 1110 m 3 s −1 (Fig. 9c). As a comparison, the Q difference between a normal year and a dry year is about 300 m 3 s −1 and the total amount of water to be transferred from the Yangtze through the South-to-North Water Diversion is aimed to be 1420 m 3 s −1 (NSBD, 2011). Seasonal and interannual flow variability is a major challenge facing the Yangtze, and future research in the seasonal LUC influence and interaction with the monsoon system is needed. Note, however, that our estimates are associated with parameter sensitivity (see Fig. S9) and large uncertainties as discussed in the Limitations.
We note that the relevance of considering TMR governance depends on future LUC. The simulated Q in this paper follows from a rather extreme LUC scenario (from potential to current land use). The current LUC in this study is 15 million km 2 cropland and 28 million km 2 pasture conversion (Ramankutty et al., 2008). As a comparison, models used in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4) estimated cropland changes from −1.2 to +12 million km 2 between 2000 and 2050 (IPCC, 2007). A more recent multi-model comparison ranged cropland conversion until 2050 from −1 to +8.5 million km 2 across different scenarios (Schmitz et al., 2014). In total, the potential land for agricultural conversion has been estimated at 17 million km 2 (Schmitz et al., 2014). Thus, future LUC can be considerable, and potential TMR impacts on Q will be dependent on the type and geographi-  Note also that the figure has two y axes, m 3 s −1 to the left and % to the right cal distribution of LUC, as well as dependent on prevailing winds, hydroclimate, and orography.

Limitations
In interpreting our results, it should be noted that our approach only accounts for the TMR effects. The frequency or intensity of P is assumed to remain unaffected by thermal layer processes or circulation perturbation, which may introduce a bias into the quantitative estimates of hydrological flows under water-limited conditions (i.e. semi-arid regions and temperate regions during summertime) (Medvigy et al., 2011). Furthermore, vegetation response to P is not simulated, such as forest dieback from increased fire risk under drying conditions. Human modification of Q through dams and climate change (Haddeland et al., 2014) are also not considered in this study. In addition, the land-use change over land may affect the above ocean processes mainly through modification of the energy balance and circulation in monsoon regions, which we do not account for. Changes in freshwater discharge to the oceans might have implications for ocean circulation and climate, as studies of for example river discharge to the Arctic Ocean showed (Peterson, 2002(Peterson, , 2006. However, moisture recycling's buffering effect (which mitigates river flow changes) should have a mitigating effect on the ocean's response to freshwater inflow. Otherwise, precipitation over the ocean can influence ocean salinity (IPCC, 2013) and precipitation patterns over land can be influenced by sea surface temperature (Xie et al., 2010), but we consider this outside the scope of our study and likely to be of minor importance for the research questions that we address. Our TMR analyses should, thus, be seen as an inquiry to better understand the relative importance of local and remote LUC effects on Q from a water balance perspective, rather than an exact prediction. Nevertheless, due to the inevitable recycling of moisture in the global hydrological cycle, uncertainties in the magnitude are unlikely to affect our key conclusions that upwind extra-basin LUC can be essential for Q. The magnitude of our estimated P (Fig. 5) and Q from LUC is conservative in comparison to the literature (Spracklen and Garcia-Carreras, 2015). For example, a meta-  analysis of 96 different general circulation models (GCMs) and regional climate model (RCM) deforestation simulations showed that under 10 % conversion of Amazon forest to pasture or soybean production, the inter-quartile range of rainfall change in the Amazon basin is 0 to −4 % (Spracklen and Garcia-Carreras, 2015). In comparison, the STEAM-WAM2layers approach with change from potential to current land-use change (i.e. 8.8 % deforestation extent in the Amazon) causes a rainfall reduction of 0.4 % in the Amazon and thus falls in the conservative range. In addition, our analyses concern mean annual Q, and can also be considered conservative in the sense that seasonal signals are expected to be much stronger.
The limitations of our methods should also been seen in the light of the strengths and limitations of alternative methods for studying hydrological LUC effects; see Table S1. The most complex and coupled modelling approaches account for the highest number of feedback processes. However, the high degree of freedom in GCMs also contributes to the high sensitivity of precipitation to initial conditions   Figure 9. Impacts of human land-use change on mean annual precipitation source (i.e. P import + P basin-recycling ) and evaporation sinks (i.e. E export + E basin-recycling ) for selected river basins (boundaries in dark yellow). and the low signal-to-noise ratios. For example, a scenario replacing natural with present-day land cover only detected a significant response in less than 5 % of all grid cells in a single model analysis (Findell et al., 2007) and less than 5 % in non-perturbed grid cells across seven different models (Pitman et al., 2009). Regional deforestation scenarios generate higher ratios of significant results near the source of change, albeit noise remains high in distant regions (Werth and Avis-sar, 2002). The challenges in simulating precipitation due to cloud formation, aerosol representation, and inherent uncertainties in circulation response (Aloysius et al., 2016;Koren et al., 2012;Shepherd, 2014), and non-closure of water balance in semi-coupled modelling approaches (Bring et al., 2015), also contribute to a high model dependence in estimates of river flow change from LUC (Kundzewicz et al., 2007). Thus, the sign, magnitude, and location of im-  Figure 10. Nation influence on river flow change depending on whether TMR is taken into account in the 15 studied basins. Without considering TMR, river flow change influence originates entirely from evaporation change within the basin. With consideration of TMR, nation influence to river flow change is considered as the sum of absolute changes in precipitation import and the sum of absolute changes in evaporation export (Methods). Single country contribution below 5 % is bundled into category "Other".
pacts vary widely among models (Aloysius et al., 2016;Pitman et al., 2009). Observation-based methods relate presence of vegetation or irrigation to precipitation or river flows using statistical methods, often in combination with moisture tracking to determine the geographical origin of rainfall (DeAngelis et al., 2010;Kustu et al., 2010Kustu et al., , 2011Spracklen et al., 2012). Limitations of this type of method include variations in data quality, challenges in isolating effects of land use from climate variability, and difficulties establishing causation from correlation (Matin and Bourque, 2015). Key elements missing in all approaches including our own are socio-economic dynamics and landscape resilience, which are complex issues currently explored in experimental model settings (Nitzbon et al., 2017;Reyer et al., 2015).

Future research outlook
Research of land-use change dynamics and moisture recycling is becoming increasingly detailed, looking at e.g. the role of transpiration for moisture recycling during dry periods Staal et al., 2018), studying the interplay with groundwater use and recharge (Keune et al., 2018), and identifying influential source areas within a basin (Weng et al., 2018;Staal et al., 2018). Nevertheless, a key challenge for considering TMR effects in water governance is the modelling uncertainties and inherent variabilities associated with landatmosphere feedback processes. The most complex modelling approaches account for the highest number of feedback processes. However, the sign, magnitude, and location of impacts vary widely even among state-of-the-art climate models (Pitman et al., 2009;Aloysius et al., 2016). Key future improvements in climate models' ability to simulate P from LUC will contribute to the governability of TMR. Indepth examination of differences in model simulation of P (e.g. the ongoing Precipitation Driver Response Model Intercomparison Project Myhre et al., 2017) is one step in this direction. Tracking moisture in coupled climate models could further help identify causes of simulated differences in atmospheric and hydrological outputs. Key elements missing in current research on LUC effects on hydrological flows include socio-economic dynamics and landscape resilience, which are complex issues currently explored in experimental model settings (Nitzbon et al., 2017;Reyer et al., 2015). In the meantime, "no-regret" policies in river basin management, where TMR objectives align with other aims, can potentially be explored in conjunction with LUC scenarios that include TMR effects.

Conclusions
We analysed the potential impact of human LUC on Q worldwide through TMR, and separately looked at the remote and local LUC effects of relevance to water governance. Despite the river basin being the standard unit in water governance and water resource management, we find that Q are ultimately dependent on the modifications in both incoming P and outflowing E. At the global scale, Q as a response to LUC is almost halved by taking TMR into account. Due to variations in moisture recycling patterns and LUC, the magnitude and spatial sources of the TMR effect on Q vary substantially among individual basins. In some basins, the remote LUC effect on Q exceeded local withinbasin effects (e.g. in the Amazon), while in others, TMR introduced considerable foreign nation influence on Q (e.g. in the Yangtze). International governance arrangements of teleconnnected LUC influence could be needed, even for river basins that today are not considered transboundary. We conclude that consideration of TMR is essential for understanding Q modifications and managing water resources in a rapidly changing and tele-coupled world (Liu et al., 2013) facing increasing pressure on both land (Schmitz et al., 2014) and water (Mekonnen and Hoekstra, 2016). Further research in both climate modelling and water governance strategies is needed to internalize land-atmosphere interactions in future water resource considerations.