1. Introduction
Nature's contributions to people (NCP), the benefits or detriments of living nature to people's quality of life (Díaz et al., Reference Díaz, Pascual, Stenseke, Martín-López, Watson, Molnár, Hill, Chan, Baste, Brauman, Polasky, Church, Lonsdale, Larigauderie, Leadley, van Oudenhoven, van der Plaat, Schröter, Lavorel and Shirayama2018), are increasingly recognized as an important concept in conservation research and sustainability governance (Brauman et al., Reference Brauman, Garibaldi, Polasky, Aumeeruddy-Thomas, Brancalion, DeClerck, Jacob, Mastrangelo, Nkongolo, Palang, Pérez-Méndez, Shannon, Shrestha, Strombom and Verma2020; Global Biodiversity Framework, 2021; IPBES, 2019; Millennium Ecosystem Assessment, Reference Millennium Ecosystem Assessment2005; Wood et al., Reference Wood, Jones, Johnson, Brauman, Chaplin-Kramer, Fremier, Girvetz, Gordon, Kappel, Mandle, Mulligan, O'Farrell, Smith, Willemen, Zhang and DeClerck2018). Well-recognized local benefits from NCP include water security, food security, hazard risk reduction, material well-being, cultural values, and climate security.
Vegetation-regulated moisture recycling (VMR), that is, the role of vegetation in generating evaporation and maintaining precipitation, is one crucial way that nature contributes to people via supporting global food security and economic prosperity (Keys et al., Reference Keys, Wang-Erlandsson and Gordon2016). Around a fifth of precipitation falling on global land areas is supplied by vegetation through VMR, and regionally, up to half of precipitation can originate from VMR (Keys et al., Reference Keys, Wang-Erlandsson and Gordon2016). Thus, through VMR, reductions in evaporation associated with land cover and land use change in upwind areas can lead to decreased downwind precipitation (see Box 1). An improved quantitative understanding is needed of the spatially explicit regions that are most essential for generating, and benefitting from, VMR.
The global hydrological cycle over land includes surface, subterranean, and atmospheric flows of water (Oki & Kanae, Reference Oki and Kanae2006), of which the latter are the largest (at 115,000 km3 yr–1 of terrestrial precipitation and 65,000 km3 yr–1 of terrestrial evaporation). In this work, we focus on the atmospheric flows of water, specifically tracking the contribution of evaporation from the land surfaces into the atmosphere, following the moisture flows across the planet, and locating its fate as precipitation elsewhere (van der Ent et al., Reference van der Ent, Savenije, Schaefli and Steele-Dunne2010; Figure 1). Moisture tracking has been studied for decades (Koster et al., Reference Koster, Jouzel, Suozzo, Russell, Broecker, Rind and Eagleson1986; Trenberth, Reference Trenberth1999), including tracking the atmospheric flows of water associated with terrestrial surfaces (Dirmeyer et al., Reference Dirmeyer, Koster and Guo2006; Dirmeyer & Brubaker, Reference Dirmeyer and Brubaker2007).
The important role of vegetation in supplying evaporation flows to the atmosphere and downwind precipitation has been demonstrated globally (Keys et al., Reference Keys, Wang-Erlandsson and Gordon2016; Wei et al., Reference Wei, Dirmeyer, Wisser, Bosilovich and Mocko2012), as well as regionally, including in the Amazon (Spracklen & Garcia-Carreras, Reference Spracklen and Garcia-Carreras2015; Staal et al., Reference Staal, Flores, Aguiar, Bosmans, Fetzer and Tuinenburg2020b), the Indo-Gangetic plain (Tuinenburg et al., Reference Tuinenburg, Hutjes, Stacke, Wiltshire and Lucas-Picher2014), Eastern China (Guo et al., Reference Guo, van der Ent, Klingaman, Demory, Vidale, Turner, Stephan and Chevuturi2019), the North American great plains (DeAngelis et al., Reference DeAngelis, Dominguez, Fan, Robock, Kustu and Robinson2010), and across Siberia (Oshima et al., Reference Oshima, Tachibana and Hiyama2015). Forests typically generate more evaporation than short vegetation types such as grasslands, as they are able to intercept more rainfall on the canopy, forest floor, stems, and epiphytes, as well as store and access more water from soil moisture, bedrock moisture, and groundwater. However, even short vegetation can generate substantial evaporation, particularly over large areas. Where soil moisture is not a constraining factor, grassland evaporation can during its active biomass growing season exceed forest evaporation (Teuling et al., Reference Teuling, Seneviratne, Stöckli, Reichstein, Moors, Ciais, Luyssaert, van den Hurk, Ammann, Bernhofer, Dellwik, Gianelle, Gielen, Grünwald, Klumpp, Montagnani, Moureaux, Sottocornola and Wohlfahrt2010). A desert planet is estimated to have a substantially less active water cycle with substantially less rainfall (Kleidon et al., Reference Kleidon, Fraedrich and Heimann2000), and human deforestation is already estimated to have reduced precipitation substantially (Wang-Erlandsson et al., Reference Wang-Erlandsson, Fetzer, Keys, van der Ent, Savenije and Gordon2018). Vegetation's role in remote precipitation is especially important since transpiration, the water flows that vegetation use for photosynthesis and biomass growth is typically more important in supporting precipitation over longer distances and during dry periods (O'Connor et al., Reference O'Connor, Dekker, Staal, Tuinenburg, Rebel and Santos2021; Staal et al., Reference Staal, Tuinenburg, Bosmans, Holmgren, van Nes, Scheffer, Zemp and Dekker2018; van der Ent et al., Reference van der Ent, Wang-Erlandsson, Keys and Savenije2014; Wang-Erlandsson et al., Reference Wang-Erlandsson, van der Ent, Gordon and Savenije2014).
Recent work emphasizes the outsized importance of densely vegetated areas, such as tropical forests, in both providing moisture for other downwind areas as well as for recycling water internally (Spracklen et al., Reference Spracklen, Baker, Garcia-Carreras and Marsham2018). For example, in the interior Amazon basin, there is evidence that the moisture recycled by the tropical forest contributes to the resilience of the forest itself, stabilizing rainfall to avoid falling below critical thresholds that could lead to (potentially irreversible) savannization (O'Connor et al., Reference O'Connor, Dekker, Staal, Tuinenburg, Rebel and Santos2021; Staal et al., Reference Staal, Tuinenburg, Bosmans, Holmgren, van Nes, Scheffer, Zemp and Dekker2018). In East Africa, the forested Kenyan highlands provide important moisture recycling both internally within the highlands themselves as well as for adjacent regions including the cross-border regions in Uganda and Rwanda (Keys et al., Reference Keys, Warrier, van der Ent, Galvin and Boone2022). Globally, tropical forests may buffer themselves against climate shocks, though this may change as global climate change proceeds into the 21st century (Staal et al., Reference Staal, Fetzer, Wang-Erlandsson, Bosmans, Dekker, van Nes, Rockström and Tuinenburg2020a).
Recent work has identified critical natural assets (CNA), the most important natural and semi-natural ecosystems for providing beneficial NCP, comprising 30–44% of the total global land area (Chaplin-Kramer et al., Reference Chaplin-Kramer, Neugarten, Sharp, Collins, Polasky, Hole, Schuster, Strimas-Mackey, Mulligan, Brandon, Diaz, Fluet-Chouinard, Gorenflo, Johnson, Kennedy, Keys, Longley-Wood, McIntyre, Noon and Watson2022). While this work included VMR among the 14 NCP considered, it only considered the benefits of VMR to rainfed working lands, rather than to the CNA themselves; the remote functional connections between downwind CNA areas and the upwind source areas whose moisture they rely on have yet to be mapped.
A CNA area is fundamentally an ecosystem that provides a critical benefit to people, and to the extent that all ecosystems depend on precipitation, or more precisely the maintenance of their current precipitation regime, disturbing that precipitation regime could result in alterations to the ecosystem that could diminish its capacity to provide benefits to people. From a policy standpoint, such interdependencies can be important to understand, because if CNA regions are maintained but the ecosystems they depend on for precipitation are not, land use in those source regions could lead to reductions in (or potential collapse of) NCP in dependent ecosystems (Staal et al., Reference Staal, Fetzer, Wang-Erlandsson, Bosmans, Dekker, van Nes, Rockström and Tuinenburg2020a). Likewise, if CNA for moisture regulation are not maintained, this could lead to reductions in precipitation in downwind areas, with implications, for example, river flows, water supply availability (Keys et al., Reference Keys, Wang-Erlandsson and Gordon2018), crop production (Bagley et al., Reference Bagley, Desai, Dirmeyer and Foley2012; Oliveira et al., Reference Oliveira, Costa, Soares-Filho and Coe2013), and self-amplified forest loss (Zemp et al., Reference Zemp, Schleussner, Barbosa, Hirota, Montade, Sampaio, Staal, Wang-Erlandsson and Rammig2017).
In this work, we explore for the first time the importance of sources of precipitation for local CNA regions globally, as well as the global areas that are reliant on CNA evaporation, specifically how global or continental-scale moisture recycling supports the ecosystems that provide more landscape-scale benefits. We further aim to explore the types and spatial distribution of anthropogenically modified landscapes, in an effort to understand the demographic and land management challenges that may be associated with sources and sinks of CNA moisture. This study will help identify regions to be prioritized for protection, restoration, and improved management, in ways that are synergistic with other conservation objectives (Keys et al., Reference Keys, Wang-Erlandsson, Gordon, Galaz and Ebbesson2017; te Wierik et al., Reference te Wierik, Gupta, Cammeraat and Artzy-Randrup2020).
2. Data and methods
This study first maps the moisture recycling associated with different CNA regions, and secondly, examines the composition of the CNA regions globally with regard to VMR (Keys et al., Reference Keys, Wang-Erlandsson and Gordon2016) and the anthropogenic biomes known as anthromes (Ellis & Ramankutty, Reference Ellis and Ramankutty2008).
2.1 Identifying the 90% CNA regions
To recognize the multitude of contributions that ecosystems provide to human societies, the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES) has underlined the importance of quantifying NCP (Díaz et al., Reference Díaz, Pascual, Stenseke, Martín-López, Watson, Molnár, Hill, Chan, Baste, Brauman, Polasky, Church, Lonsdale, Larigauderie, Leadley, van Oudenhoven, van der Plaat, Schröter, Lavorel and Shirayama2018). CNA were mapped by Chaplin-Kramer et al. (Reference Chaplin-Kramer, Neugarten, Sharp, Collins, Polasky, Hole, Schuster, Strimas-Mackey, Mulligan, Brandon, Diaz, Fluet-Chouinard, Gorenflo, Johnson, Kennedy, Keys, Longley-Wood, McIntyre, Noon and Watson2022) for 12 of nature's local contributions to people (including crop pollination, fodder for livestock, riverine fish harvest, marine fish harvest, nitrogen retention, sediment retention, flood mitigation, coastal protection, timber, fuelwood, access to nature, and marine recreation), identified through multi-objective optimization. While VMR was included in a separate prioritization to identify CNA for ‘global NCP’ (which included just VMR and carbon storage) and then overlapped with the CNA for these 12 local NCP, we focus here only on the local NCP to avoid circularity in the use of the VMR layer. Global maps of each of the 12 NCP mapped onto natural and semi-natural habitats at 300 m resolution (ESA-CCI land cover) were used to select 2 km planning units that achieve different target levels (ranging from 5 to 100%) of each NCP within each country in the minimum area. The NCP–area relationship asymptotes at around 90% of the total current value of all NCP in 30% of global land area, after which point much more area is required for each additional increment of NCP value; for this reason the 90% target was defined as ‘critical’. This analysis uses only the local CNA (optimized at the country level) rather than the global CNA (optimized globally), because the latter includes only two NCP that are relevant at only larger scales, carbon storage, and moisture recycling (from the same model used here). The local CNA regions at the original resolution of 2 km were resampled to 1.5 degrees to match the resolution of the VMR layer, using ‘near’ to produce a 0–1 mask and using ‘average’ to produce a proportional area raster for the amount of CNA area within the 1.5 degree grid cell. For this analysis, we focus on areas in which CNA (providing 90% of NCP at 2 kpm) occupy 90% of a VMR grid cell (at 1.5 degree), hereafter CNA regions. We note that we do not disaggregate the 90% CNA areas into constituent NCPs. We use the aggregated 90% CNA areas, to specifically understand the role of CNA hotspots in the broader atmospheric water cycle.
We also note that we use the lens of NCP, which is an inherently anthropocentric view of nature. These methods to delineate the functional interdependencies between ecosystems could be applied to other conservation priority maps, such as biodiversity hotspots or key biodiversity areas.
2.2 Moisture tracking with WAM-2layers
The Water Accounting Model 2-layers (WAM-2layers) is a Eulerian atmospheric moisture tracking model, which permits tracking evaporation arising from the surface of the Earth, the water vapor flowing through the atmosphere, and the eventual fate of precipitation elsewhere. We track the moisture budget in the ERA-Interim reanalysis, using six types of input data, including two-dimensional surface pressure, precipitation, and evaporation, and three-dimensional zonal and meridional winds, and specific humidity (Dee et al., Reference Dee, Uppala, Simmons, Berrisford, Poli, Kobayashi, Andrae, Balmaseda, Balsamo, Bauer, Bechtold, Beljaars, van de Berg, Bidlot, Bormann, Delsol, Dragani, Fuentes, Geer and Vitart2011). These data were downloaded at the 1.5 degree by 1.5 degree grid resolution. The two-dimensional data were downloaded at the three-hourly increment, and the three-dimensional data were downloaded at the six-hourly increment. We provide an overview of the WAM-2layers below, but other dimensions of model performance are explored and evaluated elsewhere, including the original implementation (van der Ent et al., Reference van der Ent, Savenije, Schaefli and Steele-Dunne2010), modification for two atmospheric layers (Wang-Erlandsson et al., Reference Wang-Erlandsson, van der Ent, Gordon and Savenije2014), evaluation of variability (Keys et al., Reference Keys, Barnes, van der Ent and Gordon2014), and comparison with results from regional climate models (van der Ent et al., Reference van der Ent, Tuinenburg, Knoche, Kunstmann and Savenije2013).
There are several procedures for the WAM-2layers analysis. First, the three-dimensional moisture data are integrated into two layers, representing an upper and lower level of the atmosphere. The boundary between the upper and lower layers generally corresponds to the approximate 800 hectopascal level, which can assist in separating upper and lower-level wind shear (Keys et al., Reference Keys, Barnes, van der Ent and Gordon2014). Second, the fluxes of moisture entering upward (evaporation), horizontally (water vapor transport along prevailing winds), and downward (precipitation) are calculated at the 15-minute timestep. This procedure is calculated for numerical stability (van der Ent et al., Reference van der Ent, Savenije, Schaefli and Steele-Dunne2010). Third, we use the 90% CNA data calculated above as a mask whereby we can track moisture forward in time (i.e. tracking evaporation from CNA origins through the atmosphere to its fate as precipitation elsewhere on Earth), and backward in time (i.e. tracking precipitation falling on CNA regions backward through the atmosphere to its origin as evaporation elsewhere on Earth).
The results of this analysis yield 15 years of moisture tracking information, from 2000 to 2014, aggregated at the monthly timestep. CNA moisture tracking information was calculated first for all ‘local’ CNA regions globally, and then separately for CNA regions within individual continents, including Africa, Australia, Eurasia, North America, and South America.
2.3 Vegetation-regulated moisture recycling
To distinguish the importance of vegetation for fostering current moisture recycling relationships, we employ previously calculated data on the role of vegetation in global moisture recycling (Keys et al., Reference Keys, Wang-Erlandsson and Gordon2016). These data provide an estimate of the fractional evaporation that arises globally and that can be attributed to the present composition and configuration of vegetation on terrestrial land surfaces. The data were generated, in part, by the Simple Terrestrial Evaporation to Atmosphere Model (STEAM) (Wang-Erlandsson et al., Reference Wang-Erlandsson, van der Ent, Gordon and Savenije2014), and the WAM-2layers. The data are available at the monthly timestep and employ the same ERA-Interim atmospheric water cycle data as that used in the present study, except for the evaporation data which are simulated by STEAM.
We note that the source of the three-dimensional atmospheric data (ERA-Interim) and the time period of analysis (2000–2014) for the moisture tracking data used in this analysis is the same as that used for the VMR analysis (Keys et al., Reference Keys, Wang-Erlandsson and Gordon2016). Thus, we can be confident that the three-dimensional moisture tracking (particularly the wind patterns) is identical between the two analyses.
2.4 Synthesizing source and sink data with anthropogenic biomes
Global data on anthropogenic biomes (i.e. anthromes) were originally downloaded at the 5 arc minute resolution. We perform a nearest neighbor spatial interpolation to upscale the anthrome data to the resolution of the ERA-Interim analysis, which is considerably coarser. However, because our analysis is inherently continental and global in nature, this coarsening of the data is considered acceptable. Future work that aims to explore country- or smaller-scale relationships would likely require a different approach. The global map of these anthromes is available in the supplementary information.
3. Results
Our results are broadly divided into two sections. First, we examine the results of moisture recycling associated with different CNA regions. Second, we examine the composition of the CNA regions by employing two additional datasets, VMR (Keys et al., Reference Keys, Wang-Erlandsson and Gordon2016) and the anthropogenic biomes known as anthromes (Ellis & Ramankutty, Reference Ellis and Ramankutty2008).
3.1 Substantial global moisture regulation occurs on land for both sources and sinks for CNA
Globally, CNA regions are spread across all populated continents (Figure 2, black outlines). The largest cluster of CNA regions in this analysis is located in Sub-Saharan Africa, in and around the Congo River basin. We find that source regions providing evaporation to precipitation within CNA regions (Figure 2a) are spread across all continents, and that taken collectively, land is responsible for 65% of the moisture falling as precipitation in CNA regions. Likewise, 44% of evaporation arising from CNA regions falls as precipitation on land elsewhere in the world.
3.2 CNA regions are associated with substantial internal moisture recycling
Moisture recycling dynamics reveal that CNA regions send and receive a substantial amount of moisture from and to one another (Figure 3). In this case, we find that CNA regions receive about 27% of average annual precipitation from other CNA regions. When restricted to specific continents, African and South American CNA have very high internal recycling from regional CNA, with 27 and 24%, respectively, of precipitation coming from intracontinental CNA regions. This demonstrates the outsize importance of CNA regions, particularly in those parts of the world, for providing supporting ecosystem services in the form of moisture recycling regulation.
3.3 Interpreting importance of contiguous ecosystems for moisture recycling
We explore the importance of current vegetation globally by using the VMR metric (Keys et al., Reference Keys, Wang-Erlandsson and Gordon2016). Using this, we identify the parts of the world that are sending (or receiving) a substantial amount of moisture to (or from) CNA regions. We define these source areas by the 10 mm yr−1 boundary of contribution, meaning that each year, regions enclosed in the pink boundaries contribute 10 mm yr−1 or more of evaporation to CNA precipitation (Figure 4). We see large, contiguous areas across much of the tropical world contributing both large fractions (Figure 4a) and absolute amounts (Figure 4b) of moisture to CNA precipitation.
The regions that are contributing 10 mm yr−1 or more of evaporation annually include substantial regions in tropical Central Africa and the Amazon. Also, the upper basin of the Congo River, the Sudd wetland in South Sudan, and the highlands around the Great Lakes of East Africa all provide substantial evaporation that is recycled as precipitation in CNA regions. Beyond the tropics, the key sources originate from the Himalayas, Eastern Europe, and western Russia, as well as the southeastern United States and the Pacific Northwest of North America.
A similar pattern is seen for moisture sinks as for moisture sources. Again, we see large, contiguous areas across much of the tropical world receiving both large fractions (Figure 5a) and absolute amounts (Figure 5b) of moisture from CNA evaporation. We define these as sink areas that receive 10 mm yr−1 or more of precipitation from CNA evaporation. The locations within the pink boundaries in Figure 5, that receive 10 mm yr−1 or more of precipitation from CNA evaporation, are more spatially extensive than the source regions displayed in Figure 4. Sinks extend beyond the tropics to include much of Sub-Saharan Africa, South America, substantial portions of Russia and China, and half of North America.
3.4 Anthropogenic biomes permit new perspective on CNA interconnections
Overlaying anthropogenic land uses (i.e. anthromes) with VMR, we find that as moisture sources, forests are generally above the average moisture recycling efficiency (Figure 6, top left; dashed line shows the area vs. volume relationship), meaning that they provide more moisture volume per unit area to CNA regions than do other evaporation sources, while rangelands and croplands tend to have equal or below average efficiency. For the VMR sinks (Figure 6, top right) we see a more heterogeneous result, with a mix of rangelands and forests above the average volume per area, and wild woodlands notably below average. Populated woodlands and rangelands are a dominant source of precipitation to CNA regions, and they also receive the greatest volume of CNA moisture (Figure 6). In contrast, wild woodlands are spatially the most frequent land use receiving moisture from CNA regions (Figure 6, top right: dark green dot at right side of sink graph).
In looking at the anthrome characteristics for the evaporative source regions of CNA, we see that there are substantial forested areas that provide evaporation to CNA. Tropical South America, northeastern Brazil, and the entire Guiana Shield provide substantial evaporation for CNA precipitation. Beyond the tropics, forested mountainous areas in central China and the Himalaya are critical sources of moisture for CNA regions. In Eurasia, the regions that provide substantial moisture recycling are dominated by heavily populated rainfed croplands, while sources in southeastern North America are predominantly residential woodlands and rainfed croplands (for more information see Supplementary Figure 1).
The locations that receive CNA precipitation include nearly all the tropical forests of Africa and South America, as well as tropical forests in Borneo and parts of southeast Asia. In extratropical South America and Africa, rangelands are the dominant landscape, while in North America and Eurasia, there is a heterogeneous cover of dense settlements, cultivated lands, rangelands, and woodlands.
4. Discussion
4.1 Managed lands play an integral role in CNA stewardship
This work enables the direct evaluation of the importance of different types of landscapes, and is unique in enabling the distinction between the performance of wildlands and managed lands. Wildlands, particularly tropical forests, often receive the majority of the attention with regard to climate and biodiversity concerns and are often prioritized in conservation conversations, yet wild tropical forests are not the only key sources of precipitation regulation for the ecosystems most important to people (i.e. containing CNA). We find that managed lands are also significant for modulating intra-CNA moisture recycling, serving as both sources of CNA precipitation and recipients of CNA evaporation. The importance of moisture recycling to managed lands such as populated woodlands and rangelands illustrates a key way that CNA regions support human well-being: providing moisture needed for rainfed livelihoods such as silviculture and pastoral production (Keys et al., Reference Keys, van der Ent, Gordon, Hoff, Nikoli and Savenije2012; Keys & Wang-Erlandsson, Reference Keys and Wang-Erlandsson2018). These lands require careful stewardship to ensure they are able to continue supporting CNA regions with VMR benefits. CNA regions internally recycle a large proportion of their moisture, but these areas still depend significantly on the lands beyond their boundaries for rainfall.
We find that managed woodlands are nearly as important as wild forests for consistently providing above-average volumes of water per unit area to CNA regions (Figure 6), despite natural forests' role in modulating and maintaining stable moisture recycling patterns globally (Meier et al., Reference Meier, Schwaab, Seneviratne, Sprenger, Lewis and Davin2021; Molina et al., Reference Molina, Salazar, Martínez, Villegas and Arias2019; Mu et al., Reference Mu, Biggs and De Sales2021; O'Connor et al., Reference O'Connor, Dekker, Staal, Tuinenburg, Rebel and Santos2021; Pranindita et al., Reference Pranindita, Wang-Erlandsson, Fetzer and Teuling2021; Staal et al., Reference Staal, Tuinenburg, Bosmans, Holmgren, van Nes, Scheffer, Zemp and Dekker2018; Wang-Erlandsson et al., Reference Wang-Erlandsson, Fetzer, Keys, van der Ent, Savenije and Gordon2018). The disproportionate importance of managed woodlands for CNA is likely a consequence of the spatial proximity of managed woodlands to CNA areas. Thus, both managed and wild forests are critical sources of resilience and stability to precipitation regimes for CNA regions globally, and both should be considered alongside CNA regions in conservation and management decision-making.
The importance of rangelands as a source of moisture for CNA regions, as well as a recipient of CNA moisture, is less well-appreciated. Rangelands, which are predominantly located in global arid and semi-arid regions and may comprise varying combinations of grasslands, shrublands, and partial treecover, can be highly sensitive to changes in precipitation patterns, including seasonality (Sircely et al., Reference Sircely, Conant and Boone2019), interannual variability (Wilcox et al., Reference Wilcox, Sorice and Young2011), long-term climatic changes (Boone et al., Reference Boone, Conant, Sircely, Thornton and Herrero2018; Foley et al., Reference Foley, Coe, Scheffer and Wang2003), and interactions with human societies (Galvin et al., Reference Galvin, Thornton, de Pinho, Sunderland and Boone2006). In addition to this sensitivity, dryland vegetation can play a complex role in responding to ecological disturbance, including fires, overgrazing or denuding of vegetation, soil compaction, woody encroachment, and introduction of invasive species (Miralles et al., Reference Miralles, Nieto, McDowell, Dorigo, Verhoest, Liu, Teuling, Johannes Dolman, Good and Gimeno2016; Turnbull et al., Reference Turnbull, Wilcox, Belnap, Ravi, D'Odorico, Childers, Gwenzi, Okin, Wainwright, Caylor and Sankey2012). Our findings highlight the need to further study the moisture recycling inter-relationships that rangelands can have, both with other systems (such as providing moisture to CNA regions) and in terms of their reliance on CNA regions, to better understand implications for management.
4.2 Transboundary moisture regulation requires transboundary governance
A key implication of our results is that moisture recycling connects CNA regions to other landscapes, including in adjacent and distant countries. This transboundary dimension invites a discussion about what forms of governance might be possible in the context of moisture recycling and the atmospheric water cycle, specifically in the context of forests themselves. Keys et al. (Reference Keys, Wang-Erlandsson, Gordon, Galaz and Ebbesson2017) suggest that transboundary moisture recycling governance exists along a continuum of increasing complexity, from simple bilateral relationships (e.g. proof of harm [Allum, Reference Allum, Bratspies and Miller2006]) to complex networks of overlapping and confounding connection, with less stringent but more flexible polycentric governance options (Galaz et al., Reference Galaz, Olsson, Hahn, Folke and Svedin2008). Likewise, te Wierik et al. (Reference te Wierik, Gupta, Cammeraat and Artzy-Randrup2020) explore governance as it relates to both atmospheric water and green water (i.e. the water that exists in the soil column and vegetation but that does not form surface runoff). This opens different challenges than transboundary moisture recycling governance alone since this framing also considers how land use change can modulate green water, and correspondingly, evaporation flows. Despite this, such a framing is critical for understanding transboundary moisture recycling, particularly as it relates to interlinkages with the science of NCP which is the foundation for understanding CNA regions globally.
Other work has explored the importance of forests in fostering transboundary moisture recycling resilience. Ellison et al. (Reference Ellison, Claassen, Van Noordwijk, Sullivan, Vira, Xu, Archer and Haywood2018) explore the institutional scales at which forest–water relationships are governed, including the presence or absence of international agreements that can be meaningfully brought to bear in transboundary forest- and water-related topics. Likewise, other work has suggested that forest management should be explicitly linked to the maintenance of moisture recycling patterns, specifically given the fact that terrestrial landscapes recycle moisture many times as it cascades across the land surface (Food and Agriculture Organization of the United Nations, 2019). What our analysis adds to this body of evidence is that transboundary governance of moisture recycling is important not just for water management but for the variety of other ways that natural ecosystems support human well-being. The connections between forests and transboundary moisture recycling governance should be considered in the discourse around the global conservation of wilderness and semi-natural areas (Beyer et al., Reference Beyer, Venter, Grantham and Watson2020; Riggio et al., Reference Riggio, Baillie, Brumby, Ellis, Kennedy, Oakleaf, Tait, Tepe, Theobald, Venter, Watson and Jacobson2020).
4.3 Uncertainties and research frontiers
This research connects several areas of study that have not traditionally been well-linked, specifically a largely academic research effort on moisture recycling dynamics, the conceptual and theoretical discussion of Earth system governance, and conservation planning. It also raises some key uncertainties that must be addressed to deepen understanding of how long-distance connections could be better integrated into conservation planning. We propose five new research frontiers based on our current work.
First, the scale dependence of moisture recycling phenomena is a key uncertainty. Our results demonstrate the average exchange of moisture among CNA regions for the period of time 2000–2014. This provides an important starting point for understanding the interannual patterns of the recent, observed past. However, while the data we use is suitable for global and regional analysis of patterns at the monthly or annual timestep, neither is sufficient for analyzing specific patterns, especially if being put to use for policy. Different moisture tracking models (Tuinenburg et al., Reference Tuinenburg, Theeuwen and Staal2020; van der Ent et al., Reference van der Ent, Savenije, Schaefli and Steele-Dunne2010; Vázquez et al., Reference Vázquez, Nieto, Liberato and Gimeno2021; Wei et al., Reference Wei, Dirmeyer, Wisser, Bosilovich and Mocko2012) and numerous types of climate data (Bosilovich et al., Reference Bosilovich, Robertson, Takacs, Molod and Mocko2017; e.g. Dee et al., Reference Dee, Uppala, Simmons, Berrisford, Poli, Kobayashi, Andrae, Balmaseda, Balsamo, Bauer, Bechtold, Beljaars, van de Berg, Bidlot, Bormann, Delsol, Dragani, Fuentes, Geer and Vitart2011) are available and appropriate for different scales. Global analyses may be similar at regional to global spatial scales and monthly to annual timescales (Keys et al., Reference Keys, Barnes, van der Ent and Gordon2014; van der Ent et al., Reference van der Ent, Tuinenburg, Knoche, Kunstmann and Savenije2013), but different models and datasets will yield different results at finer spatial and temporal resolutions. Multiple simulations using multiple data sources (and possibly multiple moisture tracking models) will be required for specific policy relevance in particular contexts, to ensure robust results and better ability to explore uncertainties, both of which are essential to decision-making.
Likewise, exploration of different time horizons and dynamic processes is needed to inform intra- and interannual variability and long-term change. Seasonality may drive more important changes in certain regions than mean annual characteristics. For example, the relatively larger importance of moisture recycling during dry seasons can be relevant for conservation policy in semi-arid and arid regions (Falkenmark et al., Reference Falkenmark, Wang-Erlandsson and Rockström2019; Staal et al., Reference Staal, Tuinenburg, Bosmans, Holmgren, van Nes, Scheffer, Zemp and Dekker2018; Wang-Erlandsson et al., Reference Wang-Erlandsson, van der Ent, Keys, Fetzer, Taniguchi and Gordonn.d.). Dynamic climate simulations that can capture more detailed dynamic processes could also be used to better constrain and explore different types of uncertainty. While we do not consider climate change in this analysis, other work has begun examining how future changes in precipitation and rising temperatures could modulate moisture recycling (Findell et al., Reference Findell, Keys, van der Ent, Lintner, Berg and Krasting2019). Accounting for climate change will be critical for robust long-term biodiversity conservation planning (Arneth et al., Reference Arneth, Shin, Leadley, Rondinini, Bukvareva, Kolb, Midgley, Oberdorff, Palomo and Saito2020; Hannah et al., Reference Hannah, Roehrdanz, Marquet, Enquist, Midgley, Foden, Lovett, Corlett, Corcoran, Butchart, Boyle, Feng, Maitner, Fajardo, McGill, Merow, Morueta-Holme, Newman, Park and Svenning2020; Pecl et al., Reference Pecl, Araújo, Bell, Blanchard, Bonebrake, Chen, Clark, Colwell, Danielsen, Evengård, Falconi, Ferrier, Frusher, Garcia, Griffis, Hobday, Janion-Scheepers, Jarzyna, Jennings and Williams2017). Future work should also consider interactions between climate trajectories and changing land use, for example, cascading impacts, surface water trade-offs, and geopolitical change. Additional questions related to longer-term sustainability and environmental justice, drawing on the increasing burden of climate change impacts, will also be important to consider (Menton et al., Reference Menton, Larrea, Latorre, Martinez-Alier, Peck, Temper and Walter2020).
Second, our results clearly identify ecologically critical, yet potentially fragile parts of the coupled Earth system including the tropical forests in both the Amazon and the Congo as key sources and sinks of moisture recycling. Recent advances in defining the specific interactions of the entire water cycle in terms of Earth system tipping points provide new frontiers of research and practice related to moisture recycling (Gleeson et al., Reference Gleeson, Wang-Erlandsson, Porkka, Zipper, Jaramillo, Gerten, Fetzer, Cornell, Piemontese, Gordon, Rockström, Oki, Sivapalan, Wada, Brauman, Flörke, Bierkens, Lehner, Keys and Famiglietti2020a, Reference Gleeson, Wang-Erlandsson, Zipper, Porkka, Jaramillo, Gerten, Fetzer, Cornell, Piemontese, Gordon, Rockström, Oki, Sivapalan, Wada, Brauman, Flörke, Bierkens, Lehner, Keys and Famiglietti2020b). Work by Zipper et al. (Reference Zipper, Jaramillo, Wang-Erlandsson, Cornell, Gleeson, Porkka, Häyhä, Crépin, Fetzer, Gerten, Hoff, Matthews, Ricaurte-Villota, Kummu, Wada and Gordon2020) highlights the social and ecological scales at which different aspects of the water cycle can be reasonably managed. We find that moisture recycling is simultaneously a transboundary phenomenon and that any management will likely need to be integrated with intersectoral resource governance. This is echoed by Ahlström et al. (Reference Ahlström, Hileman, Wang-Erlandsson, García, Moore, Jonas, Pranindita, Kuiper, Fetzer, Jaramillo and Svedin2021), suggesting opportunities for engagement with resource management authorities, regional scientific organizations, and other forms of transboundary decision-making (Keys et al., Reference Keys, Wang-Erlandsson, Gordon, Galaz and Ebbesson2017; te Wierik et al., Reference te Wierik, Gupta, Cammeraat and Artzy-Randrup2020). The advent of computationally efficient simulations of Earth system dynamics (which include human societies) permits a new window into understanding the practical consequences of transgressing social and ecological tipping points (Wunderling et al., Reference Wunderling, Gelbrecht, Winkelmann, Kurths and Donges2020, Reference Wunderling, Donges, Kurths and Winkelmann2021). While this remains a theoretical aspect of Earth system science, the ability to link tipping point simulations to current policy discourse permits a realistic representation of Earth system surprises alongside conservation-related policies (particularly those related to tipping-prone systems like the Amazon rainforest). Such methods could enable forecasting of the ecological consequences of specific policy decisions (e.g. relaxed enforcement of protected Amazon forest, potential resulting savannization [Hirota et al., Reference Hirota, Holmgren, Van Nes and Scheffer2011; Nepstad et al., Reference Nepstad, Stickler, Filho and Merry2008; Staal et al., Reference Staal, Fetzer, Wang-Erlandsson, Bosmans, Dekker, van Nes, Rockström and Tuinenburg2020a]), and consequences for downwind human systems via reductions in moisture recycling to Bolivia, Paraguay, and Uruguay (Keys & Wang-Erlandsson, Reference Keys and Wang-Erlandsson2018). These research questions can now be reasonably explored given recent conceptual, computational, and scientific advances.
Third, context-specific cases could be explored that enable higher-resolution modeling of coupled land use change and moisture recycling impacts, but with explicit inclusion of management criteria (e.g. different regulatory regimes in protected areas) alongside distinct political contexts (e.g. cooperative, antagonistic). Beyond addressing the resolution concerns addressed above, this sort of tailored case analysis would permit the production of management-relevant insights even more useful for addressing real-world problems. In particular, a blind spot in much of this work is the acute lack of moisture recycling analysis in the Congo basin (including both ground-based and modeling studies). Central Africa is among the most important CNA regions globally (Chaplin-Kramer et al., Reference Chaplin-Kramer, Neugarten, Sharp, Collins, Polasky, Hole, Schuster, Strimas-Mackey, Mulligan, Brandon, Diaz, Fluet-Chouinard, Gorenflo, Johnson, Kennedy, Keys, Longley-Wood, McIntyre, Noon and Watson2022), and it is clear from our analysis that it is also critical from a moisture recycling perspective. Thus, we need a much better, broader understanding of this region from an integrated land and moisture recycling perspective, as well as from a socio-economic perspective. This is particularly urgent given increased agricultural and ecological drought risks across Africa, as well as carbon sink-to-source conversion across the tropics due to tree mortality (Hubau et al., Reference Hubau, Lewis, Phillips, Affum-Baffoe, Beeckman, Cuní-Sanchez, Daniels, Ewango, Fauset, Mukinzi, Sheil, Sonké, Sullivan, Sunderland, Taedoumg, Thomas, White, Abernethy, Adu-Bredu and Zemagho2020).
Fourth, there is a fundamental need to consider systemic vulnerability and situated vulnerability, such that different socio-economic contexts make people more (or less) directly nature-dependent. Some work has attempted to explore moisture recycling vulnerability from a variety of perspectives (Keys & Wang-Erlandsson, Reference Keys and Wang-Erlandsson2018; Keys et al., Reference Keys, van der Ent, Gordon, Hoff, Nikoli and Savenije2012, Reference Keys, Wang-Erlandsson and Gordon2018), yet very little work has explicitly focused on the nexus of moisture recycling, ecosystem conservation, and human well-being. Keys and Falkenmark (Reference Keys and Falkenmark2018) highlight the importance of considering the ability to change moisture recycling (and detect such a change) in the context of Sustainable Development Goal (SDG) achievement, given that moisture recycling might drive surprising changes in ecosystem function but be delayed substantially in time. Weng (Reference Weng2020) suggests that thoughtful reforestation in the Amazon could lead to substantial improvements in the stability of the atmospheric water cycle for Bolivian farmers. Recent studies suggest that afforestation and forest restoration have substantial global potential to increase water availability both locally and regionally through moisture recycling (Cui et al., Reference Cui, Lian, Huntingford, Gimeno, Wang, Ding, He, Xu, Chen, Gentine and Piao2022; Hoek van Dijke et al., Reference Hoek van Dijke, Herold, Mallick, Benedict, Machwitz, Schlerf, Pranindita, Theeuwen, Bastin and Teuling2022).
Fifth, the results of our work should be updated regularly to reflect changes in underlying conservation efforts and to see whether those changes can be detected at national, regional, or global scales in the moisture recycling signature. This work points to several areas of future research that must be pursued to better inform national and international decision-making. Likewise, the entire analysis could be updated to reflect major policy changes, or even pre-emptively simulated to support more-informed decision-making. For example, policy implementing new targets set by the Convention on Biological Diversity (CBD) would benefit from multiple scenarios that consider moisture recycling changes to reflect the conservation-related consequences for the water cycle.
5. Conclusions
Conserving CNA regions in isolation is not sufficient to protect them and support their functioning; better conservation and management of the lands outside CNA regions is needed as well, especially in semi-populated managed rangelands and woodlands. A better understanding of the ways in which VMR connects spatially distant wild and working lands is an essential contribution to the global sustainability conversation. Transboundary water management in international river basins is already recognized as an important yet extremely difficult-to-achieve objective (Mason et al., Reference Mason, Ward, Watson, Venter and Runting2020; Zeitoun et al., Reference Zeitoun, Goulden and Tickner2013). Similarly, carbon management for climate is a global challenge, and fighting biodiversity loss while sustaining NCP requires international cooperation to succeed (Schleicher et al., Reference Schleicher, Zaehringer, Fastré, Vira, Visconti and Sandbrook2019; TWI2050-The World in, 2018). Like these objectives, VMR crosses the boundaries of watersheds and nations, yet it does not receive the same attention in international agreements such as the SDGs or the CBD.
We recommend VMR be considered as a crucial NCP and prioritized alongside biodiversity, climate change mitigation, and integrated watershed management in the international sustainable development conversation – not only among national governments, but within UN policy-setting contexts, international aid investment decision-making, and corporate business practice improvement. The best available science should underpin global action for achieving sustainable development, including the SDGs (Sachs et al., Reference Sachs, Schmidt-Traub, Kroll, Lafortune, Fuller and Woelm2020a, Reference Sachs, Schmidt-Traub and Lafortune2020b) and the CBD's Global Biodiversity Framework (Global Biodiversity Framework, 2021). Better accounting for processes like VMR in land management decision-making would bring us one step closer to achieving such ambitious goals, particularly in regard to providing greater resilience in moderating local climate.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/sus.2023.24.
Acknowledgements
None.
Author contributions
P. W. K., P. M. C., R. C.-K., and L. W.-E. all conceived of the analysis. P. W. K. performed the analysis. P. W. K., P. M. C., R. C.-K., and L. W.-E. wrote the manuscript.
Funding statement
P. W. K. acknowledges support from National Aeronautics and Space Administration (NASA) Grant number 80NSSC19K0182; L. W.-E. acknowledges financial support from the European Research Council through the ‘Earth Resilience in the Anthropocene’ project (no. ERC-2016-ADG 743080), and from Formas the Swedish Research Council (grant no. 2019-01220); P. M. C. acknowledges support from the Gordon and Betty Moore Foundation.
Competing interests
None.
Research transparency and reproducibility
All data used in this analysis are available at the Mountain Scholar permanent data repository: http://dx.doi.org/10.25675/10217/234391