Abrupt climate fluctuations in Tibet as imprints of multiple meltwater events during the early to mid-Holocene

Understanding the impact of meltwater discharge during the final stage of the Laurentide Ice Sheet (LIS) has important implications for predicting sea level rise and climate change. Here we present a high-resolution ice-core isotopic record from the central Tibetan Plateau (TP), where the climate is sensitive to the meltwater forcing, and explore possible signals of the climate response to potential LIS meltwater discharges in the early to mid-Holocene. The record shows four abrupt large fluctuations during the 7-9 ka BP (kiloannum before present), reflecting large shifts of the mid-latitude westerlies and the Indian summer monsoon (ISM) over this period, and they corresponded to possible LIS freshwater events documented in other paleoclimate records. Our study suggests that multiple rapid meltwater discharge events might have occurred during the final stage of LIS. The finding implies the possibility of rapid sea level rise and unstable climate in the transition zone between the mid-latitude westerlies and the ISM due to fast polar ice retreat under the anthropogenic global warming.

Fossil mosses are emitting methane after maritime Antarctic glacier retreat

In the extraordinary weather conditions of the austral summer of 2023, fossil mosses thawed out from under the Bellingshausen Ice Dome, King George Island, Southern Shetland Archipelago of maritime Antarctica. At the end of the austral summer, we directly measured greenhouse gas fluxes (CH4 and CO2) from the surface of fossil mosses. We showed that fossil mosses were strong emitters of CH4 and weak emitters of CO2. The real-time measured CH4 emissions reached 0.173 μmol m-2 s-1, which is comparable to CH4 efflux in water bodies or wet tundra in the Arctic.

Climate impacts on tree-ring stable isotopes across the Northern Hemispheric boreal zone

Boreal regions are changing rapidly with anthropogenic global warming. In order to assess risks and impacts of this process, it is crucial to put these observed changes into a long-term perspective. Summer air temperature variability can be well reconstructed from conifer tree rings. While the application of stable isotopes can potentially provide complementary climatic information over different seasons. In this study, we developed new triple stable isotope chronologies in tree-ring cellulose (δ¹³C_trc, δ¹⁸O_trc, δ²H_trc) from a study site in Canada. Additionally, we performed regional aggregated analysis of available stable isotope chronologies from 6 conifers' tree species across high-latitudinal (HL) and - altitudinal (HA) as well as Siberian (SIB) transects of the Northern Hemispheric boreal zone. Our results show that summer air temperature still plays an important role in determining tree-ring isotope variability at 11 out of 24 sites for δ¹³C_trc, 6 out of 18 sites for δ¹⁸O_trc and 1 out of 6 sites for δ²H_trc. Precipitation, relative humidity and vapor pressure deficit are significantly and consistently recorded in both δ¹³C_trc and δ¹⁸O_trc along HL. Summer sunshine duration is captured by all isotopes, mainly for HL and HA transects, indicating an indirect link with an increase in air and leaf temperature. A mixed temperature-precipitation signal is preserved in δ¹³C_trc and δ¹⁸O_trc along SIB transect. The δ²H_trc data obtained for HL-transect provide information not only about growing seasonal moisture and temperature, but also capture autumn, winter and spring sunshine duration signals. We conclude that a combination of triple stable isotopes in tree-ring studies can provide a comprehensive description of climate variability across the boreal forest zone and improve ecohydrological reconstructions.

Hydrogeochemical and isotopic characterization of the Gioia Tauro coastal Plain (Calabria - southern Italy): A multidisciplinary approach for a focused management of vulnerable strategic systems

This work pursues the hydro-geochemical and isotopic characterization of the complex groundwater system of the Gioia Tauro Plain, one of the most important industrialized and agricultural coastal areas of southern Italy. The anthropic pressure exposes the water resources at risk of depletion and quality degradation making the plain groundwater a system of high scientific and social interest. The plain is characterized by a shallow aquifer, mostly recharged by local rains and a deep aquifer apparently less influenced by local precipitation. Both aquifers are mainly Ca-HCO₃ waters except for localized sectors where Na-HCO₃, Na-Cl and Ca-SO₄ waters are present. In deep aquifer, both prolonged interaction with sedimentary rocks, mainly deriving from the erosion of crystalline rocks, and direct cation exchange represent the primary factors controlling the formation of Na-HCO₃ waters. Mixing processes between these waters and either connate brine and/or deep thermal waters contribute to the formation of isolated high salinity Na-Cl-rich waters. In shallow aquifer, inputs of N-rich sewage and agriculture-related contaminants, and SOx emissions in proximity of the harbor are responsible of the increasing nitrate and sulphate concentrations, respectively. The Cl/Br and NO₃/Cl ratios highlight contamination mainly linked to agricultural activities and contribution of wastewater. Along the northern boundary, the warmest groundwater (Na-Cl[SO₄]) were found close to a bend of the main strike-slip fault system, locally favouring the rising of B- and Li-rich deep waters, testifying the influence of geological-structural features on deep water circulation. Despite the high-water demand, a direct marine intrusion is localized in a very restricted area, where we observed an incipient groundwater-seawater mixing (seawater contribution ≤7 %). The qualitative and quantitative conditions of the shallow aquifer still have acceptable levels because of the relatively high recharge inflow. A reliable hydrogeochemical conceptual model, able to explain the compositional variability of the studied waters, is proposed.

Application of hydrochemical and multi-isotopic (87Sr/86Sr, δ13C-DIC, δ2H-H2O, δ18O-H2O) tools to determine contamination sources and processes in the Guadalhorce River Basin, southern Spain

The integrated use of multi-isotopic (⁸⁷Sr/⁸⁶Sr, δ¹³C-DIC, δ²H-H₂O, δ¹⁸O-H₂O) and hydrochemical data was applied in the highly anthropized Guadalhorce river basin, southern Spain, to improve the knowledge about water contamination sources and processes and to achieve improved water resource management. The results obtained highlight the importance of the use of isotopes as tracers of pollutants. DIC, δ²H-H₂O, δ¹⁸O-H₂O and δ¹³C-DIC allowed differentiating two water recharge end members: direct rainwater, infiltrated into the upper and lower detritic aquifers of the sub-basins, and the Guadalhorce dam system, which act as a source in some groundwater and surface waters of the lower sub-basin. ⁸⁷Sr/⁸⁶Sr data supported the existing conclusions in relation to pollution sources in the study area. The Triassic basement (evaporites) of the carbonate and detritic aquifers of the basin generally controls the natural ⁸⁷Sr/⁸⁶Sr composition in waters of the upper sub-basin. Only one groundwater sample reflects the influence of a human organic source (sewage) in its composition. On the other hand, mixing of human inorganic (fertilizers and detergents) strontium sources is required to explain the ⁸⁷Sr/⁸⁶Sr contents of the lower sub-basin waters. Discriminating the use of domestic detergents as another anthropogenic source of strontium and sulphate in waters is a novel finding in this research. The conclusions reached can be extrapolated to other anthropized basins.

Deciphering the origin of groundwater inflow into the Talave tunnel (SE Spain)

The Talave tunnel (TT) is an infrastructure of a major water transfer from the Tajo river basin (center Spain) to the Segura river basin (SE Spain), crossing the Júcar river basin. The tunnel was drilled between 1969 and 1978. It is 32 km long, N/NW-S/SE oriented, has a maximum depth of 320 m, intersects several aquifers, and its southern stretch follows the eastern boundary of the Alcadozo aquifer. The TT drilling perturbed groundwater flow in two river basins, and the induced groundwater inlets generated social and administrative concern lasting until today. The main objectives of this paper are understanding the historical and current tunnel-massif hydrodynamic relationships, and deciphering the origin of groundwater inflow into the tunnel. The first objective was approached analyzing the discharge flow evolution since the drilling until 2016, together with old (1970s) and recent (2014-2017) piezometric data. For the second objective, hydrochemical and isotopic data were generated between 2014 and 2018 from discrete and integrated discharge to the tunnel. Attaining both goals benefited of recent studies on groundwater recharge and functioning in the Alcadozo aquifer. Discharge flows stabilized in the early 2000s. 340 hm3 were drained between 1969 and 2016, producing a piezometric drawdown between 15 and 120 m along the tunnel. The main inflow zones correspond to tectonic fractures concentrated in the middle and southern sections of the tunnel. The existence of a hydrogeological divide between the Júcar and the Segura river basins some 3 km to the N of the watershed divide implies that all groundwater comes from the Segura hydrogeological basin. The isotopes suggest that groundwater comes mainly from regional flow lines originated to the W of the tunnel, together with some local recharge. The effects of tunneling can be used to assess the impacts of imminent groundwater development planned by the Basin Authority.

Isotope and hydrochemical systematics of groundwater from a multi-tiered aquifer in the central parts of Indo-Gangetic Plains, India - Implications for groundwater sustainability and security

The Indo-Gangetic multi-aquifer system provides water supplies to the most populous regions of the Indian subcontinent, however precise knowledge on the sources and dynamics of groundwater is still missing. Environmental isotopes (²H, ¹⁸O, ¹³C, ³H and ¹⁴C) and hydrochemical modeling tools were used in this study in the multi-tiered aquifers underlying the Middle Gangetic Plains (MGP) to investigate the source of recharge, aquifer dynamics and inter-connectivity among aquifers. Within a depth span of 300 m, three aquifers, with contrasting recharge sources and dynamics, were delineated in this Sone-Ganga-Punpun interfluve region, with limited cross-aquifer hydraulic interconnections. The chemistry evolves from Ca-HCO₃ to Na-Ca-HCO₃ in the shallow semiconfined Aquifer-I with a mean transit time of 20-23 years. The dominant recharge to Aquifer-I is from the river inflows and rainwater percolation through paleochannels. The semi-confined to confined Aquifer-II holds fresh quality groundwater with mixed water facies (Mg/Ca-Na-HCO₃). The modeled age of Aquifer-II groundwater is found to be 205-520 years, which is supported by presence of negligible tritium and minor variations in stable isotopes. Outcrop regions of Aquifer-II sediments in the marginal alluvial areas and deep-seated paleochannels in the southwestern part are the potential zones for Aquifer-II recharge. A deep confined Aquifer-III with fresh quality of groundwater is identified below 220 m. This aquifer is characterized by old age (~3.5 to 4.7 ka BP) and enriched δ¹⁸O (-5.7‰). These results along with the existing paleoclimate records of this region infer that Aquifer-III is recharged during an arid climate. The marginal alluvial plains are the probable recharge zones for Aquifer-III. This study helped in conceptualizing the groundwater flow paths in multi-tiered aquifers of MGP. The knowledge and understanding would extend crucial inputs for the sustainable development of deep aquifers not only in the MGP but also in other regions of Indo-Gangetic Plains.

Comprehending the groundwater recharge of a coastal city in humid tropical setting using stable isotopes

Unscientific urbanization in coastal cities has enforced the need for understanding groundwater recharge sources and processes for sustainable development. In this paper, stable isotope compositions of precipitation, groundwater and river water were determined to understand the significant recharge sources of phreatic aquifers in the two differently urbanized environments, viz. urban and peri urban clusters of Kozhikode District, Kerala, and southern India. The two monsoon systems viz., southwest and northeast, are the major source of groundwater recharge in the region, but due to change in landuse pattern in the last decade has significantly altered groundwater recharge. Hence with the aid of stable isotope ratios of water, estimation of point recharge of rainwater to the groundwater in different hydrgeological setting of the area was done. The monsoon rains contributed 35% in the urban alluvial aquifers, up to 39% in the urban laterite aquifer and 42% in the peri urban laterite aquifer. An attempt was made to correlate the landuse changes in the past decade in the region with the groundwater availability.

Using a multi-disciplinary approach to characterize groundwater systems in arid and semi-arid environments: The case of Biskra and Batna regions (NE Algeria)

This study presents a multi-disciplinary approach for the hydrogeological assessment and characterization of water resources in typical arid and semi-arid areas with high anthropogenic pressure, and where environmental conditions and political context prevent extensive field surveys. The use of a three-dimensional (3D) hydrogeological conceptual model, integrating hydrochemical and multi-isotope data, is proposed for the Batna and Biskra area (NE Algeria). Geological data were assembled in 3D geological software, from which a 3D hydrogeological conceptual model was constructed, which included the delineation of groundwater flow directions. The isotopic characterization, including deuterium and oxygen isotopic composition of water (δ²H and δ¹⁸O), and tritium (³H), provided information regarding recharge sources, flow pathways and residence times of groundwaters. Hydrochemical parameters, measured on the same samples, supported the interpretation of isotope data. All data were processed in a geographic information system (GIS) environment. The effectiveness of this approach was tested on a complex system of aquifers with high hydrogeological heterogeneity. Results show the important role the tectonic setting of an area can play in the hydrogeology and hydrochemistry of its principal groundwater systems. The fault network in the study region connects different aquifers, resulting in the mixing of groundwaters. The region most influenced by geological structures is the southern part of the study area, close to Biskra city. In fact, besides a limited contribution of recharge from rain and surface water derived from flood events, the recharge of the Cenozoic aquifers seems to proceed from the ascension of deeper Cretaceous groundwaters through the fault network, as indicated by temperature, bulk chemistry and in particular δ²H, δ¹⁸O and ³H results. In contrast, results suggest that the recharge of the low mineralized Maastrichtian waters is primarily influenced by local precipitation and runoff in the mountainous northern part of the study area. Tritium content, low salinity, and bulk chemistry all suggest such waters to be a mix of pre-bomb (deeper flow-lines within the aquifer) and recent water, with no contribution from the deepest Continental Intercalaire groundwaters. The proposed approach reduces ambiguity about the studied aquifer systems, greatly improves the conceptual understanding of their behaviour, and could provide insights into the vulnerability of the aquifers to different anthropogenic pollution phenomena. The methodology used appears to be a valid tool that could be applied to other geographical areas, to inform the design and implementation of efficient management strategies aimed at improving the quality and availability of water resources. Moreover, three-dimensional modelling methods are becoming increasingly applied to different aspects of groundwater management, to obtain a detailed picture of subsurface conditions.

Snowmelt as a determinant factor in the hydrogeological behaviour of high mountain karst aquifers: The Garcés karst system, Central Pyrenees (Spain)

Time series of environmental tracers (groundwater stable isotope composition, electrical conductivity and temperature) and concentration breakthrough curves of artificial tracers (uranine, eosine, amino-G and naphtionate) have been analyzed to characterize fast preferential and slow matrix in-transit recharge flows in the Paleocene-Eocene limestone aquifer of the Ordesa and Monte Perdido National Park, an alpine karst system drained by a water table cave, a rare hydrological feature in high mountain karst systems with similar characteristics. Snowmelt favors the areal recharge of the system. This process is reflected in the large proportion of groundwater flowing through the connected porosity structure of the karst aquifer, which amounts the 75% of the total system water discharge. From the perspective of water resources recovery, the water capacity of the fissured-porous zone (matrix) represents 99% of the total karst system storage. The volume associated to the karst conduits is very small. The estimated mean travel times are 9 days for conduits and 475 days for connected porosity. These short travel times reveal high vulnerability of the karst system to pollutants in broad sense and a great impact of climate change on the associated water resources.

Annual study of hydraulic characteristics in surface flow constructed wetlands using hydrogen and oxygen stable isotope technology

Complex flow patterns and hydraulic characteristics could reduce the utilization efficiency of constructed wetland (CW), and consequently, its pollutant removal performance. Thus, it is of great importance to explore the internal flow patterns of CWs. Isotopic molecules exist naturally in CWs and have special properties under liquid conditions; using hydrogen and oxygen isotope technology cannot only reduce secondary pollution but also reflect the hydraulic characteristics of CWs. In the present study, the annual variation of isotopic composition in field-scale CW was investigated to evaluate the long-term feasibility of stable isotopic technology characterizing hydraulic flow patterns. The relationship between nutrients concentration distribution and flow pattern variation in CW under different seasons was discussed as well. Results demonstrated that isotope ¹⁸O/¹⁶O distribution could be used to determine the internal flow pattern of CW throughout the year, except for preferential flow area of CWs in winter, since more hydraulic retention time is needed to ensure the change of water isotopes due to the small evaporation in winter. Lower ammonia nitrogen concentration was observed in the stagnant area, while the total phosphorus concentration of the stagnant area increased during winter. And more attention should be paid to aquatic plants during the CW design, since it has significant influence on the hydraulic flow patterns of CW.

Spatial and temporal evolution of groundwater arsenic contamination in the Red River delta, Vietnam: Interplay of mobilisation and retardation processes

Geogenic arsenic (As) contamination of groundwater poses a major threat to global health, particularly in Asia. To mitigate this exposure, groundwater is increasingly extracted from low-As Pleistocene aquifers. This, however, disturbs groundwater flow and potentially draws high-As groundwater into low-As aquifers. Here we report a detailed characterisation of the Van Phuc aquifer in the Red River Delta region, Vietnam, where high-As groundwater from a Holocene aquifer is being drawn into a low-As Pleistocene aquifer. This study includes data from eight years (2010-2017) of groundwater observations to develop an understanding of the spatial and temporal evolution of the redox status and groundwater hydrochemistry. Arsenic concentrations were highly variable (0.5-510 μg/L) over spatial scales of <200 m. Five hydro(geo)chemical zones (indicated as A to E) were identified in the aquifer, each associated with specific As mobilisation and retardation processes. At the riverbank (zone A), As is mobilised from freshly deposited sediments where Fe(III)-reducing conditions occur. Arsenic is then transported across the Holocene aquifer (zone B), where the vertical intrusion of evaporative water, likely enriched in dissolved organic matter, promotes methanogenic conditions and further release of As (zone C). In the redox transition zone at the boundary of the two aquifers (zone D), groundwater arsenic concentrations decrease by sorption and incorporations onto Fe(II) carbonates and Fe(II)/Fe(III) (oxyhydr)oxides under reducing conditions. The sorption/incorporation of As onto Fe(III) minerals at the redox transition and in the Mn(IV)-reducing Pleistocene aquifer (zone E) has consistently kept As concentrations below 10 μg/L for the studied period of 2010-2017, and the location of the redox transition zone does not appear to have propagated significantly. Yet, the largest temporal hydrochemical changes were found in the Pleistocene aquifer caused by groundwater advection from the Holocene aquifer. This is critical and calls for detailed investigations.

How water isotopes (18O, 2H, 3H) within an island freshwater lens respond to changes in rainfall

Coastal aquifers provide an important source of water globally. Understanding how groundwater responds to changes in rainfall recharge is important for sustainable development. To this end, we investigate how water isotopes (¹⁸O, ²H, ³H) and chloride (Cl) concentrations within an island freshwater lens respond under varying rainfall conditions in a region experiencing climate change. Uniquely, this study presents a three year dataset of groundwater collected seasonally between May 2013 and August 2016 from ten wells. Variation in all tracers was observed. The Cl and tritium (³H) show opposing seasonal variation in some sections of the lens, with higher Cl observed in the austral summer when less rainfall occurs and evapotranspiration is highest. The opposite occurs in the austral winter months when ³H increases from atmospheric input via rainfall recharge, and Cl is diluted. An overall decline in ³H values and enrichment in stable water isotopes over the study period was also observed. This study shows that understanding groundwater of freshwater lenses should not rely on a single sampling campaign because seasonal variability is large. The identification of a dual recharge regime, with contributions from both winter rainfall and episodic events, has important implications for understanding the future fate of the freshwater lens on Rottnest Island. The finding that episodic rainfall is a major contributor to groundwater recharge is important and can only be assessed with a multi-year isotope dataset for groundwater and rainfall.

Identification of suites of traits that explains drought resistance and phenological patterns of plants in a semi-arid grassland community

Grassland ecosystems are comprised of plants that occupy a wide array of phenological niches and vary considerably in their ability to resist the stress of seasonal soil-water deficits. Yet, the link between plant drought resistance and phenology remains unclear in perennial grassland ecosystems. To evaluate the role of soil water availability and plant drought tolerance in driving phenology, we measured leaf hydraulic conductance (K_sat), resistance to hydraulic failure (P₅₀), leaf gas exchange, plant and soil water stable isotope ratios (δ¹⁸O), and several phenology metrics on ten perennial herbaceous species in mixed-grass prairie. The interaction between P₅₀ and δ¹⁸O of xylem water explained 67% of differences in phenology, with lower P₅₀ values associated with later season activity, but only among shallow-rooted species. In addition, stomatal control and high water-use efficiency also contributed to the late flowering and late senescence strategies of plants that had low P₅₀ values and relied upon shallow soil water. Alternatively, plants with deeper roots did not possess drought-tolerant leaves, but had high hydraulic efficiency, contributing to their ability to efficiently move water longer distances while maintaining leaf water potential at relatively high values. The suites of traits that characterize these contrasting strategies provide a mechanistic link between phenology and plant-water relations; thus, these traits could help predict grassland community responses to changes in water availability, both temporally and vertically within the soil profile.

Geographic classification of U.S. Washington State wines using elemental and water isotope composition

Wine fraud leaves wineries vulnerable to damage in reputation and potential lost revenue. To reduce this risk for wines from Washington State (WA), USA, advanced analytical instrumentation and statistical methods were employed to geographically classify 133 wines from 4 major wine producing regions, including 70 wines from WA. Analyses of 37 elements and 2 water isotopes were performed with Triple Quadrupole Inductively Coupled Plasma Mass Spectrometry and Cavity Ring-Down Spectroscopy, respectively. Linear discriminant analysis resulted in 96.2% discrimination, achieved with 11 parameters (Mn, Zn, Pb, Ni, As, D/H, La, Ce, Si, Zr and Sr) that were linearly combined into 3 functions. WA wines were uniquely distinguished in large part with low D/H ratios and Mn concentrations derived from the isotopically light precipitation and volcanic loess soils encountered in this region, respectively. This study is the first of its kind to focus on the authentication of WA wines.

Arsenic-rich shallow groundwater in sandy aquifer systems buffered by rising carbonate waters: A geochemical case study from Mannar Island, Sri Lanka

Major ion, trace elements, and stable isotope analyses were performed on groundwater samples collected from Mannar Island in the northern Indian Ocean. Arsenic concentrations up to 34μg/L have been observed in groundwater samples from the island. In addition, 23% of extensively used shallow drinking water wells showed higher arsenic contents than the recommended value by the World Health Organization (10μg/L). Groundwater in the island showed pH values between 6.9 and 8.9 and was dominated by Na⁺, K⁺, Ca²⁺, Mg²⁺, HCO₃^-, Cl^- and SO₄^2-. The δ¹⁸O_H2O and δ²H_H2O composition of most groundwater plotted very close to the local meteoric waterline, however, some wells showed enriched isotope compositions that are most likely due to evaporation. Sea water intrusion in this island was likely of minor importance as indicated by the major ion composition. An approximated mass balance calculation using chloride concentrations indicated that out of the 35 investigated wells only 6 near-shore wells were influenced by sea water intrusion up to about 15%. Even though this is a sandy aquifer, groundwaters were characterized with higher contents of dissolved inorganic carbon (DIC) (2.11-10.9mmol/L). The corresponding δ¹³C_DIC values varied from -19.4‰ to -6.5‰. Except for a few samples with values approaching -20‰, these isotope values are more typical for carbonate dissolution and equilibration of CO₂ in the aquifer. This study shows that the underlying carbonate system may buffer the aqueous geochemistry of the groundwater on the island. The high arsenic content in groundwater may have been mobilized through reductive dissolution of Fe-Mn oxides and oxy-hydroxides that are coated on sandy aquifer materials. The lower content of DOC (0.2-1.5mmol/L) provides evidence for the restricted formation of pyrite in the aquifer.

Seasonal variation of δ18O and δ2H in leaf water of Fagus sylvatica L. and related water compartments

The study aims to assess variability in leaf water isotopic enrichment occurring in the field under natural conditions. We focused on seasonal variation and difference between sun-exposed and shaded leaves. Isotopic composition (δ¹⁸O, δ²H) of leaf water was monitored in a beech tree (Fagus sylvatica L.) growing in the forest-meadow ecotone together with δ¹⁸O (²H) of water compartments which are in close relation to this signal, namely twig and soil water. The sampling was carried out in approximately two-week intervals during five consecutive vegetation seasons. The δ¹⁸O (²H) data showed a distinct seasonal pattern and a consistency in relative differences between the seasons and sample categories. Leaf water was the most isotopically enriched water compartment. The leaf water enrichment decreased toward the autumn reflecting the change in δ¹⁸O (²H) of source water and evaporative demands. The soil and twig water isotopic signal was depleted against current precipitation as it partly retained the isotopic signature from winter precipitation however the seasonal pattern of soil and twig water followed that of precipitation. No significant differences between sun-exposed and shaded samples were detected. Nevertheless, the observed strong seasonal pattern of isotope composition of leaf, twig and soil water should be taken into account when using leaf water enrichment for further calculations or modeling.

Vertical variation in the amplitude of the seasonal isotopic content of rainfall as a tool to jointly estimate the groundwater recharge zone and transit times in the Ordesa and Monte Perdido National Park aquifer system, north-eastern Spain

The time series of stable water isotope composition relative to meteorological stations and springs located in the high mountainous zone of the Ordesa and Monte Perdido National Park are analyzed in order to study how the seasonal isotopic content of precipitation propagates through the hydrogeological system in terms of the aquifer recharge zone elevation and transit time. The amplitude of the seasonal isotopic composition of precipitation and the mean isotopic content in rainfall vary along a vertical transect, with altitudinal slopes for δ¹⁸O of 0.9‰/km for seasonal amplitude and -2.2‰/km for isotopic content. The main recharge zone elevation for the sampled springs is between 1950 and 2600m·a.s.l. The water transit time for the sampled springs ranges from 1.1 to 4.5yr, with an average value of 1.85yr and a standard deviation of 0.8yr. The hydrological system tends to behave as a mixing reservoir.