# Decline in Iran's groundwater recharge

Roohollah Noori, Mohsen Maghrebi, Søren Jessen, Sayed Bateni, Essam Heggy, Saman Javadi, Mojtaba Noury, Séverin Pistre, Soroush Abolfathi, Amir Aghakouchak

► To cite this version:

Roohollah Noori, Mohsen Maghrebi, Søren Jessen, Sayed Bateni, Essam Heggy, et al., Decline in Iran's groundwater recharge. Nature Communications, 2023, 14 (1), pp.6674. 10.1038/s41467-023-42411-2. hal-04564738

HAL Id: hal-04564738

https://hal.science/hal-04564738

Submitted on 5 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Decline in Iran's groundwater recharge

Roohollah Noori, Mohsen Maghrebi, Søren Jessen, Sayed Bateni, Essam Heggy, Saman Javadi, Mojtaba Noury, Séverin Pistre, Soroush Abolfathi, Amir Aghakouchak

► To cite this version:

Roohollah Noori, Mohsen Maghrebi, Søren Jessen, Sayed Bateni, Essam Heggy, et al., Decline in Iran's groundwater recharge. Nature Communications, 2023, 14 (1), pp.6674. 10.1038/s41467-023-42411-2. hal-04564738

HAL Id: hal-04564738

https://hal.science/hal-04564738v1

Submitted on 5 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

# Decline in Iran's groundwater recharge

Received: 20 February 2023

Accepted: 11 October 2023

Published online: 21 October 2023

<figure><img src="image_3.png" /></figure>

Check for updates

Roohollah Noori Ⓓ1,2,*, Mohsen Maghrebi<fnref n="1" />, Søren Jessen Ⓒ3, Sayed M. Bateni⁴,
Essam Heggy Ⓒ5,6, Saman Javadi⁷, Mojtaba Noury Ⓒ8, Severin Pistre Ⓒ9,
Soroush Abolfathi Ⓒ10 & Amir AghaKouchak Ⓒ11,12,13

Groundwater recharge feeds aquifers supplying fresh-water to a population over 80 million in Iran—a global hotspot for groundwater depletion. Using an extended database comprising abstractions from over one million ground-water wells, springs, and qanats, from 2002 to 2017, here we show a significant decline of around -3.8 mm/yr in the nationwide groundwater recharge. This decline is primarily attributed to unsustainable water and environmental resources management, exacerbated by decadal changes in climatic conditions. However, it is important to note that the former's contribution outweighs the latter. Our results show the average annual amount of nationwide groundwater recharge (i.e., ~40 mm/yr) is more than the reported average annual runoff in Iran (i.e., ~32 mm/yr), suggesting the surface water is the main contributor to groundwater recharge. Such a decline in groundwater recharge could further exacerbate the already dire aquifer depletion situation in Iran, with devastating consequences for the country's natural environment and socio-economic development.

Securing groundwater is essential since it globally supplies approximately half of the irrigation water and most of the domestic water demands<fnref n="1" />. Given that 60% of the global population resides in regions that are expected to experience water scarcity and reduction in clean surface water availability by 2050², reliance on groundwater resources will inevitably further increase in the future. However, despite the importance of groundwater as one of the pillars of food and water security, extractions from aquifer systems exceed their natural recharge, leading to the rapid depletion of aquifers across the world³⁻⁵. Therefore, re-establishing the balance between the amount of groundwater withdrawal and recharge is essential to sustainably use groundwater resources, especially in (semi)arid regional scales such as Iran, which is known as a global hotspot for aquifers depletion⁶⁻⁸.

The imbalance between groundwater withdrawal and recharge (i.e., groundwater depletion) was first reported in some aquifers in Iran

in the 1970s. However, unpublished national observations revealed groundwater depletion in some plains from as early as the 1950s. This coincided with the gradual replacement of Persian qanats, which were sustainable groundwater extraction systems and UNESCO World Cultural Heritage sites⁹, with (semi)deep wells. This shift was a consequence of the modernization in the agricultural sector⁸. Currently, Iran records the volume of groundwater withdrawal from over one million extraction points, including wells, springs, and qanats, through its national groundwater monitoring network that covers all 30 national hydrological basins (Fig. S1). However, obtaining ground-truth data on nationwide groundwater recharge remains challenging due to the site-specific nature of measurements and the complex interaction between various hydroclimatic and physiographic factors. Factors such as precipitation, evapotranspiration, topography, depth to the groundwater table, geology, and land use/cover (LUC) play a crucial

<footnote n="1">¹Graduate Faculty of Environment, University of Tehran, Tehran, Iran. ²Faculty of Governance, University of Tehran, Tehran, Iran. ³Department of Geosciences and Natural Resource Management, University of Copenhagen, Copenhagen, Denmark. ⁴Department of Civil and Environmental Engineering and Water Resources Research Center, University of Hawaii at Manoa, Honolulu, HI, USA. ⁵Viterbi School of Engineering, University of Southern California, Los Angeles, CA, USA. ⁶Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA. ⁷Department of Water Engineering, College of Abouraihan, University of Tehran, Tehran, Iran. ⁸Iran Water Resources Management Company, Ministry of Energy, Tehran, Iran. ⁹HydroSciences Montpellier, University of Montpellier, CNRS, IRD, Montpellier, France. ¹⁰School of Engineering, University of Warwick, Coventry, UK. ¹¹Department of Civil and Environmental Engineering, University of California, Irvine, USA. ¹²Department of Earth System Science, University of California, Irvine, USA. ¹³Institute for Water, Environment and Health, United Nations University, Hamilton, ON, Canada. ✉e-mail: noor@ut.ac.ir</footnote>
<figure><img src="image_2.png" /></figure>

# Decline in Iran's groundwater recharge

Received: 20 February 2023

Accepted: 11 October 2023

Published online: 21 October 2023

<figure><img src="image_3.png" /></figure>

Check for updates

Roohollah Noori Ⓓ1,2,*, Mohsen Maghrebi<fnref n="1" />, Søren Jessen Ⓒ3, Sayed M. Bateni⁴,
Essam Heggy Ⓒ5,6, Saman Javadi⁷, Mojtaba Noury Ⓒ8, Severin Pistre Ⓒ9,
Soroush Abolfathi Ⓒ10 & Amir AghaKouchak Ⓒ11,12,13

Groundwater recharge feeds aquifers supplying fresh-water to a population over 80 million in Iran—a global hotspot for groundwater depletion. Using an extended database comprising abstractions from over one million ground-water wells, springs, and qanats, from 2002 to 2017, here we show a significant decline of around -3.8 mm/yr in the nationwide groundwater recharge. This decline is primarily attributed to unsustainable water and environmental resources management, exacerbated by decadal changes in climatic conditions. However, it is important to note that the former's contribution outweighs the latter. Our results show the average annual amount of nationwide groundwater recharge (i.e., ~40 mm/yr) is more than the reported average annual runoff in Iran (i.e., ~32 mm/yr), suggesting the surface water is the main contributor to groundwater recharge. Such a decline in groundwater recharge could further exacerbate the already dire aquifer depletion situation in Iran, with devastating consequences for the country's natural environment and socio-economic development.

Securing groundwater is essential since it globally supplies approximately half of the irrigation water and most of the domestic water demands<fnref n="1" />. Given that 60% of the global population resides in regions that are expected to experience water scarcity and reduction in clean surface water availability by 2050², reliance on groundwater resources will inevitably further increase in the future. However, despite the importance of groundwater as one of the pillars of food and water security, extractions from aquifer systems exceed their natural recharge, leading to the rapid depletion of aquifers across the world³⁻⁵. Therefore, re-establishing the balance between the amount of groundwater withdrawal and recharge is essential to sustainably use groundwater resources, especially in (semi)arid regional scales such as Iran, which is known as a global hotspot for aquifers depletion⁶⁻⁸.

The imbalance between groundwater withdrawal and recharge (i.e., groundwater depletion) was first reported in some aquifers in Iran

in the 1970s. However, unpublished national observations revealed groundwater depletion in some plains from as early as the 1950s. This coincided with the gradual replacement of Persian qanats, which were sustainable groundwater extraction systems and UNESCO World Cultural Heritage sites⁹, with (semi)deep wells. This shift was a consequence of the modernization in the agricultural sector⁸. Currently, Iran records the volume of groundwater withdrawal from over one million extraction points, including wells, springs, and qanats, through its national groundwater monitoring network that covers all 30 national hydrological basins (Fig. S1). However, obtaining ground-truth data on nationwide groundwater recharge remains challenging due to the site-specific nature of measurements and the complex interaction between various hydroclimatic and physiographic factors. Factors such as precipitation, evapotranspiration, topography, depth to the groundwater table, geology, and land use/cover (LUC) play a crucial

<footnote n="1">¹Graduate Faculty of Environment, University of Tehran, Tehran, Iran. ²Faculty of Governance, University of Tehran, Tehran, Iran. ³Department of Geosciences and Natural Resource Management, University of Copenhagen, Copenhagen, Denmark. ⁴Department of Civil and Environmental Engineering and Water Resources Research Center, University of Hawaii at Manoa, Honolulu, HI, USA. ⁵Viterbi School of Engineering, University of Southern California, Los Angeles, CA, USA. ⁶Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA. ⁷Department of Water Engineering, College of Abouraihan, University of Tehran, Tehran, Iran. ⁸Iran Water Resources Management Company, Ministry of Energy, Tehran, Iran. ⁹HydroSciences Montpellier, University of Montpellier, CNRS, IRD, Montpellier, France. ¹⁰School of Engineering, University of Warwick, Coventry, UK. ¹¹Department of Civil and Environmental Engineering, University of California, Irvine, USA. ¹²Department of Earth System Science, University of California, Irvine, USA. ¹³Institute for Water, Environment and Health, United Nations University, Hamilton, ON, Canada. ✉e-mail: noor@ut.ac.ir</footnote>
role in this process¹⁰,¹¹. In fact, spatial and temporal variability of groundwater recharge, as an index of groundwater availability¹², is not accessible even locally in most parts of the world, including Iran. This scarcity of data indicates how rarely this variable is monitored. A possible way to obtain estimates of groundwater recharge over a large spatial scale is using global/regional hydrologic models. However, the recharge estimated by these models is not certain, especially in semi-arid regions such as Iran, where uncertainties in evapotranspiration and precipitation may lead to large errors in the results¹³,¹⁴, hindering the detection of regional imbalances between groundwater recharge and human withdrawals⁸.

Following the study conducted by Noori et al.⁸ on Iran's groundwater depletion, which raised concerns about the lack of the country's groundwater recharge data, we estimate its recharge using an extensive, nationwide database in all 30 hydrological basins of Iran from 2002 to 2017. In this study period, the country's groundwater resources were largely depleted due to anthropogenic and natural drivers. Our results suggest a significant decline in the nationwide groundwater recharge, primarily attributed to unsustainable development practices. By elucidating the anthropogenic and natural drivers contributing to this decline in Iran's groundwater recharge as a representative of a (semi)arid region in the world, we aim to foster a more comprehensive understanding of the regional and global impacts of these large-scale depletions. Our findings advocate for the urgent need for transformative changes in groundwater management and policies, especially in regions experiencing increasing hydroclimatic fluctuations.

# Results and discussion

## Estimated groundwater recharge

The estimated groundwater recharges for each aquifer, totaling 666 aquifers, aggregated in Iran's hydrological basins from 2002 to 2017, are depicted in Fig. 1. The Method's section details how these values were obtained. Our findings reveal a statistically significant decrease in nationwide groundwater recharge of about 35% (p value ≤ 0.001) during the study period, which has resulted in a groundwater imbalance in many basins (Fig. 1). The annual change in groundwater recharge varies from -10.3 mm/yr in the Great Karoon basin to +1.9 mm/yr in the Anzali basin, averaging about -3.8 mm/yr across the country (Table 1). The average groundwater recharge across the country and for the study period is estimated at 39.6 mm/yr, which is slightly above the upper end of the range globally reported for (semi) arid areas (i.e., 0.2 to 35 mm/yr)¹⁵. However, although Iran is mainly under a (semi)arid climate, some areas are wet with annual mean precipitation exceeding 1000 mm/yr (Fig. S1). As such, the nationwide recharge rate is higher than the global range reported for (semi) arid areas.

The mean amounts of estimated groundwater recharge lie within -44% to +21% of the recharge reported previously for some individual aquifers in Iran¹⁶⁻²⁰. Differences between the methods used in the previous studies and the estimation method developed for this study could explain inconsistencies in the results. The observed progressive decline in groundwater recharge could account for the reduction in the studied aquifers in 2022 compared to their previous levels. Moreover, previous studies have also expressed concerns regarding the diminished recharge in various aquifers across the country²¹⁻²⁶. Overall, the existing literature, alongside the records of Iran's water authorities (see Note S1), unanimously agree that groundwater resource overdrafts and declining recharge have not only pushed the country into a state of water bankruptcy but also posed a significant threat to the country's socioeconomic development.

## Groundwater recharge ratio

The countrywide groundwater recharge ratio (Fig. 2), defined as the fraction of recharge to precipitation, declined from 21% in 2006 to 14%

in 2017, averaging about 17% during our study period, closely approximating the average groundwater recharge ratio at the global scale (i.e., 16%)³. The combination of declined recharge and overexploitation of aquifers has resulted in approximately 422 out of 609 Iranian plains being in a critical/prohibited state. In this context, the term critical/prohibited signifies a condition where digging new wells is prohibited, except for drinking water use. Meanwhile, many of the free plains (i.e., where drilling new wells is legal) are located in desert or mountainous regions, where the aquifers have poor potential for groundwater usage. However, on a smaller geographical scale, i.e., in basins, the average recharge ratio during our study period varied from 3.3% (in Anzali basin) to 39.8% (in Bakhtegan basin) (Fig. 3). In general, the average recharge ratio is lower in the basins with the lowest groundwater depletion such as Anzali and Haraz basins⁸, which are located in the wet region in the north of Iran. On the contrary, the average recharge ratio is higher in the basins with the highest decrease in groundwater storage, such as Bakhtegan and Salt Lake basins²⁷, which are located in relatively dry zones in the central and eastern parts of Iran, where the groundwater resources dominantly contribute to the supply of water for human and environment (up to 90%). Interestingly, both Anzali and Haraz basins, which had the lowest recharge ratio, were the only basins that demonstrated a decreasing trend in groundwater electrical conductivity during the study period⁸. This decline in electrical conductivity could be attributed to the reduced input of contaminants through infiltration into the groundwater²⁸.

## Main drivers of declined groundwater recharge

While Noori et al.⁸ and Ashraf et al.²⁷ showed that the current extraction rates are unsustainable, other hydroclimatic and anthropogenic drivers (e.g., precipitation, evaporation, and changes in the LUC processes at the surface) governing the decline in aquifer recharge need to be further investigated. To better understand the impact of the natural and anthropogenic drivers on the decline in Iran's groundwater recharge, the changes in the annual mean precipitation across the country and each basin were investigated during the 2000 to 2018 period (Fig. S2). Our results suggest a non-significant change in annual mean precipitation in Iran (p value ≤ 0.05), in line with the results reported by Moshir Panahi et al.²⁹ for the period from 1986 to 2016. No significant change was also detected in the annual mean precipitation in 26 out of the 30 basins investigated in this study (p value ≤ 0.05). Similarly, no significant change was observed for the evapotranspiration reported in Iran, for the period from 1986 to 2016²⁹. Precipitation and (potential) evapotranspiration are the main natural drivers of groundwater recharge. Other natural-based controls on recharge include, among other factors, geology and topography³⁰; however, these factors remained relatively stable in Iran over the observation period. This suggests that human interventions have dominantly impacted the decline in Iran's groundwater recharge. It is noteworthy that insignificant changes in natural drivers do not necessarily mean that recharge generation mechanisms remain unchanged. For example, snow contributes more than rainfall to groundwater recharge³¹,³². The study conducted by Safarianzengir et al.³³ highlighted a strong correlation between the decline in groundwater and a decreasing trend in snowfall in Iran (2000–2019). Recent observations also reported a decline in snowfall in different parts of Iran during the last decades³⁴⁻³⁶.

However, our results reveal the average amount of nationwide groundwater recharge was 39.6 mm/yr during the study period, which is more than the reported average annual runoff in Iran, i.e., about 32 mm/yr, for the same period (2001–2016)²⁹, suggesting the national surface water resources are likely seeping into the groundwater basins. This finding is further supported by Saemian et al.³⁷ and Maghrebi et al.³⁸ studies, which reported that Iran's surface water resources, as the main contributor of groundwater recharge, have largely declined during the
role in this process¹⁰,¹¹. In fact, spatial and temporal variability of groundwater recharge, as an index of groundwater availability¹², is not accessible even locally in most parts of the world, including Iran. This scarcity of data indicates how rarely this variable is monitored. A possible way to obtain estimates of groundwater recharge over a large spatial scale is using global/regional hydrologic models. However, the recharge estimated by these models is not certain, especially in semi-arid regions such as Iran, where uncertainties in evapotranspiration and precipitation may lead to large errors in the results¹³,¹⁴, hindering the detection of regional imbalances between groundwater recharge and human withdrawals⁸.

Following the study conducted by Noori et al.⁸ on Iran's groundwater depletion, which raised concerns about the lack of the country's groundwater recharge data, we estimate its recharge using an extensive, nationwide database in all 30 hydrological basins of Iran from 2002 to 2017. In this study period, the country's groundwater resources were largely depleted due to anthropogenic and natural drivers. Our results suggest a significant decline in the nationwide groundwater recharge, primarily attributed to unsustainable development practices. By elucidating the anthropogenic and natural drivers contributing to this decline in Iran's groundwater recharge as a representative of a (semi)arid region in the world, we aim to foster a more comprehensive understanding of the regional and global impacts of these large-scale depletions. Our findings advocate for the urgent need for transformative changes in groundwater management and policies, especially in regions experiencing increasing hydroclimatic fluctuations.

# Results and discussion

## Estimated groundwater recharge

The estimated groundwater recharges for each aquifer, totaling 666 aquifers, aggregated in Iran's hydrological basins from 2002 to 2017, are depicted in Fig. 1. The Method's section details how these values were obtained. Our findings reveal a statistically significant decrease in nationwide groundwater recharge of about 35% (p value ≤ 0.001) during the study period, which has resulted in a groundwater imbalance in many basins (Fig. 1). The annual change in groundwater recharge varies from -10.3 mm/yr in the Great Karoon basin to +1.9 mm/yr in the Anzali basin, averaging about -3.8 mm/yr across the country (Table 1). The average groundwater recharge across the country and for the study period is estimated at 39.6 mm/yr, which is slightly above the upper end of the range globally reported for (semi) arid areas (i.e., 0.2 to 35 mm/yr)¹⁵. However, although Iran is mainly under a (semi)arid climate, some areas are wet with annual mean precipitation exceeding 1000 mm/yr (Fig. S1). As such, the nationwide recharge rate is higher than the global range reported for (semi) arid areas.

The mean amounts of estimated groundwater recharge lie within -44% to +21% of the recharge reported previously for some individual aquifers in Iran¹⁶⁻²⁰. Differences between the methods used in the previous studies and the estimation method developed for this study could explain inconsistencies in the results. The observed progressive decline in groundwater recharge could account for the reduction in the studied aquifers in 2022 compared to their previous levels. Moreover, previous studies have also expressed concerns regarding the diminished recharge in various aquifers across the country²¹⁻²⁶. Overall, the existing literature, alongside the records of Iran's water authorities (see Note S1), unanimously agree that groundwater resource overdrafts and declining recharge have not only pushed the country into a state of water bankruptcy but also posed a significant threat to the country's socioeconomic development.

## Groundwater recharge ratio

The countrywide groundwater recharge ratio (Fig. 2), defined as the fraction of recharge to precipitation, declined from 21% in 2006 to 14%

in 2017, averaging about 17% during our study period, closely approximating the average groundwater recharge ratio at the global scale (i.e., 16%)³. The combination of declined recharge and overexploitation of aquifers has resulted in approximately 422 out of 609 Iranian plains being in a critical/prohibited state. In this context, the term critical/prohibited signifies a condition where digging new wells is prohibited, except for drinking water use. Meanwhile, many of the free plains (i.e., where drilling new wells is legal) are located in desert or mountainous regions, where the aquifers have poor potential for groundwater usage. However, on a smaller geographical scale, i.e., in basins, the average recharge ratio during our study period varied from 3.3% (in Anzali basin) to 39.8% (in Bakhtegan basin) (Fig. 3). In general, the average recharge ratio is lower in the basins with the lowest groundwater depletion such as Anzali and Haraz basins⁸, which are located in the wet region in the north of Iran. On the contrary, the average recharge ratio is higher in the basins with the highest decrease in groundwater storage, such as Bakhtegan and Salt Lake basins²⁷, which are located in relatively dry zones in the central and eastern parts of Iran, where the groundwater resources dominantly contribute to the supply of water for human and environment (up to 90%). Interestingly, both Anzali and Haraz basins, which had the lowest recharge ratio, were the only basins that demonstrated a decreasing trend in groundwater electrical conductivity during the study period⁸. This decline in electrical conductivity could be attributed to the reduced input of contaminants through infiltration into the groundwater²⁸.

## Main drivers of declined groundwater recharge

While Noori et al.⁸ and Ashraf et al.²⁷ showed that the current extraction rates are unsustainable, other hydroclimatic and anthropogenic drivers (e.g., precipitation, evaporation, and changes in the LUC processes at the surface) governing the decline in aquifer recharge need to be further investigated. To better understand the impact of the natural and anthropogenic drivers on the decline in Iran's groundwater recharge, the changes in the annual mean precipitation across the country and each basin were investigated during the 2000 to 2018 period (Fig. S2). Our results suggest a non-significant change in annual mean precipitation in Iran (p value ≤ 0.05), in line with the results reported by Moshir Panahi et al.²⁹ for the period from 1986 to 2016. No significant change was also detected in the annual mean precipitation in 26 out of the 30 basins investigated in this study (p value ≤ 0.05). Similarly, no significant change was observed for the evapotranspiration reported in Iran, for the period from 1986 to 2016²⁹. Precipitation and (potential) evapotranspiration are the main natural drivers of groundwater recharge. Other natural-based controls on recharge include, among other factors, geology and topography³⁰; however, these factors remained relatively stable in Iran over the observation period. This suggests that human interventions have dominantly impacted the decline in Iran's groundwater recharge. It is noteworthy that insignificant changes in natural drivers do not necessarily mean that recharge generation mechanisms remain unchanged. For example, snow contributes more than rainfall to groundwater recharge³¹,³². The study conducted by Safarianzengir et al.³³ highlighted a strong correlation between the decline in groundwater and a decreasing trend in snowfall in Iran (2000–2019). Recent observations also reported a decline in snowfall in different parts of Iran during the last decades³⁴⁻³⁶.

However, our results reveal the average amount of nationwide groundwater recharge was 39.6 mm/yr during the study period, which is more than the reported average annual runoff in Iran, i.e., about 32 mm/yr, for the same period (2001–2016)²⁹, suggesting the national surface water resources are likely seeping into the groundwater basins. This finding is further supported by Saemian et al.³⁷ and Maghrebi et al.³⁸ studies, which reported that Iran's surface water resources, as the main contributor of groundwater recharge, have largely declined during the
<figure><img src="image_4.png" /></figure>

**Fig. 1 |** Temporal trend of annual mean groundwater recharge (mm/yr) in Iran (dark color line with circle marker) and in each of the 30 basins from 2002 to 2017. This figure shows a statistically significant decrease in nationwide ground-water recharge of about 35% (p value ≤ 0.001) during the study period. The annual change in groundwater recharge varies from -10.3 mm/yr in the Great Karoon basin to +1.9 mm/yr in the Anzali basin, averaging about -3.8 mm/yr across the country. Hd (light brown color dash line with square marker): Hirmand basin, SB (light gray color dash line with square marker): South Baluchestan basin, LD (blue color dash line with square marker): Lut Desert basin, HM (red color dash line with square marker): Hamun Mashkil basin, Pn (green color dash line with square marker): Patargan basin, Ak (light gray color dash line with diamond marker): Atrak basin, SkD (light brown color dash line with diamond marker): Siahkooh Desert basin, CD (blue color dash line with square marker): Central Desert basin, Sn (red color dash line with diamond marker): Sirjan basin, As (green color dash line with diamond marker): Aras basin, SgD (blue color dash line with plus marker): Saghand Desert

basin, Sj (light gray color dash line with plus marker): Sedij basin, Mn (light brown color dash line with plus marker): Mehran basin, Hh (red color dash line with plus marker): Helleh basin, Hn (green color dash line with plus marker): Hamun basin, Sd (light gray color dash line with cross marker): Sefid-roud basin, LU (blue color dash line with cross marker): Lake Urmia basin, Ji (light brown color dash line with cross marker): Jarrahi basin, Qm (red color dash line with cross marker): Qareghom basin, Kn (green color dash line with cross marker): Karian basin, WBR (light brown color dash line with triangle marker): West Boundary River basin, Ai (light gray color dash line with triangle marker): Anzali basin, Kh (blue color dash line with triangle marker): Karkheh basin, Hz (red color dash line with triangle marker): Haraz basin, Qu (green color dash line with triangle marker): Qaresou basin, Gi (light gray color dash line with circle marker): Gavkhouni basin, Gn (light brown color dash line with circle marker): Gorgan basin, SL (blue color dash line with circle marker): Salt Lake basin, Bn (red color dash line with circle marker): Bakhtegan basin, and GK (green color dash line with circle marker): Great Karoon basin.

last decades, mainly due to human interventions. Approximately 56% of Iranian rivers have undergone a decline in stream-flow, which is around 2.5 times higher than values reported for the world's large rivers. Around 20% of the country's permanent rivers have transformed into seasonal rivers, and former seasonal rivers have dried up and/or become narrow streams, leading to the development of new river and fluvial regimes³⁸. Anthropogenically dwindling renewable water has also

resulted in shrinkage of almost all of the national lakes, wetlands, marshes, and ponds⁸,³⁹. This, in turn, is likely to have resulted in a reduced aquifer recharge from surface water bodies. A decline in irrigation return flows that can be attributed to the technological progresses of agricultural practices and improvements in the efficiency of irrigation water-use during the last two decades, as suggested by Noori et al.⁸, could also contribute to the decline in the country's groundwater
<figure><img src="image_4.png" /></figure>

**Fig. 1 |** Temporal trend of annual mean groundwater recharge (mm/yr) in Iran (dark color line with circle marker) and in each of the 30 basins from 2002 to 2017. This figure shows a statistically significant decrease in nationwide ground-water recharge of about 35% (p value ≤ 0.001) during the study period. The annual change in groundwater recharge varies from -10.3 mm/yr in the Great Karoon basin to +1.9 mm/yr in the Anzali basin, averaging about -3.8 mm/yr across the country. Hd (light brown color dash line with square marker): Hirmand basin, SB (light gray color dash line with square marker): South Baluchestan basin, LD (blue color dash line with square marker): Lut Desert basin, HM (red color dash line with square marker): Hamun Mashkil basin, Pn (green color dash line with square marker): Patargan basin, Ak (light gray color dash line with diamond marker): Atrak basin, SkD (light brown color dash line with diamond marker): Siahkooh Desert basin, CD (blue color dash line with square marker): Central Desert basin, Sn (red color dash line with diamond marker): Sirjan basin, As (green color dash line with diamond marker): Aras basin, SgD (blue color dash line with plus marker): Saghand Desert

basin, Sj (light gray color dash line with plus marker): Sedij basin, Mn (light brown color dash line with plus marker): Mehran basin, Hh (red color dash line with plus marker): Helleh basin, Hn (green color dash line with plus marker): Hamun basin, Sd (light gray color dash line with cross marker): Sefid-roud basin, LU (blue color dash line with cross marker): Lake Urmia basin, Ji (light brown color dash line with cross marker): Jarrahi basin, Qm (red color dash line with cross marker): Qareghom basin, Kn (green color dash line with cross marker): Karian basin, WBR (light brown color dash line with triangle marker): West Boundary River basin, Ai (light gray color dash line with triangle marker): Anzali basin, Kh (blue color dash line with triangle marker): Karkheh basin, Hz (red color dash line with triangle marker): Haraz basin, Qu (green color dash line with triangle marker): Qaresou basin, Gi (light gray color dash line with circle marker): Gavkhouni basin, Gn (light brown color dash line with circle marker): Gorgan basin, SL (blue color dash line with circle marker): Salt Lake basin, Bn (red color dash line with circle marker): Bakhtegan basin, and GK (green color dash line with circle marker): Great Karoon basin.

last decades, mainly due to human interventions. Approximately 56% of Iranian rivers have undergone a decline in stream-flow, which is around 2.5 times higher than values reported for the world's large rivers. Around 20% of the country's permanent rivers have transformed into seasonal rivers, and former seasonal rivers have dried up and/or become narrow streams, leading to the development of new river and fluvial regimes³⁸. Anthropogenically dwindling renewable water has also

resulted in shrinkage of almost all of the national lakes, wetlands, marshes, and ponds⁸,³⁹. This, in turn, is likely to have resulted in a reduced aquifer recharge from surface water bodies. A decline in irrigation return flows that can be attributed to the technological progresses of agricultural practices and improvements in the efficiency of irrigation water-use during the last two decades, as suggested by Noori et al.⁸, could also contribute to the decline in the country's groundwater
<table><caption>**Table 1 | Trend analysis results and of recharge and the lower and upper limits of the 95% confidence interval in different Iranian basins and across the country**</caption><thead><tr><th>Basin</th><th>Sen's slope (L)</th><th>Significant level (p value)</th><th>L<sub>min95</sub></th><th>L<sub>max95</sub></th></tr></thead><tbody><tr><td>Sirjan</td><td>-0.485</td><td>≤0.01</td><td>-0.682</td><td>-0.185</td></tr><tr><td>Atrak</td><td>0.027</td><td>>0.05</td><td>-0.075</td><td>0.164</td></tr><tr><td>Aras</td><td>-0.334</td><td>≤0.01</td><td>-0.531</td><td>-0.131</td></tr><tr><td>South Baluchestan</td><td>0.188</td><td>≤0.01</td><td>0.004</td><td>0.264</td></tr><tr><td>Sedij</td><td>0.558</td><td>≤0.05</td><td>0.102</td><td>0.895</td></tr><tr><td>Jarrahi</td><td>0.295</td><td>>0.05</td><td>-0.272</td><td>0.873</td></tr><tr><td>Helleh</td><td>-0.496</td><td>≤0.05</td><td>-0.804</td><td>-0.036</td></tr><tr><td>Lake Urmia</td><td>-0.183</td><td>>0.05</td><td>-0.459</td><td>0.057</td></tr><tr><td>Salt Lake</td><td>-0.804</td><td>≤0.05</td><td>-1.708</td><td>-0.136</td></tr><tr><td>Haraz</td><td>-2.423</td><td>≤0.001</td><td>-2.880</td><td>-1.721</td></tr><tr><td>Anzali</td><td>1.931</td><td>≤0.001</td><td>1.058</td><td>2.433</td></tr><tr><td>Sefid-roud</td><td>-0.455</td><td>>0.05</td><td>-1.414</td><td>0.083</td></tr><tr><td>Bakhtegan</td><td>-2.907</td><td>≤0.01</td><td>-4.349</td><td>-1.466</td></tr><tr><td>Gorgan</td><td>-3.003</td><td>≤0.05</td><td>-6.193</td><td>-0.628</td></tr><tr><td>Qareghom</td><td>0.252</td><td>>0.05</td><td>-0.324</td><td>0.604</td></tr><tr><td>Great Karoon</td><td>-10.297</td><td>≤0.001</td><td>-13.520</td><td>-8.088</td></tr><tr><td>Karkheh</td><td>-2.135</td><td>≤0.01</td><td>-3.513</td><td>-1.039</td></tr><tr><td>Mehran</td><td>-0.162</td><td>>0.05</td><td>-0.333</td><td>0.209</td></tr><tr><td>Saghand Desert</td><td>-0.115</td><td>>0.05</td><td>-0.485</td><td>0.127</td></tr><tr><td>Siahkooh Desert</td><td>-0.510</td><td>≤0.001</td><td>-0.699</td><td>-0.343</td></tr><tr><td>Lut Desert</td><td>-0.051</td><td>≤0.01</td><td>-0.077</td><td>-0.018</td></tr><tr><td>Central Desert</td><td>-0.660</td><td>≤0.001</td><td>-0.894</td><td>-0.384</td></tr><tr><td>Gavkhouni</td><td>-3.955</td><td>≤0.001</td><td>-5.148</td><td>-2.768</td></tr><tr><td>West Boundary River</td><td>0.490</td><td>>0.05</td><td>-0.459</td><td>1.165</td></tr><tr><td>Karian</td><td>-0.99</td><td>≤0.05</td><td>-1.360</td><td>-0.224</td></tr><tr><td>Patargan</td><td>-0.01</td><td>>0.05</td><td>-0.082</td><td>0.223</td></tr><tr><td>Hamun</td><td>0.23</td><td>>0.05</td><td>-0.111</td><td>0.745</td></tr><tr><td>Hamun Mashkil</td><td>0.10</td><td>>0.05</td><td>-0.074</td><td>0.266</td></tr><tr><td>Hirmand</td><td>0.03</td><td>>0.05</td><td>-0.004</td><td>0.054</td></tr><tr><td>Qaresou</td><td>-1.56</td><td>≤0.05</td><td>-2.957</td><td>-0.354</td></tr><tr><td>Iran</td><td>-3.83</td><td>≤0.001</td><td>-4.240</td><td>-2.682</td></tr></tbody></table>

$L_{min95}$ and $L_{max95}$ are the lower and upper limits of the 95% confidence interval of calculated Sen's slope (L), respectively.

recharge—a concept known as pendulum swing⁴⁰. Maghrebi et al.⁴¹ proposed that one of the key factors contributing to the failure of government plans for aquifers restoration through modernization of the irrigation systems was the decline in groundwater recharge, primarily caused by reductions in irrigation return flows. This reduction in recharge led to a significant increase in the nationwide groundwater deficit, escalating from 106 km³ to more than 140 km³ during a 12-year period. Intensive and progressive land subsidence across the country during the last two decades has also contributed to reduced soil permeability due to compaction, thereby impacting the capacity for groundwater recharge⁴²⁻⁴⁵. In support of an effect of land subsidence, Iran is ranked among the countries with highest rate of land subsidence globally²⁵. As a result, aquifers' capacities in Iran are declining due to countrywide land subsidence. In this scenario, the soil compaction reduces its permeability and hence the aquifers' ability to be recharged, leading to a decline in both the rate and the total volume of recharge. Another factor influencing groundwater recharge is the changes in LUC, as several regions in Iran have experienced extensive deforestation, desertification, and rapid urbanization⁴⁶⁻⁵¹.

## Implications of reduced groundwater recharge

Reduction in groundwater recharge combined with the nonrenewable groundwater discharge from the country's aquifers⁸,²⁷,⁵² further contributes to decline in groundwater storage and shrink in groundwater table. The current declining trend in groundwater recharge and the reducing trend in groundwater table reported by Noori et al.⁸ and Ashraf et al.²⁷ persist through this century under status of shrinking renewable water and the quo water and environmental resources management in Iran. It has resulted in a gradual exacerbation of Iran's water and food security and makes the country's landscapes prone to a wider and more severe range of already experienced disasters, such as desertification, dust storms, landslides, land subsidence, sinkholes, droughts, floods, and fires. Additionally, it leads to a decrease in soil fertility, biodiversity, and an increase in greenhouse gas emission. From a regional perspective, such a decline in groundwater recharge has devastating consequences for the country's exports of agricultural products to the neighboring countries, compromising regional food security, given that Iran is one of the major food producers in the Persian Gulf region.

Artificial recharge is undoubtedly a promising engineering solution to recover depleted groundwater resources globally. Iran is also making efforts to complete a national plan that aims to artificially recharge groundwater up to about 1 km³. However, it should be noted that this amount is considerably less than the deficit in the country's

<figure><img src="image_5.png" /><figcaption>**Fig. 2 | Temporal trend of groundwater recharge ratio (%), i.e., ratio of recharge to precipitation, in Iran from 2002 to 2017.** The countrywide groundwater recharge ratio declined from 21% in 2006 to 14% in 2017, averaging -17% during the</figcaption></figure>

study period, closely approximating the average groundwater recharge ratio at the global scale (i.e., 16%).
<table><caption>**Table 1 | Trend analysis results and of recharge and the lower and upper limits of the 95% confidence interval in different Iranian basins and across the country**</caption><thead><tr><th>Basin</th><th>Sen's slope (L)</th><th>Significant level (p value)</th><th>L<sub>min95</sub></th><th>L<sub>max95</sub></th></tr></thead><tbody><tr><td>Sirjan</td><td>-0.485</td><td>≤0.01</td><td>-0.682</td><td>-0.185</td></tr><tr><td>Atrak</td><td>0.027</td><td>>0.05</td><td>-0.075</td><td>0.164</td></tr><tr><td>Aras</td><td>-0.334</td><td>≤0.01</td><td>-0.531</td><td>-0.131</td></tr><tr><td>South Baluchestan</td><td>0.188</td><td>≤0.01</td><td>0.004</td><td>0.264</td></tr><tr><td>Sedij</td><td>0.558</td><td>≤0.05</td><td>0.102</td><td>0.895</td></tr><tr><td>Jarrahi</td><td>0.295</td><td>>0.05</td><td>-0.272</td><td>0.873</td></tr><tr><td>Helleh</td><td>-0.496</td><td>≤0.05</td><td>-0.804</td><td>-0.036</td></tr><tr><td>Lake Urmia</td><td>-0.183</td><td>>0.05</td><td>-0.459</td><td>0.057</td></tr><tr><td>Salt Lake</td><td>-0.804</td><td>≤0.05</td><td>-1.708</td><td>-0.136</td></tr><tr><td>Haraz</td><td>-2.423</td><td>≤0.001</td><td>-2.880</td><td>-1.721</td></tr><tr><td>Anzali</td><td>1.931</td><td>≤0.001</td><td>1.058</td><td>2.433</td></tr><tr><td>Sefid-roud</td><td>-0.455</td><td>>0.05</td><td>-1.414</td><td>0.083</td></tr><tr><td>Bakhtegan</td><td>-2.907</td><td>≤0.01</td><td>-4.349</td><td>-1.466</td></tr><tr><td>Gorgan</td><td>-3.003</td><td>≤0.05</td><td>-6.193</td><td>-0.628</td></tr><tr><td>Qareghom</td><td>0.252</td><td>>0.05</td><td>-0.324</td><td>0.604</td></tr><tr><td>Great Karoon</td><td>-10.297</td><td>≤0.001</td><td>-13.520</td><td>-8.088</td></tr><tr><td>Karkheh</td><td>-2.135</td><td>≤0.01</td><td>-3.513</td><td>-1.039</td></tr><tr><td>Mehran</td><td>-0.162</td><td>>0.05</td><td>-0.333</td><td>0.209</td></tr><tr><td>Saghand Desert</td><td>-0.115</td><td>>0.05</td><td>-0.485</td><td>0.127</td></tr><tr><td>Siahkooh Desert</td><td>-0.510</td><td>≤0.001</td><td>-0.699</td><td>-0.343</td></tr><tr><td>Lut Desert</td><td>-0.051</td><td>≤0.01</td><td>-0.077</td><td>-0.018</td></tr><tr><td>Central Desert</td><td>-0.660</td><td>≤0.001</td><td>-0.894</td><td>-0.384</td></tr><tr><td>Gavkhouni</td><td>-3.955</td><td>≤0.001</td><td>-5.148</td><td>-2.768</td></tr><tr><td>West Boundary River</td><td>0.490</td><td>>0.05</td><td>-0.459</td><td>1.165</td></tr><tr><td>Karian</td><td>-0.99</td><td>≤0.05</td><td>-1.360</td><td>-0.224</td></tr><tr><td>Patargan</td><td>-0.01</td><td>>0.05</td><td>-0.082</td><td>0.223</td></tr><tr><td>Hamun</td><td>0.23</td><td>>0.05</td><td>-0.111</td><td>0.745</td></tr><tr><td>Hamun Mashkil</td><td>0.10</td><td>>0.05</td><td>-0.074</td><td>0.266</td></tr><tr><td>Hirmand</td><td>0.03</td><td>>0.05</td><td>-0.004</td><td>0.054</td></tr><tr><td>Qaresou</td><td>-1.56</td><td>≤0.05</td><td>-2.957</td><td>-0.354</td></tr><tr><td>Iran</td><td>-3.83</td><td>≤0.001</td><td>-4.240</td><td>-2.682</td></tr></tbody></table>

$L_{min95}$ and $L_{max95}$ are the lower and upper limits of the 95% confidence interval of calculated Sen's slope (L), respectively.

recharge—a concept known as pendulum swing⁴⁰. Maghrebi et al.⁴¹ proposed that one of the key factors contributing to the failure of government plans for aquifers restoration through modernization of the irrigation systems was the decline in groundwater recharge, primarily caused by reductions in irrigation return flows. This reduction in recharge led to a significant increase in the nationwide groundwater deficit, escalating from 106 km³ to more than 140 km³ during a 12-year period. Intensive and progressive land subsidence across the country during the last two decades has also contributed to reduced soil permeability due to compaction, thereby impacting the capacity for groundwater recharge⁴²⁻⁴⁵. In support of an effect of land subsidence, Iran is ranked among the countries with highest rate of land subsidence globally²⁵. As a result, aquifers' capacities in Iran are declining due to countrywide land subsidence. In this scenario, the soil compaction reduces its permeability and hence the aquifers' ability to be recharged, leading to a decline in both the rate and the total volume of recharge. Another factor influencing groundwater recharge is the changes in LUC, as several regions in Iran have experienced extensive deforestation, desertification, and rapid urbanization⁴⁶⁻⁵¹.

## Implications of reduced groundwater recharge

Reduction in groundwater recharge combined with the nonrenewable groundwater discharge from the country's aquifers⁸,²⁷,⁵² further contributes to decline in groundwater storage and shrink in groundwater table. The current declining trend in groundwater recharge and the reducing trend in groundwater table reported by Noori et al.⁸ and Ashraf et al.²⁷ persist through this century under status of shrinking renewable water and the quo water and environmental resources management in Iran. It has resulted in a gradual exacerbation of Iran's water and food security and makes the country's landscapes prone to a wider and more severe range of already experienced disasters, such as desertification, dust storms, landslides, land subsidence, sinkholes, droughts, floods, and fires. Additionally, it leads to a decrease in soil fertility, biodiversity, and an increase in greenhouse gas emission. From a regional perspective, such a decline in groundwater recharge has devastating consequences for the country's exports of agricultural products to the neighboring countries, compromising regional food security, given that Iran is one of the major food producers in the Persian Gulf region.

Artificial recharge is undoubtedly a promising engineering solution to recover depleted groundwater resources globally. Iran is also making efforts to complete a national plan that aims to artificially recharge groundwater up to about 1 km³. However, it should be noted that this amount is considerably less than the deficit in the country's

<figure><img src="image_5.png" /><figcaption>**Fig. 2 | Temporal trend of groundwater recharge ratio (%), i.e., ratio of recharge to precipitation, in Iran from 2002 to 2017.** The countrywide groundwater recharge ratio declined from 21% in 2006 to 14% in 2017, averaging -17% during the</figcaption></figure>

study period, closely approximating the average groundwater recharge ratio at the global scale (i.e., 16%).
<figure><img src="image_6.png" /><figcaption>**Fig. 3** | Mean groundwater recharge ratio (%) in each Iranian basin from 2002 to 2017. The average recharge ratio during the study period varied from 3.3% (in Anzali basin) to 39.8% (in Bakhtegan basin).</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 4** | Evolution of groundwater monitoring network in Iran. The blue line shows the onset of our study when the monitoring network covers around 78% of the national aquifers' area and >50% of the number of aquifers in Iran.</figcaption></figure>

groundwater resources (see Note S1). Iran's water resource management experiences highlight that engineering solutions alone cannot prevent the alarming depletion of groundwater resources for the country⁴¹, particularly when other factors such as countrywide land subsidence can hinder the effectiveness of recharge efforts. Our primary suggestion for facilitating the restoration of the depleted aquifers is prioritizing and reducing the water demand instead of relying on nonrenewable groundwater use. This recommendation requires adequate water governance backed by efficient conservation policies that are not currently implemented in Iran. The poor governance of hydro-environmental resources, characterized by unsustainable LUC planning, inequitable allocation rules, institutional ineffectiveness,

and poor economic conditions, makes it challenging to rehabilitate groundwater resources through reduction in abstraction and implementation of efficient artificial recharge strategies. To address this concern, the existing top-down governance structure approach should be replaced with a cooperative and bottom-up approach that considers the interests of local stakeholders, particularly farmers, as the main consumers of groundwater resources.

## Methods

### Study area and data

Iran's groundwater monitoring network has been gradually established since the early 1970s (Fig. 4 and Note S2). The network covers 492 out
<figure><img src="image_6.png" /><figcaption>**Fig. 3** | Mean groundwater recharge ratio (%) in each Iranian basin from 2002 to 2017. The average recharge ratio during the study period varied from 3.3% (in Anzali basin) to 39.8% (in Bakhtegan basin).</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 4** | Evolution of groundwater monitoring network in Iran. The blue line shows the onset of our study when the monitoring network covers around 78% of the national aquifers' area and >50% of the number of aquifers in Iran.</figcaption></figure>

groundwater resources (see Note S1). Iran's water resource management experiences highlight that engineering solutions alone cannot prevent the alarming depletion of groundwater resources for the country⁴¹, particularly when other factors such as countrywide land subsidence can hinder the effectiveness of recharge efforts. Our primary suggestion for facilitating the restoration of the depleted aquifers is prioritizing and reducing the water demand instead of relying on nonrenewable groundwater use. This recommendation requires adequate water governance backed by efficient conservation policies that are not currently implemented in Iran. The poor governance of hydro-environmental resources, characterized by unsustainable LUC planning, inequitable allocation rules, institutional ineffectiveness,

and poor economic conditions, makes it challenging to rehabilitate groundwater resources through reduction in abstraction and implementation of efficient artificial recharge strategies. To address this concern, the existing top-down governance structure approach should be replaced with a cooperative and bottom-up approach that considers the interests of local stakeholders, particularly farmers, as the main consumers of groundwater resources.

## Methods

### Study area and data

Iran's groundwater monitoring network has been gradually established since the early 1970s (Fig. 4 and Note S2). The network covers 492 out
of 609 plains (Fig. S3), including 666 aquifers (Fig. S4), and contains 12,293 active piezometers (Fig. S4) that measure the groundwater level at monthly intervals. The total area of aquifers is around 272,000 km² and the monitoring network covered more than 96% of this area (i.e., 261,000 km²) by the end of 2017. The remaining area (i.e., -11,000 km²), includes 133 relatively small and less important aquifers distributed across the country³.

Although the groundwater levels were measured in active piezo-meters from the early 1960s, the comprehensive database of Iran's groundwater resources was established from the water year of 2001-2002, when -78% of the aquifer areas became covered by the groundwater monitoring network³³. However, this database does not include groundwater recharge. To address this knowledge gap, we utilized the national groundwater database to calculate the groundwater recharge across Iran in each of the 30 basins during 2002 to 2017 period. This time interval was selected as our study period due to the data availability and adequate coverage of the country's aquifers using the groundwater monitoring network (Note S2).

## Groundwater recharge calculation

Groundwater recharge includes any penetrated water into the saturated zone³⁴. We used the water budget method, i.e., Eq. (1)¹⁵, to estimate groundwater recharge in each aquifer across Iran.

$$ \text{Recharge} = Q_{off} + \Delta S \quad (1) $$

where, $Q_{off}$ equals egressed water from the aquifer (positive flux out of aquifer), and $\Delta S$ is the groundwater storage change (negative when storage decreases) (see, Fig. S5).

In this study, the representative hydrograph was determined for each aquifer. Then the annual change in groundwater storage was estimated by considering the aquifer area, specific yield, and the annual average change in the groundwater level. Given the large number of aquifers investigated in this study (i.e., 666), we aggregated the calculated recharge to obtain values for each hydrological basin and for the entire country. Details about the aggregated data are given in Note S3 and also supplemented in https://zenodo.org/record/8382150.

Egressed water from the aquifers (i.e., $Q_{off}$) includes all anthropogenic groundwater withdrawals and natural groundwater losses. Anthropogenic groundwater withdrawals are fully measured based on sampling campaigns and field inventories in over a million wells and qanats on approximately a five-year cycle during wet and dry periods in Iran. The same measurements as the 5-year field campaigns are done twice a year during wet and dry periods at a small but representative set of selected extraction points in different aquifers. Natural groundwater losses can take place through different routes such as natural springs, evapotranspiration, and outgoing lateral fluxes from the aquifers. Fortunately, the national groundwater database includes the annual groundwater discharge through natural springs, which is measured as frequently as the anthropogenic groundwater withdrawals and thus also based on sampling campaigns and field inventories. Evaporation through the saturated zone is usually negligible when the depth to the groundwater table is large (>3 m depth). We supposed that the groundwater loss through evaporation, including transpiration, is negligible because the groundwater depth was more than 3 m in almost all the aquifers⁸. Groundwater losses through aquifer discharge to surface waters are also negligible due to the severe decline in the countrywide groundwater table during the last five decades (>20 m)⁵⁵. This hypothesis gains further support when data show both a widespread decline in Iran's river flows³⁸ and drying out the national wetlands and inland lakes⁸,⁵⁶. The natural groundwater losses through outgoing lateral groundwater fluxes are not considered due to the lack of data, which introduces uncertainty to the recharge estimate for individual aquifers. However, for the

nationwide recharge estimate, these losses can be reasonably ignored by disregarding cross-border groundwater fluxes. It should be noted that incoming lateral fluxes from the adjacent aquifers are represented in the recharge itself. Another issue is that surface water basins are determined by topography, and they do not necessarily match with regional aquifers, for which the territory is identified by underground geology. As a result, depending on how aquifers and surface drainage basins are situated in a given region, there might be additional fluxes, both incoming and outgoing to a given aquifer, which represent how the aquifer interacts with adjacent drainage basins. Lack of data in this regard limited our ability to consider such detailed information in our study. However, only 6 out of 666 studied aquifers were shared in two adjacent drainage basins. Therefore, we expected that this issue would have a minimal impact on our results. In summary, the simplifications implemented in this study may introduce uncertainties in the results. To address this issue, we conducted an uncertainty analysis that accounts for the impacts of such simplifications in our study.

## Data analysis

Groundwater data processing, including aggregation, division, multiplication, normalization, and statistical analysis, for 666 national aquifers during the study period (2002-2017), was performed using Pivot Table in Microsoft Excel environment, ArcGIS, and MATLAB. All figures were drawn in ArcGIS and MATLAB environments, Microsoft Office software, and Grapher program.

For a description of the change in groundwater recharge and precipitation, the magnitude of temporal trends was calculated by Sen's slope estimator method⁵⁷. A Mann-Kendall test⁵⁸,⁵⁹ was also applied to differentiate between significant and insignificant trends in both groundwater recharge and precipitation. These statistical analyses were done by MAKESENS 1.0 software, introduced by the Finnish Meteorological Institute⁶⁰.

## Data uncertainty

The results obtained in this study heavily rely on the data quality of the national groundwater database. Factors including poor information on aquifer type (e.g., confined or unconfined aquifer types), specific yield, and hypotheses made can potentially introduce errors/uncertainties in the calculated groundwater recharge for Iran and in each of the 30 basins. Detailed information on the data and the corresponding sources of errors/uncertainties is given by Noori et al.⁸.

An uncertainty analysis was conducted by determining the lower and upper limits of the 95% confidence interval for magnitude of the trends in groundwater recharge calculated using Sen's slope estimator⁵³ (see the Methods section for details). Table 1 shows the results of the uncertainty analysis for the calculated groundwater recharge across 30 basins and the countrywide values. By considering a significant level of 0.05 (i.e., $\alpha = 0.05$), the lower and upper limits of the confidence interval show the acceptable deviation from the magnitude of the trends for groundwater recharge across different basins and for the entire country.

## Data availability

The raw data used in this study are publicly available via Data Archive of the Iran Water Resources Management Company (IWRMC) at https://stu.wrm.ir/register.asp. Alternatively, the processed data presented in this paper are supplemented in https://zenodo.org/record/8382150.

## References

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

2. Boretti, A. & Rosa, L. Reassessing the projections of the world water development report. npj Clean. Water 2, 15 (2019).
of 609 plains (Fig. S3), including 666 aquifers (Fig. S4), and contains 12,293 active piezometers (Fig. S4) that measure the groundwater level at monthly intervals. The total area of aquifers is around 272,000 km² and the monitoring network covered more than 96% of this area (i.e., 261,000 km²) by the end of 2017. The remaining area (i.e., -11,000 km²), includes 133 relatively small and less important aquifers distributed across the country³.

Although the groundwater levels were measured in active piezo-meters from the early 1960s, the comprehensive database of Iran's groundwater resources was established from the water year of 2001–2002, when -78% of the aquifer areas became covered by the groundwater monitoring network³³. However, this database does not include groundwater recharge. To address this knowledge gap, we utilized the national groundwater database to calculate the groundwater recharge across Iran in each of the 30 basins during 2002 to 2017 period. This time interval was selected as our study period due to the data availability and adequate coverage of the country's aquifers using the groundwater monitoring network (Note S2).

## Groundwater recharge calculation

Groundwater recharge includes any penetrated water into the saturated zone³⁴. We used the water budget method, i.e., Eq. (1)¹⁵, to estimate groundwater recharge in each aquifer across Iran.

$$ \text{Recharge} = Q_{off} + \Delta S \quad (1) $$

where, $Q_{off}$ equals egressed water from the aquifer (positive flux out of aquifer), and $\Delta S$ is the groundwater storage change (negative when storage decreases) (see, Fig. S5).

In this study, the representative hydrograph was determined for each aquifer. Then the annual change in groundwater storage was estimated by considering the aquifer area, specific yield, and the annual average change in the groundwater level. Given the large number of aquifers investigated in this study (i.e., 666), we aggregated the calculated recharge to obtain values for each hydrological basin and for the entire country. Details about the aggregated data are given in Note S3 and also supplemented in https://zenodo.org/record/8382150.

Egressed water from the aquifers (i.e., $Q_{off}$) includes all anthropogenic groundwater withdrawals and natural groundwater losses. Anthropogenic groundwater withdrawals are fully measured based on sampling campaigns and field inventories in over a million wells and qanats on approximately a five-year cycle during wet and dry periods in Iran. The same measurements as the 5-year field campaigns are done twice a year during wet and dry periods at a small but representative set of selected extraction points in different aquifers. Natural groundwater losses can take place through different routes such as natural springs, evapotranspiration, and outgoing lateral fluxes from the aquifers. Fortunately, the national groundwater database includes the annual groundwater discharge through natural springs, which is measured as frequently as the anthropogenic groundwater withdrawals and thus also based on sampling campaigns and field inventories. Evaporation through the saturated zone is usually negligible when the depth to the groundwater table is large (>3 m depth). We supposed that the groundwater loss through evaporation, including transpiration, is negligible because the groundwater depth was more than 3 m in almost all the aquifers⁸. Groundwater losses through aquifer discharge to surface waters are also negligible due to the severe decline in the countrywide groundwater table during the last five decades (>20 m)⁵⁵. This hypothesis gains further support when data show both a widespread decline in Iran's river flows³⁸ and drying out the national wetlands and inland lakes⁸,⁵⁶. The natural groundwater losses through outgoing lateral groundwater fluxes are not considered due to the lack of data, which introduces uncertainty to the recharge estimate for individual aquifers. However, for the

nationwide recharge estimate, these losses can be reasonably ignored by disregarding cross-border groundwater fluxes. It should be noted that incoming lateral fluxes from the adjacent aquifers are represented in the recharge itself. Another issue is that surface water basins are determined by topography, and they do not necessarily match with regional aquifers, for which the territory is identified by underground geology. As a result, depending on how aquifers and surface drainage basins are situated in a given region, there might be additional fluxes, both incoming and outgoing to a given aquifer, which represent how the aquifer interacts with adjacent drainage basins. Lack of data in this regard limited our ability to consider such detailed information in our study. However, only 6 out of 666 studied aquifers were shared in two adjacent drainage basins. Therefore, we expected that this issue would have a minimal impact on our results. In summary, the simplifications implemented in this study may introduce uncertainties in the results. To address this issue, we conducted an uncertainty analysis that accounts for the impacts of such simplifications in our study.

## Data analysis

Groundwater data processing, including aggregation, division, multiplication, normalization, and statistical analysis, for 666 national aquifers during the study period (2002–2017), was performed using Pivot Table in Microsoft Excel environment, ArcGIS, and MATLAB. All figures were drawn in ArcGIS and MATLAB environments, Microsoft Office software, and Grapher program.

For a description of the change in groundwater recharge and precipitation, the magnitude of temporal trends was calculated by Sen's slope estimator method⁵⁷. A Mann-Kendall test⁵⁸,⁵⁹ was also applied to differentiate between significant and insignificant trends in both groundwater recharge and precipitation. These statistical analyses were done by MAKESENS 1.0 software, introduced by the Finnish Meteorological Institute⁶⁰.

## Data uncertainty

The results obtained in this study heavily rely on the data quality of the national groundwater database. Factors including poor information on aquifer type (e.g., confined or unconfined aquifer types), specific yield, and hypotheses made can potentially introduce errors/uncertainties in the calculated groundwater recharge for Iran and in each of the 30 basins. Detailed information on the data and the corresponding sources of errors/uncertainties is given by Noori et al.⁸.

An uncertainty analysis was conducted by determining the lower and upper limits of the 95% confidence interval for magnitude of the trends in groundwater recharge calculated using Sen's slope estimator⁵³ (see the Methods section for details). Table 1 shows the results of the uncertainty analysis for the calculated groundwater recharge across 30 basins and the countrywide values. By considering a significant level of 0.05 (i.e., $\alpha = 0.05$), the lower and upper limits of the confidence interval show the acceptable deviation from the magnitude of the trends for groundwater recharge across different basins and for the entire country.

## Data availability

The raw data used in this study are publicly available via Data Archive of the Iran Water Resources Management Company (IWRMC) at https://stu.wrm.ir/register.asp. Alternatively, the processed data presented in this paper are supplemented in https://zenodo.org/record/8382150.

## References

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

2. Boretti, A. & Rosa, L. Reassessing the projections of the world water development report. npj Clean. Water 2, 15 (2019).
3. Wada, Y. et al. Global depletion of groundwater resources. *Geo-phys. Res. Lett.* **37**, L20402 (2010).

4. Gleeson, T., Wada, Y., Bierkens, M. F. & Van Beek, L. P. Water balance of global aquifers revealed by groundwater footprint. *Nature* **488**, 197–200 (2012).

5. Dalin, C., Wada, Y., Kastner, T. & Puma, M. J. Groundwater depletion embedded in international food trade. *Nature* **543**, 700–706 (2017).

6. Döll, P. et al. Global-scale assessment of groundwater depletion and related groundwater abstractions: combining hydrological modeling with information from well observations and GRACE satellites. *Water Resour. Res.* **50**, 5698–5720 (2014).

7. Maghrebi, M. et al. Iran's groundwater hydrochemistry. *Earth Space Sci.* **8**, e2021EA001793 (2021).

8. Noori, R. et al. Anthropogenic depletion of Iran's aquifers. *Proc. Natl Acad. Sci. USA* **118**, e2024221118 (2021).

9. Maghrebi, M. et al. Anthropogenic decline of ancient, sustainable water systems: qanats. *Groundwater* **61**, 139–146 (2023).

10. De Vries, J. J. & Simmers, I. Groundwater recharge: an overview of processes and challenges. *Hydrogeol. J.* **10**, 5–17 (2002).

11. Jasechko, S. et al. The pronounced seasonality of global ground-water recharge. *Water Resour. Res.* **50**, 8845–8867 (2014).

12. Herbert, C. & Döll, P. Global assessment of current and future groundwater stress with a focus on transboundary aquifers. *Water Resour. Res.* **55**, 4760–4784 (2019).

13. Wartenburger, R. et al. Evapotranspiration simulations in ISIMIP2a–evaluation of spatio-temporal characteristics with a comprehensive ensemble of independent datasets. *Environ. Res. Lett.* **13**, 75001 (2018).

14. Reinecke, R. et al. Uncertainty of simulated groundwater recharge at different global warming levels: a global-scale multi-model ensemble study. *Hydrol. Earth Syst. Sci.* **25**, 787–810 (2021).

15. Scanlon, B. R., Healy, R. W. & Cook, P. G. Choosing appropriate techniques for quantifying groundwater recharge. *Hydrogeol. J.* **10**, 18–39 (2002).

16. Ahmadi, T. et al. Estimation of groundwater recharge using various methods in Neishaboor Plain, Iran. *Groundwater Modeling and Management under Uncertainty* (pp. 9–15. CRC Press, Taylor & Francis Group, London, UK, 2012).

17. Arabi Javanmard, Z. & Jafari, H. Estimating recharge to the Alehtar aquifer in Lorestan province. *J. Water Soil* **33**, 37–49 (2019).

18. Pahlevani Majdabady, M., Rasoulzadeh, A., Kanooni, A. & Ahmadzadeh, G. Estimation of groundwater recharge originating from agricultural irrigation and rainfall in Shyramyn plain, Iran. *Irrig. Drain.* **69**, 107–120 (2020).

19. Babaei, M. & Ketabchi, H. Estimation of groundwater recharge rate using a distributed model (case study of Rafsanjan aquifer, Kerman province). *Iran. J. Soil Water Res.* **51**, 1457–1468 (2020).

20. Moosavi, V. & Hayatzadeh, M. Groundwater recharge modeling using semi-distributed SWAT Model, case study: Marvast Plain. *Watershed Eng. Manag.* **14**, 282–298 (2022).

21. Mohammadi, Z., Salimi, M. & Faghih, A. Assessment of groundwater recharge in a semi-arid groundwater system using water balance equation, southern Iran. *J. Afr. Earth Sci.* **95**, 1–8 (2014).

22. Zamani, T., Karimi, H., Tavakoli, M. & Alimoradi, S. Investigation the factors affecting the groundwater drawdown in Mehran Plain. *Ilam. Hydrogeol.* **2**, 17–28 (2018).

23. Aghazadeh, N., Chitsazan, M. & Mirzaee, Y. The potential and actual urban aquifer recharge and site selection for artificial recharge using GIS and AHP methods (case study: Urmia urban aquifer). *J. Adv. Appl. Geol.* **32**, 168–179 (2019).

24. Zareian, M. J. & Eslamian, S. Groundwater withdrawal adjustment based on changes in groundwater balance components (a case study: an arid region in central Iran). *Arab. J. Geosci.* **14**, 1822 (2021).

25. IPRC (Iranian Parliament Research Center). Challenges of unsustainable groundwater withdrawal in Iran: current situation and land subsidence crisis (accessed 08 June 2023) https://rc.majlis.ir/fa/report/show/1773166 (2023).

26. Maghsoudi, R., Javadi, S., Shourian, M. & Golmohammadi, G. Determining the optimal aquifer exploitation under artificial recharge using the combination of numerical models and particle swarm optimization. *Hydrology* **10**, 100 (2023).

27. Ashraf, S., Nazemi, A. & AghaKouchak, A. Anthropogenic drought dominates groundwater depletion in Iran. *Sci. Rep.* **11**, 9135 (2021).

28. Uhl, A. et al. Making waves: pulling the plug—climate change effects will turn gaining into losing streams with detrimental effects on groundwater quality. *Water Res.* **220**, 118649 (2022).

29. Moshir Panahi, D. et al. Variability and change in the hydro-climate and water resources of Iran over a recent 30-year period. *Sci. Rep.* **10**, 7450 (2020).

30. Dripps, W. R. & Bradbury, K. R. The spatial and temporal variability of groundwater recharge in a forested basin in northern Wisconsin. *Hydrol. Process.* **24**, 383–392 (2010).

31. Wu, W. Y. et al. Divergent effects of climate change on future groundwater availability in key mid-latitude aquifers. *Nat. Commun.* **11**, 3710 (2020).

32. Lone, S. A. et al. Meltwaters dominate groundwater recharge in cold arid desert of Upper Indus River Basin (UIRB), western Himalayas. *Sci. Total Environ.* **786**, 147514 (2021).

33. Safarianzengir, V. et al. Monitoring and analysis of changes in the depth and surface area snow of the Mountains in Iran using remote sensing data. *J. Indian Soc. Remote Sens.* **48**, 1479–1494 (2020).

34. Araghi, A. & Mousavi-Baygi, M. Variability in snowfall/total precipitation-day ratio in Iran. *Theor. Appl. Climatol.* **140**, 547–558 (2020).

35. Emami, F. & Koch, M. Modeling the impact of climate change on water availability in the Zarrine River Basin and inflow to the Boukan Dam. *Iran. Clim.* **7**, 51 (2019).

36. Habibi, M., Babaeian, I. & Schöner, W. Changing causes of drought in the Urmia Lake Basin–increasing influence of evaporation and disappearing snow cover. *Water* **13**, 3273 (2021).

37. Saemian, P. et al. How much water did Iran lose over the last two decades? *J. Hydrol. Reg. Stud.* **41**, 101095 (2022).

38. Maghrebi, M. et al. Spatiotemporal changes in Iranian rivers' discharge. *Elem. Sci. Anth* **11**, 00002 (2023).

39. Malekmohammadi, B., Uvo, C. B., Moghadam, N. T., Noori, R. & Abolfathi, S. Environmental risk assessment of wetland ecosystems using Bayesian belief networks. *Hydrology* **10**, 16 (2023).

40. Maghrebi, M., Noori, R. & AghaKouchak, A. Iran: Renovated irrigation network deepens water crisis. *Nature* **618**, 238 (2023).

41. Van Emmerik, T. H. M. et al. Socio-hydrologic modeling to understand and mediate the competition for water between agriculture development and environmental health: Murrumbidgee River basin, Australia. *Hydrol. Earth Syst. Sci.* **18**, 4239–4259 (2014).

42. Motagh, M. et al. Land subsidence in Iran caused by widespread water reservoir overexploitation. *Geophys. Res. Lett.* **35**, L16403 (2008).

43. Haghighi, M. H. & Motagh, M. Ground surface response to continuous compaction of aquifer system in Tehran, Iran: results from a long-term multi-sensor InSAR analysis. *Remote Sens. Environ.* **221**, 534–550 (2019).

44. Negahdary, M. Shrinking aquifers and land subsidence in Iran. *Science* **376**, 1279–1279 (2022).

45. Sadeghi, H. et al. Assessing the vulnerability of Iran to subsidence hazard using a hierarchical FUCOM-GIS framework. *Remote Sens. Appl. Soc. Environ.* **31**, 100989 (2023).

46. Madani, K., AghaKouchak, A. & Mirchi, A. Iran's socio-economic drought: challenges of a water-bankrupt nation. *Iran. Stud.* **49**, 997–1016 (2016).
3. Wada, Y. et al. Global depletion of groundwater resources. *Geo-phys. Res. Lett.* **37**, L20402 (2010).

4. Gleeson, T., Wada, Y., Bierkens, M. F. & Van Beek, L. P. Water balance of global aquifers revealed by groundwater footprint. *Nature* **488**, 197–200 (2012).

5. Dalin, C., Wada, Y., Kastner, T. & Puma, M. J. Groundwater depletion embedded in international food trade. *Nature* **543**, 700–706 (2017).

6. Döll, P. et al. Global-scale assessment of groundwater depletion and related groundwater abstractions: combining hydrological modeling with information from well observations and GRACE satellites. *Water Resour. Res.* **50**, 5698–5720 (2014).

7. Maghrebi, M. et al. Iran's groundwater hydrochemistry. *Earth Space Sci.* **8**, e2021EA001793 (2021).

8. Noori, R. et al. Anthropogenic depletion of Iran's aquifers. *Proc. Natl Acad. Sci. USA* **118**, e2024221118 (2021).

9. Maghrebi, M. et al. Anthropogenic decline of ancient, sustainable water systems: qanats. *Groundwater* **61**, 139–146 (2023).

10. De Vries, J. J. & Simmers, I. Groundwater recharge: an overview of processes and challenges. *Hydrogeol. J.* **10**, 5–17 (2002).

11. Jasechko, S. et al. The pronounced seasonality of global ground-water recharge. *Water Resour. Res.* **50**, 8845–8867 (2014).

12. Herbert, C. & Döll, P. Global assessment of current and future groundwater stress with a focus on transboundary aquifers. *Water Resour. Res.* **55**, 4760–4784 (2019).

13. Wartenburger, R. et al. Evapotranspiration simulations in ISIMIP2a–evaluation of spatio-temporal characteristics with a comprehensive ensemble of independent datasets. *Environ. Res. Lett.* **13**, 75001 (2018).

14. Reinecke, R. et al. Uncertainty of simulated groundwater recharge at different global warming levels: a global-scale multi-model ensemble study. *Hydrol. Earth Syst. Sci.* **25**, 787–810 (2021).

15. Scanlon, B. R., Healy, R. W. & Cook, P. G. Choosing appropriate techniques for quantifying groundwater recharge. *Hydrogeol. J.* **10**, 18–39 (2002).

16. Ahmadi, T. et al. Estimation of groundwater recharge using various methods in Neishaboor Plain, Iran. *Groundwater Modeling and Management under Uncertainty* (pp. 9–15. CRC Press, Taylor & Francis Group, London, UK, 2012).

17. Arabi Javanmard, Z. & Jafari, H. Estimating recharge to the Alehtar aquifer in Lorestan province. *J. Water Soil* **33**, 37–49 (2019).

18. Pahlevani Majdabady, M., Rasoulzadeh, A., Kanooni, A. & Ahmadzadeh, G. Estimation of groundwater recharge originating from agricultural irrigation and rainfall in Shyramyn plain, Iran. *Irrig. Drain.* **69**, 107–120 (2020).

19. Babaei, M. & Ketabchi, H. Estimation of groundwater recharge rate using a distributed model (case study of Rafsanjan aquifer, Kerman province). *Iran. J. Soil Water Res.* **51**, 1457–1468 (2020).

20. Moosavi, V. & Hayatzadeh, M. Groundwater recharge modeling using semi-distributed SWAT Model, case study: Marvast Plain. *Watershed Eng. Manag.* **14**, 282–298 (2022).

21. Mohammadi, Z., Salimi, M. & Faghih, A. Assessment of groundwater recharge in a semi-arid groundwater system using water balance equation, southern Iran. *J. Afr. Earth Sci.* **95**, 1–8 (2014).

22. Zamani, T., Karimi, H., Tavakoli, M. & Alimoradi, S. Investigation the factors affecting the groundwater drawdown in Mehran Plain. *Ilam. Hydrogeol.* **2**, 17–28 (2018).

23. Aghazadeh, N., Chitsazan, M. & Mirzaee, Y. The potential and actual urban aquifer recharge and site selection for artificial recharge using GIS and AHP methods (case study: Urmia urban aquifer). *J. Adv. Appl. Geol.* **32**, 168–179 (2019).

24. Zareian, M. J. & Eslamian, S. Groundwater withdrawal adjustment based on changes in groundwater balance components (a case study: an arid region in central Iran). *Arab. J. Geosci.* **14**, 1822 (2021).

25. IPRC (Iranian Parliament Research Center). Challenges of unsustainable groundwater withdrawal in Iran: current situation and land subsidence crisis (accessed 08 June 2023) https://rc.majlis.ir/fa/report/show/1773166 (2023).

26. Maghsoudi, R., Javadi, S., Shourian, M. & Golmohammadi, G. Determining the optimal aquifer exploitation under artificial recharge using the combination of numerical models and particle swarm optimization. *Hydrology* **10**, 100 (2023).

27. Ashraf, S., Nazemi, A. & AghaKouchak, A. Anthropogenic drought dominates groundwater depletion in Iran. *Sci. Rep.* **11**, 9135 (2021).

28. Uhl, A. et al. Making waves: pulling the plug—climate change effects will turn gaining into losing streams with detrimental effects on groundwater quality. *Water Res.* **220**, 118649 (2022).

29. Moshir Panahi, D. et al. Variability and change in the hydro-climate and water resources of Iran over a recent 30-year period. *Sci. Rep.* **10**, 7450 (2020).

30. Dripps, W. R. & Bradbury, K. R. The spatial and temporal variability of groundwater recharge in a forested basin in northern Wisconsin. *Hydrol. Process.* **24**, 383–392 (2010).

31. Wu, W. Y. et al. Divergent effects of climate change on future groundwater availability in key mid-latitude aquifers. *Nat. Commun.* **11**, 3710 (2020).

32. Lone, S. A. et al. Meltwaters dominate groundwater recharge in cold arid desert of Upper Indus River Basin (UIRB), western Himalayas. *Sci. Total Environ.* **786**, 147514 (2021).

33. Safarianzengir, V. et al. Monitoring and analysis of changes in the depth and surface area snow of the Mountains in Iran using remote sensing data. *J. Indian Soc. Remote Sens.* **48**, 1479–1494 (2020).

34. Araghi, A. & Mousavi-Baygi, M. Variability in snowfall/total precipitation-day ratio in Iran. *Theor. Appl. Climatol.* **140**, 547–558 (2020).

35. Emami, F. & Koch, M. Modeling the impact of climate change on water availability in the Zarrine River Basin and inflow to the Boukan Dam. *Iran. Clim.* **7**, 51 (2019).

36. Habibi, M., Babaeian, I. & Schöner, W. Changing causes of drought in the Urmia Lake Basin–increasing influence of evaporation and disappearing snow cover. *Water* **13**, 3273 (2021).

37. Saemian, P. et al. How much water did Iran lose over the last two decades? *J. Hydrol. Reg. Stud.* **41**, 101095 (2022).

38. Maghrebi, M. et al. Spatiotemporal changes in Iranian rivers' discharge. *Elem. Sci. Anth* **11**, 00002 (2023).

39. Malekmohammadi, B., Uvo, C. B., Moghadam, N. T., Noori, R. & Abolfathi, S. Environmental risk assessment of wetland ecosystems using Bayesian belief networks. *Hydrology* **10**, 16 (2023).

40. Maghrebi, M., Noori, R. & AghaKouchak, A. Iran: Renovated irrigation network deepens water crisis. *Nature* **618**, 238 (2023).

41. Van Emmerik, T. H. M. et al. Socio-hydrologic modeling to understand and mediate the competition for water between agriculture development and environmental health: Murrumbidgee River basin, Australia. *Hydrol. Earth Syst. Sci.* **18**, 4239–4259 (2014).

42. Motagh, M. et al. Land subsidence in Iran caused by widespread water reservoir overexploitation. *Geophys. Res. Lett.* **35**, L16403 (2008).

43. Haghighi, M. H. & Motagh, M. Ground surface response to continuous compaction of aquifer system in Tehran, Iran: results from a long-term multi-sensor InSAR analysis. *Remote Sens. Environ.* **221**, 534–550 (2019).

44. Negahdary, M. Shrinking aquifers and land subsidence in Iran. *Science* **376**, 1279–1279 (2022).

45. Sadeghi, H. et al. Assessing the vulnerability of Iran to subsidence hazard using a hierarchical FUCOM-GIS framework. *Remote Sens. Appl. Soc. Environ.* **31**, 100989 (2023).

46. Madani, K., AghaKouchak, A. & Mirchi, A. Iran's socio-economic drought: challenges of a water-bankrupt nation. *Iran. Stud.* **49**, 997–1016 (2016).
47. Shirvani, Z., Abdi, O., Buchroithner, M. F. & Pradhan, B. Analysing spatial and statistical dependencies of deforestation affected by residential growth: Gorganrood Basin, Northeast Iran. *Land Degrad. Dev.* **28**, 2176–2190 (2017).

48. Habibi, S., Basirat, M., & Razavi, M.H. Towards national urban policy in I.R. of Iran. In: Developing National Urban Policies, 231–254, https://doi.org/10.1007/978-981-15-3738-7_9 (Springer, Singapore, 2020).

49. Sadeghi, R., Abbasi-Shavazi, M.J., & Shahbazin, S. Internal migration in Iran. In: Internal Migration in the Countries of Asia, 295–317, https://doi.org/10.1007/978-3-030-44010-7_15 (Springer, Cham, 2020).

50. Eskandari Dameneh, H. et al. Desertification of Iran in the early twenty-first century: assessment using climate and vegetation indices. *Sci. Rep.* **11**, 20548 (2021).

51. Emadodin, I., Reinsch, T. & Taube, F. Drought and desertification in Iran. *Hydrology* **6**, 66 (2019).

52. Olen, S. M. Water bankruptcy in Iran. *Nat. Sustain.* **4**, 924–924 (2021).

53. IWRMC (Iran Water Resources Management Company). Water budget report of the Iran study areas (accessed 15 October 2022). Ministry of Energy, Tehran, Iran (2018).

54. Schicht, R.J. & Walton, W.C. Hydrologic budgets for three small watersheds in Illinois. *Illinois State Water Survey, Report of Investigation* **40**, https://www.isws.illinois.edu/pubdoc/RI/ISWSRI-40.pdf (1961).

55. Nabavi, E. Failed policies, falling aquifers: unpacking groundwater overabstraction in Iran. *Water Altern.* **11**, 699–724 (2018).

56. AghaKouchak, A. et al. Aral Sea syndrome desiccates Lake Urmia: call for action. *J. Great Lakes Res.* **41**, 307–311 (2015).

57. Sen, P. K. Estimates of the regression coefficient based on Kendall's tau. *J. Am. Stat. Assoc.* **63**, 1379–1389 (1968).

58. Mann, H. B. Nonparametric tests against trend. *Econometrica* **13**, 245–259 (1945).

59. Kendall, M. G. Rank correlation methods (Oxford University Press, 1975).

60. MAKESENS 1. Mann-Kendall test and Sen's slope estimates for trend of annual data. Finnish Meteorological Institute. https://en.ilmatieteenlaitos.fi/makesens (2002).

## Acknowledgements

We thank the Iran Water Resources Management Company for supporting us with groundwater data.

## Author contributions

R.N. and M.M. conceived the study. R.N. and M.M. carried out the analyses. R.N., M.M., S.Je., S.M.B., S.Ja. and M.N. contributed to the data. R.N., S.Je., S.M.B., E.H., S.P., S.A. and A.A. contributed discussions and modeling insights. R.N. wrote the article. All authors reviewed, edited, and approved the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-023-42411-2.

**Correspondence** and requests for materials should be addressed to Roohollah Noori.

**Peer review information** *Nature Communications* thanks Marc Bierkens and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

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

**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/.

© The Author(s) 2023
47. Shirvani, Z., Abdi, O., Buchroithner, M. F. & Pradhan, B. Analysing spatial and statistical dependencies of deforestation affected by residential growth: Gorganrood Basin, Northeast Iran. *Land Degrad. Dev.* **28**, 2176–2190 (2017).

48. Habibi, S., Basirat, M., & Razavi, M.H. Towards national urban policy in I.R. of Iran. In: Developing National Urban Policies, 231–254, https://doi.org/10.1007/978-981-15-3738-7_9 (Springer, Singapore, 2020).

49. Sadeghi, R., Abbasi-Shavazi, M.J., & Shahbazin, S. Internal migration in Iran. In: Internal Migration in the Countries of Asia, 295–317, https://doi.org/10.1007/978-3-030-44010-7_15 (Springer, Cham, 2020).

50. Eskandari Dameneh, H. et al. Desertification of Iran in the early twenty-first century: assessment using climate and vegetation indices. *Sci. Rep.* **11**, 20548 (2021).

51. Emadodin, I., Reinsch, T. & Taube, F. Drought and desertification in Iran. *Hydrology* **6**, 66 (2019).

52. Olen, S. M. Water bankruptcy in Iran. *Nat. Sustain.* **4**, 924–924 (2021).

53. IWRMC (Iran Water Resources Management Company). Water budget report of the Iran study areas (accessed 15 October 2022). Ministry of Energy, Tehran, Iran (2018).

54. Schicht, R.J. & Walton, W.C. Hydrologic budgets for three small watersheds in Illinois. *Illinois State Water Survey, Report of Investigation* **40**, https://www.isws.illinois.edu/pubdoc/RI/ISWSRI-40.pdf (1961).

55. Nabavi, E. Failed policies, falling aquifers: unpacking groundwater overabstraction in Iran. *Water Altern.* **11**, 699–724 (2018).

56. AghaKouchak, A. et al. Aral Sea syndrome desiccates Lake Urmia: call for action. *J. Great Lakes Res.* **41**, 307–311 (2015).

57. Sen, P. K. Estimates of the regression coefficient based on Kendall's tau. *J. Am. Stat. Assoc.* **63**, 1379–1389 (1968).

58. Mann, H. B. Nonparametric tests against trend. *Econometrica* **13**, 245–259 (1945).

59. Kendall, M. G. Rank correlation methods (Oxford University Press, 1975).

60. MAKESENS 1. Mann-Kendall test and Sen's slope estimates for trend of annual data. Finnish Meteorological Institute. https://en.ilmatieteenlaitos.fi/makesens (2002).

## Acknowledgements

We thank the Iran Water Resources Management Company for supporting us with groundwater data.

## Author contributions

R.N. and M.M. conceived the study. R.N. and M.M. carried out the analyses. R.N., M.M., S.Je., S.M.B., S.Ja. and M.N. contributed to the data. R.N., S.Je., S.M.B., E.H., S.P., S.A. and A.A. contributed discussions and modeling insights. R.N. wrote the article. All authors reviewed, edited, and approved the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-023-42411-2.

**Correspondence** and requests for materials should be addressed to Roohollah Noori.

**Peer review information** *Nature Communications* thanks Marc Bierkens and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

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

**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/.

© The Author(s) 2023