water stress research articles

An official website of the United States government

The .gov means it’s official. Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

The site is secure. The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

• Publications
• Account settings

Preview improvements coming to the PMC website in October 2024. Learn More or Try it out now .

• Advanced Search
• Journal List
• Am J Public Health
• v.101(8); Aug 2011

Water Stress and Water Scarcity: A Global Problem

Water scarcity is a result of an imbalance between the supply of and demand for water sources in a geographical area. A limited or inadequate water supply is one of the toughest issues facing the world's poorest countries, but the populations of these countries are not the only populations to endure the burden of water scarcity. 1 , 2 As Levy and Sidel correctly point out, distribution in water-strained settings is the source of numerous conflicts globally. 3 Most dialogues about this topic, either explicitly or implicitly, talk about water stress and conflict as an issue limited to the low- and middle-income tropical or subtropical countries, as most of their examples affirm. 3 As a result, limited literature and discussion exists on the public health effects of water scarcity in wealthier countries or regions such as the United States, 4 , 5 Europe, 6 Australia, 7 , 8 or the Asia-Pacific where water scarcity threatens regional food security through a number of ecological, political, and economic pathways. 9 Perhaps this paucity of discourse stems from a false notion that these countries and regions do not currently and will not contend with the further consequences of population growth, urbanization, and anthropogenic and natural climate change. However, such an assumption would be far from the truth 10 as emerging empirical literature from these regions confirms. 5 – 9

A jar of water sits on a counter in the new science building at the Bertschi School in Seattle, WA, which is one of several dozen nationwide that recycle and manage rainwater and wastewater on site. Photograph by Ted S. Warren. Printed with permission of AP Wide World.

Insufficient access to water and ensuing water stress is undoubtedly a major global public health threat and challenge for our generation. As Levy and Sidel discuss, many roles exist for public health professionals in shaping and “ensuring equitable access to freshwater.” 3 (p779) Where public health workers devote their efforts and pick their battles is often determined by the literature and surrounding debate. In this instance, I pose that a broader debate and research agenda is needed. The majority of collective public health efforts ought to be directed to those countries striving to attain Millennium Development Goal (MDG) Target 7C: “Halve, by 2015, the proportion of the population without sustainable access to safe drinking water and basic sanitation” (available at: http://www.un.org/millenniumgoals ). However, our dialogue and concern should not be limited to only those countries striving to attain MDG 7.C. By doing so we, as public health workers, risk not only overlooking a number of populations with imminent water issues in our own backyards 4 , 5 , 7 , 8 but also missing some easy successes in populations with enormous per capita water consumption compared with those populations struggling with MDG Target 7C.

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

• View all journals
• My Account Login
• Explore content
• About the journal
• Publish with us
• Sign up for alerts
• Open access
• Published: 21 January 2022

Hotspots for social and ecological impacts from freshwater stress and storage loss

• Xander Huggins   ORCID: orcid.org/0000-0002-6313-8299 1 , 2 ,
• Tom Gleeson   ORCID: orcid.org/0000-0001-9493-7707 1 , 3 ,
• Matti Kummu   ORCID: orcid.org/0000-0001-5096-0163 4 ,
• Samuel C. Zipper   ORCID: orcid.org/0000-0002-8735-5757 5 ,
• Yoshihide Wada   ORCID: orcid.org/0000-0003-4770-2539 6 ,
• Tara J. Troy   ORCID: orcid.org/0000-0001-5366-0633 1 &
• James S. Famiglietti 2 , 7  

Nature Communications volume  13 , Article number:  439 ( 2022 ) Cite this article

24k Accesses

42 Citations

167 Altmetric

Metrics details

• Water resources

Humans and ecosystems are deeply connected to, and through, the hydrological cycle. However, impacts of hydrological change on social and ecological systems are infrequently evaluated together at the global scale. Here, we focus on the potential for social and ecological impacts from freshwater stress and storage loss. We find basins with existing freshwater stress are drying (losing storage) disproportionately, exacerbating the challenges facing the water stressed versus non-stressed basins of the world. We map the global gradient in social-ecological vulnerability to freshwater stress and storage loss and identify hotspot basins for prioritization ( n  = 168). These most-vulnerable basins encompass over 1.5 billion people, 17% of global food crop production, 13% of global gross domestic product, and hundreds of significant wetlands. There are thus substantial social and ecological benefits to reducing vulnerability in hotspot basins, which can be achieved through hydro-diplomacy, social adaptive capacity building, and integrated water resources management practices.

Similar content being viewed by others

Diagnosing challenges and setting priorities for sustainable water resource management under climate change

Ibrahim Nourein Mohammed, John D. Bolten, … Derek Vollmer

Prefectures vulnerable to water scarcity are not evenly distributed across China

Fubo Zhao, Xi Wang, … Shailesh Kumar Singh

Importance and vulnerability of the world’s water towers

W. W. Immerzeel, A. F. Lutz, … J. E. M. Baillie

Introduction

Humans and ecosystems, alongside hydrological and biogeochemical cycles, are a deeply coupled and global social-ecological system 1 , 2 , 3 . Social-ecological systems are complex adaptive systems formed by interactions and feedbacks between biophysical and social processes 4 , 5 , and freshwater is fundamental for flourishing and resilient ecosystems, societies, and the larger Earth System 3 , 6 , 7 , 8 . Global freshwater storage and flows are dynamic with oscillations and persistent trends driven by the combined influence of human activity, climate change, and natural variability on sub-seasonal to multi-decadal timescales 9 , 10 . Human dominance over the water cycle is increasingly recognized 11 and the continued development of global hydrological models enables higher fidelity representations of human and climate change impacts on water resources 12 . However, the science of understanding the reciprocal impacts of freshwater stress and storage trends on humans and ecosystems remains in its infancy at the global scale. We argue that studying both directions of this coupled social-ecological system (i.e., social-ecological activity impacts on freshwater and freshwater impacts on social-ecological activity) is crucial to confronting global freshwater challenges, yet the latter has received considerably less attention. In this paper, we consider the potential for freshwater stress and storage loss to impact humans and ecosystems. We do this by synthesizing a subset of the few but critical global ecohydrological and sociohydrological datasets with freshwater storage, freshwater withdrawal, and streamflow datasets (see Supplementary Table  1 ).

We seek to build on the existing literature on global freshwater scarcity and security topics, which broadly address social and ecological impacts of freshwater-related stresses and hazards. We refer here to freshwater scarcity studies as those which evaluate the ratios of water use to streamflow and streamflow per capita, typically at the basin scale, e.g., 13 , 14 , 15 , 16 and to freshwater security studies as those which integrate multidimensional indicators of physical, chemical, socioeconomic, and institutional factors and aggregate using grid-based, basin, or administrative discretization schemes, e.g., 17 , 18 , 19 . While water scarcity assessments exclusively focus on freshwater stress, it is only one element of physical water security and thus only one component of water security assessments. Both approaches, however, have important limitations that constrain their ability to support specific conclusions and drive policy implementation regarding the impacts of freshwater-related hazards such as stress and storage loss.

For instance, freshwater scarcity assessments typically apply globally-consistent classification schemes which do not represent important spatial variations in social and ecological sensitivities and responses. A few holistic derivatives of freshwater scarcity assessments, such as the social water stress 20 and water poverty 21 indices, have only been evaluated at the national level. Alternatively, water security assessments consider water scarcity as just one of many input variables. These assessments typically aggregate multidimensional indicators of different aspects of water security, which can lead to similar water security outcomes with different input indicator combinations. As a result, water scarcity impacts become challenging to isolate from final water security assessment results. Furthermore, this aggregation approach does not consider interactions or relationships between elements of water security which are critical determinants of social-ecological system behavior 22 .

In this paper, we combine the strengths of water scarcity and water security research, and address their limitations by integrating concepts from social-ecological systems research. We combine concepts from these fields to address the following core objectives of this study: (1) Asses the global co-occurrence of freshwater stress and freshwater storage trends at the basin scale. (2) Analyze the relationship between social adaptive capacity and ecological sensitivity indicators with freshwater stress and storage trends. (3) Derive the global gradient in social-ecological vulnerability to freshwater stress and storage trends by considering all indicators listed above, and identify hotspot basins as those with high vulnerability values with respect to the global distribution. (4) Evaluate current levels of integrated water resources management within hotspot basins. Basins, at various scales, are an increasingly used and particularly suitable geospatial unit of analysis for hydrologically-based social-ecological systems analysis 23 . In this study, all analyses are performed at a large basin scale ( n  = 1204, median area ~70,000 km 2 ). Input data align to the year 2015 as best as possible, and data are summarized to the basin scale by computing the area-weighted basin average or within-basin sum, depending on the intensive or extensive nature of each dataset (see “Methods” section and Supplementary Information ). See Box  1 for definitions of key terminology used in this paper.

The global co-occurrence of freshwater stress and freshwater storage trends

We mapped freshwater stress and trends in freshwater storage at the basin scale and analyzed the co-occurrence of these phenomena (Fig.  1 ).

a Freshwater stress, derived from freshwater withdrawal and streamflow datasets (see “Methods” section). b Freshwater storage trend per basin. c Combinations of freshwater stress and storage trend per basin, which together derive basin freshwater status (shown in Fig.  2b ). Values overlaying the legend indicate the number of basins satisfying each set of conditions. For categorical plotting purposes only, ±3 mm year −1 is used as the threshold denoting a clear directional storage trend, based on the error level of the underlying observations 25 . d – g The exposure of social-ecological activity to freshwater stress and storage trends. Each plot represents storage trends as the x -axis coordinate, and log-transformed freshwater stress as the y-axis coordinate with the size of each circle based on the basin’s value respective to each plotting dimension.

Freshwater stress represents the state of demand-driven water scarcity 15 and is defined as the ratio of freshwater withdrawal to streamflow (Fig.  1a ). Trends in freshwater storage, conversely, represent the evolution of total storage, defined as the vertical sum of groundwater, soil moisture, surface water, and snow water equivalent storages (Fig.  1b ). Freshwater stress and storage are linked, as freshwater storage becomes a required source of water during periods when demands exceed supply. As climate change intensifies hydrological extremes globally, the strategic importance of the world’s largest store of liquid freshwater, groundwater, will only continue to increase 24 . Though studies have focussed on global assessments of freshwater stress 13 , 14 , 15 and trends in freshwater storage 9 , no study to date has mapped these two variables against one another. Doing so provides important context to differentiate basins of equal freshwater stress, as drying trends are likely to exacerbate challenges derived from freshwater stress, while wetting trends may yield offsetting effects. However, as freshwater stress calculations do not differentiate between withdrawals sourced from streamflow or storage, the two variables are not necessarily independent.

We found that 201 (42%) of the 478 currently stressed basins (withdrawal/streamflow > 0.10) are simultaneously losing freshwater storage (Fig.  1c ). These basins are located in south and southwestern USA, northeastern Brazil, central Argentina, Algeria, and concentrate throughout the Middle East, the Caucasus, northern India, and northern China. Predominantly, these regions are agriculturally significant and heavily irrigated 9 , with the exception of a few basins in South America whose trends are likely the product of natural variability 9 . Conversely, 98 (21%) of the currently stressed basins are gaining freshwater storage. The storage trends in these basins have largely been attributed to natural variability with the exception of central India, whose trends are partially attributed to groundwater recovery following groundwater policy change 9 . The remaining 179 stressed basins have freshwater storage trends that are smaller than can be definitively interpreted from the satellites monitoring these trends 25 . This skew towards negative storage trends (i.e., drying) in the world’s water-stressed basins dissipates and even reverses in the non-stressed basins, where drying and wetting trends are found in 23% and 32% of the 726 non-stressed basins, respectively. While previous work has shown that the world’s dry regions are becoming drier while the wet regions are becoming wetter 26 , this work reveals that the stressed regions of the world are becoming drier while the non-stressed regions of the world have no clear overall trend in freshwater storage.

The encompassed human population, food crop production, gross domestic product (GDP), biodiversity, and wetlands enumerate the potential social-ecological impacts from the current state of global freshwater stress and storage trends. Around 2.2 billion people, 27% of global food crop production, and 28% of global GDP live, grow, and situate in freshwater stressed basins that are drying (Fig.  1d–f ). These totals represent an upper limit as not all social and ecological activity within these basins will be affected by freshwater stress and storage loss, which will depend on local levels of adaptive capacity and ecological sensitivity 22 (our focus in the subsequent sections). Conversely, 1.2 billion people, 24% of global food crop production, and 19% of global GDP are found in stressed basins that are wetting. We find less taxonomic biodiversity in the freshwater stressed and drying basins, and greater biodiversity in unstressed and wetting basins. Roughly the same number of wetlands of international importance are found in stressed and drying basins as in stressed and wetting basins. While these totals represent the magnitude of potentially affected biodiversity and wetlands, taxonomic biodiversity is only one of many critical facets of biodiversity 27 , and freshwater stress and storage trends are but two of many variables impacting global biodiversity 28 . Thus, we urge caution in interpreting the role of freshwater stress and storage in driving differences in these biodiversity distributions.

The most vulnerable populations to freshwater stress and storage loss

To better characterize social vulnerability, freshwater stress and storage loss must be placed in the context of social adaptability. We mapped and analyzed the co-occurrence of freshwater stress and storage trends with an existing global dataset of social adaptive capacity 23 summarized at the basin scale (Fig.  2 ). Social adaptive capacity (Fig.  2a ), or adaptability, represents “the ability of the system to respond to disturbances” 29 and is derived based on input indicators of governance, economic strength, and human development. This consideration of social adaptability enables more representative estimates of social, agricultural, and economic activity that are vulnerable to the co-occurrence of freshwater stress and storage loss. To consider freshwater stress and storage loss together, we developed the basin freshwater status indicator (Box  1 ) where higher values indicate co-occurring freshwater stress and storage loss (Fig.  2b , see “Methods” section).

a Social adaptive capacity, or adaptability, per basin. b Basin freshwater status, representing the combination of freshwater stress and storage trend per basin (see “Methods” section). c Combinations of basin freshwater status and social adaptability. Values overlaying the legend indicate the number of basins satisfying each set of conditions. d – g The exposure of social-ecological activity to basin freshwater status ( x -axis coordinate) and social adaptive capacity ( y -axis coordinate), with the size of each circle scaled based on the basin’s value respective to each plotting dimension. These distributions are summarized below each plot. P notation represents the percentile distribution.

We found 73 basins to possess low levels of social adaptability and severe basin freshwater status (Fig.  2c ). These basins concentrate in Northern, and Eastern Africa, the Arabian Peninsula, and Western, Central, and Southern Asia; although vulnerable basins are also found in northeast Brazil, Southern Africa, and northern China. These basins encompass approximately 1.2 billion people, 12% of global food crop production, and 6% of global GDP (Fig.  2d–f ). Conversely, 119 and 49 basins are found to have similarly severe basin freshwater status yet have moderate or high levels of social adaptability, respectively. These basins are located in southwestern USA and Mexico, Chile and Argentina, the Arabian Peninsula, regions surrounding the Caspian Sea, western Australia, and the North China Plain.

These differences in social adaptability across basins with severe freshwater status (i.e., co-occurring freshwater stress and storage loss) raise important economic considerations. First, greater social adaptability likely coincides with greater technological and economic capacity to pursue development. This development may consume greater volumes of freshwater and drive basins towards greater levels of freshwater stress or storage loss, while simultaneously increasing institutional and technical capacity to cope with limited water resources. Furthermore, freshwater stress and storage loss are not certain to induce negative economic impacts on basins, and can lead to positive impacts if a region is able to leverage its comparative advantages (e.g., irrigation efficiency) among other stressed regions 30 . Second, the divergent economic situation facing basins with severe freshwater status is particularly evident on a per-capita basis. In severe freshwater status, low adaptability basins, there resides 17% of the global population yet only 6% of global GDP. Conversely, in severe freshwater status basins with moderate-and-greater social adaptability, there resides 14% of the global population and an outsized 18% of global GDP (Fig.  2d, f ). It is thus paramount that global initiatives prioritize and link economic inequality with freshwater goals. One such example is Sustainable Development Goal (SDG) 6.4 (“reduce the number of people suffering from water scarcity”), which we argue should increasingly be linked to targets of SDG 10 (“reduce inequality within and among countries”).

Box 1 Key terminology as used in this paper. See Methods for further information

Freshwater stress: The ratio of annual freshwater withdrawal (W) to annual streamflow (Q). We refer to basins with W/Q ≥ 10% as stressed basins and those with W/Q ≥ 40% as highly stressed basins .

Freshwater storage trends: Year-over-year trends in total freshwater storage based on satellite observations over the 2002–2016 time period. Total freshwater storage is a vertically aggregated measure of water storage that includes groundwater, soil water, surface water, canopy water, and ice and snow water equivalents where present. For simplicity, we refer to negative freshwater storage trends as drying trends or storage loss and positive trends as wetting trends or storage gain .

Basin freshwater status: An integrated indicator that combines normalized freshwater stress and normalized freshwater storage trends at the basin scale. High indicator scores are assigned to basins with co-occurring freshwater stress and drying trends. We refer to high freshwater status scores through status severity .

Vulnerability: The likelihood of society and ecosystems to experience harms due to exposure to freshwater stress and storage loss when considered together as a basin’s freshwater status. This vulnerability definition is an application of Turner et al.’s generic definition 29 . Vulnerability is quantified using social adaptability, ecological sensitivity, and basin freshwater status indicators. Social adaptability and ecological sensitivity indicators are described in the text and Methods.

Hotspot basin: Highlighted basins that possess the greatest vulnerability scores. We identify hotspot basins to support their prioritization in global water resources and integrated management initiatives. Basins are considered hotspots if sorted into “high” and “very high” vulnerability classes following a categorical classification of the numerical vulnerability results.

Hotspot basins found on all continents

We mapped the global gradient in social-ecological vulnerability to freshwater stress and storage loss at the basin scale and, from this, identified those with the greatest vulnerability as hotspot basins (Fig.  3 ). Hotspot mapping has been a successful endeavor within the field of conservation biogeography 31 , 32 , and many global hydrology studies have identified regions of exceptional water scarcity and security challenges e.g., 13 , 14 , 15 , 17 , 18 , 19 . Here, we seek to combine and apply these concepts in an integrated global social-ecological vulnerability context. As a useful reference, biodiversity hotspots aim to “maximize the number of species “saved” given available resources” by asking “where are places rich in species and under threat?” 33 . For comparison, the aim of our hotspot mapping is to ‘minimize the social and ecological impacts of freshwater stress and storage loss given available resources’ by asking “what basins with sensitive ecosystems and limited social adaptive capacity are exposed to freshwater stress and storage loss?”

a – d Social-ecological vulnerability results. a Hotspot basins of social-ecological vulnerability to freshwater stress and storage loss. b Vulnerability classification, based on the product of basin freshwater status and social-ecological sensitivity to freshwater stress and storage loss (see “Methods” section). c Histograms of the global distribution of vulnerability classes by basin count and surface area. d Summarized social-ecological activity within transitional and hotspot basins. e Ecological vulnerability results, presented as vulnerability classes. f Social vulnerability results, presented as vulnerability classes. Vulnerability classes for e and f are derived using the same methods as shown for social-ecological vulnerability in b .

We conceptualize vulnerability as the product of (i) ecological sensitivity, (ii) social adaptive capacity, and (iii) basin freshwater status. To represent ecological sensitivity, we derived an indicator using data products from two global ecohydrological studies that assess broad ecosystem sensitivity to freshwater storage and use (see “Methods” section). To represent social adaptability, we utilized the same adaptive capacity dataset as used in the previous section (Fig.  2a ). To classify the derived global vulnerability results into hotspot basins, we implemented a simple classification algorithm developed for heavy-tailed distributions 34 , which appropriately describes the global vulnerability distribution.

The most vulnerable basins are constrained to regions confronting co-occurring freshwater stress and storage loss. When considering social and ecological vulnerability individually (Fig.  3e, f ), we find spatial variation between ecological vulnerability (Fig.  3e ) and social vulnerability (Fig.  3f ). For instance, several basins in affluent nations with sensitive ecosystems reveal high ecological vulnerability but low social vulnerability (southwestern USA; western Australia). Conversely, several basins in Eastern Africa and northeastern India possess high social vulnerability but low to moderate ecological vulnerability. While these differences are notable and could impact regional strategies, it remains essential in most, if not all, regions that social and ecological vulnerabilities be confronted simultaneously 4 . For this purpose, we combined ecological sensitivity and adaptive capacity indicators into a combined social-ecological sensitivity indicator (see “Methods” section) to map combined social-ecological vulnerability (Fig.  3a ).

We identify 168 basins, representing 14% of all basins and 11% of the global land area considered in our study, as vulnerability hotspots (Fig.  3a–c ). These hotspot basins consist of basins receiving “high” and “very high” vulnerability scores through our classification procedure. Of the 168 basins, 78 (6% of all basins) are classified in the most-severe “very high” vulnerability class, while 90 (7% of all basins) are classified in the “high” vulnerability class. We also identified 232 basins (19% of all basins) as “transitional” basins, which are not classified alongside basins with null vulnerability yet also do not possess extreme values within the global vulnerability distribution. The 78 hotspot basins with “very high” vulnerability represent the multiple epicenters for potential social and ecosystem impacts from freshwater stress and storage loss. These basins are found in Argentina, northeastern Brazil, the American southwest, Mexico, Northern, Eastern, and Southern Africa, the Middle East and Arabian Peninsula, the Caucasus, West Asia, northern India and Pakistan, Southeastern Asia, and northern China.

A total of over 1.5 billion people, 17% of global food crop production, and 13% of global GDP are found within hotspot basins (Fig.  3d ). Of these, ~300 million people, 4% of global food crop production, and 4% of global GDP situate within the 78 “very high” vulnerability basins. Consistent with the relationship between biodiversity and basin freshwater status, we find the most vulnerable basins to be less taxonomically biodiverse than less vulnerable basins. While it is possible that these lower biodiversity levels may have eroded due to freshwater stress and storage loss, a proper investigation is outside the scope of this study and would require a wider array of pressures to be considered. The hotspot basins encompass 157 wetlands of international importance, which we highlight to prioritize their conservation in these vulnerable environments (Supplementary Table  2 ).

While the degree of social-ecological activity within hotspot basins is substantial, the global proportion of each dimension found in hotspot basins is roughly proportional to the fraction of basins within each vulnerability class. Thus, as the hotspot basins do not contribute disproportionately to global totals of social-ecological activity, we find it important to restate and clarify the motivating purpose of this hotspot mapping. The hotspot basins do not identify the greatest contributors to global social-ecological activity that face severe freshwater challenges. Rather, the hotspot basins are those with sensitive ecosystems and adaptability-limited societies exposed to the co-occurrence of freshwater stress and storage loss, and thus are the basins most likely to suffer social and ecological harms due to these freshwater conditions.

The identification of hotspot basins shows high levels of consistency across two uncertainty analyses and a sensitivity analysis focused on the impacts of subjective methodological decisions (Supplementary Section  4 ). We consider individually the impacts of (i) uniform over-estimation and under-estimation of each data input (spatially uniform uncertainty) and (ii) heterogeneous uncertainty in each data input (spatially variable uncertainty) on our hotspot basin results. Performing 10,000 realizations for each uncertainty analysis, we find that 98% of the identified transitional and hotspot basins are identified as at least transitional basins in over 50% the realizations considering spatially uniform uncertainty, and 96% when considering spatially variable uncertainty (Supplementary Figs.  8 and 9 ). The subjectivity-focused sensitivity analysis considered 24 alternative methodological configurations, and revealed that our identified transitional and hotspot basins are consistently identified across the majority of configurations (Supplementary Fig.  10 ).

Implementation of integrated water resources management is inconsistent across hotspot basins

We compared national implementation levels of integrated water resources management (IWRM) with our global vulnerability results (Fig.  4 ). For IWRM implementation data, we rely on the IWRM Data Portal 35 which tracks progress on SDG 6.5.1 (“IWRM implementation at the national scale”).

a Map of IWRM implementation overlaid by hotspot basin results. b Scatterplot of individual basin values of social-ecological vulnerability ( x -axis) and IWRM implementation ( y -axis). Transboundary basins are represented by concentric red circles, with the number of circles representing the number of nations present within each basin. See text for interpretation of labels 1, 2, and 3.

IWRM is defined as “a process which promotes the co-ordinated development and management of water, land and related resources, in order to maximize the resultant economic and social welfare in an equitable manner without compromising the sustainability of vital ecosystems” 36 , while the SDG framework notes that IWRM implementation “supports all Goals across the 2030 Agenda” 37 . Thus, as the IWRM paradigm seeks to guide management of water resources to minimize trade-offs between human well-being, ecological health, and water resources sustainability, assessing implementation levels of IWRM against our vulnerability results provides insight regarding the performance of IWRM globally while simultaneously emphasizing the broad sustainability implications within hotspot basins.

Globally, we find no direct relationship between vulnerability and IWRM implementation at the basin scale. There is thus a wide range of IWRM implementation across all levels of social-ecological vulnerability to freshwater stress and storage loss, and there is no indication that IWRM implementation levels are greatest where they are most needed. This finding likely derives from variations in proactive versus reactive governance and management approaches to freshwater challenges across the globe. As our analysis is conducted at a snapshot in time (input data align to ~2015), we can only generate hypotheses about the performance of IWRM globally. For example, basins with high levels of IWRM implementation and low vulnerability (label 1 in Fig.  4b ) have either proactively implemented IWRM, have effectively reduced their vulnerability through IWRM implementation, or simply benefit from a favorable intersection of regional climate and economy.

Alternatively, basins with low levels of IWRM and low vulnerability can be categorized as non-proactive in their IWRM implementation (label 2 in Fig.  4b ). We place particular emphasis here on basins with low levels of IWRM where vulnerability is high (label 3 in Fig.  4b ), which we argue should be the priority basins and regions of SDG 6.5-focused initiatives. Identified nations with low levels of IWRM implementation and very high vulnerability include Afghanistan, Algeria, Argentina, Egypt, India, Iraq, Kazakhstan, Mexico, Somalia, Ukraine, Uzbekistan, and Yemen. As one-third (36%) of all hotspot basins are transboundary (Fig.  4b ), improving basin-level IWRM implementation will require multilateralism and hydro-diplomacy and cannot be left to individual nations acting alone. Furthermore, we observe a lower level of IWRM implementation across hotspot basins that are transboundary versus non-transboundary hotspot basins (mean basin IWRM Data Portal score = 50 vs. 56), suggesting greater multilateralism and cooperation are needed in transboundary basins.

Implications for hotspot basins

There are many possible social and ecological implications for hotspot basins, however these depend on basin-specific relationships among freshwater stress, storage trends, ecological sensitivity, and adaptive capacity. Ecologically, hotspot basins are more likely to suffer from transgressed environmental flows 38 with ensuing impacts on freshwater and riparian ecosystems 39 . Freshwater storage loss is linked to increased drought frequency 40 , and falling water tables simultaneously decrease ecosystem resilience to drought by limiting root water-uptake 41 , harming groundwater-dependent ecosystems 42 , and perturbing the land surface energy balance 43 . Where freshwater stress and depletion are driven by irrigation, the ecological implications are not limited to within the basin, as irrigation can modify moisture recycling and precipitation patterns across precipitationsheds 44 , with cascading impacts on potentially distant ecosystems 45 . Socially, hotspot basins are forced to confront shrinking water resources to satisfy domestic, industrial, and agricultural demands 46 , 47 . Where water tables are dropping and wells are at risk of running dry 48 , groundwater access may intensify social inequalities as only the wealthy may be able to afford to drill deeper wells 49 . Such conditions may trigger conflict and have implications on international security 28 . Of the near-700 water conflicts documented since 2000 by The Water Conflict Chronology 50 , two-thirds (68%) are found within either transitional or hotspot basins from this analysis. As this study considers only a subset of water related stresses, this serves to indicate the potentially leading role freshwater stress and storage loss may play in instigating, contextualizing, and sustaining social conflict.

Opportunities for global ecohydrology and sociohydrology

This analysis is a representation of the best-available data for addressing global social-ecological impacts from freshwater stress and storage loss. However, it simultaneously highlights important limitations of existing ecohydrology and sociohydrology research at the global scale. For instance, we were required to develop an indicator to represent ecological sensitivity to freshwater from two existing ecohydrological datasets. However, these two datasets do not represent an exhaustive set of ecosystem processes and functions which may be impacted by changes in freshwater. We also do not address sub-grid variability considerations in either ecological sensitivity or social adaptability, which would independently benefit from specific study. Further, owing to a lack of alternatives, we adopted a rather reductionist and general derivation of social adaptive capacity as a proxy for social ability to respond and adapt to freshwater stress and storage loss. Future work that targets societal responses and relationships to specific freshwater stresses (e.g., Di Baldassarre et al. 51 but applied at the global scale and for various hydrological hazards) would be highly relevant for subsequent studies of this kind. We also only considered direct impacts of freshwater stress and storage loss, and thus did not consider indirect and non-local impacts such as the water and food security impacts of virtual water trade 52 , 53 , which does not yet have data available globally at a sub-national scale.

In closing, we seek to establish a thematic connection between this work and Abbott et al. 11 , who found human activity to be largely absent in water cycle diagrams and suggested this inaccuracy contributes to a “misunderstanding of global hydrology by policymakers, researchers and the public”. In a similar spirit, we observe that social-ecological system impacts of hydrological change are under-considered in global hydrological studies and we argue this underrepresentation contributes to a lack of awareness or misunderstanding of ecohydrological and sociohydrological connections at the global scale. Addressing hydrological phenomena simultaneously as products and drivers of change within the global social-ecological system can only elevate the consideration freshwater will receive in complex, multi-objective, multi-disciplinary decision making. Freshwater stress and storage trends are only two of several critical aspects of freshwater with broad social-ecological sustainability and resilience implications 3 . Pending data availability, similar analyses can be performed for seasonal and inter-annual variability in water storage and several quality considerations. Developing such a network of studies will support a more comprehensive understanding of the social-ecological resilience implications of global hydrological change.

Study approach

The overall study approach is summarized and illustrated in Supplementary Section  1 and Supplementary Figs.  1 and 2 . Below, we focus on the specific methods performed in our analysis. A flow chart of our methodology is shown in Supplementary Fig.  3 .

Data selection

All underlying data to this study were retrieved from pre-existing, published, and open datasets. We used 12 geospatial datasets and one basin scheme. We aligned our input data to the year 2015 as best as possible. We summarize input dataset selection, justification, and relevant preprocessing in Supplementary Table  1 (Supplementary Section  2 ).

Geospatial methods

We performed all analyses using the basin scheme of HydroBASINS level 4 54 and at the spatial resolution of 0.5 degrees using the World Geodetic Reference System 1984 ellipsoid (WGS 84). All raw input datasets were harmonized to 0.5 degrees using methods outlined in Supplementary Table  1 . All data were summarized to the basin scale by (i) calculating the area-weighted average value of intensive properties within each basin using an algorithm to calculate cell area on a geographic grid 55 , or (ii) calculating the within-basin sum of extensive properties. We masked basins in Greenland, northern Canada, and several small islands from our analysis due to inconsistent data coverage. This masking reduced the HydroBASINS level 4 discretization scheme to 1204 basins from an original set of 1341 (90% retention by count, 97% by surface area).

Hotspot basins

For a more extensive summary of the theory and justification of our approach, see Supplementary Section  1 . In brief, we base our analysis on the vulnerability definition of Turner et al. 29 although there is considerable consensus around the general principles of vulnerability 56 . Our approach is similar to that of Varis et al. 23 , a recent global assessment of river basin resilience using social-ecological principles. However, our study is of a narrower, more specific scope directed at the potential for social and ecological impacts from freshwater stress and storage loss rather than evaluating broad basin resilience to a wide range of ecological vulnerabilities.

We analyze social-ecological system vulnerability as the product of (i) exposure to freshwater stress and storage loss, represented by the basin freshwater status indicator, and (ii) the combination of ecological sensitivity and social adaptive capacity, which we represent using derived indicators.

Basin freshwater status is derived to represent exposure to the spatial co-occurrence and severity of freshwater stress and freshwater storage trends. This approach enables trends in freshwater storage to differentiate basins of equal freshwater stress as storage trends can aggravate or offset existing stress levels depending on the direction of the storage trend and the existing stress level. The derivation of this indicator is summarized in the “Basin freshwater status” section.

The social-ecological sensitivity indicator represents the ecological ability to either adapt or absorb freshwater stress or storage perturbations and social ability to adapt proactively and reactively to generic stresses 23 . The derivation of this social-ecological sensitivity indicator is summarized in the section below entitled “Ecosystem sensitivity and social adaptive capacity”. We note that other common quantitative vulnerability approaches incorporate additional considerations such as scaling the sensitivity term by the proximity of the system state to a critical threshold, and undertaking a probabilistic approach to exposure 57 . Such thresholds, however, are highly uncertain, spatially-variable, or unknown in most complex adaptive systems 58 , including human-water systems at the global scale 59 . Thus, we did not perform this threshold proximity scaling of the sensitivity term. Furthermore, as we evaluated the current state of freshwater stress and the existing trend in freshwater storage, there was no probabilistic component to our analysis. Vulnerability was thus represented in our analysis through Eq. ( 1 ).

where V represents vulnerability of the social-ecological system, S represents the combined social-ecological sensitivity indicator, and B represents the basin freshwater status, per basin i .

Basin freshwater status

We derived basin freshwater status to compress the bivariate relationship of freshwater stress and storage trends into a single indicator. The indicator is a composite of both freshwater stress and storage trend inputs which we individually normalized using the value of 0.4 times annual streamflow ( Q ) per basin. The composite indicator (Eq.  4 ) is the arithmetic mean of the normalized freshwater stress indicator (Eq.  2 ) and the normalized storage trend indicator (Eq.  3 ). Our freshwater stress calculation differs from other freshwater stress studies which consider water sharing rules between upstream and downstream basins 15 and treat arid regions separately 60 . As these existing freshwater stress dataset are not available at our operating scale of HydroBASINS level 4, we calculated the ratio of withdrawal to streamflow within each basin and interpret this ratio as an approximate, relative measure of freshwater demand to within-basin generated renewable freshwater. We normalized both freshwater stress and freshwater storage trends by 0.4Q as this threshold is used throughout the freshwater stress literature to denote high basin stress levels 13 , 15 , 61 , 62 . Following common approaches in other global indicator-based assessments, we bound both normalized indicators to a maximum magnitude of 1, i.e., we set an upper limit of 1 for the normalized freshwater stress results by setting all values >1 to 1, and set upper and lower limits of 1 and −1 for the normalized storage trends by setting all vales <−1 to −1 and >1 to 1. We also “flipped” the normalized storage trend indicator (i.e., multiplied the indicator by −1) so that drying trends drying trends correspond to positive indicator scores for consistency with the freshwater stress indicator for which greater (more positive) values correspond with greater levels of stress. Basin freshwater status was calculated as the arithmetic mean of these two indicators, with a minimum value set to 0 as negative values are possible where wetting trends offset existing freshwater stress. Where large earthquakes interfered with storage trend observations (i.e., the 2011 Tohoku earthquake and the 2004 Sumatra-Andaman earthquake) 9 , basin freshwater status was set to the independent normalized freshwater stress indicator alone. Input data for these indicators are shown in Supplementary Fig.  4 .

where F is the freshwater stress indicator, T is the normalized freshwater storage trend indicator, and B is basin freshwater status. W represents annual freshwater withdrawals (mm year −1 ), Q represents annual streamflow (mm year −1 ), and \(\frac{{{{{{\rm{d}}}}}}TWS}{{{{{{\rm{d}}}}}}t}\) represents year-over-year trends in freshwater storage (mm year −1 ), per basin i .

Ecosystem sensitivity and social adaptive capacity

We derived an indicator to represent general ecological sensitivity to freshwater stress and storage loss as no existing dataset fit this use. We combined data products from de Graaf et al. 63 and Seddon et al. 64 , which represent the most relevant global ecohydrological studies. de Graaf et al. used a global surface water-groundwater model to estimate the groundwater head decline at which environmental flow limits are transgressed for all basins in which there is currently groundwater pumping. Seddon et al. developed the Vegetation Sensitivity Index to quantify the sensitivity in vegetation productivity to anomalies in three climate variables: water, temperature, and cloudiness, where water anomalies were represented by the ratio of actual evapotranspiration to potential transpiration. We incorporated only the water-specific component of the Vegetation Sensitivity Index. In brief, the de Graaf et al. dataset represents the sensitivity of environmental flows to changes in groundwater storage in basins where groundwater is currently being withdrawn, while the Seddon et al. dataset represents vegetation sensitivity to anomalies in soil moisture and shallow groundwater storage.

To combine these ecohydrological datasets of different dimensions which simultaneously contribute to our overall understanding of ecological sensitivity to changes in freshwater storage, we performed a purely statistical, percentile-based approach. At the gridded resolution of 0.5 degrees, we converted both input datasets into area-weighted percentile datasets, where a value of 1 represents the grid cells with the global maximum values (99th–100th percentile) per dataset and values of 0.01 represent the grid cells with the global minimum values (0th–1st percentile) per dataset and where all values in-between apply to an equal proportion of the land surface. We then (i) computed the average value of both percentile-transformed datasets within each basin, (ii) averaged the two independent averages from (i) to combine the two datasets, and (iii) then normalized all basins by the global maximum basin value. This approach ensured that the most-sensitive basin received an ecosystem sensitivity indicator score of 1 and the least-sensitive basins received scores near 0. It is readily acknowledged, however, that a process-based derivation of ecosystem sensitivity that integrates groundwater, streamflow, and soil moisture considerations would be a superior alternative however none exist to our knowledge. The input data and derivation of the ecosystem sensitivity indicator are shown in Supplementary Fig.  5 .

To reconcile the ecological sensitivity indicator (where higher values correspond with more sensitive basins) with the concept of social adaptive capacity, we inverted the dataset of social adaptive capacity so that greater values corresponded with lower adaptability. We sourced the adaptive capacity dataset from Varis et al. 23 , as outlined in the main text. Social adaptive capacity is shown in Supplementary Fig.  6 .

We combined ecosystem sensitivity and adaptive capacity indicators using the fuzzy sum operation (Eq.  4 ). The fuzzy sum is an increasing linear combination operator that ensures the sum is no less than the largest input value, yet the contribution of subsequent inputs decrease as inputs overlap. If all inputs are normalized [0, 1], the fuzzy sum converges on the upper limit (i.e., 1). The fuzzy sum is increasingly used in geospatial applications, such as in landscape integrity mapping 65 and we used the fuzzy sum as it eliminates the need to weight the ecosystem sensitivity and adaptive capacity inputs, which would introduce an additional degree of subjectivity. The fuzzy sum operation to derive the social-ecological sensitivity indicator is shown in Eq. ( 5 ). The resulting social-ecological sensitivity indicator and its input datasets are shown in Supplementary Fig.  7 .

where S represents the social-ecological sensitivity indicator, E represents ecological sensitivity, and A represents social adaptive capacity, all per basin i . Note that the (1–A) term represents the inversion of the adaptive capacity dataset.

Hotspot identification

As our vulnerability analysis yields a relative gradient in global social-ecological vulnerability, classifying individual basins into vulnerability categories presents a particular challenge as no process-based thresholds are identified in the literature to differentiate the results. For example, while other prominent hotspot mapping initiatives such as the biodiversity hotspots are based on strict criteria (e.g., the biodiversity hotspots must contain \({\ge}\!\) 0.5% or \({\ge}\!\) 1500 of the world’s plant species and must have lost \({\ge}\!\) 70% or more of its primary vegetation 32 ), our social-ecological vulnerability hotspots are deeply interdisciplinary and multivariate, and thus make such explicit criteria challenging to identify. As vulnerability is not directly observable, the pragmatic approach is often to measure and classify relative vulnerability 57 .

Given the heavy-tailed distribution of vulnerability results, we selected a relative classification algorithm, the Head/Tail Breaks method 34 , to categorize the basins into vulnerability classes and subsequently hotspot basins. The Head/Tail Breaks classification scheme was developed to better represent the hierarchical structure of heavy-tailed distributions compared to other common methods such as Jenks natural breaks optimization. The classification scheme partitions the data into “head” and “tail” classes based on the arithmetic mean of the distribution and recursively re-partitions the “head” class based on the arithmetic mean of the “head” values until a skewness threshold is reached. For consistency between dimensions, we applied three iterations of the algorithm to all vulnerability distributions (i.e., social, ecological, and social-ecological) to partition the data consistently into four classes rather than use a skew-based stopping criterion. We classified the lowest-level class as non-hotspots (low vulnerability), the 2nd-level class as “transitional” basins (moderate vulnerability), and the 3rd-level and 4th-level classes as hotspots (high and very high vulnerability). As vulnerability is derived as the product of basin freshwater status and social ecological sensitivity, the class breaks can be represented by reciprocal function curves as shown in Fig.  3b (labeled as “Curves of equal vulnerability”). While implementing theory-based thresholds of vulnerability would be preferred, it is not realistic given current data availability and process knowledge of ecohydrological and sociohydrological systems at the global scale. We justify our approach as a data-driven, natural classification that characterizes the global gradient in social-ecological vulnerability to freshwater stress and storage loss.

Uncertainty and sensitivity analyses

We performed two uncertainty analyses to explore the impact on our hotspot mapping from potential input data uncertainty, and we performed one sensitivity analysis to identify the impact of subjective decision making in our methodology. These analyses are presented and described in Supplementary Section  4 and Supplementary Figs.  8 – 10 .

Data availability

The basin vulnerability and hotspot basin classification results from this study have been deposited in the University of Victoria’s Scholars Portal Dataverse ( https://doi.org/10.5683/SP3/SLR3GF ). The raw data used in this study are accessible through several data and code repositories. We provide data sources, persistent web links, and descriptions of all data used in this study in Supplementary Table  1 .

Code availability

All analyses were conducted using the R project for statistical computing 66 . R packages necessary for analysis and visualization include: raster 67 , sf 68 , gdalUtils 69 , spatstat.geom 70 , ggplot2 71 , tmap 72 , and scico 73 , 74 , 75 . Composite figures were assembled in Affinity Designer ( https://affinity.serif.com/en-us/designer/ ). The code generated to conduct this study is available in the GitHub repository https://github.com/XanderHuggins/HotspotBasins and is archived on Zenodo ( https://doi.org/10.5281/zenodo.5728475 ) 76 .

Falkenmark, M. Water and mankind: a complex system of mutual interaction. Ambio 6 , 3–9 (1977).

Google Scholar  

Falkenmark, M., Wang-Erlandsson, L. & Rockström, J. Understanding of water resilience in the anthropocene. J. Hydrol. X 2 , 100009 (2019).

Article   Google Scholar  

Falkenmark, M. & Wang-Erlandsson, L. A water-function-based framework for understanding and governing water resilience in the Anthropocene. One Earth 4 , 213–225 (2021).

Folke, C., Biggs, R., Norström, A. V., Reyers, B. & Rockström, J. Social-ecological resilience and biosphere-based sustainability science. Ecol. Soc. 21 , art41 (2016).

Ostrom, E. A general framework for analyzing sustainability of socio-ecological systems. Science 325 , 419–422 (2009).

Article   ADS   MathSciNet   CAS   PubMed   MATH   Google Scholar  

Steffen, W. et al. Planetary boundaries: guiding human development on a changing planet. Science 347 , 1259855 (2015).

Article   PubMed   Google Scholar  

Gleeson, T. et al. The water planetary boundary: interrogation and revision. One Earth 2 , 223–234 (2020).

Rockström, J. et al. The unfolding water drama in the anthropocene: towards a resilience‐based perspective on water for global sustainability. Ecohydrology 7 , 1249–1261 (2014).

Rodell, M. et al. Emerging trends in global freshwater availability. Nature 557 , 651–659 (2018).

Article   ADS   CAS   PubMed   PubMed Central   Google Scholar  

Haddeland, I. et al. Global water resources affected by human interventions and climate change. Proc. Natl Acad. Sci. 111 , 3251–3256 (2014).

Article   ADS   CAS   PubMed   Google Scholar  

Abbott, B. W. et al. Human domination of the global water cycle absent from depictions and perceptions. Nature Geoscience 12 , 533–540 (2019).

Article   ADS   CAS   Google Scholar  

Bierkens, M. F. P. Global hydrology 2015: state, trends, and directions. Water Resour. Res. 51 , 4923–4947 (2015).

Article   ADS   Google Scholar  

Vörösmarty, C. J. Global water resources: vulnerability from climate change and population growth. Science 289 , 284–288 (2000).

Article   ADS   PubMed   Google Scholar  

Kummu, M., Ward, P. J., de Moel, H. & Varis, O. Is physical water scarcity a new phenomenon? Global assessment of water shortage over the last two millennia. Environ. Res. Lett. 5 , 034006 (2010).

Kummu, M. et al. The world’s road to water scarcity: shortage and stress in the 20th century and pathways towards sustainability. Sci. Rep. 6 , 38495 (2016).

Wada, Y., Wisser, D. & Bierkens, M. F. P. Global modeling of withdrawal, allocation and consumptive use of surface water and groundwater resources. Earth Syst. Dyn. 5 , 15–40 (2014).

Padowski, J. C., Gorelick, S. M., Thompson, B. H., Rozelle, S. & Fendorf, S. Assessment of human-natural system characteristics influencing global freshwater supply vulnerability. Environ. Res. Lett. 10 , 104014 (2015).

Gain, A. K., Giupponi, C. & Wada, Y. Measuring global water security towards sustainable development goals. Environ. Res. Lett. 11 , 124015 (2016).

Vörösmarty, C. J. et al. Global threats to human water security and river biodiversity. Nature 467 , 555–561 (2010).

Ohlsson, L. Water conflicts and social resource scarcity. Phys. Chem. Earth Part B 25 , 213–220 (2000).

Sullivan, C. A. Calculating a water poverty index. World Dev. 30 , 1195–1210 (2002).

Chapin, F. S., Folke, C. & Kofinas, G. P. In Principles of Ecosystem Stewardship: Resilience-based Natural Resource Management in a Changing World (eds. Folke, C. et al.) 3–28 (Springer, 2009).

Varis, O., Taka, M. & Kummu, M. The planet’s stressed river basins: too much pressure or too little adaptive capacity? Earth’s Future 7 , 1118–1135 (2019).

Article   ADS   PubMed   PubMed Central   Google Scholar  

Taylor, R. G. et al. Ground water and climate change. Nat. Clim. Change 3 , 322–329 (2013).

Vishwakarma, B. D., Devaraju, B. & Sneeuw, N. What is the spatial resolution of GRACE satellite products for hydrology? Remote Sens. 10 , 852 (2018).

Famiglietti, J. S. The global groundwater crisis. Nat. Clim. Change 4 , 945–948 (2014).

Su, G. et al. Human impacts on global freshwater fish biodiversity. Science 371 , 835–838 (2021).

Blowes, S. A. et al. The geography of biodiversity change in marine and terrestrial assemblages. Science 366 , 339–345 (2019).

Turner, B. L. et al. A framework for vulnerability analysis in sustainability science. Proc. Natl Acad. Sci. USA 100 , 8074–8079 (2003).

Dolan, F. et al. Evaluating the economic impact of water scarcity in a changing world. Nat. Commun. 12 , 1915 (2021).

Possingham, H. P. & Wilson, K. A. Turning up the heat on hotspots. Nature 436 , 919–920 (2005).

Myers, N., Mittermeier, R. A., Mittermeier, C. G., da Fonseca, G. A. B. & Kent, J. Biodiversity hotspots for conservation priorities. Nature 403 , 853–858 (2000).

Whittaker, R. J. et al. Conservation biogeography: assessment and prospect. Divers. Distrib. 11 , 3–23 (2005).

Jiang, B. Head/tail breaks: a new classification scheme for data with a heavy-tailed distribution. Prof. Geogr. 65 , 482–494 (2013).

UNEP & UNEP-DHI. IWRM Data Portal. http://iwrmdataportal.unepdhi.org/ (2021).

Hassing, J., Ipsen, N., Clausen, T. J., Larsen, H. & Lindgaard-Jørgensen, P. Integrated Water Resources Management in Action (UNESCO, 2009).

UNEP. Progress on integrated water resources management. Global baseline for SDG 6Indicator 6.5.1: degree of IWRM implementation . https://www.unwater.org/app/uploads/2018/11/SDG6_Indicator_Report_651_Progress-on-Integrated-Water-Resources-Management_ENGLISH_2018.pdf (2018).

Richter, B. D., Davis, M. M., Apse, C. & Konrad, C. A presumptive standard for environmental flow protection. River Res. Applic. 28 , 1312–1321 (2012).

Poff, N. L. & Zimmerman, J. K. H. Ecological responses to altered flow regimes: a literature review to inform the science and management of environmental flows. Freshw. Biol. 55 , 194–205 (2010).

Pokhrel, Y. et al. Global terrestrial water storage and drought severity under climate change. Nat. Clim. Change 11 , 226–233 (2021).

Fan, Y., Miguez-Macho, G., Jobbágy, E. G., Jackson, R. B. & Otero-Casal, C. Hydrologic regulation of plant rooting depth. PNAS 114 , 10572–10577 (2017).

Eamus, D., Zolfaghar, S., Villalobos-Vega, R., Cleverly, J. & Huete, A. Groundwater-dependent ecosystems: recent insights from satellite and field-based studies. Hydrol. Earth Syst. Sci. 19 , 4229–4256 (2015).

Kollet, S. J. & Maxwell, R. M. Capturing the influence of groundwater dynamics on land surface processes using an integrated, distributed watershed model. Water Resour. Res. (2008).

Keys, P. W. et al. Invisible water security: moisture recycling and water resilience. Water Security 8 , 100046 (2019).

Article   PubMed   PubMed Central   Google Scholar  

Keys, P. W. & Wang-Erlandsson, L. On the social dynamics of moisture recycling. Earth System. Earth Syst. Dyn. 9 , 829–847 (2018).

Falkenmark, M., Lundqvist, J. & Widstrand, C. Macro-scale water scarcity requires micro-scale approaches. Nat. Resour. Forum 13 , 258–267 (1989).

Article   CAS   PubMed   Google Scholar  

Gleick, P. H. Basic water requirements for human activities: meeting basic needs. Water Int. 21 , 83–92 (1996).

Jasechko, S. & Perrone, D. Global groundwater wells at risk of running dry. Science 372 , 418–421 (2021).

Perrone, D. Groundwater overreliance leaves farmers and households high and dry. One Earth 2 , 214–217 (2020).

Pacific Institute. Water Conflict Chronology (Island Press, 2019).

Di Baldassarre, G. et al. Socio-hydrology: conceptualising human-flood interactions. Hydrol. Earth Syst. Sci. 17 , 3295–3303 (2013).

Dalin, C., Konar, M., Hanasaki, N., Rinaldo, A. & Rodriguez-Iturbe, I. Evolution of the global virtual water trade network. PNAS 109 , 5989–5994 (2012).

Rosa, L., Chiarelli, D. D., Tu, C., Rulli, M. C. & D’Odorico, P. Global unsustainable virtual water flows in agricultural trade. Environ. Res. Lett. 14 , 114001 (2019).

Lehner, B. & Grill, G. Global river hydrography and network routing: baseline data and new approaches to study the world’s large river systems. Hydrol. Process. 27 , 2171–2186 (2013).

Santini, M., Taramelli, A. & Sorichetta, A. ASPHAA: a GIS-based algorithm to calculate cell area on a latitude-longitude (Geographic) regular grid. Trans. GIS 14 , 351–377 (2010).

Lei, Y., Wang, J., Yue, Y., Zhou, H. & Yin, W. Rethinking the relationships of vulnerability, resilience, and adaptation from a disaster risk perspective. Nat. Hazards 70 , 609–627 (2014).

Luers, A. L., Lobell, D. B., Sklar, L. S., Addams, C. L. & Matson, P. A. A method for quantifying vulnerability, applied to the agricultural system of the Yaqui Valley, Mexico. Glob. Environ. Change 13 , 255–267 (2003).

Scheffer, M. et al. Early-warning signals for critical transitions. Nature 461 , 53–59 (2009).

Xu, L., Mao, F., Famiglietti, J. S., Pomeroy, J. W. & Pahl-Wostl, C. Conceptualizing Cascading Effects of Resilience in Human–Water Systems. in Multisystemic resilience: Adaptation and Transformation in Contexts of Change 744–767 (Oxford University Press, 2021).

Hofste, R. W. et al. Aqueduct 3.0: Updated Decision-Relevant Global Water Risk Indicators (World Resources Institute, 2019).

Falkenmark, M. & Lundqvist, J. Towards water security: political determination and human adaptation crucial. Nat. Resour. Forum 22 , 37–51 (1998).

Smakhtin, V., Revenga, C. & Döll, P. A pilot global assessment of environmental water requirements and scarcity. Water Int. 29 , 307–317 (2004).

Article   CAS   Google Scholar  

de Graaf, I. E. M., Gleeson, T., van Beek, L. P. H. R.Sutanudjaja, E. H. & Bierkens, M. F. P. Environmental flow limits to global groundwater pumping. Nature 574 , 90–94 (2019).

Seddon, A. W. R., Macias-Fauria, M., Long, P. R., Benz, D. & Willis, K. J. Sensitivity of global terrestrial ecosystems to climate variability. Nature 531 , 229–232 (2016).

Kennedy, C. M., Oakleaf, J. R., Theobald, D. M., Baruch-Mordo, S. & Kiesecker, J. Managing the middle: a shift in conservation priorities based on the global human modification gradient. Glob. Change Biol. 25 , 811–826 (2019).

R Core Team. R: a language and environment for statistical computing. Version 4.1.1 https://www.r-project.org/ (2021).

Hijmans, R. J. raster: geographic data analysis and modeling. R package version 3.5-2. https://cran.r-project.org/package=raster (2021).

Pebesma, E. et al. sf: simple features for R. R package version 1.0–3. https://cran.r-project.org/package=sf (2021).

Greenberg, J. A. & Mattiuzzi, M. gdalUtils: wrappers for the geospatial data abstraction library (GDAL) utilities. R package version 2.0–3. https://cran.r-project.org/package=gdalUtils (2020).

Baddeley, A., Turner, R. & Rubak, E. spatstat.geom: geometric functionality of the ‘spatstat’ family. R package version 2.3-0. https://cran.r-project.org/package=spatstat.geom (2021).

Wickham, H. ggplot2: elegant graphics for data analysis. R package version 3.3–5. https://cran.r-project.org/package=ggplot2 (2021).

Tennekes, M. et al. tmap: thematic maps. R package version 3.3-2. https://cran.r-project.org/package=tmap (2021).

Crameri, F., Shephard, G. E. & Heron, P. J. The misuse of colour in science communication. Nat. Commun. 11 , 5444 (2020).

Crameri, F. Scientific colour maps (Version 7.0.0). Zenodo https://doi.org/10.5281/zenodo.4491293 (2021).

Pedersen, T. L. scico: colour palettes based on the scientific colour-maps. R package version 1.2-0. https://cran.r-project.org/package=scico (2020).

Huggins, X. et al. Code for: hotspots for social and ecological impacts from freshwater stress and storage loss. Zenodo. https://doi.org/10.5281/zenodo.5728475 (2021).

Download references

Acknowledgements

We would like to thank colleagues at the Groundwater Science and Sustainability Research Group at the University of Victoria, and the Global Institute for Water Security at the University of Saskatchewan for reading and discussing earlier versions of this manuscript. X.H. was supported by the Natural Sciences and Engineering Research Council (NSERC) of Canada through an Alexander Graham Bell Canada Graduate Scholarship.

Author information

Authors and affiliations.

Department of Civil Engineering, University of Victoria, Victoria, Canada

Xander Huggins, Tom Gleeson & Tara J. Troy

Global Institute for Water Security, University of Saskatchewan, Saskatoon, Canada

Xander Huggins & James S. Famiglietti

School of Earth and Ocean Sciences, University of Victoria, Victoria, Canada

Tom Gleeson

Water and Development Research Group, Aalto University, Espoo, Finland

Matti Kummu

Kansas Geological Survey, University of Kansas, Lawrence, KS, USA

Samuel C. Zipper

International Institute for Applied Systems Analysis, Laxenburg, Austria

Yoshihide Wada

School of Environment and Sustainability, University of Saskatchewan, Saskatoon, Canada

James S. Famiglietti

You can also search for this author in PubMed   Google Scholar

Contributions

X.H. conceived the study with advice from T.G., J.S.F., M.K., S.C.Z., Y.W. and T.J.T. X.H. assembled input data, performed the analysis, and wrote the manuscript with greatest input from T.G. and J.S.F. All authors, X.H., T.G., J.S.F., M.K., S.C.Z., Y.W. and T.J.T., discussed results and edited the manuscript at all stages.

Corresponding author

Correspondence to Tom Gleeson .

Ethics declarations

Competing interests.

The authors declare no competing interests.

Peer review

Peer review information.

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

Additional information

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

Supplementary information

Supplementary information, peer review file, rights and permissions.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/ .

Reprints and permissions

About this article

Cite this article.

Huggins, X., Gleeson, T., Kummu, M. et al. Hotspots for social and ecological impacts from freshwater stress and storage loss. Nat Commun 13 , 439 (2022). https://doi.org/10.1038/s41467-022-28029-w

Download citation

Received : 21 May 2021

Accepted : 22 December 2021

Published : 21 January 2022

DOI : https://doi.org/10.1038/s41467-022-28029-w

Share this article

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

This article is cited by

A new high-resolution global topographic factor dataset calculated based on srtm.

• Hongming Zhang
• Linlin Yuan

Scientific Data (2024)

Mapping the diversity of land uses following deforestation across Africa

• Robert N. Masolele
• Diego Marcos
• Martin Herold

Scientific Reports (2024)

• Shailesh Kumar Singh

Communications Earth & Environment (2023)

Adaptive Water Management: On the Need for Using the Post-WWII Science in Water Governance

• Peder Hjorth
• Kaveh Madani

Water Resources Management (2023)

UiO-66 octahedrons for adsorptive removal of direct blue-6: process optimization, interaction mechanism, and phytotoxicity assessment

• Harish Kumar Rajendran
• Lingaraj Sahoo

Environmental Science and Pollution Research (2023)

By submitting a comment you agree to abide by our Terms and Community Guidelines . If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Quick links

• Explore articles by subject
• Guide to authors
• Editorial policies

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.