ABSTRACT Moisture recycling in precipitation is an important hydrological process, accounting for ~ 67% globally. North-east India, home to the world's wettest place, boasts vast

wetlands and forest-cover. Despite its proximity to the coast, we find locally recycled moisture to be the primary annual source of rainfall (~ 45%). During the pre-monsoon season, the

enriched δ18O (~ − 0.7 ‰) and high d-excess (~ 14 ‰) are ascribed to enhanced transpiration, owing to atmospheric instability which causes Nor’westers. During the Monsoon season, oceanic

flux provides increased surficial moisture, enabling deep-localised convection via evaporation. Significant localised recycling, even during the Monsoon season is estimated (~ 38%), with

forest-cover, our study highlights their indispensable role in governing regional hydro-meteorology and water availability. SIMILAR CONTENT BEING VIEWED BY OTHERS PROTO-MONSOON RAINFALL AND

GREENING IN CENTRAL ASIA DUE TO EXTREME EARLY EOCENE WARMTH Article 29 January 2024 HUMAN-INDUCED WARMING ACCELERATES LOCAL EVAPOTRANSPIRATION AND PRECIPITATION RECYCLING OVER THE TIBETAN

PLATEAU Article Open access 20 July 2024 GREENING OF INDIA AND REVIVAL OF THE SOUTH ASIAN SUMMER MONSOON IN A WARMER WORLD Article Open access 09 November 2024 INTRODUCTION The enormous

evaporation (4.60 × 105 bcm/year) from the oceans is the most important component of the global hydrological cycle. However, ~ 90% of the evaporated vapour precipitates back into the oceans,

leaving only a small fraction (4.9 × 104 bcm/year) of it reaching the continents. The continental precipitation received globally amounts to 1.20 × 105 bcm/year, thus necessitating 8.1 × 

104 bcm/year (~ 67%) from continental recycling of moisture, in order to balance the global water budget1. The distribution of global continental recycling is extremely inhomogeneous, with

places like Amazon basin (high) and Australia (low) lying on both extremes of the spectrum2. Similarly, India having very large variability in precipitation pattern is also expected to have

a large degree of spatial inhomogeneity in continental recycling. The Indian water budget dictates ~ 40% of continentally recycled moisture in precipitation3. North-East India, home to the

wettest place on the planet (Mawsynram), is a major biodiversity hotspot4. This region experiences dynamic weather governed by complex orography, with the Himalayas to the north and the

Shillong Plateau further south, coupled with presence of the largest (by volume) Indian river (Brahmaputra)5, as well as vast wetlands and forest cover6,7, akin to the Amazon. Thus, we

expect a maximum (>> 40%) role of localized recycling in this region. The existing literature also identifies this region with maximum localized recycling in India. However, the

maximum estimated contribution of recycled moisture from the existing studies is limited to ~ 25%2,8,9. The average annual rainfall of NE India is ~ 2000 mm, and ~ 25% of it rains during the

pre-monsoon (March–May) period7,10,11. The pre-monsoon heavy rainfall is associated with Nor’westers (atmospheric instabilities caused when warm and moist southwest winds are overrun by the

cold and dry westerlies12). During the ISM (Indian Summer Monsoon, June–September) season, NE India receives ~ 1500 mm of rainfall. This hefty rainfall together with glacial melting during

summers causes severe floods in Brahmaputra River5,13, and observed through the GRACE (Gravity Recovery and Climate Experiment) TWSA (Terrestrial water storage anomaly) seasonal

fluctuations, indicating increased water availability14 (Fig. 1). These floods refill the water in low lying flood plains and widespread wetlands of the region. The stable isotopes of

oxygen(δ18O) and hydrogen(δ2H) in water are efficient tracers of hydrological processes, as their isotopic values vary predictably during physical processes such as evaporation,

condensation, sublimation and deposition15. Several isotopic studies in the Amazon have already identified the significant role of the abundant wetlands and forest cover in

recycling16,17,18. Declining forest cover and wetlands in the Amazon have been reported to affect the precipitation pattern remarkably19,20. Like the Amazon, several studies have also

documented worrying trends of declining forest cover and wetlands along with studies suggesting the change in regional precipitation patterns for North-east India21,22. Previous studies have

attempted to explain the variability of the ISM rainfall with the help of speleothem-based oxygen isotope reconstruction in the region coupled with use of ISOGSM223,24. Such approaches

however are hindered by coarse spatial resolution of ISOGSM2 model, sparse availability of speleothem records as well as complex orography and moisture dynamics. A continuous time-series

measurement of stable water isotopes in rain as well as surface/river water can help shed light on the moisture transport pathways as well as identify major sources of moisture in the

region. However, the isotopic studies in NE India are very limited8,25,26. With this backdrop, the current study aims to decipher the role of Shillong Plateau and Brahmaputra in influencing

the regional hydro-meteorology, and to estimate seasonal contribution of continentally recycled moisture in regional precipitation using stable isotopes (δ18O, δ2H), satellite-observations,

reanalysis data as well as HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory) model simulations. RESULTS AND DISCUSSION PRE-MONSOON During the pre-monsoon period, NE India

receives a significant amount of rainfall (500 mm) owing to Nor’westers. Nor’westers carry moisture from the BOB further inland, therefore its effect on the regional hydro-meteorology is

more pronounced south of the Shillong Plateau as observed from relatively high values of CAPE (Convective Available Potential Energy), Q (specific humidity at surface), ET

(Evapotranspiration) flux (Fig. 2a,c,e,g) and is accompanied by frequent thunderstorms along with heavy rainfall10,27 when compared to the rest of India28,29. This contrasts with

observations directly north of the Shillong Plateau. The pre-monsoonal rainfall received over this region renders the atmosphere moist and conditionally unstable, hence promoting

convection26,30. During the Pre-monsoon season, we observe highly enriched values of δ18O (~ − 0.7‰) and high d-excess (~ 14‰) in daily precipitation samples collected from Jorhat, located

to the north of the Shillong Plateau (Fig. 3a). Chakraborty et al.8 have tried to attribute higher values of δ18O in precipitation to the increased role of localized recycling dominated by

transpiration. They have further explained this by suggesting uptake of progressively evaporated (isotopically enriched) soil moisture by the plants, so that the transpired vapour is

isotopically enriched. While this argument can explain the enriched δ18O, it fails to explain the increased d-excess because evaporation decreases the d-excess of residual soil

moisture15,31,32,33. This high d-excess can only be explained if the source d-excess values are significantly high. We hypothesise that the flooding in Brahmaputra [major source of high

d-excess (~ 12.8‰)] during the previous years’ ISM period, recharges the surface as well as subsurface water, which in turn is transpired by the plants5,13. Our hypothesis is further

substantiated by similar d-excess values in precipitation and river water (Fig. 2b) and higher localized recycling (~ 53%), estimated with the help of a simple empirical model based on

HYSPLIT backward wind trajectories (Fig. 4a). In the HYSPLIT simulations, we have identified four major sources for the origin of vapour- Bay of Bengal, Arabian Sea, Local recycling, and

continental which refers to the rest of recycled moisture but not of local origin (refer Section "Hysplit clustered trajectories" for more details). The increasing trend in monthly

estimated localized recycling (inset; bottom right Fig. 4 a) concurs with the increasing trend of evapo-transpiration given by Chakraborty et al.8. The Shillong Plateau serves to demarcate

this region into two distinct halves, with the enriched isotopic values in precipitation (δ18O ~ − 0.7‰) observed to the north, which is in contrast to the relatively depleted values

reported by Breitenbach et al.25, Laskar et al.34, and Tanoue et al.35 (δ18O ~ − 2.4 ‰), to the south of the Plateau. The data reported in IAEA/WMO global network of isotopes in

precipitation (GNIP) for the stations at Dhaka (2009–2018) and Sylhet (2009–2016), which are in Bangladesh to the south of the Shillong Plateau, also show similar relatively depleted trends

during the pre-monsoon season. Although, based on Rayleigh distillation, one would have expected a successive decreasing trend in both δ18O and d-excess with longer rainout history15,32.

This suggests a more dominant role of isotopically depleted BOB moisture to the south of the plateau, whereas its role on the north is diminished due to the orographic barrier of Shillong

Plateau. In the south of the plateau the enhanced values of CAPE, ET and Q indicate a greater role of localized recycling compared to the north (Fig. 2a,c,e,g). However, the huge influx of

BOB moisture dilutes the isotopic footprint of recycled moisture in the south of plateau. On the other hand, these meteorological parameters indicate lesser role of localized recycling to

the north, but the unavailability of significant BOB moisture gives room to the recycled moisture to express its isotopic signature in the precipitation. The highest values of δ18O (~ 0.8‰)

are observed during March along with lower d-excess (~ 11.7‰) compared to April–May (δ18O ~ − 1.2‰, d-excess ~ 14.7). This observation indicates significant evaporation from the falling

raindrop due to high temperature and low rainfall36. The same can also be inferred from the lower value of LMWL (Local Meteoric Water Line) slope (6.96) for March, when compared to the GMWL

(Global Meteoric Water Line)37. Towards the end of the pre-monsoon season, a systematic slight depletion in both δ18O and dexcess is observed. This can be ascribed to the small influx of

isotopically depleted BOB moisture25,38circumventing the Shillong plateau, as seen from the wind patterns as well (Fig. 4d). INDIAN SUMMER MONSOON North-east India receives maximum (~ 1500 

mm) rainfall during the ISM season (Jun-Sep). ISM is set about with the northward migration of the ITCZ (Intertropical Convergence Zone), thus driving a huge influx of oceanic [BOB and AS

(Arabian Sea)] moisture inland. Unlike Pre-monsoon, the regional hydro-meteorology to the north and south of the Shillong Plateau are analogous (Fig. 2b,d,f,h). The huge influx of oceanic

moisture drives a sharp increasing trend of available surface water both to the north and south of the Plateau. This in turn results in lowering the BLH (Boundary Layer Height, Supplementary

Fig. 1) and CIN (Convective Inhibition) and increase in ET, CAPE, and Q values (Fig. 2b,d,f,h). Breitenbach et al.25 and Laskar et al.34 have reported an overall decreasing trend in δ18O at

Mawlong and Hailakandi respectively, located immediately south of the Shillong Plateau. A similar trend is reported by Tanoue et al.35 at three locations in Bangladesh (Dhaka, Sylhet,

Chittagong), which is further corroborated with GNIP data from Dhaka and Sylhet. We also observe a similar depleting trend at Jorhat in the first half of monsoon (June-July). Within the

depleting trend, a baseline shift in δ18O can be observed in mid of June coinciding with the onset of the ISM, which is marked by a depletion in δ18O by ~ 4–5 ‰ (Fig. 3a) as well as an

increase in the GRACE TWSA values. Following this baseline shift, there is a continuous trend of depletion in both δ18O and d-excess. This can be accounted for by the huge influx of oceanic

moisture driven by northward migration of organized convection35, carrying isotopically depleted δ18O and lower d-excess. The oceanic contribution in precipitation estimated using HYSPLIT

for June-July is ~ 57%. Even with the ISM and huge oceanic influx, the locally recycled contribution is significant (~ 37%). This is because unlike Pre-monsoon, the moisture laden winds

brought by the oceanic influx approach the Shillong Plateau almost perpendicularly (Fig. 4b,e). This serves to limit the oceanic contribution to precipitation to the north of the Plateau and

results in forced orographic ascent which is evident from Fig. 4b (inset, top left), with marked ascent of air parcel near 25°N coinciding with the start of the Plateau. Furthermore, owing

to increased moisture availability and cloud cover, the transpiration rate is significantly reduced26,39, thus limiting the contribution of enriched transpired moisture to precipitation8,40.

As a result, we do not observe an enriching trend in δ18O associated with significant moisture recycling. Despite the homogeneity in hydro-meteorological parameters to the north and south

of the Plateau, we observe a disparity in our observed d-excess trend when compared to the decreasing trend reported by Breitenbach et al.25 and Laskar et al.34 in the latter half of the ISM

season (Aug–Sep). We would have expected a decreasing trend in d-excess accompanied by a depleting trend in δ18O at Jorhat on account of successive Rayleigh distillation process15,32.

However, we observe a depleting trend in δ18O without any corresponding trend in d-excess (Fig. 3a). A similar conspicuous isotopic pattern akin to Jorhat is also observed at Tezpur, also

located to the south of the Shillong Plateau38,41. One of the causal factors for the observed disparity could be the weakening of the ISM due to southward propagation of the ITCZ leading to

lowering of oceanic influx. This is augmented by the orographic barrier posed by the Shillong Plateau, further hindering the oceanic contribution to the north of the Plateau42. Another

causal mechanism is that in the latter half of the ISM, this region experiences frequent floods in Brahmaputra5,13, which is known to bring glacial fed water having depleted δ18O (~ − 10.3

‰) and high d-excess (~ 13.1‰)43 (Fig. 3b). The increased surface water availability leads to favourable conditions promoting local recycling of moisture via evaporation2. The admixture of

Brahmaputra moisture (high d-excess, low δ18O) with BOB influx (low d-excess, low δ18O) could explain the depleting trend in δ18O without any discernible accompanying trend in the d-excess.

The increased proportion of locally recycled moisture is also corroborated with the higher recycled contribution (~ 39% local, ~ 11% continental), accompanied by a drop in oceanic

contribution (from ~ 57% to ~ 50%), estimated using HYSPLIT (inset; bottom right Fig. 4b). From the foregoing, even during the ISM season, the locally recycled moisture plays an important

role to the north of the Shillong Plateau. The overall estimated contribution from localized recycling in precipitation is ~ 38%. This result indicates a relatively higher proportion of

localized recycling compared to the previous studies in this region2,8,9,26. The high proportion of recycling reported in our study is strongly backed up by the underlying regional

hydro-meteorology. As with the increased surficial moisture during the ISM season, there is a significant lowering of viscosity in the boundary layer, hence lowering BLH44,45 (Supplementary

Fig. 1). This is further accompanied by a reduction in the CIN and an increasing trend in spatial extent of higher CAPE values coinciding with the regions having high Q at the surface, which

further gets accentuated post-floods. Also, from Fig. 5, we observe significant updraft beginning nearer to the surface and extending to higher altitudes, indicating deeper localized

convective cells drawing moisture from the surface. A combination of all the above-mentioned factors together results in an environment conducive to stronger local recycling. POST-MONSOON

North-east India receives very scanty rainfall during the Post-monsoon (Oct–Dec) due to the retreat of the ITCZ46, resulting in systematic decline of oceanic contribution along with

weakening of monsoonal winds (Fig. 4f). The backward wind trajectory map also suggests that the distance traversed by the winds progressively reduces from Oct–Dec, thus drawing moisture from

further inland. The estimated contribution of locally recycled moisture also shows an increasing trend from Oct (~ 58%)–Dec (~ 93%). We observe a progressively increasing trend in both δ18O

and d-excess in precipitation (Fig. 3a) with shortened rainout history (Fig. 4c). IMPLICATIONS OF HIGH MOISTURE RECYCLING Our study has brought to light the vital role of locally recycled

moisture to precipitation throughout the year (~ 45%), at least to the north of the Shillong Plateau. During the Pre-monsoon season, we highlight the increased role of transpired moisture to

precipitation, hence, signifying the importance of vegetation and forest cover. During the ISM season, influx of oceanic moisture and floods in Brahmaputra greatly increases the surficial

water availability in low lying floodplains and wetlands, promoting enhanced recycling via evaporation. Several studies have reported dwindling wetlands and forest cover due to

overexploitation, urbanisation, siltation and deforestation6,47,48. This puts an increased stress on ecosystems and biodiversity as well as loss of livelihood and economy. Owing to climate

change, a notable shift in precipitation pattern has been reported49,50. On top of this, our study highlights the vital role of forest cover and wetlands in the regional hydrology and

precipitation patterns. Further anthropogenic activities such as rapid urbanization, building of dams or barrages could greatly endanger the regional hydrology and water availability.

METHODS ISOTOPIC ANALYSIS We collected daily rainwater samples at Jorhat (26.72°N, 94.18°E) situated in the state of Assam (marked in Fig. 1a(ii)), for the period 2010–11. We also collected

samples from the Brahmaputra River during 2008–2015 with most samples collected from two sites Pandu (26.17°N, 91.67°E) and Dibrugarh (27.50° N, 94.84°E). The oxygen and hydrogen isotopic

analyses (δ18O and δD) were done by standard equilibration method in which water samples are equilibrated respectively with CO2 and H2. The equilibrated CO2 and H2 gases were analysed in

isotope ratio mass spectrometer Delta V Plus in continuous flow mode using Gas bench II at Physical Research Laboratory (PRL), Ahmedabad, to measure the 18O/16O and D/H ratios to compute the

δ18O and δD values51. The four secondary laboratory standards prepared in bulk and stored at PRL were analysed by Isotope Hydrology Laboratory at IAEA, Vienna on special request from IWIN

National Program to obtain its authentic δ18O and δD values. Based on repeated analyses of multiple aliquots of these four secondary laboratory standards, the reproducibility of measurement

was found to be better than 0.1 ‰ for δ18O and 1 ‰ for δD. The oxygen and hydrogen isotopic composition are expressed in terms of abundance ratios of heavy to light isotopes (R = 18O/16O or

D/H) and reported as δ in per mil (‰) notation [δ18O or δD = (Rsample/Rstd − 1) × 1000]. Rstd is the ratio in the Vienna Standard Mean Ocean Water (VSMOW). Another derived parameter used in

this study is d-excess (= δD − 8 × δ18O), defined by Dansgaard, 1964 to study kinetic effects associated with evaporation of water. The temporal variation in the δ18O, δD and d-excess rain,

in conjunction with meteorological parameters has been used to discern various processes and factors such as vapor source, rainout history, localized recycling and role of Shillong

Plateau36,52,53,54. REANALYSIS DATA Both long-term as well as daily hydro-meteorological parameters were derived from ERA 5 dataset

(https://www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era5). In order to monitor the behaviour of CAPE, CIN, ET, and other meteorological parameters over a broader scale, we have

prepared monthly climatology for the study area (23.72°N–28°N, and 89°E–96.18°E) from 1979 to 2018 with resolution of 0.25° × 0.25°. The long-term climatological records are used to validate

whether the understanding of possible causal rain-forming mechanisms obtained with the help of isotopic records is consistent with the underlying hydro-meteorology of the region. It is

necessary to compute the height at which maximum convergence and condensation in clouds occurs and its variation with season in order to accurately represent the rain forming pathways with

the help of HYSPLIT backward wind trajectories. In order to resolve this, we estimate the pressure and geopotential height above msl (mean sea level) for which the CRWC (Cloud Specific Rain

Water Content) peaks and assume it to best represent the approximate cloud condensation height. This exercise is repeated for those days during 2010–2011 when we observe rainfall at Jorhat.

Monthly mean pressure for CRWC maxima is computed from the daily estimates. HYSPLIT backward wind trajectories are then generated with their respective cloud condensation height for each

month. HYSPLIT BACKWARD WIND TRAJECTORIES In the present study wind, trajectory analyses were done using Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) Model from

https://www.ready.noaa.gov/HYSPLIT.php55. The backward trajectories were generated with Jorhat (26.72°N, 94.18°E) as the starting location for the HYSPLIT run. Ensembles of four 120 h back

trajectories starting at every 6-h interval (00:00, 06:00, 12:00, 18:00) were obtained for every daily rain event. Thus, for each daily rain event, four trajectories were monitored to

understand moisture source location. HYSPLIT MOISTURE SOURCE ESTIMATION For the daily rain samples, we obtain the pressure/geopotential height in a 1° × 1° grid centred around Jorhat for

which the CRWC is maximum. Monthly estimates for pressure/altitude at which CRWC maxima occurs is obtained by taking the mean over rainy days for the particular month. This estimate is used

to initialise the height at which the backward trajectories converge over Jorhat. Specific humidity is used as a proxy for moisture in air parcels, with decrease attributed to precipitation

and increase mapped to moisture pickup from the source region based on Sodemann et al.56. We identify four distinct sources of moisture (a) Bay of Bengal (b) Arabian Sea (c) Local (d)

Continental. In order to ascertain the component labelled as ‘Local’, we consider the region covering North-eastern part of Indian subcontinent including part of north-western Myanmar. The

component labelled as ‘Continental’ accounts for the rest of the recycled moisture that is not derived locally. We observe that the moisture derived from Bay of Bengal, Arabian Sea, and

Local components account for ~ 91% of the total contribution and are sufficient to identify the major moisture regions. The relative contributions from the four sources are calculated at

every hourly interval for each trajectory. The four trajectories considered for any given day have different percentage contributions from each source. Hence, their contributions are amount

weighted wrt the final specific humidity reported at Jorhat for the respective trajectories. In order to obtain monthly/annual values, we amount weight the daily estimates with the amount of

daily rainfall. The detailed methodology is described in Oza et al.57. HYSPLIT CLUSTERED TRAJECTORIES HYSPLIT backward wind trajectories provide an efficient tool in identifying major

moisture source pathways and source locations. Despite this, when studying over hundreds of backward wind trajectories together, representing them graphically, in a concise and coherent

fashion could be a challenge. Clustering provides an efficient way to obtain grouping of similar trajectories, hence providing a clear and concise graphical representation of different

moisture transport pathways. The conventional HYSPLIT trajectory clustering technique attempts to group trajectories based on similarity between their respective endpoints with the help of a

spatial variance method55,56,57,58. Hence, we end up neglecting the altitude-based information which would be otherwise invaluable in identifying convection as well as understanding the

role of orographic barriers. Furthermore, the spatial variance method only includes information about the trajectory endpoints, thus discounting any variations that might have occurred

further along the trajectory pathway. Hence, we have invoked an improved clustering algorithm incorporating latitude, longitude as well as altitude at every 6 hourly intervals along the

trajectories’ journey. This approach can capture the role played by orographic barriers as well as highlight localised convection which cannot be conveyed from just a 2-dimensional

representation. For each trajectory we consider the latitude, longitude, altitude above msl at every 6-h interval along the 120 h of its journey to reduce the volume of data. Furthermore,

Principal Component Analysis (PCA) is performed to further reduce the dimensionality of the dataset. The first eight principal components are considered since they preserve > 95% of the

explained variance. Kmeans algorithm is used along with ‘means++’ initialisation for the purpose of clustering59. Elbow method is used to determine the ideal number of clusters for different

seasons. Specific humidity is plotted with the help of a colour bar along with the trajectory to map the moisture pickup regions. GRAVITY RECOVERY AND CLIMATE EXPERIMENT (GRACE) GRACE level

3 gridded data was used to observe TWSA with 2008–2013 as the baseline wrt which the anomalies are computed since this covers the period of both our as well as some past studies14. The

optional scaling factors/ gridded gain factors have been multiplied to the original product to make it comparable to model data such as GLDAS. The TWSA data was used to infer the period

during which we observe a sudden increase in moisture availability in the region, indicative of onset of the ISM as well as advent of floods. DATA AVAILABILITY The authors declare that the

reanalysis data can be obtained from https://www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era5. The HYSPLIT backward wind trajectories can be generated from

https://www.ready.noaa.gov/HYSPLIT.php. The GRACE TWSA data can be found at https://grace.jpl.nasa.gov/data/get-data/. The GNIP data used in the study can be accessed at

http://www-naweb.iaea.org/napc/ih/IHS_resources_gnip.html. The stable water isotopic data as well as code used to run the models are available from the corresponding author upon reasonable

request. REFERENCES * van der Ent, R. J. & Tuinenburg, O. A. The residence time of water in the atmosphere revisited. _Hydrol. Earth Syst. Sci._ 21, 779–790 (2017). Article  ADS  Google

Scholar  * Risi, C., Noone, D., Frankenberg, C. & Worden, J. Role of continental recycling in intraseasonal variations of continental moisture as deduced from model simulations and water

vapor isotopic measurements. _Water Resour. Res._ 49, 4136–4156 (2013). Article  ADS  CAS  Google Scholar  * Gupta, S. K. & Deshpande, R. D. Water for India in 2050: First-order

assessment of available options. _Curr. Sci._ 86, 9 (2004). Google Scholar  * Kuttippurath, J. _et al._ Observed rainfall changes in the past century (1901–2019) over the wettest place on

Earth. _Environ. Res. Lett._ 16, 024018 (2021). Article  ADS  Google Scholar  * Rao, M. P. _et al._ Seven centuries of reconstructed Brahmaputra River discharge demonstrate underestimated

high discharge and flood hazard frequency. _Nat. Commun._ 11, 6017 (2020). Article  ADS  CAS  Google Scholar  * Bassi, N., Kumar, M. D., Sharma, A. & Pardha-Saradhi, P. Status of

wetlands in India: A review of extent, ecosystem benefits, threats and management strategies. _J. Hydrol. Regional Stud._ 2, 1–19 (2014). Article  Google Scholar  * Dikshit, K. R. &

Dikshit, J. K. Natural vegetation: Forests and grasslands of North-East India. in _North-East India: Land, People and Economy_ (eds. Dikshit, K. R. & Dikshit, J. K.) 213–255 (Springer,

2014). https://doi.org/10.1007/978-94-007-7055-3_9. * Chakraborty, S. _et al._ Linkage between precipitation isotopes and biosphere-atmosphere interaction observed in northeast India. _Npj

Clim. Atmos. Sci._ 5, 1–11 (2022). Article  Google Scholar  * Pathak, A., Ghosh, S. & Kumar, P. Precipitation recycling in the Indian subcontinent during summer monsoon. _J.

Hydrometeorol._ 15, 2050–2066 (2014). Article  ADS  Google Scholar  * Mahanta, R., Sarma, D. & Choudhury, A. Heavy rainfall occurrences in northeast India. _Int. J. Climatol._ 33,

1456–1469 (2013). Article  Google Scholar  * Murata, F. _et al._ dominant synoptic disturbance in the extreme rainfall at cherrapunji, northeast India, based on 104 years of rainfall data

(1902–2005). _J. Clim._ 30, 8237–8251 (2017). Article  ADS  Google Scholar  * Roy, S. C. & Chatterji, G. Origin of nor’westers. _Nature_ 124, 481–481 (1929). Article  ADS  Google Scholar

  * Dhar, O. N. & Nandargi, S. A study of floods in the Brahmaputra basin in India. _Int. J. Climatol._ 20, 771–781 (2000). Article  Google Scholar  * Reager, J. T., Thomas, B. F. &

Famiglietti, J. S. River basin flood potential inferred using GRACE gravity observations at several months lead time. _Nat. Geosci._ 7, 588–592 (2014). Article  ADS  CAS  Google Scholar  *

Gat, J. R. _Isotope Hydrology: A Study of the Water Cycle_ (World Scientific, 2010). Book  Google Scholar  * Salati, E., DallOlio, A., Matsui, E. & Gat, J. R. Recycling of water in the

Amazon basin: An isotopic study. _Water Resourc. Res._ 15, 1250–1258 (1979). Article  ADS  CAS  Google Scholar  * Victoria, R. L., Martinelli, L. A., Mortatti, J. & Richey, J. Mechanisms

of water recycling in the Amazon basin: Isotopic insights. _Ambio_ 20, 384–387 (1991). Google Scholar  * Wright, J. S. _et al._ Rainforest-initiated wet season onset over the southern

Amazon. _PNAS_ 114, 8481–8486 (2017). Article  ADS  CAS  Google Scholar  * Leite-Filho, A. T., Soares-Filho, B. S., Davis, J. L., Abrahão, G. M. & Börner, J. Deforestation reduces

rainfall and agricultural revenues in the Brazilian Amazon. _Nat. Commun._ 12, 2591 (2021). Article  ADS  CAS  Google Scholar  * Spracklen, D. V. & Garcia-Carreras, L. The impact of

Amazonian deforestation on Amazon basin rainfall. _Geophys. Res. Lett._ 42, 9546–9552 (2015). Article  ADS  Google Scholar  * Lele, N. & Joshi, P. K. Analyzing deforestation rates,

spatial forest cover changes and identifying critical areas of forest cover changes in North-East India during 1972–1999. _Environ. Monit. Assess._ 156, 159 (2008). Article  Google Scholar 

* Sudhakar Reddy, C. _et al._ Quantification and monitoring of deforestation in India over eight decades (1930–2013). _Biodivers. Conserv._ 25, 93–116 (2016). Article  Google Scholar  *

Kathayat, G. _et al._ Interannual oxygen isotope variability in Indian summer monsoon precipitation reflects changes in moisture sources. _Commun. Earth Environ._ 2, 1–10 (2021). Article 

Google Scholar  * Kathayat, G. _et al._ Protracted Indian monsoon droughts of the past millennium and their societal impacts. _Proc. Natl. Acad. Sci._ 119, e2207487119 (2022). Article  CAS 

Google Scholar  * Breitenbach, S. F. M. _et al._ Strong influence of water vapor source dynamics on stable isotopes in precipitation observed in Southern Meghalaya, NE India. _Earth Planet.

Sci. Lett._ 292, 212–220 (2010). Article  ADS  CAS  Google Scholar  * Pradhan, R., Singh, N. & Singh, R. P. Onset of summer monsoon in Northeast India is preceded by enhanced

transpiration. _Sci. Rep._ 9, 18646 (2019). Article  ADS  CAS  Google Scholar  * Das, S., Sarkar, A., Das, M. K., Rahman, Md. M. & Islam, Md. N. Composite characteristics of Nor’westers

based on observations and simulations. _Atmos. Res._ 158–159, 158–178 (2015). Article  Google Scholar  * Vinay Kumar, P. & Venkateswara Naidu, C. Is pre-monsoon rainfall activity over

India increasing in the recent era of global warming?. _Pure Appl. Geophys._ 177, 4423–4442 (2020). Article  ADS  Google Scholar  * Narayanan, P., Sarkar, S., Basistha, A. & Sachdeva, K.

Trend analysis and forecast of pre-monsoon rainfall over India. _Weather_ 71, 94–99 (2016). Article  ADS  Google Scholar  * Stevens, B. Atmospheric moist convection. _Annu. Rev. Earth

Planet. Sci._ 33, 605–643 (2005). Article  ADS  MathSciNet  CAS  MATH  Google Scholar  * Bershaw, J. Controls on deuterium excess across Asia. _Geosciences_ 8, 257 (2018). Article  ADS 

Google Scholar  * Clark, I. D. & Fritz, P. _Environmental Isotopes in Hydrogeology_ (CRC Press, 1997). https://doi.org/10.1201/9781482242911. Book  Google Scholar  * Cui, J., Tian, L.,

Biggs, T. W. & Wen, R. Deuterium-excess determination of evaporation to inflow ratios of an alpine lake: Implications for water balance and modeling. _Hydrol. Process._ 31, 1034–1046

(2017). Article  ADS  Google Scholar  * Laskar, A. H. _et al._ Stable isotopic characterization of Nor’westers of southern Assam, NE India. _J. Clim. Chang._ 1, 75–87 (2015). Article  Google

Scholar  * Tanoue, M. _et al._ Seasonal variation in isotopic composition and the origin of precipitation over Bangladesh. _Prog. Earth Planet Sci._ 5, 77 (2018). Article  ADS  Google

Scholar  * Oza, H., Ganguly, A., Padhya, V. & Deshpande, R. D. Hydrometeorological processes and evaporation from falling rain in Indian sub-continent: Insights from stable isotopes and

meteorological parameters. _J. Hydrol._ 591, 125601 (2020). Article  CAS  Google Scholar  * Chen, F. _et al._ Relationship between sub-cloud secondary evaporation and stable isotopes in

precipitation of Lanzhou and surrounding area. _Quatern. Int._ 380–381, 68–74 (2015). Article  Google Scholar  * Sinha, N. _et al._ Isotopic investigation of the moisture transport processes

over the Bay of Bengal. _J. Hydrol. X_ 2, 100021 (2019). Article  CAS  Google Scholar  * Rawson, H. M., Begg, J. E. & Woodward, R. G. The effect of atmospheric humidity on

photosynthesis, transpiration and water use efficiency of leaves of several plant species. _Planta_ 134, 5–10 (1977). Article  CAS  Google Scholar  * Chakraborty, S., Belekar, A. R., Datye,

A. & Sinha, N. Isotopic study of intraseasonal variations of plant transpiration: An alternative means to characterise the dry phases of monsoon. _Sci. Rep._ 8, 8647 (2018). Article  ADS

  CAS  Google Scholar  * Sinha, N., Chakraborty, S. & Mohan, P. M. Modern rain-isotope data from Indian island and the mainland on the daily scale for the summer monsoon season. _Data

Brief_ 23, 103793 (2019). Article  Google Scholar  * Grujic, D. _et al._ Formation of a Rain Shadow: O and H stable isotope records in authigenic clays from the Siwalik group in eastern

Bhutan. _Geochem. Geophys. Geosyst._ 19, 3430–3447 (2018). Article  CAS  Google Scholar  * Lambs, L., Balakrishna, K., Brunet, F. & Probst, J. L. Oxygen and hydrogen isotopic composition

of major Indian rivers: A first global assessment. _Hydrol. Process._ 19, 3345–3355 (2005). Article  ADS  CAS  Google Scholar  * Ek, M. B. & Holtslag, A. A. M. Influence of soil

moisture on boundary layer cloud development. _J. Hydrometeorol._ 5, 86–99 (2004). Article  ADS  Google Scholar  * Findell, K. L. & Eltahir, E. A. B. Atmospheric controls on soil

moisture-boundary layer interactions. Part I: Framework development. _J. Hydrometeorol._ 4, 552–569 (2003). Article  ADS  Google Scholar  * Syroka, J. & Toumi, R. On the withdrawal of

the Indian summer monsoon. _Q. J. R. Meteorol. Soc._ 130, 989–1008 (2004). Article  ADS  Google Scholar  * Bhatta, L. D. _et al._ Ecosystem service changes and livelihood impacts in the

maguri-motapung wetlands of Assam. _India. Land_ 5, 15 (2016). Article  Google Scholar  * Choudhury, B. A., Saha, S. K., Konwar, M., Sujith, K. & Deshamukhya, A. Rapid drying of

northeast India in the last three decades: Climate change or natural variability?. _J. Geophys. Res. Atmos._ 124, 227–237 (2019). Article  ADS  Google Scholar  * Das, D. Changing climate and

its impacts on Assam, Northeast India. _Bandung J. Glob. South_ 2, 26 (2016). Article  Google Scholar  * Deka, R. L., Mahanta, C., Pathak, H., Nath, K. K. & Das, S. Trends and

fluctuations of rainfall regime in the Brahmaputra and Barak basins of Assam, India. _Theor. Appl. Climatol._ 114, 61–71 (2013). Article  ADS  Google Scholar  * Maurya, A. S., Shah, M.,

Deshpande, R. D. & Gupta, S. K. Protocol for δ18O and δD analyses of water sample using Delta V plus IRMS in CF Mode with Gas Bench II for IWIN National Programme at PRL, Ahmedabad. in

_11th ISMAS Triennial Conference of Indian Society for Mass Spectrometry_ vol. 314, 314–317 (Indian Society for Mass Spectrometry Hyderabad, 2009). * Deshpande, R. D. & Gupta, S. K.

Oxygen and hydrogen isotopes in hydrological cycle: new data from IWIN national programme. _Proc. Indian Natl. Sci. Acad._ 78, 321–331 (2012). CAS  Google Scholar  * Deshpande, R. D. &

Gupta, S. K. National programme on isotope fingerprinting of waters of India (IWIN). _Glimpses of Geosciences Research in India, the Indian Report to IUGS, Indian National Science Academy

_(eds Singhvi, AK, Bhattacharya, A. & Guha, S.), 10–16 (2008). * Oza, H., Padhya, V., Ganguly, A. & Deshpande, R. D. Investigating hydrometeorology of the Western Himalayas: Insights

from stable isotopes of water and meteorological parameters. _Atmos. Res._ 268, 105997 (2022). Article  CAS  Google Scholar  * Stein, A. F. _et al._ NOAA’s HYSPLIT atmospheric transport and

dispersion modeling system. _Bull. Am. Meteor. Soc._ 96, 2059–2077 (2015). Article  ADS  Google Scholar  * Sodemann, H., Schwierz, C. & Wernli, H. Interannual variability of Greenland

winter precipitation sources: Lagrangian moisture diagnostic and North Atlantic Oscillation influence. _J. Geophys. Res. Atmos._ 113, D3 (2008). Google Scholar  * Oza, H. _et al._

Hydrometeorological processes in semi-arid western India: insights from long term isotope record of daily precipitation. _Clim. Dyn._ 54, 2745–2757 (2020). Article  Google Scholar  * Su, L.,

Yuan, Z., Fung, J. C. H. & Lau, A. K. H. A comparison of HYSPLIT backward trajectories generated from two GDAS datasets. _Sci. Total Environ._ 506–507, 527–537 (2015). Article  ADS 

Google Scholar  * Ahmed, M., Seraj, R. & Islam, S. M. S. The k-means algorithm: A comprehensive survey and performance evaluation. _Electronics_ 9, 1295 (2020). Article  Google Scholar 

Download references ACKNOWLEDGEMENTS The work reported here is carried out under the aegis of a National Programme on Isotope Fingerprinting of Waters of India (IWIN) which was initially

funded jointly by the Department of Science and Technology (DST), Government of India, vide Grant No. IR/S4/ESF-05/2004, and the Physical Research Laboratory (PRL), a Unit under the

Department of Space, Government of India. The IWIN National Programme is currently sustained exclusively by the PRL. AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * Geosciences Division,

Physical Research Laboratory, Navrangpura, Ahmedabad, 380009, India Akash Ganguly, Harsh Oza, Virendra Padhya, Amit Pandey, Swagatika Chakra & R. D. Deshpande * Indian Institute of

Technology Gandhinagar, Gandhinagar, 382355, India Akash Ganguly Authors * Akash Ganguly View author publications You can also search for this author inPubMed Google Scholar * Harsh Oza View

author publications You can also search for this author inPubMed Google Scholar * Virendra Padhya View author publications You can also search for this author inPubMed Google Scholar * Amit

Pandey View author publications You can also search for this author inPubMed Google Scholar * Swagatika Chakra View author publications You can also search for this author inPubMed Google

Scholar * R. D. Deshpande View author publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS R.D.D. conceptualized the study. V.P. was responsible for data

curation of stable water isotopes. R.D.D., A.G. and H.O. interpreted the isotopic data. A.G. devised the methodology and was responsible for analysis of reanalysis data, HYSPLIT model

simulations and statistical analysis of the data. A.G., H.O. and R.D.D. drafted the manuscript. A.G., H.O., A.P. and S.C. helped in preparation of the graphics. CORRESPONDING AUTHOR

Correspondence to Akash Ganguly. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing interests. ADDITIONAL INFORMATION PUBLISHER'S NOTE Springer Nature remains

neutral with regard to jurisdictional claims in published maps and institutional affiliations. SUPPLEMENTARY INFORMATION SUPPLEMENTARY FIGURE 1. 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 licence, and indicate if changes were made. The images or other third party

material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the

article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Ganguly, A., Oza, H.,

Padhya, V. _et al._ Extreme local recycling of moisture via wetlands and forests in North-East Indian subcontinent: a Mini-Amazon. _Sci Rep_ 13, 521 (2023).

https://doi.org/10.1038/s41598-023-27577-5 Download citation * Received: 29 August 2022 * Accepted: 04 January 2023 * Published: 10 January 2023 * DOI:

https://doi.org/10.1038/s41598-023-27577-5 SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not

currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative