Artificial Tracer Tests Research Articles

Multi-scale hydrostructural approach for karst environment. Application to the Arcier hydrosystem (eastern France)

Based on a multi-scale and hydrostructural approach, this study presents the most relevant methodology to be applied to a karst hydrosystem in order to get a full understanding of underground water flow. It implies a complete structural analysis, from the hydrosystem scale to the outcrop scale, including the intermediate scale of the major geological structures. We illustrate the method in the Arcier hydrosystem, in the northwestern border of the Jura fold-and-thrust belt (Eastern France).Field mapping and structural analysis allow to update the geological vision of the hydrosystem with two kink-type fault propagation folds, including a trishear kinematic model, on either side of a plateau presenting a hollow-and-dome configuration. Fracturing analysis reveals a fault-fracture network that we infer governs the entire hydrosystem. A Riedel pattern is highlighted, characterized by a N–S-striking (N355° ± 5), sinistral strike-slip, regional shear zone. Then, two 3D geological models, at different scales, constructed with MOVE and Visual Karsys softwares are combined with water levels and artificial tracer tests. It reveals a multilayer aquifer and a redefinition of groundwater circulations for the Arcier hydrosystem.The results demonstrate a strong geological control of karstic hydrosystems on groundwater circulations, proving that classical hydrogeological methods, such as natural and/or artificial tracers, must be combined with rigorous geological analysis. Moreover, the multi-scale approach provides an explanation of groundwater circulation based on the intersection between 3D geometry of impervious layers delimiting the aquifers and their base water level, instead of the 2D view (section or map) requiring systematic recourse to inferred vertical faults to cross permeability barriers vertically or laterally. This study also brings a new vision to the local protection of the water resource.

Recharge dynamic and flow-path geometry controls of solute transport in karst aquifer

Solute transport is complex in karst aquifer due to its heterogeneous structure and polytropic hydrodynamic condition. Hydrological monitoring and artificial tracer test were conducted in two flow paths under different recharge conditions, and the Advection–Dispersion Equation and Dual Region Advection–Dispersion models were successfully calibrated to simulate the breakthrough curves. Intermittent concentrated recharge at sinkholes and quick hydrological response at spring outlet often occur after rainstorms. Single symmetrical shape and bimodal shape with long tail were characterized in the breakthrough curves. Stronger recharge hydrodynamic condition drives a larger transport velocity with higher solute transport efficiency. The tracer injection time during a concentrated recharge event is found to significantly affect the solute transport process. Constantly pushed by the follow-up recharge water, the piston transport process of solute in conduit enables higher recovery rate for the tracer injected at the initial recharge process. Solute transient storage in fissures occurs while tracer is injected at the peak flow. A negative correlation between dispersivity and flow velocity is found in karst conduits. The reduction in cross-section area of karst conduits causes high dispersivity and obvious tailing in breakthrough curves due to the blocking function of solute by the uneven bottom surface in karst conduits. The solute transport process is significantly controlled by the changes in hydrodynamic and flow cross-section due to concentrated recharge events, which deepen our understanding of solute transport in karst groundwater.

Assessment of urban landslide groundwater characteristics and origin using artificial tracers, hydro-chemical and stable isotope approaches

In the framework of landslides, the hydrogeological features play an essential role in slope stability, governing water movement and thus resulting in modification of the effective stress in the soil. In this framework, the hydrogeological conceptualization of landslide areas and the identification of groundwater origin are key points to developing risk mitigation measures. In fact, groundwater recharge cannot always be attributed to local precipitation alone. Mixing processes between water derived from local infiltration and deep water upflow along tectonic lineaments or anthropogenic water can affect the groundwater balance on a local scale. This study aims to define the potential groundwater origin of one of the highest risk urban landslides in central Italy and to define a hydrogeological conceptual model by exploiting its existing drainage system network. This research is based on a multiple-techniques approach based on hydrological water balance, artificial tracer tests during recharge period, seasonal monitoring of the water stable-isotope content, hydro-chemical survey during low-flow periods, and analysis of the piezometric level fluctuation. All these analyses are coupled with a detailed reconstruction of the geology of the area depicted from boreholes and drill holes. Two groundwater bodies have been evidenced from the study. The shallower one is located in the landslide unstable zone and is hydraulically connected to a deeper groundwater body hosted in the underlying bedrock. Results highlighted that the local rainfall regime could not fully explain the hydro-chemical facies. Local water contributions to the landslide area coming from leakage of the urban sewerage system have been evidenced, excluding deep groundwater upflow from the fault system.

3D Discrete Fracture Network Modelling from UAV Imagery Coupled with Tracer Tests to Assess Fracture Conductivity in an Unstable Rock Slope: Implications for Rockfall Phenomena

The stability of a rock slope is strongly influenced by the pattern of groundwater flow through the fracture system, which may lead to an increase in the water pressure in partly open joints and the consequent decrease in the rock wall strength. The comprehension of the fracture pattern is a challenging but vital aspect in engineering geology since the fractures’ spatial distribution, connectivity, and aperture guide both the water movement and flow quantity within the rock volume. In the literature, the most accepted methods to hydraulically characterise fractured rocks in situ are the single borehole packer test, the high-resolution flow meters for fractures, and the artificial tracer tests performed in boreholes. However, due to the high cost a borehole requires and the general absence of wells along coastal cliffs, these methods may not be appropriate in rockfall-prone areas. In this study, an unsaturated rocky cliff, strongly affected by rockfalls, was investigated by combining kinematic analysis, Discrete Fracture Network (DFN) modelling, and artificial tracer tests. The DFN model and potential rock block failure mechanisms were derived from high-resolution 3D virtual outcrop models via the Structure from Motion (SfM) photogrammetry technique. An artificial tracer was injected using a double ring infiltrometer atop the recharge zone of the slope to determine the infiltration rate and validate the DFN results. The DFN and tracer test methods are frequently used at different spatial scales and for different disciplines. However, the integration of digital photogrammetry, DFN, and tracer tests may represent a new step in rockfall and landslide studies. This approach made possible the identification of groundwater flow patterns within the fracture system and revealed about a 10-day tracer transit time from the injection area and the monitored slope, with similar conductivity values gathered from both the DFN and tracer test. Planar and wedge failures with volumes ranging from 0.1 and 1 m3 are the most probable failure mechanisms in the areas. The results were consistent with the delay between the intense rainfall and the slope failures previously documented in the study area and with their mechanisms.

Reaction of the carbonate Sibillini Mountains Basal aquifer (Central Italy) to the extensional 2016–2017 seismic sequence

Hydrogeological perturbations in response to earthquakes are widely described worldwide. In carbonate aquifers, a post-seismic discharge increase is often attributed to an increase of bulk permeability due to co-seismic fracturing and the attention on the role of faults to explain the diversion of groundwater is increasing. We focus on the reaction of carbonate hydrogeological basins to extensional seismicity, taking as an example the effects of the Central Italy 2016–2017 seismic sequence, on the Basal aquifer of the Sibillini Mountains area. Geo-structural, seismological and ground deformation data were collected and merged with artificial tracer tests results and with a 4-years discharge and geochemical monitoring campaign. The main NNW-directed groundwater flow was diverted to the west and a discharge deficit was observed at the foot-wall of the activated fault system with a relevant discharge increase, accompanied by geochemical variations, at the fault system hanging-wall. The observed variations are consistent with the combined action of a permeability increase along the activated fault systems, which modified the predominant pre-seismic along-strike regional flow, and with hydraulic conductivity increase due to fracturing, determining a fast aquifers emptying. We show that the prevailing mechanism depends on the aquifer systems position with respect to the activated faults.

Transit Time index (TTi) as an adaptation of the humification index to illustrate transit time differences in karst hydrosystems: application to the karst springs of the Fontaine de Vaucluse system (southeastern France)

Abstract. Transit time can be estimated thanks to natural tracers, but few of them are usable in the 0–6-month range. The main purpose of this work is to analyze the potential of the ratio of heavy- to light-weight organic compounds (the humification index (HIX); Ohno, 2002; Zsolnay et al., 1999) as a natural tracer of short transit time (Blondel et al., 2012). Critical analysis of former studies shows that although the link between HIX and transit time seems consistent, the whole methodological approach needs to be consolidated. Natural organic matter fluorescence from 289 groundwater samples from four springs and 10 flow points located in the unsaturated zone of the Vaucluse karst system is characterized by parallel factor analysis (PARAFAC) thanks to the excitation–emission matrix (EEM), thus (i) allowing for the identification of main fluorescent compounds of sampled groundwater and (ii) evidencing the inadequacy of HIX 2D emission windows to characterize groundwater organic matter. We then propose a new humification index called the Transit Time index (TTi) based on the Ohno (2002) formula but using PARAFAC components of heavy and light organic matter from our samples instead of 2D windows. Finally, we evaluate TTi relevance as a transit time tracer by (i) performing a detailed analysis of its dynamics on a selected spring (Millet) and (ii) comparing its mean value over karst springs of the Vaucluse karst system. Principal component analysis (PCA) of TTi and other hydrochemical parameters monitored at Millet spring put in relief the different ranges of transit time associated with the different organic matter compounds. PCA results also provide evidence that TTi can detect a small proportion of fast infiltration water within a mix, while other natural tracers of transit time provide no or less sensitive information. TTi distributions at monitored karst springs are consistent with relative transit times expected for the small-scale, short average transit time systems. TTi thus appears as a relevant qualitative tracer of transit time in the 0–6-month range where existing tracers fail and may remain applicable, even in the case of anthropic contamination thanks to PARAFAC modeling. Transforming it into quantitative information is a challenging task which may be possible thanks to intensive studies of organic matter degradation kinetics in natural waters with the help of radiogenic isotope usage or an artificial tracer test.

Tracing flow directions of the Pentung Allogenic River in Gunungsewu Karst Area, Gunungkidul, using artificial tracer test

Previous studies in the Gunungsewu Karst Area categorized the Pentung Allogenic River catchment that drains to the Sawah Ombo Pond as the upper part of the Bribin Karst Drainage Basin (KDB). However, more detailed results from natural isotope studies after 2010 have provided reasons to postulate that the groundwater does not flow from the catchment to Bribin KDB. Therefore, the research aimed to trace the direction of Pentung’s groundwater flows in Gunungsewu’s karstic hydrogeological system. To achieve this, an artificial tracer test was used. There were two fluorescent dye injection namely uranine and tinopal. The results confirmed that the allogenic river’s catchment supplied groundwater to Beton KDB instead of Bribin KDB. While the artificial tracers were not detected in the water at Gilap Cave and Bribin Underground River, they appeared in Jomblangan Cave and Sriti Spring, parts of Beton KDB.

Groundwater Resources in a Complex Karst Environment Involved by Wind Power Farm Construction

The need to produce energy from clean energy sources has caused public administrations and private companies to look for suitable places. The windiness detected in the eastern area of the Matese karst massif (southern Italy) has favored the construction of wind farms to produce electricity from clean energy sources. During the installation of the first wind turbines, some alterations in the supply of drinking water, fed by the springs of this area, were attributed by the population to this installation. Therefore, in order to assess whether there has been an impact produced by the wind farms on the quality of groundwater, a detailed hydrogeological study was developed. Karst hydrogeological features of the area were mapped, focusing on endorheic areas, sinkholes and karst springs. Artificial tracer tests were then carried out to investigate groundwater flow circulation and connection between surface karst landforms and springs. Chemical and physical characteristics of the groundwater were monitored during the construction of the wind farms and, for the following months, by infield measurements and laboratory analysis of spring water samples. This study highlights that wind farms mainly develop along the boundary of endorheic areas, which are important recharge zones for groundwater resources, and are directly connected to the major karst springs through sinkholes and a dense network of karst conduits. The results of the monitoring did not reveal any anomalies in the quality of the water and, therefore, any alterations cannot be attributed to the wind farms. Our investigation appears useful for a better understanding of the possible actual and future effects of the wind farms on both groundwater circulation and spring water quality in this karst area.

Characterizing the hierarchical groundwater flow systems in Karstic Xiangxi River Basin, West Hubei, Central China

Karstification degree is a key factor in controlling the characteristics of karst GFP. In this study, hydrogeology, karst development, hydrograph, artificial and natural tracers and hydrochemistry were integrated to characterize the groundwater flow patterns (GFP) in the Cambrian-Ordovician karst aquifer in Xiangxi River Basin. Groundwater samples were collected from three karst groundwater points: Wulongdong spring (WLD), Shuimoxi spring (SMX), and groundwater exposed by Borehole ZK03 (ZK03), which have different discharge variations. The hydrogeological analysis and the artificial tracer tests showed that groundwater runoff channel for WLD, SMX and ZK03 are dominated by conduits, fracture networks and micro fissures-matrix, respectively. The natural responses (discharge, electrical conductivity and temperature) of WLD to rainfall indicated that flow in this system is discharged by conduit, fracture and matrix, respectively. Stable isotopes (2H, 18O) of precipitation and groundwater suggested that the WLD and SMX recharge at the adjacent zone with mean elevations of 878 m and 676 m, respectively; while ZK03 recharge at mean elevations of 1338 m from thousands of meters away. Mean residence times calculated by artificial and natural tracers (18O, CFC) of conduit flow in WLD (7–19 h), fracture flow in WLD and SMX (243 and 600 days), and matrix flow in ZK03 (41 years) were obviously increased gradually. Hydro-chemical analysis shows that ZK03 has higher concentrations of ions that related to carbonate dissolution, indicating groundwater in ZK03 has longer time to interact with rocks. The results show that GFP in study aquifer are characterized by local fast conduit flow, local intermediate fracture networks flow and regional slow matrix flow with progressively slower groundwater refresh rate. The conceptual model was also proposed to enable a better understanding of groundwater circulation, water resources and vulnerability of the karst system in Central China.

A Holistic Approach to Study Groundwater-Surface Water Modifications Induced by Strong Earthquakes: The Case of Campiano Catchment (Central Italy)

Carbonate aquifers are characterised by strong heterogeneities and their modelling is often a challenging aspect in hydrological studies. Understanding carbonate aquifers can be more complicated in the case of strong seismic events which have been widely demonstrated to influence groundwater flow over wide areas or on a local scale. The 2016–2017 seismic sequence of Central Italy is a paradigmatic example of how earthquakes play an important role in groundwater and surface water modifications. The Campiano catchment, which experienced significant discharge modifications immediately after the mainshocks of the 2016–2017 seismic sequence (Mmax = 6.5) has been analysed in this study. The study area is within an Italian national park (Sibillini Mts.) and thus has importance from a naturalistic and socio-economic standpoint. The research strategy coupled long-period artificial tracer tests (conducted both before and after the main earthquakes), geochemical and discharge analyses and isotope hydrology with hydrogeological cross-sections. This study highlights how the seismic sequence temporarily changed the behaviour of the normal faults which act predominantly as barriers to flow in the inter-seismic period, with water flow being normally favoured along the fault strikes. On the contrary, during earthquakes, groundwater flow can be significantly diverted perpendicularly to fault-strikes due to co-seismic fracturing and a consequent permeability increase. The interaction between groundwater and surface water is not only important from the point of view of scientific research but also has significant implications at an economic and social level.

Detection and Quantification of Dam Leakages Based on Tracer Tests: A Field Case Study

Leakage is a common phenomenon in dams, and its early detection is critical for dam safety. In the present study, a new method based on tracer tests is applied to detect and quantify leakage in the Wanyao Dam, Jiangshan City, China. The objective is to detect the leakage zone of a dam wall by combining the natural tracer test and the artificial tracer test. Temperature, electrical well-logging tests with nature tracers, and the artificial tracer test with salt (NaCl) were conducted using 48 and 5 pre-existing boreholes, respectively. Using natural tracer tests, the 48 boreholes are categorized into 4 leakage classes: (1) Class 1, high connectivity within whole borehole; (2) Class 2 high connectivity at lower depths; (3) Class 3, weaken connectivity; and (4) Class 4, safe boreholes with no connectivity. Using the proposed method, specific leakage rates of some boreholes were estimated. The results of the new method are validated by comparison with those from natural tracer tests, site-investigation, and historical observation data. Overall, the new tracer test has the following merits: (1) low cost, (2) environment friendliness, and (3) is simple to apply. Moreover, the proposed method improves the accuracy of traditional tracer tests for detecting leakage zones.

Transport-based source tracking of contaminants in a karst aquifer: Model implementation, proof of concept, and application to event-based field data

Identification and location of contamination sources is crucial for water resource protection — especially in karst aquifers which provide 25% of the world´s population with water but are highly vulnerable to contamination. Transport-based source tracking is proposed and verified here as a complementary approach to microbial and chemical source tracking in karst aquifers for identifying and locating such sources of contamination and for avoiding ambiguities that might arise from using one method alone. The transport distance is inversely modelled from contaminant breakthrough curves (BTC), based on analytical solutions of the 1D two-region non-equilibrium advection dispersion equation using GNU Octave. Besides the BTC, the model requires reliable estimates of transport velocity and input time. The model is shown to be robust, allows scripted based, automated 2D sensitivity analyses (interplay of two parameters), and can be favourable when distributed numerical models are inappropriate due to insufficient data. Sensitivity analyses illustrate that the model is highly sensitive to the input time, the flow velocity, and the fraction of the mobile fluid region. A conclusive verification approach was performed by applying the method to synthetic data, tracer tests, and event-based field data. Transport distances were correctly modelled for a set of artificial tracer tests using a discharge-velocity relationship that could be established for the respective karst catchment. For the first time such an approach was shown to be applicable to estimate the maximum distance to the contamination source for coliform bacteria in karst spring water combined with microbial source tracking. However, prediction intervals for the transport distance can be large even in well-studied karst catchments mainly related to uncertainties in the flow velocity and the input time. Using a maximum transport distance is proposed to account for less permeable, “slower” pathways. In general, transport-based source tracking might be used wherever transport can be described by the 1D two-region non-equilibrium model, e.g. rivers and fractured or porous aquifers.

The Role of Faults in Groundwater Circulation before and after Seismic Events: Insights from Tracers, Water Isotopes and Geochemistry

The interaction between fluids and tectonic structures such as fault systems is a much-discussed issue. Many scientific works are aimed at understanding what the role of fault systems in the displacement of deep fluids is, by investigating the interaction between the upper mantle, the lower crustal portion and the upraising of gasses carried by liquids. Many other scientific works try to explore the interaction between the recharge processes, i.e., precipitation, and the fault zones, aiming to recognize the function of the abovementioned structures and their capability to direct groundwater flow towards preferential drainage areas. Understanding the role of faults in the recharge processes of punctual and linear springs, meant as gaining streams, is a key point in hydrogeology, as it is known that faults can act either as flow barriers or as preferential flow paths. In this work an investigation of a fault system located in the Nera River catchment (Italy), based on geo-structural investigations, tracer tests, geochemical and isotopic recharge modelling, allows to identify the role of the normal fault system before and after the 2016–2017 central Italy seismic sequence (Mmax = 6.5). The outcome was achieved by an integrated approach consisting of a structural geology field work, combined with GIS-based analysis, and of a hydrogeological investigation based on artificial tracer tests and geochemical and isotopic analyses.

Tracer tests to infer the drainage of the multiple porosity aquifer of Luxembourg Sandstone (Grand-Duchy of Luxembourg): implications for drinking water protection

The Luxembourg Sandstone is the main strategic aquifer of the Grand-Duchy of Luxembourg, providing one third of the total drinking water resource in the country. In the unconfined part, the aquifer discharges through more than 400 springs, which are used for supply. The long-lasting diffuse contamination by pesticides (more recently Metolachor ESA) has been previously simulated assuming the aquifer as an equivalent porous medium, and using environmental tracers and baseflow characteristics. A statistical analysis of the whole data set of artificial tracer tests made for this aquifer (112 injections involving 102 springs) emphasized the role of the main fissures (NE–SW) for preferential groundwater flow and the importance of including them in a conceptual model involving a multiple-porosity aquifer. The results show that these fissures are influencing the whole flow system. The average peak velocity is 16 m/h and the average first arrival (maximum) velocity is 52 m/h (32 m/h for injection-well to spring traces alone), with a recorded maximum velocity of 350 m/h. The threats of infiltrating streams within catchments have also been characterized, suggesting the usefulness of tracer tests for optimizing protective measures. The main benefit of the study is the improvement in the ability to delineate safeguard zones accurately, according to more representative velocities and the strong anisotropy.

Identification of Karst Underground River Catchment Areas with Artificial Tracer Tests and Water Balance in Banteng Cave Springs (Karst Gombong Selatan, Central Java)

The karst hills of Gombong Selatan have abundant potential water resources, especially in locations that have underground springs and rivers. The connectivity between the subsurface passageways that is difficult to know can threaten the potential of water resources, one of which is due to pollution caused by uncontrolled human activities. Therefore, identification of catchment systems and boundaries of water catchment areas in karst aquifers is needed that can contribute to sustainable water resources management policies. This catchment identification needs to be conducted because previous studies have never explored the eastern side of this karst area. This study aims to (1) define the underground river flow connectivity (upstream-downstream) of Banteng Cave; and (2) limiting the water catchment area of Banteng Cave. The method used to determine the subsurface connectivity system was carried out through an artificial tracer test, while the catchment area was delineated using a water balance approach. The results showed that the underground river of Banteng Cave has connectivity with Lake Blembeng, as evidenced by a change in watercolour after the tracer test and breakthrough curve (BTC) analysis of the tracing test results. BTC analysis shows that the Banteng Cave passageway has one main passage and does not have a tunnel branch. The estimated area of the Banteng Cave catchment used a water balance approach, which is 141.73 hectares. The Banteng Cave karst catchment conditions are dominated by the formation of valleys and karst cones accompanied by the appearance of valleys and karst hills that are quite evenly distributed, indicating that the Banteng Cave karst catchment is included in the advanced karst development phase. Furthermore, this research contributes significantly to increase knowledge regarding the characteristics of void karst development in aquifers which in the future are very important for determining water resources management policies.

Hydrogeological Characterization of the Gunungsewu Karst Area in the Upper Reach of Gremeng Karst Drainage Basin, Indonesia

The definition of a karst drainage basin in the aquifer of the karst area is an absolute thing to do as a water resource management unit. This research aimed to characterize the hydrogeology of the Gunungsewu Karst Area in the upper reach of the Gremeng Karst Drainage Basin. For this purpose, it was divided into three stages, namely geological survey, speleological survey, and artificial tracer test. The results indicate that the area observed lies in two or more geological formations: Semilir Formation, composed of sandstone and tuff from the ancient volcano Wonodadi eruptions, and Wonosari Formation, in which carbonate rocks consisting of massive coral limestone and bedded chalky limestone predominate. In volcanic rocks, the surface rivers have developed into allogenic streams flowing from underlying beds exposed updip. The contact between the two different formations results in the formation of ponors and springs. The artificial tracer test revealed two underground river systems with single conduits controlled by the allogenic rivers and large slope gradients in the study area.

An integrated hydrogeological approach to evaluate the leakage potential from a complex and fractured karst aquifer, example of Abolabbas Dam (Iran)

The design and construction of dams are risky in karst regions with active tectonic settings and strong heterogeneity due to their leakage potential. Hydrogeological characterization by karst-specific methods can be used to reduce the risk of improper dam construction. Numerous karstic dams have been constructed in the world in simple tectonic settings. In this study, a combination of geological and hydrogeological investigations is used with extensive data collection to evaluate the leakage potential of the Abolabbas Dam construction project at the Malagha karst aquifer, a tectonized highland karst region in the southwest of Iran. Two thrust faults which bounded the Malagha karst aquifer (overthrusting) caused extensive crushed zones, continuous lengths with high permeabilities in most of the boreholes, extensive open joints, and small-size solution cavities. The results of artificial tracer tests imply that the flow regime is diffuse. The flat water table, diffuse flow regime, extensive crushed zones, and continuous high-permeabilities most probably imply that the groundwater flows through an extensive network of joints without main conduits to drain the system. Despite a diffuse flow regime, an extensive grout curtain both in length and depth is required to prevent leakage because the loose materials of the crushed zones will be washed out without it under high hydraulic gradients after dam impoundment. The impermeable formation is too deep due to overthrusting. Therefore, a huge leakage may occur under the hanging grout curtain. The results allow a clearer understanding of the location and properties of fracture zones, erodible materials, high permeability zones, groundwater flow regime, and zones of preferential flow which are essential for reducing the failure of dam construction in such environments. Overall, it is demonstrated that integrating different complementary methods can help to minimize the uncertainties inherent in engineering problems in similar karst scarce-data regions.

Groundwater circulation and earthquake-related changes in hydrogeological karst environments: a case study of the Sibillini Mountains (central Italy) involving artificial tracers

Using artificial tracer tests, this study confirms the presence of a single basal aquifer feeding the springs in the wide and complex hydrogeological boundary of the Sibillini Mountains of central Italy. The tracer was introduced into the sinkhole of the Castelluccio di Norcia plain. The tracer test results, observed at the studied springs, highlight the changes induced by the 2016/2017 earthquake in the water circulation of the aquifer system. In particular, the seismic events increased the hydraulic conductivity of the basal aquifer, with a consequent increase in the springs’ flow rates at the western hydrogeological boundary and a decrease in the flow rates at the eastern Adriatic hydrogeological boundary. This phenomenon is in accordance with the hydro-structural framework of the area. The study also investigated the relation between groundwater circulation and tracer behaviour in the springs during the pre- and post-earthquake periods. The tracer test results led to the formulation of hypotheses about water circulation of the area. The trend of the tracer breakthrough curves demonstrates that the upper portion of the basal aquifer is characterised by fast water circulation similar to that in the epiphreatic area of karst aquifers, while a slow circulation due to fissures with interconnected drains occurs in the deepest portion of the aquifers. The obtained results highlight how tectonics and karst can affect the hydrogeological setting of the Apennine carbonate chain; in particular, the seismicity of central Italy may alter groundwater circulation for a long period of time when great magnitude earthquakes occur.

Variability in Groundwater Flow and Chemistry in the Houzhai Karst Basin, Guizhou Province, China

ABSTRACT Understanding how karst aquifers store and transmit water and contaminants is an ongoing problem in hydrogeology. Multiple flow paths and recharge heterogeneity contribute to the complexity of these systems. This study explored karst-conduit connectivity and water-chemistry variability within the Houzhai catchment in Guizhou Province, China. Artificial tracer tests were conducted during both the monsoon and dry seasons to understand temporal variability in connectivity and water velocity between karst features. Multiple flow paths through the catchment were activated during the monsoon season and partially abandoned during the dry season. Additionally, gradient reversals during monsoonal high-flow events and as a result of pumping were observed. Synoptic water samples from several karst features taken during both monsoon and dry seasons elucidated spatial and temporal variability within the catchment. Water residence time was generally longer during the dry season, and flow within the Houzhai catchment was determined to be temporally dependent. Time-series sampling at the outlet spring following a monsoonal storm event captured chemical variability and identified multiple flow paths. Overall, this study refines widely applicable methods for studying karst systems to this catchment and provides a foundation for future studies in similar settings.

Coupling SKS and SWMM to Solve the Inverse Problem Based on Artificial Tracer Tests in Karstic Aquifers

Artificial tracer tests constitute one of the most powerful tools to investigate solute transport in conduit-dominated karstic aquifers. One can retrieve information about the internal structure of the aquifer directly by a careful analysis of the residence time distribution (RTD). Moreover, recent studies have shown the strong dependence of solute transport in karstic aquifers on boundary conditions. Information from artificial tracer tests leads us to propose a hypothesis about the internal structure of the aquifers and the effect of the boundary conditions (mainly high or low water level). So, a multi-tracer test calibration of a model appeared to be more consistent in identifying the effects of changes to the boundary conditions and to take into consideration their effects on solute transport. In this study, we proposed to run the inverse problem based on artificial tracer tests with a numerical procedure composed of the following three main steps: (1) conduit network geometries were simulated using a pseudo-genetic algorithm; (2) the hypothesis about boundary conditions was imposed in the simulated conduit networks; and (3) flow and solute transport were simulated. Then, using a trial-and-error procedure, the simulated RTDs were compared to the observed RTD on a large range of simulations, allowing identification of the conduit geometries and boundary conditions that better honor the field data. This constitutes a new approach to better constrain inverse problems using a multi-tracer test calibration including transient flow.