Determination of groundwater solute transport parameters in finite element modelling using tracer injection and withdrawal testing data

The groundwater tracer injection and withdrawal tests are often carried out for the determination of aquifer solute transport parameters. However, the parameter analyses encounter a great difficulty due to the radial flow nature and the variability of the temporal boundary conditions. An adaptive methodology for the determination of groundwater solute transport parameters using tracer injection and withdrawal test data had been developed and illustrated through an actual case. The methodology includes the treatment of the tracer boundary condition at the tracer injection well, the normalization of tracer concentration, the groundwater solute transport finite element modelling and the method of least squares to optimize the parameters. An application of this methodology was carried out in a field test in the South of Hanoi city. The tested aquifer is Pleistocene aquifer, which is a main aquifer and has been providing domestic water supply to the city since the French time. Effective porosity of 0.31, longitudinal dispersivity of 2.2 m, and hydrodynamic dispersion coefficients from D = 220 m²/d right outside the pumping well screen to D =15.8 m²/d right outside the tracer injection well screen have been obtained for the aquifer at the test site. The minimal sum of squares of the differences between the observed and model normalized tracer concentration is 0.00119, which is corresponding to the average absolute difference between observed and model normalized concentrations of 0.035 5 (while 1 is the worst and 0 is the best fit).     

Hydrochemical evidence for mixing of river water and groundwater during high-flow conditions, lower Suwannee River basin, Florida, USA

 Karstic aquifers are highly susceptible to rapid infiltration of river water, particularly during periods of high flow. Following a period of sustained rainfall in the Suwannee River basin, Florida, USA, the stage of the Suwannee River rose from 3.0 to 5.88 m above mean sea level in April 1996 and discharge peaked at 360 m³/s. During these high-flow conditions, water from the Suwannee River migrated directly into the karstic Upper Floridan aquifer, the main source of water supply for the area. Changes in the chemical composition of groundwater were quantified using naturally occurring geochemical tracers and mass-balance modeling techniques. Mixing of river water with groundwater was indicated by a decrease in the concentrations of calcium, silica, and ²²²Rn; and by an increase in dissolved organic carbon (DOC), tannic acid, and chloride, compared to low-flow conditions in water from a nearby monitoring well, Wingate Sink, and Little River Springs. The proportion (fraction) of river water in groundwater ranged from 0.13 to 0.65 at Wingate Sink and from 0.5 to 0.99 at well W-17258, based on binary mixing models using various tracers. The effectiveness of a natural tracer in quantifying mixing of river water and groundwater was related to differences in tracer concentration of the two end members and how conservatively the tracer reacted in the mixed water. Solutes with similar concentrations in the two end-member waters (Na, Mg, K, Cl, SO₄, SiO₂) were not as effective tracers for quantifying mixing of river water and groundwater as those with larger differences in end-member concentrations (Ca, tannic acid, DOC, ²²²Rn, HCO₃). Received, March 1999 / Revised, July 1999 / Accepted, July 1999     

Simulation of a reactive tracer experiment using stochastic hydraulic conductivity fields

Reactive tracer tests are performed to derive flow, transport and in situ biodegradation parameters. This paper describes the 3D simulation of a reactive tracer test using the transition probability geostatistical approach. Fifty different equally probable aquifer realizations were generated based on the geological information of 107 boreholes in an area of 62,500 m². One realization was chosen for the reactive transport simulation based on the results of groundwater flow modeling and on particle tracking calculations for the site. Field velocities at the site vary between 0.4 and 3 m/d. The transport of the reactive tracers deuterium ring labeled toluene-d₅ and fully deuterated toluene-d₈ was simulated and first-order biodegradation rates of 0.017 d⁻¹ for toluene-d₅ and 0.012 d⁻¹ for toluene-d₈ were determined.     

Groundwater age,mixing and flow rates in the vicinity of large open pit mines,Pilbara region,northwestern Australia

Determining groundwater ages from environmental tracer concentrations measured on samples obtained from open bores or long-screened intervals is fraught with difficulty because the sampled water represents a variety of ages. A multi-tracer technique (Cl, ¹⁴C, ³H, CFC-11, CFC-12, CFC-113 and SF₆) was used to decipher the groundwater ages sampled from long-screened production bores in a regional aquifer around an open pit mine in the Pilbara region of northwest Australia. The changes in tracer concentrations due to continuous dewatering over 7 years (2008–2014) were examined, and the tracer methods were compared. Tracer concentrations suggest that groundwater samples are a mixture of young and old water; the former is inferred to represent localised recharge from an adjacent creek, and the latter to be diffuse recharge. An increase in ¹⁴C activity with time in wells closest to the creek suggests that dewatering of the open pit to achieve dry mining conditions has resulted in change in flow direction, so that localised recharge from the creek now forms a larger proportion of the pumped groundwater. The recharge rate prior to development, calculated from a steady-state Cl mass balance, is 6 mm/y, and is consistent with calculations based on the ¹⁴C activity. Changes in CFC-12 concentrations with time may be related to the change in water-table position relative to the depth of the well screen.     

A conceptual framework of groundwater flow in some crystalline aquifers in Southeastern Ghana

A conceptual groundwater flow model was developed for the crystalline aquifers in southeastern part of the Eastern region, Ghana. The objective was to determine approximate levels of groundwater recharge, estimate aquifer hydraulic parameters, and then test various scenarios of groundwater extraction under the current conditions of recharge. A steady state groundwater flow model has been calibrated against measured water levels of 19 wells in the area. The resulting recharge is estimated to range from 8.97 × 10^?5 m/d to 7.14 × 10^?4 m/d resulting in a basin wide average recharge of about 9.6% of total annual precipitation, which results in a basin wide quantitative recharge of about 2.4 million m³/d in the area. This compares to recharge estimated from the chloride mass balance of 7.6% of precipitation determined in this study. The general groundwater flow in the area has also been determined to conform to the general northeast–southwest structural grain of the country. The implication is that the general hydrogeology is controlled by post genetic structural entities imposed on the rocks to create ingresses for sufficient groundwater storage and transport. Calibrated aquifer hydraulic conductivities range between 0.99 m/d and over 19.4 m/d. There is a significant contribution of groundwater discharge to stream flow in the study area. Increasing groundwater extraction will have an effect on stream flow. This study finds that the current groundwater extraction levels represent only 0.17% of the annual recharge from precipitation, and that groundwater can sustain future increased groundwater demands from population growth and industrialization.     

Groundwater flow and solute transport at the Mourquong saline-water disposal basin, Murray Basin, southeastern Australia

Saline groundwater and drainage effluent from irrigation are commonly stored in some 200 natural and artificial saline-water disposal basins throughout the Murray-Darling Basin of Australia. Their impact on underlying aquifers and the River Murray, one of Australia's major water supplies, is of serious concern. In one such scheme, saline groundwater is pumped into Lake Mourquong, a natural groundwater discharge complex. The disposal basin is hydrodynamically restricted by low-permeability lacustrine clays, but there are vulnerable areas in the southeast where the clay is apparently missing. The extent of vertical and lateral leakage of basin brines and the processes controlling their migration are examined using (1) analyses of chloride and stable isotopes of water (²H/¹H and ¹⁸O/¹⁶O) to infer mixing between regional groundwater and lake water, and (2) the variable-density groundwater flow and solute-transport code SUTRA. Hydrochemical results indicate that evaporated disposal water has moved at least 100 m in an easterly direction and that there is negligible movement of brines in a southerly direction towards the River Murray. The model is used to consider various management scenarios. Salt-load movement to the River Murray was highest in a "worst-case" scenario with irrigation employed between the basin and the River Murray. Present-day operating conditions lead to little, if any, direct movement of brine from the basin into the river. Electronic Publication     

Integrated use of geophysical, hydrological and geographic information system (GIS) methods in enhancing the groundwater quality in a fluoride-endemic terrain (Andhra Pradesh, India)

The concept of groundwater recharge and quality improvement is often implemented in arid and semi-arid areas with depleted aquifers. Nalgonda district in Andhra Pradesh, India, has endemic fluoride, with concentrations in drinking water varying between 3 and 8?mg/l. Numerous techniques adopted in the recent past for defluoridizing groundwater proved to have limitations. The integrated approach of a geographic information system (GIS) and an analytic hierarchy process (AHP), to identify suitable sites for recharge structures over an area of ??115?km², is highlighted. Further, to validate the delineated sites, a micro-watershed basin (2?km²) was selected for detailed recharge assessment and site feasibility studies through geophysical and tracer tests. Groundwater velocity (7?m/day) and flow direction through fractures in the shallow horizon were established through tracer experiments. The efficacy of the recommended recharge structures and their impact on groundwater quality were assessed over a period of 5?years, from 2002 to 2007, and the mean groundwater fluoride concentration of?>?3.5?mg/l over the study area was brought down to?<?1.5?mg/l.     

Use of geochemical tracers for estimating groundwater influxes to the Big Sioux River,eastern South Dakota,USA

Understanding the spatial distribution and variability of geochemical tracers is crucial for estimating groundwater influxes into a river and can contribute to better future water management strategies. Because of the much higher radon (²²²Rn) activities in groundwater compared to river water, ²²²Rn was used as the main tracer to estimate groundwater influxes to river discharge over a 323-km distance of the Big Sioux River, eastern South Dakota, USA; these influx estimates were compared to the estimates using Cl^? concentrations. In the reaches overall, groundwater influxes using the ²²²Rn activity approach ranged between 0.3 and 6.4 m³/m/day (mean 1.8 m³/m/day) and the cumulative groundwater influx estimated during the study period was 3,982–146,594 m³/day (mean 40,568 m³/day), accounting for 0.2–41.9% (mean 12.5%) of the total river flow rate. The mean groundwater influx derived using the ²²²Rn activity approach was lower than that calculated based on Cl^? concentration (35.6 m³/m/day) for most of the reaches. Based on the Cl^? approach, groundwater accounted for 37.3% of the total river flow rate. The difference between the method estimates may be associated with minimal differences between groundwater and river Cl^? concentrations. These assessments will provide a better understanding of estimates used for the allocation of water resources to sustain agricultural productivity in the basin. However, a more detailed sampling program is necessary for accurate influx estimation, and also to understand the influence of seasonal variation on groundwater influxes into the basin.     

Hydrogeology of a coal-seam gas exploration area, southeastern British Columbia, Canada: Part 2. Modeling potential hydrogeological impacts associated with depressurizing

A three-dimensional, finite-element flow model was used to assess the hydrogeological effects of depressurizing coalbeds lying in the Weary Creek exploration block, Elk River valley, southeastern British Columbia, Canada. The simulation results permit, at an early stage, assessment of the environmental and economic implications of how the flow system may respond to depressurization. Estimated reservoir conditions for the coal-seam gas targets lying within the Upper Jurassic–Lower Cretaceous Mist Mountain Formation indicate that the coalbeds must be depressurized by up to 350 m to attain the critical gas desorption pressure. The simulations suggest that depressurizing has little effect on groundwater flux to the Elk River. Simulated water production for three depressurizing wells operating under steady-state, single-phase flow for initial reservoir conditions of 13 and 16.5 cm³/g is 645 m³/d (4,057 barrels/d) and 355 m³/d (2,233 barrels/d), respectively. Groundwaters collected from monitoring wells have relatively low salinity, ranging from about 250–1,300 mg/L. The groundwater is supersaturated with respect to Ca–Mg–Fe carbonates (calcite, dolomite, and siderite) and Al-bearing silicates, including kaolinite and illite. Dissolved trace-metal concentrations are low; only Fe, Cd, Cr, and Zn exceed Canadian water-quality guidelines for aquatic life. Groundwaters were devoid of the more soluble monocyclic aromatic organic compounds, including benzene, toluene, ethylbenzene, and polycyclic aromatic compounds, including naphthalene. Electronic Publication     

Biodegradability of trimethylbenzene isomers under the strictly anaerobic conditions in groundwater contaminated by gasoline

Trimethylbenzene (TMB), as a constituent of gasoline, is often expected to be used as a conservative tracer in anaerobic BTEX-contaminated groundwater site to correct for attenuation due to dispersion, dilution and sorption along a flow path. To evaluate the suitability of using TMB as a tracer and to better understand biodegradability of TMB in contaminated groundwater by gasoline under anaerobic conditions, laboratory microcosms were conducted with mixed nitrate/iron/sulfate electron-acceptor amendments, using aquifer materials collected from Canadian Forces Base (CFB), Borden, Ontario, Canada. The results showed that under denitrifying conditions, biodegradation of 1,3,5-TMB, 1,2,4-TMB and 1,2,3-TMB were relatively slow and after 204 days of incubation approximately 27, 24, and 16% of the initial concentrations, respectively, were degraded in the microcosms. Under sulfate-reducing conditions, TMB isomers were recalcitrant. In contrast, significant biodegradation of TMB was observed under iron-reducing conditions. 1,3,5-TMB, 1,2,4-TMB and 1,2,3-TMB were degraded to 44, 47, and 24% of initial concentrations with first-order biodegradation rate constants of 0.003, 0.006 and 0.013 d⁻¹, respectively. This study indicates that TMB biodegradation is insignificant under nitrate and sulfate-reducing conditions but significant under iron-reducing conditions. Therefore, the use of TMB as a tracer for interpreting removal of other biodegradable gasoline constituents such as BTEX requires caution, especially in the presence of iron-reducing conditions.     

Aquifer characteristics and its modeling around an industrial complex,Tuticorin, Tamil Nadu,India: A case study

Anthropogenic pollution of shallow groundwater resources due to industrial activities is becoming a cause of concern in the east coastal belt of the state of Tamil Nadu, India. Integrated hydrogeological, geophysical and tracer studies were carried out in the coastal region encompassing an industrial complex. The objective has been to gain knowledge of aquifer characteristics, ascertaining groundwater movement and its flow direction, which would in turn reveal the possibility of contamination of groundwater regime and its better management. The results of multi-parameters and model study indicate that the velocity of groundwater flow ranges from 0.013 m/d to 0.22m/d in and around the industrial complex in upstream western part of the catchment and 0.026 m/d to 0.054m/d in the downstream eastern part, near the coast. These parameters are vital for the development of groundwater management scheme.     

A regional groundwater-flow model for sustainable groundwater-resource management in the south Asian megacity of Dhaka,Bangladesh

The water resources that supply most of the megacities in the world are under increased pressure because of land transformation, population growth, rapid urbanization, and climate-change impacts. Dhaka, in Bangladesh, is one of the largest of 22 growing megacities in the world, and it depends on mainly groundwater for all kinds of water needs. The regional groundwater-flow model MODFLOW-2005 was used to simulate the interaction between aquifers and rivers in steady-state and transient conditions during the period 1981–2013, to assess the impact of development and climate change on the regional groundwater resources. Detailed hydro-stratigraphic units are described according to 150 lithology logs, and a three-dimensional model of the upper 400 m of the Greater Dhaka area was constructed. The results explain how the total abstraction (2.9 million m³/d) in the Dhaka megacity, which has caused regional cones of depression, is balanced by recharge and induced river leakage. The simulated outcome shows the general trend of groundwater flow in the sedimentary Holocene aquifers under a variety of hydrogeological conditions, which will assist in the future development of a rational and sustainable management approach.

     

桂林东区峰林平原岩溶地下水示踪试验与分析

为查明桂林东区峰林平原地下水系统的结构特征,在岩溶水文地质调查的基础上,通过人工取样并采用高分辨率野外荧光仪,对桂林东区峰林平原区投放的示踪剂荧光素钠的接收情况进行检测。结果表明：峰林平原区地下水的出露比较多,主要以沟渠、湖塘、民井等为主,其中1号点初现时间为第22天,最大视流速为19.5 m/d,平均视流速4.9 m/d,说明本试验段内地下径流为典型的面状流态,岩溶含水介质比较均匀。峰林平原区岩溶含水层结构岩溶化程度很高,呈网络状,存在统一的地下水位,岩溶发育较均一,同时也存在较大的裂隙甚至较大规模的管道和地下河系统,容易形成岩溶地下水优先流。     

Changes of groundwater flow field of Luanhe River Delta under the human activities and its impact on the ecological environment in the past 30 years

The Luanhe River Delta is located in the center of the Circum-Bohai Sea Economic Zone. It enjoys rapid economic and social development while suffering relatively water scarcity. The overexploitation of groundwater in the Luanhe River Delta in recent years has caused the continuous drop of groundwater level and serious environmental and geological problems. This study systematically analyzes the evolution characteristics of the population, economy, and groundwater exploitation in the Luanhe River Delta and summarizes the change patterns of the groundwater flow regime in different aquifers in the Luanhe River Delta according to previous water resource assessment data as well as the latest groundwater survey results. Through comparison of major source/sink terms and groundwater resources, the study reveals the impacts of human activities on the groundwater resources and ecological environment in the study area over the past 30 years from 1990 to 2020. The results are as follows. The average annual drop rate of shallow groundwater and the deep groundwater in the centers of depression cones is 0.4 m and 1.64 m, respectively in the Luanhe River Delta in the past 30 years. The depression cones of shallow and deep groundwater in the study area cover an area of 545.32 km² and 548.79 km², respectively, accounting for more than 10% of the total area of the Luanhe River Delta. Overexploitation of groundwater has further aggravated land subsidence. As a result, two large-scale subsidence centers have formed, with a maximum subsidence rate of up to 120 mm/a. The drop of groundwater level has induced some ecological problems in the Luanhe River Delta area, such as the zero flow and water quality deterioration of rivers and continuous shrinkage of natural wetlands and water. Meanwhile, the proportion of natural wetland area to the total wetland area has been decreased from 99% to 8% and the water area from 1776 km² to 263 km². These results will provide data for groundwater overexploitation control, land subsidence prevention, and ecological restoration in plains and provide services for water resources management and national land space planning.© 2021 China Geology Editorial Office.     

Use of geothermal methods in outlining deep groundwater flow systems in Paleozoic interior basins of Brazil

Results of regional-scale geothermal studies are presented, providing new insights into the characteristics of deep groundwater flow systems in the Paleozoic sedimentary basins in the Amazon, Paraná and Parnaíba regions of Brazil. The study makes use mainly of bottom-hole temperature data sets for oil wells, the depths of which vary from 1,000 to 4,000 m. The techniques employed in data analysis have allowed identification of non-linear features in vertical distributions of temperature, produced by deep groundwater flows in the study area. According to the results obtained, vertical velocities of subsurface flows are found to fall in the range 10^?10 to 10^?9 m/s, while the horizontal velocities are significantly higher, of the order 10^?8 m/s. Identification of large-scale down flows has allowed inferences as to the existence of lateral movements of groundwater. The basins in the Amazon region are found to be characterized by widespread down flow of groundwater, implying the existence of distributed recharge systems operating on regional scales. There is a systematic decrease in horizontal velocities along the direction from west to east. This feature is considered indicative of gravity driven flows induced by episodes of uplift, since Miocene times, in the Andean region.     

Estimating groundwater recharge beneath irrigated farmland using environmental tracers fluoride, chloride and sulfate

Accurate recharge estimation is essential for effective groundwater management, especially in the North China Plain, where irrigation return flow is significant to vertical recharge but brings difficulty for recharge estimation. Three environmental tracers (F^?, Cl^? and SO₄ ^2?) were used to estimate vertical recharge based on the mass balance and cumulative methods. Four boreholes were dry-drilled to 5–25 m depth beneath irrigated farmland and one was drilled to 5 m beneath non-irrigated woodland; soil samples were collected in all boreholes at set depths. The results indicated that F^?, Cl^? and SO₄ ^2?were suitable tracers beneath the non-irrigated woodland, yielding recharge rates of 16.9, 18.8 and 19.4 mm/year, respectively. Recharge estimation was not straightforward when taking account of crop type, irrigation and/or fertilizer use. After comparing with previous research, conclusions were drawn: Cl^? was an appropriate tracer for irrigated farmland when taking account of Cl^? input from irrigation and absorption by crops; recharge rates were 65.9–126.8 mm/year. However, F^? was a more suitable tracer for irrigated regions where account is made of the proportion of precipitation to irrigation return flow, provided low F^? concentrations can be measured reliably.     

青海省格尔木冲洪积扇地下水资源评价及其开发利用价值

【研究目的】格尔木市地处柴达木盆地南缘，是青藏高原上正在崛起的新兴城市，随着各个工业园区的建设及中国最大的钾肥生产基地的建成，地下水资源的重要作用更加凸显，查明其地下水资源量对工矿企业及钾肥生产基地的运行至关重要。【研究方法】基于资料收集、动态长观、钻孔抽水试验等手段，采用补给量总和法以及排泄量总和法结合数值模拟等方法，评价了区域地下水天然资源量及允许开采量。【研究结果】格尔木河冲洪积扇地下水天然资源量为199.8×10⁴ m³/d，允许开采量为100×10⁴ m³/d。【结论】格尔木地区存在百万吨级的地下水源，6座水源（已建/在建/规划）规划开采量达百万吨时，地下水有充足的补给保证，可持续稳定开采。最终提交了B+C级地下水储量100×10⁴ m³/d的百万吨级地下水源。创新点：基于最新一次的地下水详查工作，完善了区域地下水流模型；在设定模型入境水量场景时，采用“丰枯年组合”的方法，充分考验了含水层的调节能力和连枯年持续供水能力，评价的结果可信度高。     

Tracing of residual multiple DNAPL sources in the subsurface using 222Rn as a natural tracer at an industrial complex in Wonju, Korea

Depth-discrete tracing of residual dense non-aqueous phase liquid (DNAPL) sources in the subsurface is of great importance in making decisions related to contaminated groundwater remediation. Temporal variations in the natural tracer ²²²Rn and contaminant concentrations in groundwater contaminated with multiple chlorinated contaminants, such as trichloroethene, carbon tetrachloride, and chloroform, were examined to trace residual multiple DNAPL contaminants at an industrial complex in Wonju, Korea. The ²²²Rn activities and multiple DNAPL concentrations in the groundwater fluctuated irregularly according to the groundwater recharge. The natural tracer ²²²Rn in groundwater present in the soil layer, originating from the underlying crystalline biotite granite, had a wide range from 29,000 to 179,000 Bq/m³, and total concentrations of chlorinated solvents ranged from 0.06 to 17.77 mg/l, indicating the ambiguous results of ²²²Rn for tracing the residual DNAPL sources. In this paper, a method is presented to locate zones with a high probability of containing depth-discrete residual multiple DNAPL sources using ²²²Rn and considering relative contaminant concentrations. The results demonstrate that the combination of the ²²²Rn activities as a natural tracer and the relative contaminant concentrations is able to be used as a useful tool for tracing residual DNAPLs.     

The study of hydrodynamic behaviour of a complex karst system under low-flow conditions using natural and artificial tracers (the catchment of the Unica River,SW Slovenia)

In this study a multi-tracer test with fluorescent tracers was combined with time series analyses of natural tracers to characterize the dynamics of the solute transport through different recharge pathways and to study hydraulic behaviour of a binary karst system under low-flow conditions. Fluorescent tracer testing included the introduction of uranine, amidorhodamine G, or naphthionate at three injection points. Sampling and monitoring took place at two karst springs (Malenščica, Unica) and at two underground rivers (Pivka, Rak) recharging the Unica River at the Polje of Planina, SW Slovenia. Other monitored parameters included precipitation, spring or underground river discharge, water temperature, and electrical conductivity. Water samples were collected and analyzed for total organic carbon, Mg²⁺, SO₄ ²⁻, and NO₃ ⁻ in the laboratory. In the study area, results of the tracer test suggest that contaminant transport in karst may be retarded for several weeks during low-flow conditions followed by increases in contaminant concentrations after subsequent rainfall events. Based on interpretation of tracer concentration breakthrough curves, low apparent dominant flow velocities (i.e., between 5.8 and 22.8 m/h through the well developed karst conduits, and 3.6 m/h through the prevailing vadose zone with a dominant influence of a diffuse recharge) were obtained. Together with analyses of hydro-chemographs the artificial tracing identified different origins of water recharging the studied aquifer. During prolonged low-water conditions the Malenščica spring is mainly recharged from the karst aquifer and the Unica spring by the sinking Pivka River. After more intensive rainfall events allogenic recharge from Cerknica prevails in the Malenščica spring, while the Unica spring drains mainly the allogenic water from the Pivka Valley. These findings of alternating hydraulic connections and drainage areas due to respective hydrological conditions are important and should be considered when monitoring water quality, implementing groundwater protection measures, and optimizing future water exploitation.     

Dynamics and anthropogenic impacts of multiple karst flow systems in a mountainous area of South China

The Xiangxi River basin, South China, is a steep terrane with well-developed karst features and an important Cambrian-Ordovician aquifer. Meteoric water in this mountainous area features a mean δ¹⁸O elevation gradient of –2.4?‰/km. This gradient was used to estimate mean recharge elevations of 760 m for Shuimoxi (SMX) spring, 1,060 m for Xiangshuidong (XSD) spring, and 1,430 m for drill hole ZK03, indicating multiple flow paths in the Cambrian-Ordovician karst aquifer. Mean residence times of 230 and 320 days and ～2 years were estimated for these features, respectively, using the damped running average model that predicts the isotopic variations in groundwater from those in precipitation. Groundwater in the regional karst flow system has the longest residence time, the highest recharge elevation, the longest flow paths, the lowest addition of anthropogenic components, and the greatest amount of water–rock interaction as indicated by its higher dissolved solids, Mg²⁺ concentrations and Mg/Ca ratios than the springs. In contrast, the local and shallow karst flow systems respond rapidly to recharge events. Artificial tracer tests prove that these shallow karst systems can also quickly transmit anthropogenic contaminants, indicating that they are highly vulnerable to human impacts, which include the enrichment of NO₃ ^–. The intensity of water–rock interaction and groundwater vulnerability are mainly determined by the structure and dynamics of the multiple karst flow systems.