Isotopic and geochemical characterization of salinization in the shallow aquifers of a reclaimed subsiding zone: The southern Venice Lagoon coastland

The coastal plain bordering the southern Venice Lagoon is a reclaimed lowland characterized by high subsidence rate, and ground level and water-table depth below sea level. In this agricultural region, where the surface hydrologic network is entirely artificially controlled by irrigation/drainage canals, salinization problems have long been encountered in soils and groundwaters. Here we use isotopic and geochemical tracers to improve our understanding of the origin of salinization and mineralization of the semi-confined aquifer (0–40 m), and the freshwater inputs to this hydrological system. Water samples have been collected at different seasons in the coastal Adriatic Sea, lagoon, rivers and irrigation canals, as well as in the semi-confined aquifer at depths between 12 and 35 m (14 boreholes), and in the first confined aquifer (three boreholes drilled between 40 and 80 m depth). Stable isotopes (δ¹⁸O and δD) and conductivity profiles show that direct saline intrusion from the sea or the lagoon is observed only in a restricted coastal strip, while brackish groundwaters are found over the entire topographic and piezometric depression in the centre of the study area. Fresh groundwaters are found only in the most western zone. The sharp isotopic contrast between the western and central regions suggests disconnected hydrological circulations between these two parts of the shallow aquifer. The border between these two regions also corresponds to the limits of the most strongly subsiding zone.Our results can be interpreted in terms of a four end-member mixing scheme, involving (1) marine water from the lagoon or the open sea, (2) alpine and pre-alpine regional recharge waters carried either by the main rivers Adige, Bacchiglione and Brenta (irrigation waters) or by the regional groundwater circulation, (3) local precipitation, and (4) evaporated waters infiltrated from the surface. Infiltration from the surface is also revealed by the stratification of the electrical conductivity profiles, showing that the brackish groundwaters are overlain by a shallow layer of less saline water all over the central depression. In the first confined aquifer, the groundwaters have isotopic compositions similar to the deep groundwaters of the Venetian confined aquifers (40–400 m depth). The isotopic data and the Br/Cl ratio show that the origin of the salinization of the phreatic aquifer can be ascribed to seawater intrusion alone, with no indication of the involvement of deep brines (identified at 450 m depth) in the process.The chemical composition of the saline and brackish groundwaters is characterized by an excess of sodium and a deficit of calcium compared to conservative mixing between fresh groundwaters and seawater. This suggests that the phreatic aquifer is progressively freshening, as a consequence of the beneficial influence of the extensive irrigation/drainage network, including raised canals acting as a hydraulic barrier along the coast. This freshening tendency may have been lasting since the reclamation in the mid-twentieth century, and has probably been accelerated by the ban on groundwater abstraction since the 1970s.     

Deciphering groundwater flow between the Complex Terminal and Plio-Quaternary aquifers in Chott Gharsa plain (southwestern Tunisia) using isotopic and chemical tools

Abstract

The Complex Terminal (CT) and Plio-Quaternary (P-Q) aquifers in the Chott Gharsa plain in southwestern Tunisia have been investigated with the aid of chemical and isotopic tools. It has been demonstrated that groundwater from the CT is mainly of palaeo-origin, especially in the western and central parts of the plain where the most negative values of δ¹⁸O and δ²H were observed (between??8.1 and??7.6‰ for δ¹⁸O, and??60 to??57‰ for δ²H), combined with low concentrations of radiocarbon (6.8–7.5 pmc) and absence of tritium. Modern recharge of the aquifer occurs only in the eastern part of the system where younger waters were observed, as indicated by their stable isotope composition, relatively high radiocarbon content and presence of tritium. Groundwater from the P-Q multi-layer aquifer represents mixtures of ascending deep CT waters and modern water recharging the P-Q aquifer system. Isotope mass balance was used to quantify mixing proportions. The calculations showed that the contribution of deep CT groundwater to the P-Q aquifer system reaches about 75% in the western and central parts of the plain where the CT aquifer remains strongly artesian. This contribution decreases to about 15% towards the eastern part of the plain, as a consequence of significant reduction of artesian pressure in this area of the CT aquifer. Chemical data suggest that mineralization of the studied groundwater systems is controlled mainly by dissolution of evaporative minerals (halite, anhydrite and gypsum) and cation exchange reactions with the matrix, possibly enhanced by recent anthropogenic disturbance of the system caused by lowering of the water table due to heavy exploitation and return flow of saline irrigation water into the P-Q aquifer.

Editor D. Koutsoyiannis; Associate editor E. Custodio

Citation Yangui, H., Abidi, I., Zouari, K., and Rozanski, K., 2012. Deciphering groundwater flow between the Complex Terminal and Plio-Quaternary aquifers in Chott Gharsa plain (southwestern Tunisia) using isotopic and chemical tools. Hydrological Sciences Journal, 57 (5), 967–984.     

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.     

Recharge sources and hydrogeological effects of irrigation and an influent river identified by stable isotopes in the Motril‐Salobreña aquifer (Southern Spain)

Measurement of the stable isotopes oxygen‐18 and deuterium in water is an important tool to characterize aquifer recharge sources. In the driest areas of the Mediterranean, this application is of special interest due to the scarcity of water and the resulting common incidence of human influence on natural hydrological systems. The Motril‐Salobreña detrital aquifer (southern Spain) is a clear example of such an impact as inhabitants have designed irrigation systems and a dam was recently built across the course of the Guadalfeo River, which feeds the aquifer. The sampling of (river or ground) water has allowed the determination of stable isotope contents (oxygen‐18 and deuterium), both temporally and spatially, and the relative importance of the main recharge sources in certain sectors. In addition, we were able to infer seasonal trends and to improve existing knowledge of the main flow paths and the position of a seasonal groundwater divide. Data analysis shows evaporation plays a minor role (despite the high temperatures in the zone), scarce rainwater influence, and the overwhelming contribution of recharge from the Guadalfeo River and from the carbonate aquifer (Escalate aquifer) in contact with the Motril‐Salobreña aquifer. Irrigation return flow during the summer months comprises the main recharge due to the significant volumes of water that infiltrate. The construction of the dam will almost certainly entail great changes in the current dynamics of the hydrogeology of the Motril‐Salobreña aquifer; therefore, knowledge of its behaviour is crucial in order to carry out sustainable use of its groundwater resources. Copyright © 2011 John Wiley & Sons, Ltd.     

Origin and type of rainfall for recharge of a karstic aquifer in the western Mediterranean: a case study from the Sierra de Gador–Campo de Dalias (southeast Spain)

Isotope signatures in precipitation from the Global Network for Isotopes in Precipitation around the Mediterranean basin and literature data are compared with isotopic data from a large karstic aquifer in southeast Spain to explain the origin and type of the precipitation events dominating recharge. Analysis of the deuterium excess d at the scale of the Mediterranean basin and at the regional scale allows us to understand the isotopic context of the study area: Campo de Dalias and the Sierra de Gador (Almería province). The origin of precipitation can be determined from its d value. The d value changes as a function of the initial evaporation condition. It depends on the relative humidity and temperature during the evaporation producing the water vapour of the clouds. The water vapour, which dominates the study area, is generated in two areas: the Atlantic Ocean (d = 10‰) and the western Mediterranean basin (d = 15‰). With increasing precipitation volume, the western Mediterranean character dominates. These heavier storms contribute mainly to recharge, as illustrated by the d value of 13·6‰ in deep groundwater of the Campo de Dalias. Weighted d values increase with the volume of precipitation, giving a significant relationship for the southern and eastern coasts of the Iberian Peninsula. This selectivity of d to monthly precipitation was used to estimate the return period of precipitation leading to aquifer recharge at 0·9–4·9 years. Moderate rainfall, which occurs more frequently, still represents ～60–90% of the total precipitation. One of the challenges to meet ever‐growing water demands is to increase recharge from moderate events yielding intermediate quantities per event, but forming the bulk of the annual precipitation. Copyright © 2006 John Wiley & Sons, Ltd.     

Factors controlling stable oxygen,hydrogen and carbon isotope ratios in regional groundwater of the eastern United Arab Emirates (UAE)

Increasing groundwater salinity and depletion of the aquifers are major concerns in the UAE. Isotopes of oxygen, hydrogen, and carbon concentrations in groundwater were used to estimate evaporation loss using the isotopes of oxygen and hydrogen, and using a carbon isotope to trace inorganic carbon cycling in two main aquifers in the eastern part of the United Arab Emirates. The δD‐δ¹⁸O of groundwater samples plotted on a line given by: δD = 4 δ¹⁸O + 4 ·4 (r² = 0·4). In comparison, the local meteoric water line (LMWL) has been defined by the line: δD = 8 δ¹⁸O + 15. In order to better understand the system investigated, samples were separated into two groups based on the δD‐δ¹⁸O relationship. These are (1) samples that plot above the LMWL (δD = 6·1 δ¹⁸O + 12·4, r² = 0·8) and which are located predominantly in the north of the study area, and (2) samples that plot below the LMWL (δD = 5·6 δ¹⁸O + 6·2, r² = 0·8) and which are mostly distributed in the south. Slopes for both the groups are similar and lower than that for LMWL indicating potential evaporation of recharging water. However, the y‐intercept, which differs between the two groups, suggests evaporation of return flow and evapotranspiration in the unsaturated zone to be more significant in the south. This is attributed to intense agricultural activities in the region. Samples from the eastern Gravel Plain aquifer have δ¹³C and dissolved inorganic carbon (DIC) values in the range from ? 10 to 17‰, and 12–100 mg C/l, respectively, while the range for those from the Ophiolite aquifer is from ? 11 to ? 16.4‰, and 16–114 mg C/l respectively. This suggests the control of C‐3 and C‐4 plants on DIC formation, an observation supported by the range δ¹³C of soil organic matter (from ? 18·5 to ? 22·1‰.) Copyright © 2007 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.     

Evaluation of historical nitrate sources in groundwater and impact of current irrigation practices on groundwater quality

Abstract

Large-scale agricultural activities cause deterioration of groundwater resources throughout the world. This study focuses on the irrigated agricultural district of Jinghuiqu in Shaanxi, China. The objectives are to identify the main sources and processes that control nitrate transport, and to assess the impact of long-term irrigation practices on nitrate leaching. The hydrogeochemistry of major solute constituents and stable isotope ratios of NO₃ ^? in groundwater were used to identify historical sources of contamination and transformation processes occurring in the shallow groundwater of the Jinghuiqu irrigation district. The N-NO₃ ^? concentrations in groundwater ranged from 0.0 to more than 100 mg L^?1 in April 1990, and from 0.47 to 42.0 mg L^?1 in December 2009. Our measurements of N-NO₃ ^? show that the main reasons for this dramatic decline in N-NO₃ ^? concentrations from 1990 to 2009 are denitrification, which controls nitrogen types, together with the sharp decline of the groundwater table. The isotope ratios in collected samples showed that the source of nitrate was mainly manure, a result that corresponds with historical data showing that manure was the main nitrate source before the 1980s, and that fertilizers became dominant thereafter. A numerical model was then used to evaluate future impacts of current irrigation practices on groundwater sources. The HYDRUS-1D model was employed first to evaluate the water movement and the N-NO₃ ^? solute transport for a single irrigation pulse experiment carried out under field conditions, and then in evaluating the long-term impact of irrigation on N-NO₃ ^? leaching. Simulations showed that, after introducing irrigation, the downward drainage fluxes of N-NO₃ ^?, and the N-NO₃ ^? and N-tot concentrations increased at all depths within the 2-m soil profile. This indicates that N concentrations below 1-m depth and in the groundwater could become even more elevated with increased future irrigation.

Editor D. Koutsoyiannis

Citation Liu, X.-H., Sun, S.-J., Ji, P., and ?im?nek, J., 2013. Evaluation of historical nitrate sources in groundwater and impact of current irrigation practices on groundwater quality. Hydrological Sciences Journal, 58 (1), 1–15.     

Isotopic tracing for conceptual models of groundwater hydrodynamics in multilayer aquifer systems of central and southern Tunisia

Abstract

New and previously published data sets including stable and radiogenic isotope measurements (¹⁸O, ²H, ³H, ¹³C and ¹⁴C) were used to investigate, conceptualize and compare groundwater hydrodynamics within three major multilayer aquifer systems located in central and southern Tunisia. It has been demonstrated that the investigated aquifer systems contain modern and palaeoclimatic waters. Modern groundwaters, which refer to contemporaneous and post-nuclear recharge waters, are characterized by enriched stable isotope contents, high carbon-14 activities and high to moderate tritium concentration. While palaeoclimatic groundwaters, which refer to Late Pleistocene and Early Holocene recharge waters, are distinguished by their depleted stable isotope contents, low carbon-14 activities and insignificant tritium concentrations. Established conceptual models have elucidated the groundwater hydrodynamics within the studied aquifer systems. They show that groundwater mixing occurs between end-members from the shallow and deep aquifers that migrate by downward and upward leakage towards the intermediate aquifer.

Editor D. Koutsoyiannis; Associate editor S. Faye

Citation Dassi, L. and Tarki, M., 2014. Isotopic tracing for conceptual models of groundwater hydrodynamics in multilayer aquifer systems of central and southern Tunisia. Hydrological Sciences Journal, 59 (6), 1240–1258. http://dx.doi.org/10.1080/02626667.2014.892206     

Hydrogeochemistry and environmental isotopes of ground water in Jeju volcanic island,Korea: implications for nitrate contamination

Ground water from springs and public supply wells was investigated for hydrochemistry and environmental isotopes of ³H, ¹⁸O and D in Jeju volcanic island, Korea. The wells are completed in a basaltic aquifer and the upper part of hydrovolcanic sedimentary formation. Nitrate contamination is conspicuous in the coastal area where most of the samples have nitrate concentrations well above 1 mg NO₃ N/l. Agricultural land use seems to have a strong influence on the distribution of nitrate in ground water. Comparison of stable isotopic compositions of precipitation and ground water show that ground water mostly originates from rainy season precipitation without significant secondary modification and that local recharge is dominant. ³H concentration of ground water ranged from nearly zero to 5 TU and is poorly correlated with vertical location of well screens. The occurrence of the ³H‐free, old ground water is due to the presence of low permeability layers near the boundary of the basaltic aquifer and the hydrovolcanic sedimentary formation, which significantly limits ground water flow from the upper basaltic aquifer. The old ground water exhibited background‐level nitrate concentrations despite high nitrate loadings, whereas young ground water had considerably higher nitrate concentrations. This correlation of ³H and nitrate concentration may be ascribed to the history of fertilizer use that has increased dramatically since the early 1960s in the island. This suggests that ³H can be used as a qualitative indicator for aquifer vulnerability to nitrate contamination. Copyright © 2005 John Wiley & Sons, Ltd.     

A comparison of peak and plate electrodes in electrical resistivity tomography: application to the chalky groundwater of the Beauvais aquifer (northern part of the Paris basin,France)

Electrical resistivity tomography was used in order to explore an experimental site of the LaSalle Beauvais Polytechnic Institute (France). The test was conducted along a profile line of 315 m length, using 64 electrodes deployed at an inter‐electrode spacing of 5 m, and the data were recorded using gradient, Wenner and pole–dipole arrays. The performance of plate electrodes (non‐conventional flat‐based) is compared with the performance of peak electrodes (conventional spike). The hydrogeophysical investigation of the chalk aquifer system of Beauvais shows that the performance of plate electrodes is satisfactory, leading to inversions of small root‐mean‐square errors. Peak and the plate electrodes were tested before and after injection of a salt tracer in the piezometer of the experimental site. The study demonstrates the usefulness of plate electrodes (efficient, less time consuming) and the possibility of aquifer characterization by a salt tracer. Copyright © 2010 John Wiley & Sons, Ltd.     

Characterizing spatial and temporal variation in 18O and 2H content of New Zealand river water for better understanding of hydrologic processes

Time series of hydrogen and oxygen stable isotope ratios (δ²H and δ¹⁸O) in rivers can be used to quantify groundwater contributions to streamflow, and timescales of catchment storage. However, these isotope hydrology techniques rely on distinct spatial or temporal patterns of δ²H and δ¹⁸O within the hydrologic cycle. In New Zealand, lack of understanding of spatial and temporal patterns of δ²H and δ¹⁸O of river water hinders development of regional and national-scale hydrological models. We measured δ²H and δ¹⁸O monthly, together with river flow rates at 58 locations across New Zealand over a two-year period. Results show: (a) general patterns of decreasing δ²H and δ¹⁸O with increasing latitude were altered by New Zealand's major mountain ranges; δ²H and δ¹⁸O were distinctly lower in rivers fed from higher elevation catchments, and in eastern rain-shadow areas of both islands; (b) river water δ²H and δ¹⁸O values were partly controlled by local catchment characteristics (catchment slope, PET, catchment elevation, and upstream lake area) that influence evaporation processes; (c) regional differences in evaporation caused the slope of the river water line (i.e., the relationship between δ²H and δ¹⁸O in river water) for the (warmer) North Island to be lower than that of the (cooler, mountain-dominated) South Island; (d) δ²H seasonal offsets (i.e., the difference between seasonal peak and mean values) for individual sites ranged from 0.50‰ to 5.07‰. Peak values of δ¹⁸O and δ²H were in late summer, but values peaked 1 month later at the South Island sites, likely due to greater snow-melt contributions to streamflow. Strong spatial differences in river water δ²H and δ¹⁸O caused by orographic rainfall effects and evaporation may inform studies of water mixing across landscapes. Generally distinct seasonal isotope cycles, despite the large catchment sizes of rivers studied, are encouraging for transit time analysis applications.     

Temporal variation of stable isotopes in a precipitation–groundwater system: implications for determining the mechanism of groundwater recharge in high mountain–hills of the Loess Plateau,China

The groundwater in shallow loess aquifers in high mountain–hills in the western Loess Plateau in China is almost the sole water resource for local residents. However, the question about how the loess groundwater naturally circulates in these high mountain–hills, characterized by low precipitation and high potential evaporation, remains unclear. The objectives of this study are to evaluate the application of hydrogen and oxygen isotopes to (1) examine temporal variations of the isotopic composition of precipitation and shallow groundwater and (2) uncover the mechanism of groundwater recharge in high mountain–hills. Results from 2 years of monitoring data show a difference in the stable isotopes for groundwater and local precipitation between the winter and summer periods. Similar to precipitation, stable isotopes in groundwater are observed to be depleted in winter and enriched in summer, particularly in oxygen isotope. A prominent characteristic is that H and O isotopes of groundwater show a very clear response to strong precipitation in the rainy season in 2013. The results highlight that local precipitation is the likely recharge source for groundwater in shallow loess aquifers. Annual recharge from local precipitation maintains the groundwater resource in the shallower loess aquifer. The mechanisms governing shallow loess groundwater recharge in high mountain–hills were evaluated. In addition to possible vertical slow percolation of soil water through the unsaturated zone, rapid groundwater recharge mechanisms have been identified as temporal preferential infiltration through sinkholes, slip surface or landslide surface and through the interface of loess layer and palaeo‐soils. Most groundwater can be recharged after a heavy rainy season. Copyright © 2015 John Wiley & Sons, Ltd.     

云南4个湖泊浮游生物碳、氮稳定同位素的季节变化及其影响因子

浮游生物是湖泊食物网的重要组成,其碳、氮稳定同位素能够反映元素地球化学循环和食物来源的波动,是了解水域生态系统结构变化的重要手段之一.本文选取云南4个不同类型湖泊,开展浮游生物碳、氮稳定同位素组成(δ13C、δ15N)的季节变化与湖泊对比研究.大型深水湖泊(抚仙湖和阳宗海)中,浮游植物δ13C值在夏、秋季(-20.34...     

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.