Isotope hydrology of the Allt a' Mharcaidh catchment,Cairngorms, Scotland: implications for hydrological pathways and residence times

The hydrology of oxygen‐18 (¹⁸O) isotopes was monitored between 1995 and 1998 in the Allt a' Mharcaidh catchment in the Cairngorm Mountains, Scotland. Precipitation (mean δ¹⁸O=−7·69‰) exhibited strong seasonal variation in δ¹⁸O values over the study period, ranging from −2·47‰ in the summer to −20·93‰ in the winter months. As expected, such variation was substantially damped in stream waters, which had a mean and range of δ¹⁸O of −9·56‰ and −8·45 to −10·44‰, respectively. Despite this, oxygen‐18 proved a useful tracer and streamwater δ¹⁸O variations could be explained in terms of a two‐component mixing model, involving a seasonally variable δ¹⁸O signature in storm runoff, mixing with groundwater characterized by relatively stable δ¹⁸O levels. Variations in soil water δ¹⁸O implied the routing of depleted spring snowmelt and enriched summer rainfall into streamwaters, probably by near‐surface hydrological pathways in peaty soils. The relatively stable isotope composition of baseflows is consistent with effective mixing processes in shallow aquifers at the catchment scale. Examination of the seasonal variation in δ¹⁸O levels in various catchment waters provided a first approximation of mean residence times in the major hydrological stores. Preliminary estimates are 0·2–0·8 years for near‐surface soil water that contributes to storm runoff and 2 and >5 years for shallow and deeper groundwater, respectively. These ¹⁸O data sets provide further evidence that the influence of groundwater on the hydrology and hydrochemistry of upland catchments has been underestimated. Copyright © 2000 John Wiley & Sons, Ltd.     

Estimation of mean residence times of subsurface waters using seasonal variation in deuterium excess in a small headwater catchment in Japan

We measured deuterium excess (d = δD ? 8δ¹⁸O) in throughfall, groundwater, soil water, spring water, and stream water for 3 years in a small headwater catchment (Matsuzawa, 0·68 ha) in the Kiryu Experimental Watershed in Japan. The d value represents a kinetic effect produced when water evaporates. The d value of the throughfall showed a sinusoidal change (amplitude: 6·9‰ relative to Vienna standard mean ocean water (V‐SMOW)) derived from seasonal changes in the source of water vapour. The amplitude of this sinusoidal change was attenuated to 1·3–6·9‰ V‐SMOW in soil water, groundwater, spring water, and stream water. It is thought that these attenuations derive from hydrodynamic transport processes in the subsurface and mixing processes at an outflow point (stream or spring) or a well. The mean residence time (MRT) of water was estimated from d value variations using an exponential‐piston flow model and a dispersion model. MRTs for soil water were 0–5 months and were not necessarily proportional to the depth. This may imply the existence of bypass flow in the soil. Groundwater in the hillslope zone had short residence times, similar to those of the soil water. For groundwater in the saturated zone near the spring outflow point, the MRTs differed between shallow and deeper groundwater; shallow groundwater had a shorter residence time (5–8 months) than deeper groundwater (more than 9 months). The MRT of stream water (8–9 months) was between that of shallow groundwater near the spring and deeper groundwater near the spring. The seasonal variation in the d value of precipitation arises from changes in isotopic water vapour composition associated with seasonal activity of the Asian monsoon mechanism. The d value is probably an effective tracer for estimating the MRT of subsurface water not only in Japan, but also in other East Asian countries influenced by the Asian monsoon. Copyright © 2006 John Wiley & Sons, Ltd.     

Temporal and spatial variation of water stable isotopes (18O and 2H) in the Kaidu River basin,Northwestern China

Water resources are the most critical factors to ecology and society in arid basins, such as Kaidu River basin. Isotope technique was convenient to trace this process and reveal the influence from the environment. In this paper, we try to investigate the temporal and spatial characteristics in stable isotope (¹⁸O and ²H) of surface water and groundwater in Kaidu River. Through the water stable isotope composition measurement, spatial and temporal characteristics of deuterium (δ²H) and oxygen 18 (δ¹⁸O) were analysed. It is revealed that (1) comparing the stream water line with the groundwater line and local meteorological water line of Urumqi City, it is found that the contribution of precipitation to surface water in stream runoff is the main source, whereas the surface water is the main source of groundwater. Groundwater is mainly drainage of surface runoff in the river; (2) in the main stream of Kaidu River, the spatial variability of river water showed a ‘heavier‐lighter‐heavier’ change along with the main stream for δ¹⁸O, and temporal variability showed higher in summer and lower in winter; (3) the δ¹⁸O and δ²H values of groundwater samples ranged from ?11.36 to ?7.97‰ and ?73.45 to ?60.05‰, respectively. There is an increasing trend of isotopic values along the groundwater flow path. The seasonal fluctuation of δ¹⁸O is not clear in most samples. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydrological processes controlling ground and surface water flow from a hypermaritime forest–peatland complex,Diana Lake Provincial Park,British Columbia,Canada

The proposed harvesting of previously undeveloped forests in north coastal British Columbia requires an understanding of hydrological responses. Hydrometric and isotopic techniques were used to examine the hydrological linkages between meteoric inputs to the surface‐groundwater system and runoff response patterns of a forest‐peatland complex. Quickflow accounted for 72–91% of peak storm discharge. The runoff ratio was lowest for open peatland areas with thick organic horizons (0·02–0·05) due to low topographic gradients and many surface depressions capable of retaining surface water. Runoff ratio increased comparatively for ephemeral surface seep flows (0·06–0·40) and was greatest in steeply sloping forest communities with more permeable soils (0·33–0·69). The dominant mechanism for runoff generation was saturated shallow subsurface flow. Groundwater fluxes from the organic horizon of seeps (1·70–1·72 m³ day^?1 m^?1) were an important component of quickflow. The homogeneous δ²H? δ¹⁸O composition of groundwater indicated attenuation of the seasonal rainfall signal by mixing during recharge. The positive correlation (r² = 0·64 and 0·38, α = 0·05) between slope index and δ¹⁸O values in groundwater suggests that the spatial pattern in the δ¹⁸O composition along the forest‐peatland complex is influenced by topography and provides evidence that topographic indices may be used to predict groundwater residence time. Copyright © 2006 John Wiley & Sons, Ltd.     

Estimating groundwater recharge from water isotope (δ2H, δ18O) depth profiles in the Densu River basin,Ghana

Abstract

Accurate estimation of groundwater recharge is essential for the proper management of aquifers. A study of water isotope (δ²H, δ¹⁸O) depth profiles was carried out to estimate groundwater recharge in the Densu River basin in Ghana, at three chosen observation sites that differ in their altitude, geology, climate and vegetation. Water isotopes and water contents were analysed with depth to determine water flow in the unsaturated zone. The measured data showed isotope enrichment in the pore water near the soil surface due to evaporation. Seasonal variations in the isotope signal of the pore water were also observed to a depth of 2.75 m. Below that depth, the seasonal variation of the isotope signal was attenuated due to diffusion/dispersion and low water flow velocities. Groundwater recharge rates were determined by numerical modelling of the unsaturated water flow and water isotope transport. Different groundwater recharge rates were computed at the three observation sites and were found to vary between 94 and 182 mm/year (± max. 7%). Further, the approximate peak-shift method was applied to give information about groundwater recharge rates. Although this simple method neglects variations in flow conditions and only considers advective transport, it yielded mean groundwater recharge rates of 110–250 mm/year (± max. 30%), which were in the same order of magnitude as computed numerical modelling values. Integrating these site-specific groundwater recharge rates to the whole catchment indicates that more water is potentially renewed than consumed nowadays. With increases in population and irrigation, more clean water is required, and knowledge about groundwater recharge rates – essential for improving the groundwater management in the Densu River basin – can be easily obtained by measuring water isotope depth profiles and applying a simple peak-shift approach.

Citation Adomako, D., Maloszewski, P., Stumpp, C., Osae, S. & Akiti, T. T. (2010) Estimating groundwater recharge from water isotope (δ²H, δ¹⁸O) depth profiles in the Densu River basin, Ghana. Hydrol. Sci. J. 55(8), 1405–1416.     

Isotopic composition (δ18O and δD) of precipitation and groundwater in a semi‐arid,mountainous area (Guadiana Menor basin,Southeast Spain)

We characterize the precipitation and groundwater in a mountainous (peaks slightly above 3000 m a.s.l.), semi‐arid river basin in SE Spain in terms of the isotopes ¹⁸O and ²H. This basin, with an extension of about 7000 km², is an ideal site for such a study because fronts from the Atlantic and the Mediterranean converge here. Much of the land is farmed and irrigated both by groundwater and runoff water collected in reservoirs. A total of approximately 100 water samples from precipitation and 300 from groundwater have been analysed. To sample precipitation we set up a network of 39 stations at different altitudes (800–1700 m a.s.l.), with which we were able to collect the rain and snowfall from 29 separate events between July 2005 and April 2007 and take monthly samples during the periods of maximum recharge of the aquifers. To characterize the groundwater we set up a control network of 43 points (23 springs and 20 wells) to sample every 3 months the main aquifers and both the thermal and non‐thermal groundwater. We also sampled two shallow‐water sites (a reservoir and a river). The isotope composition of the precipitation forms a local meteoric water line (LMWL) characterized by the equation δD = 7·72δ¹⁸O + 9·90, with mean values for δ¹⁸O and δD of − 10·28‰ and − 69·33‰, respectively, and 12·9‰ for the d‐excess value. To correlate the isotope composition of the rainfall water with groundwater we calculated the weighted local meteoric water line (WLMWL), characterized by the equation δD = 7·40δ¹⁸O + 7·24, which takes into account the quantity of water precipitated during each event. These values of (dδD/dδ¹⁸O)< 8 and d‐excess (δD–8δ¹⁸O)< 10 in each curve bear witness to the ‘amount effect’, an effect which is more manifest between May and September, when the ground temperature is higher. Other effects noted in the basin were those of altitude and the continental influence. The isotopic compositions of the groundwater are represented by the equation δD = 4·79δ¹⁸O − 18·64. The groundwater is richer in heavy isotopes than the rainfall, with mean values of − 8·48‰ for δ¹⁸O and − 59·27‰ for δD. The isotope enrichment processes detected include a higher rate of evaporation from detrital aquifers than from carbonate ones, the effects of recharging aquifers from irrigation return flow and/or from reservoirs' leakage and enrichment in δ¹⁸O from thermal water. Copyright © 2010 John Wiley & Sons, Ltd.     

Continuous monitoring of stream δ18O and δ2H and stormflow hydrograph separation using laser spectrometry in an agricultural catchment

A portable Wavelength Scanned‐Cavity Ring‐Down Spectrometer (Picarro L2120) fitted with a diffusion sampler (DS‐CRDS) was used for the first time to continuously measure δ¹⁸O and δ²H of stream water. The experiment took place during a storm event in a wet tropical agricultural catchment in north‐eastern Australia. At a temporal resolution of one minute, the DS‐CRDS measured 2160 δ¹⁸O and δ²H values continuously over a period of 36 h with a precision of ±0.08 and 0.5‰ for δ¹⁸O and δ²H, respectively. Four main advantages in using high temporal resolution stream δ¹⁸O and δ²H data during a storm event are highlighted from this study. First, they enabled us to separate components of the hydrograph, which was not possible using high temporal resolution electrical conductivity data that represented changes in solute transfers during the storm event rather than physical hydrological processes. The results from the hydrograph separation confirm fast groundwater contribution to the stream, with the first 5 h of increases in stream discharge comprising over 70% pre‐event water. Second, the high temporal resolution stream δ¹⁸O and δ²H data allowed us to detect a short‐lived reversal in stream isotopic values (δ¹⁸O increase by 0.4‰ over 9 min), which was observed immediately after the heavy rainfall period. Third, δ¹⁸O values were used to calculate a time lag of 20 min between the physical and chemical stream responses during the storm event. Finally, the hydrograph separation highlights the role of event waters in the runoff transfers of herbicides and nutrients from this heavily cultivated catchment to the Great Barrier Reef. Copyright © 2015 John Wiley & Sons, Ltd.     

Stable isotope tracers as diagnostic tools in upscaling flow path understanding and residence time estimates in a mountainous mesoscale catchment

δ¹⁸O measurements of precipitation and stream waters were used as a natural tracer to investigate hydrological pathways and residence times in the River Feshie, a complex mesoscale (231 km²) catchment in the Cairngorm Mountains of Scotland. Precipitation δ¹⁸O exhibited strong seasonal variation over the 2001–02 hydrological year, ranging from −6·9‰ in the summer, to −12·0‰ during winter snowfalls (mean δ¹⁸O −9·59‰). Although damped, this seasonality was reflected in stream water outputs at seven sampling sites in the catchment, allowing δ¹⁸O variations to be used to infer hydrological source areas. Thus, stream water δ¹⁸O was generally controlled by a seasonally variable storm flow end member, mixing with groundwater of more constant isotopic composition. Periodic regression analysis allowed the differences in this mixing process between monitoring subcatchments to be assessed more quantitatively to provide a preliminary estimate of mean stream water residence time. This demonstrated the importance of responsive hydrological pathways associated with peat and shallow alpine soils in the headwater subcatchments in producing seasonally variable runoff with short mean residence times (33–113 days). In contrast, other tributaries with more freely draining soils and larger groundwater storage in shallow aquifers provided more effective mixing of variable precipitation inputs, resulting in longer residence time estimates (178–445 days). The mean residence time of runoff leaving the Feshie catchment reflected an integration of these contrasting influences (110–200 days). These insights from δ¹⁸O measurements extend the hydrological understanding of the Feshie catchment gained from other hydrochemical tracers, and demonstrate the utility of isotope tracers in investigating hydrological processes at the mesoscale. Copyright © 2005 John Wiley & Sons, Ltd.     

Determination of groundwater discharge rates and water residence time of groundwater‐fed lakes by stable isotopes of water (18O, 2H) and radon (222Rn) mass balances

Lacustrine groundwater discharge (LGD) and the related water residence time are crucial parameters for quantifying lake matter budgets and assessing its vulnerability to contaminant input. Our approach utilizes the stable isotopes of water (δ¹⁸O, δ²H) and the radioisotope radon (²²²Rn) for determining long‐term average and short‐term snapshots in LGD. We conducted isotope balances for the 0.5‐km² Lake Ammelshainer See (Germany) based on measurements of lake isotope inventories and groundwater composition accompanied by good quality and comprehensive long‐term meteorological and isotopic data (precipitation) from nearby monitoring stations. The results from the steady‐state annual isotope balances that rely on only two sampling campaigns are consistent for both δ¹⁸O and δ²H and suggested an overall long‐term average LGD rate that was used to infer the water residence time of the lake. These findings were supported by the good agreement of the simulated LGD‐driven annual cycles of δ¹⁸O and δ²H lake inventories with the observed lake isotope inventories. However, radon mass balances revealed lower values that might be the result of seasonal LGD variability. For obtaining further insights into possible seasonal variability of groundwater–lake interaction, stable water isotope and radon mass balances could be conducted more frequently (e.g., monthly) in order to use the derived groundwater discharge rates as input for time‐variant isotope balances.     

Insights into groundwater salinization from hydrogeochemical and isotopic evidence in an arid inland basin

In the Manas River basin (MRB), groundwater salinization has become a major concern, impeding groundwater use considerably. Isotopic and hydrogeochemical characteristics of 73 groundwater and 11 surface water samples from the basin were analysed to determine the salinization process and potential sources of salinity. Groundwater salinity ranged from 0.2 to 11.91 g/L, and high salinities were generally located in the discharge area, arable land irrigated by groundwater, and depression cone area. The quantitative contributions of the evaporation effect were calculated, and the various groundwater contributions of transpiration, mineral dissolution, and agricultural irrigation were identified using hydrogeochemical diagrams and δD and δ¹⁸O compositions of the groundwater and surface water samples. The average evaporation contribution ratios to salinity were 5.87% and 32.7% in groundwater and surface water, respectively. From the piedmont plain to the desert plain, the average groundwater loss by evaporation increased from 7% to 29%. However, the increases in salinity by evaporation were small according to the deuterium excess signals. Mineral dissolution, transpiration, and agricultural irrigation activities were the major causes of groundwater salinization. Isotopic information revealed that river leakage quickly infiltrated into aquifers in the piedmont area with weak evaporation effects. The recharge water interacted with the sediments and dissolved minerals and subsequently increased the salinity along the flow path. In the irrigation land, shallow groundwater salinity and Cl^? concentrations increased but not δ¹⁸O, suggesting that both the leaching of soil salts due to irrigation and transpiration effect dominated in controlling the hydrogeochemistry. Depleted δ¹⁸O and high Cl^? concentrations in the middle and deep groundwater revealed the combined effects of mixing with paleo‐water and mineral dissolution with a long residence time. These results could contribute to the management of groundwater sources and future utilization programs in the MRB and similar areas.     

Groundwater mixing between different aquifer types in a complex structural setting discerned by elemental and stable isotope geochemistry

Karst aquifers are well known for their intricate stratigraphy and geologic structures, which make groundwater characterization challenging because flowpaths and recharge sources are complex and difficult to evaluate. Geochemical data, collected from ten closely spaced production wells constructed in two karst aquifers (Bangor Limestone (Mb) and Tuscumbia Limestone/Fort Payne Chert (Mftp)) in Trussville, north‐central Alabama, illustrate two distinctive groundwater end‐members: (1) higher major ion, dissolved inorganic carbon, conductivity, alkalinity concentrations, heavier δ¹³C ratios (max: −10.2 ± 0.2‰ Vienna Pee Dee Belemnite (PDB)) and lower residence times (mean: 19.5 ± 2 years, n = 2) of groundwater in the Mb aquifer and (2) lower constituent concentrations, lighter δ¹³C ratios (min: −13.4 ± 0.2‰ PDB) and longer residence times of groundwater (mean: 23.6 ± 2 years, n = 4) in the Mftp aquifer. Summer and fall data and the binary mixing model show aquifer inter‐flow mixing along solution fractures and confirms the distinctive groundwater geochemistry of the two aquifers. Lowering of static water levels over the summer (drawdown from 2 to 5.2 m) leads to more reducing groundwater conditions (lower Eh values) and slightly enriched δ¹⁸O and δD ratios during the fall [δ¹⁸O: −4.8 ± 0.1 to −5.4 ± 0.1‰ Vienna Standard Mean Oceanic Water (VSMOW), n = 9; δD: −25.4 ± 1 to −27.4 ± 1‰ VSMOW, n = 9] when compared with summer season samples (δ¹⁸O: −5.1 ± 0.1 to −5.7 ± 0.1‰ VSMOW, n = 11; δD: −25.0 ± 1 to −30.6 ± 1‰ VSMOW, n = 11). GIS analyses confirm the localized origin of recharge to the investigated aquifers. The combination of GIS, field parameters and geochemistry analyses can be successfully used to identify recharge sources, evaluate groundwater flow and transport pathways and to improve understanding of how groundwater withdrawals impact the sustainability and susceptibility to contamination of karst aquifers. Copyright © 2015 John Wiley & Sons, Ltd.     

Groundwater and surface water connectivity within the recharge area of Guarani aquifer system during El Niño 2014–2016

Recharge areas of the Guarani Aquifer System (GAS) are particularly sensitive and vulnerable to climate variability; therefore, the understanding of infiltration mechanisms for aquifer recharge and surface run‐off generation represent a relevant issue for water resources management in the southeastern portion of the Brazilian territory, particularly in the Jacaré‐Pepira River watershed. The main purpose of this study is to understand the interactions between precipitation, surface water, and groundwater using stable isotopes during the strong 2014–2016 El Niño Southern Oscillation event. The large variation in the isotopic composition of precipitation (from ?9.26‰ to +0.02‰ for δ¹⁸O and from ?63.3‰ to +17.6‰ for δ²H), mainly associated with regional climatic features, was not reflected in the isotopic composition of surface water (from ?7.84‰ to ?5.83‰ for δ¹⁸O and from ?49.7‰ to +33.6‰ for δ²H), mainly due to the monthly sampling frequency, and groundwater (from ?7.04‰ to ?7.76‰ for δ¹⁸O and from ?49.5‰ to ?44.7‰ for δ²H), which exhibited less variation throughout the year. However, variations in deuterium excess (d‐excess) in groundwater and surface water suggest the occurrence of strong secondary evaporation during the infiltration process, corresponding with groundwater level recovery. Similar isotopic composition in groundwater and surface water, as well as the same temporal variations in d‐excess and line‐conditioned excess denote the strong connectivity between these two reservoirs during baseflow recession periods. Isotopic mass balance modelling and hydrograph separation estimate that the groundwater contribution varied between 70% and 80%, however, during peak flows, the isotopic mass balance tends to overestimate the groundwater contribution when compared with the other hydrograph separation methods. Our findings indicate that the application of isotopic mass balance methods for ungauged rivers draining large groundwater reservoirs, such as the GAS outcrop, could provide a powerful tool for hydrological studies in the future, helping in the identification of flow contributions to river discharge draining these areas.     

Multiple‐indicator study of groundwater flow and chemistry and the impacts of river and paddy water on groundwater in the alluvial fan of the Tedori River,Japan

To investigate the source, flow paths, and chemistry of rich resources of high‐quality, shallow groundwater in the alluvial fan between the Tedori and Sai rivers in central Japan, we analysed stable isotope ratios of H, O, and Sr and concentrations of major dissolved ions and trace elements in groundwater, river water, and paddy water. The ⁸⁷Sr/⁸⁶Sr ratios of the groundwater are related to near‐surface geology: groundwater in sediment from the Tedori River has high ⁸⁷Sr/⁸⁶Sr ratios (>0.711), whereas that from the Sai River in the north of the fan has low ⁸⁷Sr/⁸⁶Sr ratios (<0.711). δ²H and δ¹⁸O values and ⁸⁷Sr/⁸⁶Sr ratios indicate that groundwater in the central and southern fans is recharged by the Tedori River, whereas recharge in the north is from the Sai River. Mg²⁺, Ca²⁺, Sr²⁺, HCO₃^?, and SO₄^2? concentrations and δ²H and δ¹⁸O values in the groundwater are high in the central fan and, except for the northern area, tend to increase with distance from the Tedori River. There are linear relationships between ⁸⁷Sr/⁸⁶Sr ratio and the reciprocal concentrations of Sr²⁺, Mg²⁺, and Ca²⁺. These geochemical characteristics suggest that as groundwater recharged from the Tedori River flows towards the central fan, it mixes with waters derived from precipitation and paddy water that have become enriched in these components during downward infiltration. These results are consistent with our hydrological analysis and numerical simulation of groundwater flow, thus verifying the validity of the model we used in our simulation of groundwater flow. Copyright © 2016 John Wiley & Sons, Ltd.     

Identification of groundwater recharge sources and processes in a heterogeneous alluvial aquifer: results from multi‐level monitoring of hydrochemistry and environmental isotopes in a riverside agricultural area in Korea

We evaluated sources and pathways of groundwater recharge for a heterogeneous alluvial aquifer beneath an agricultural field, based on multi‐level monitoring of hydrochemistry and environmental isotopes of a riverside groundwater system at Buyeo, Korea. Two distinct groundwater zones were identified with depth: (1) a shallow oxic groundwater zone, characterized by elevated concentrations of NO₃^? and (2) a deeper (>10–14 m from the ground surface) sub‐oxic groundwater zone with high concentrations of dissolved Fe, silica, and HCO₃^?, but little nitrate. The change of redox zones occurred at a depth where the aquifer sediments change from an upper sandy stratum to a silty stratum with mud caps. The δ¹⁸O and δ²H values of groundwater were also different between the two zones. Hydrochemical and δ¹⁸O? δ²H data of oxic groundwater are similar to those of soil water. This illustrates that recharge of oxic groundwater mainly occurs through direct infiltration of rain and irrigation water in the sandy soil area where vegetable cropping with abundant fertilizer use is predominant. Oxic groundwater is therefore severely contaminated by agrochemical pollutants such as nitrate. In contrast, deeper sub‐oxic groundwater contains only small amounts of dissolved oxygen (DO) and NO₃^?. The ³H contents and elevated silica concentrations in sub‐oxic groundwater indicate a somewhat longer mean residence time of groundwater within this part of the aquifer. Sub‐oxic groundwater was also characterized by higher δ¹⁸O and δ²H values and lower d‐excess values, indicating significant evaporation during recharge. We suggest that recharge of sub‐oxic groundwater occurs in the areas of paddy rice fields where standing irrigation and rain water are affected by strong evaporation, and that reducing conditions develop during subsequent sub‐surface infiltration. This study illustrates the existence of two groundwater bodies with different recharge processes within an alluvial aquifer. Copyright © 2009 John Wiley & Sons, Ltd.     

Oxygen‐18 dynamics in precipitation and streamflow in a semi‐arid agricultural watershed,Eastern Washington,USA

Understanding flow pathways and mechanisms that generate streamflow is important to understanding agrochemical contamination in surface waters in agricultural watersheds. Two environmental tracers, δ¹⁸O and electrical conductivity (EC), were monitored in tile drainage (draining 12 ha) and stream water (draining nested catchments of 6‐5700 ha) from 2000 to 2008 in the semi‐arid agricultural Missouri Flat Creek (MFC) watershed, near Pullman Washington, USA. Tile drainage and streamflow generated in the watershed were found to have baseline δ¹⁸O value of ?14·7‰ (VSMOW) year round. Winter precipitation accounted for 67% of total annual precipitation and was found to dominate streamflow, tile drainage, and groundwater recharge. ‘Old’ and ‘new’ water partitioning in streamflow were not identifiable using δ¹⁸O, but seasonal shifts of nitrate‐corrected EC suggest that deep soil pathways primarily generated summer streamflow (mean EC 250 µS/cm) while shallow soil pathways dominated streamflow generation during winter (EC declining as low as 100 µS/cm). Using summer isotopic and EC excursions from tile drainage in larger catchment (4700‐5700 ha) stream waters, summer in‐stream evaporation fractions were estimated to be from 20% to 40%, with the greatest evaporation occurring from August to October. Seasonal watershed and environmental tracer dynamics in the MFC watershed appeared to be similar to those at larger watershed scales in the Palouse River basin. A 0·9‰ enrichment, in shallow groundwater drained to streams (tile drainage and soil seepage), of δ¹⁸O values from 2000 to 2008 may be evidence of altered precipitation conditions due to the Pacific Decadal Oscillation (PDO) in the Inland Northwest. Copyright © 2009 John Wiley & Sons, Ltd.     

Using stable isotopes to understand seasonal and interannual dynamics in moisture sources and atmospheric circulation in precipitation

Stable isotopes in precipitation are useful tracers to strengthen understanding of climate change and hydrological processes. In this study, the moisture sources of 190 precipitation events in Beijing were analysed using the Hybrid Single‐particle Lagrangian Integrated Trajectory model, based on which we studied the relation between variations in precipitation δ¹⁸O and dynamics in moisture sources and atmospheric circulation in seasonal and interannual timescales. Categorization of 7 groups of moisture sources was performed, among which oceanic moisture sources presented lower δ¹⁸O in precipitation than continental moisture sources. The results show that seasonal variations of precipitation δ¹⁸O were caused by changes of moisture sources. In summer, moisture from proximal oceans dominated vapour transport to Beijing due to increasing monsoon strength and resulted in a relatively small variation of precipitation δ¹⁸O. At the interannual timescale, the variations of δ¹⁸O in summer precipitation were related to dynamics in oceanic moistures, showing depleted values when the contribution of oceanic moistures, especially the proportion of long‐distance oceanic moisture, was high. Further analysis indicated that changes of oceanic moisture sources were controlled by the strength of summer monsoons. These findings address the complexity of moisture sources in midlatitude monsoon areas and suggest that isotopic signals in precipitation have the potential to deduce changes in moisture sources and atmospheric circulation and can therefore serve for palaeoclimate reconstruction.     

Using stable isotopes of surface water and groundwater to quantify moisture sources across the Yellow River source region

Characterization of stable isotope compositions (δ²H and δ¹⁸O) of surface water and groundwater in a catchment is critical for refining moisture sources and establishing modern isotope–elevation relationships for paleoelevation reconstructions. There is no consensus on the moisture sources of precipitation in the Yellow River source region during summer season. This study presents δ²H and δ¹⁸O data from 111 water samples collected from tributaries, mainstream, lakes, and groundwater across the Yellow River source region during summertime. Measured δ¹⁸O values of the tributaries range from ?13.5‰ to ?5.8‰ with an average of ?11.0‰. Measured δ¹⁸O values of the groundwater samples range from ?12.7‰ to ?10.5‰ with an average of ?11.9‰. The δ¹⁸O data of tributary waters display a northward increase of 1.66‰ per degree latitude. The δ¹⁸O data and d‐excess values imply that moisture sources of the Yellow River source region during summertime are mainly from the mixing of the Indian Summer Monsoon and the Westerlies, local water recycling, and subcloud evaporation. Analysis of tributary δ¹⁸O data from the Yellow River source region and streamwater and precipitation δ¹⁸O data from its surrounding areas leads to a best‐fit second‐order polynomial relationship between δ¹⁸O and elevation over a 4,600 m elevation range. A δ¹⁸O elevation gradient of ?1.6‰/km is also established using these data, and the gradient is in consistence with the δ¹⁸O elevation gradient of north and eastern plateau. Such relationships can be used for paleoelevation reconstructions in the Yellow River source region.     

Water isotopic variability in Mallorca: a path to understanding past changes in hydroclimate

This paper reports the first results on δ¹⁸O and δ²H analysis of precipitations, cave drip waters, and groundwaters from sites in Mallorca (Balearic Islands, western Mediterranean), a key region for paleoclimate studies. Understanding the isotopic variability and the sources of moisture in modern climate systems is required to develop speleothem isotope‐based climate reconstructions. The stable isotopic composition of precipitation was analysed in samples collected between March 2012 and March 2013. The values are in the range reported by GNIP Palma station. Based on these results, the local meteoric water line (LMWL) δ²H = 7.9 (±0.3) δ¹⁸O + 10.8 (±2.5) was derived, with slightly lower slope than Global Meteoric Water Line. The results help tracking two main sources of air masses affecting the study sites: rain events with the highest δ¹⁸O values (> ?5‰) originate over the Mediterranean Sea, whereas the more depleted samples (< ?8‰) are sourced in the North Atlantic region. The back trajectory analysis and deuterium excess values, ranging from 0.4 to 18.4‰, further support our findings. To assess the isotopic variation across the island, water samples from eight caves were collected. The δ¹⁸O values range between ?6.9 and ?1.6‰. With one exception (Artà), the isotopic composition of waters in caves located along the coast (Drac, Vallgornera, Cala Varques, Tancada, and Son Sant Martí) indicates Mediterranean‐sourced moisture masses. By contrast, the drip water δ¹⁸O values for inland caves (Campanet, ses Rates Pinyades) or developed under a thick (>50 m) limestone cap (Artà) exhibit more negative values. A well‐homogenized aquifer supplied by rainwaters of both origins is clearly indicated by groundwater δ¹⁸O values, which show to be within 2.4‰ of the unweighted arithmetic mean of ?7.4‰. Although limited, the isotopic data presented here constitute the baseline for future studies using speleothem δ¹⁸O records for western Mediterranean paleoclimate reconstructions. Copyright © 2016 John Wiley & Sons, Ltd.     

Identification of groundwater mean transit times of precipitation and riverbank infiltration by two‐component lumped parameter models

Groundwater transit time is an essential hydrologic metric for groundwater resources management. However, especially in tropical environments, studies on the transit time distribution (TTD) of groundwater infiltration and its corresponding mean transit time (mTT) have been extremely limited due to data sparsity. In this study, we primarily use stable isotopes to examine the TTDs and their mTTs of both vertical and horizontal infiltration at a riverbank infiltration area in the Vietnamese Mekong Delta (VMD), representative of the tropical climate in Asian monsoon regions. Precipitation, river water, groundwater, and local ponding surface water were sampled for 3 to 9 years and analysed for stable isotopes (δ¹⁸O and δ²H), providing a unique data set of stable isotope records for a tropical region. We quantified the contribution that the two sources contributed to the local shallow groundwater by a novel concept of two‐component lumped parameter models (LPMs) that are solved using δ¹⁸O records. The study illustrates that two‐component LPMs, in conjunction with hydrological and isotopic measurements, are able to identify subsurface flow conditions and water mixing at riverbank infiltration systems. However, the predictive skill and the reliability of the models decrease for locations farther from the river, where recharge by precipitation dominates, and a low‐permeable aquitard layer above the highly permeable aquifer is present. This specific setting impairs the identifiability of model parameters. For river infiltration, short mTTs (<40 weeks) were determined for sites closer to the river (<200 m), whereas for the precipitation infiltration, the mTTs were longer (>80 weeks) and independent of the distance to the river. The results not only enhance the understanding of the groundwater recharge dynamics in the VMD but also suggest that the highly complex mechanisms of surface–groundwater interaction can be conceptualized by exploiting two‐component LPMs in general. The model concept could thus be a powerful tool for better understanding both the hydrological functioning of mixing processes and the movement of different water components in riverbank infiltration systems.     

Hydrologic response to valley‐scale structure in alpine headwaters

Few systematic studies of valley‐scale geomorphic drivers of streamflow regimes in complex alpine headwaters have compared response between catchments. As a result, little guidance is available for regional‐scale hydrological research and monitoring efforts that include assessments of ecosystem function. Physical parameters such as slope, elevation range, drainage area and bedrock geology are often used to stratify differences in streamflow response between sampling sites within an ecoregion. However, these metrics do not take into account geomorphic controls on streamflow specific to glaciated mountain headwaters. The coarse‐grained nature of depositional features in alpine catchments suggests that these landforms have little water storage capacity because hillslope runoff moves rapidly just beneath the rock mantle before emerging in fluvial networks. However, recent studies show that a range of depositional features, including talus slopes, protalus ramparts and ‘rock‐ice’ features may have more storage capacity than previously thought. To better evaluate potential differences in streamflow response among basins with extensive coarse depositional features and those without, we examined the relationships between streamflow discharge, stable isotopes, water temperature and the amplitude of the diurnal signal at five basin outlets. We also quantified the percentages of colluvial channel length measured along the stepped longitudinal profile. Colluvial channels, characterized by the presence of surficial, coarse‐grained depositional features, presented sediment‐rich, transport‐limited morphologies that appeared to have a cumulative effect on the timing and volume of flow downstream. Measurements taken from colluvial channels flowing through depositional landforms showed median recession constants (K_r) of 0.9–0.95, δ¹⁸O values of ≥?14.5 and summer diurnal amplitudes ≤0.8 as compared with more typical surface water recession constant values of 0.7, δ¹⁸O ≤ ?13.5 and diurnal amplitudes >2.0. Our results demonstrated strong associations between the percentage of colluvial channel length within a catchment and moderated streamflow regimes, water temperatures, diurnal signals and depleted δ¹⁸O related to groundwater influx. Copyright © 2014 John Wiley & Sons, Ltd.