Fractured Rock Aquifers Research Articles

GROUNDWATER MANAGEMENT: INTEGRATING GEOLOGICAL AND HYDROLOGICAL DATA FOR EFFECTIVE DECISION MAKING

The most valuable natural resource for all of humanity, not just for a state or nation, is water. A country's ability to prosper largely rests on how wisely it uses this resource. Thus, water, which flows in rivers and streams, can be said to be a nation's primary asset. This proves rivers are important, and further justification is needed to emphasise how important they are. The global acceptance of river basins as planning and management domains stems from the fact that water transcends national boundaries. India's rivers are one of its most defining characteristics; the people of India place great religious significance on them. Since groundwater is the most abundant freshwater resource, it is essential to pinpoint the groundwater potential zones for the ongoing and sustainable growth of socioeconomic and agricultural advances. Hard, fractured rock aquifers underlie most of south India's densely inhabited areas. The meagre amount of freshwater that these aquifers provide is the only source that people may use for cultivation and drinking. Because of the diverse nature of the aquifer owing to varied composition, degree of weathering, and density of fracturing, groundwater conditions in hard rock terrain are multivariate. Therefore, a thorough understanding of the distribution and features of the aquifer under the various hydrogeologic circumstances of the research region is necessary for the long-term utilisation of groundwater resources. This study's primary goal is to evaluate the possibility for groundwater management utilising hydrological and geological investigations.

Quantification of immiscible displacement in a rough-walled fracture through direct visualization

Subsurface substance migration in the fractured rock aquifer is mainly controlled by fractures, and immiscible fluid-fluid displacement in fractures is important to many geophysical processes and engineering activities. Using a fracture-visualization system, we present the qualitative and quantitative assessment of fracture geometry associated with fluid movement and distribution in the rough fracture. Based on fracture geometry and statistical analysis, we first conducted a quantitative study of fracture surface roughness and aperture distribution. Then, fractal dimensions of displacement front and residual oil distribution were determined by image processing procedures. Influenced by wettability and micro-scale roughness, at the end of the displacement, residual oil saturation of molded sample is lower (6.45%–25.74%), and displacement pattern is more uniform, indicating that displacement effect is better. Due to smaller differences in residual oil saturation (9.08%) under different injection directions, the impact of wettability on the displacement process is greater than that of anisotropic roughness. Additionally, the fractal dimension of the displacement front increased under low injection rates initially but decreased when the rate was increased later. Overall, visualized temporal monitoring of experimental images enabled us to provide a preliminary assessment of the impact of anisotropic roughness and the material constituting the fracture wall on invading fluid saturation and the fractal dimension of the displacement front under various injection rates.

Anisotropic effect of shear zone on groundwater potentiality in crystalline hard rock terrain

Assessment of potential groundwater recharge sites and sustainable water resource management in semi-arid crystalline rock terrain is a challenging task. Globally, analysis of remote sensing satellite imagery data for delineation of groundwater potential zones over sheared crystalline hard rock terrains has been fairly successful. But there is no existing study present at our disposal which discusses the factors controlling the inconsistent groundwater potentiality that exits along the shear zones. This study attempts to analyse the major geological factors controlling the irregular groundwater potentiality of shear zones within older crystalline rock terrain. Therefore, the study area selected for this analysis is the Purulia district of West Bengal, NE India, composed mostly of Precambrian metamorphic rocks i.e., quartzite, granite gneisses, porphyroclastic granite-gneiss, quartzo-feldspathic-granite-gneiss, mylonitic granites, quartz-biotite-granite gneiss, quartzites, carbonatites and phyllites. Satellite imagery study of IRS-P6 LISS IV standard FCC image reveals the presence of two bifurcating shear zones namely North Purulia Shear Zone (NPSZ) and South Purulia Shear Zone (SPSZ) over the study area. Careful analysis of rock structure, different lithotypes, soil thickness, electrical resistivity tomography data and water table data with an emphasis on high water table fluctuation, shows a strong spatial relation between the potentially good groundwater recharge zones and the branching/confluence sites of shear zones present in the study area. The study constructs an attempt to demonstrate the relationship between shear zone conjunctions and significant groundwater recharge sites in Precambrian crystalline fractured-rock aquifer system.

An integrated approach for characterization of a fractured-rock carbonate aquifer in the Zagros Region of Iran

Karst aquifers typically exhibit a wide range of dissolution effects that include end-members of matrix, fractured-rock, and conduit-dominated types. This study employs an integrated approach involving geological studies, hydrogeological assessments, hydrochemical and isotopic analysis, and dye tracer tests to investigate the Sarvak limestone aquifer (SLA) in the Zagros Region, southwest Iran. Key characteristics of the SLA include numerous stratification boundaries, thin limestone layers with intercalated siliceous and marl impurities, extensive joint and fracture networks, predominant autogenic recharge, and absence of notable point recharge features (e.g., sinkholes, dolines), and the exchange flow with the Bakhtiari River (B-R) has made it a unique aquifer. Numerous pieces of evidence, such as low spatial and temporal changes in groundwater level, insignificant seasonal variations in the hydrochemical and isotopic composition of water samples, and the supersaturation state of groundwater with calcite and dolomite minerals, suggested a slow flow regime in many parts of the SLA. This regime is characterized by low velocity and long residence time of groundwater. The results reveal that despite the high solubility of carbonate rocks, extensive joint networks can limit significant karst development, leading to fluid flow behavior similar to that of fractured-rock aquifers. Therefore, SLA can be considered a fractured-rock and karst aquifer. The identified paths with fast flow in the dam area are unremarkable and cannot provide a large portion of the Pelle Khan Spring discharge.

Impact of the model structure and calibration strategy on baseflow modeling in the German low mountain range

ABSTRACT Baseflow is a vital component of the water balance. The fractured hard rock aquifers of the German low mountain range are in danger of increased water stress due to climate change because they react rapidly to deficits in precipitation and groundwater tables decline sharply. Therefore, simulation software must be able to model baseflow accurately. Three soil moisture simulation and two monthly factor-based baseflow models are evaluated using two calibration strategies. Models were calibrated to total flow (S1) or stepwise to baseflow and then total flow (S2). Results were not significantly different for total flow. Regarding baseflow, S2 proved significantly better with median values (S1 calibration, validation | S2 calibration, validation) of SSE (20.3, 20.3 | 13.5, 13.8), LnNSE (0.15, 0.17 | 0.47, 0.34), and PBIAS (27.8, 21.6 | 2.5, −0.8). Parallel linear reservoir proved best at modeling baseflow with a median SSE (S2: 6.1, 5.9), LnNSE (S2: 0.64, 0.71), and PBIAS (S2: 3.8, 3.8). The new modified monthly factor approach is a simple and robust alternative with SSE (13.0, 13.3), LnNSE (0.61, 0.61), and PBIAS (9.8, −8.6). The results are useful regarding selection of baseflow model structure and calibration strategy in low mountain ranges with fractured hard rock aquifers.

Hydrochemistry and Contamination in Sedimentary and Fractured-Rock Aquifers of The Upper Rio Negro Basin – Amazonas/Brazil

The Amazon region has the largest hydrographic network on the planet. However, cities and communities along rivers obtain water from tube wells without studying the chemical quality and, consequently, the potability of their water. In specific geological contexts, there is a significant probability of natural contamination by arsenic in groundwater in the region, already attested in previous works in the Solimões, Madeira and Amazonas rivers basins, but not yet found in the Negro river basin. As the Negro is the second largest river in the Amazon Basin, it has become necessary to study the chemical quality of its groundwater, relating it to the geology of the region. In this work, water samples were collected from 21 wells and 4 water springs used for human consumption, following the protocols proposed by the US – EPA. In the field, the physicochemical parameters of the in-situ water were measured, and alkalinity tests were carried out on the samples. The cations were analyzed at INPA and the anions at the Geological Survey of Brazil (SGB). From the 25 samples, 16 had analises results within acceptable ionic balance values, and 12 of the 16 had at least one of the following elements in concentration above World Health Organization and Brazilian Guideline Values: Al, As, Fe, Mn and Pb. The presence of values below the GV, but above what is considered natural, of NO3- in 66% of the contaminated points, associated with the shallow depth of most of the wells, installed in areas of recent sedimentary deposits and that do not show any or few signs of care demonstrate that the contamination in these localities is of anthropic origin. Two samples (16% of the contaminated samples) are from deep wells installed in fractured aquifers, which in this case indicate contamination of geogenic origin. Arsenic was found in association with iron and aluminum in a well installed in alluvial sediments, which indicates that more studies are needed for this specific geologic formation to know the quality of this water. For other locations and cities, short, medium and long-term measures are needed to prevent anthropogenic contamination of aquifers.

Evaluation of Passive Vapor Diffusion Samplers to Quantify Acetylene, Ethene, and Ethane in Groundwater

AbstractAcetylene, ethene, and ethane are products from the degradation of trichloroethene (TCE) in contaminated aquifers. However, the volatility of these gases makes it challenging to avoid excessive losses during sampling. The objective of this study was to compare the quantification of acetylene, ethene, and ethane using passive vapor diffusion (PVD) samplers vs. conventional low‐flow groundwater collection. Samples were obtained from 8 to 13 monitoring wells at three sites that show potential for biotic and abiotic degradation of TCE in fractured rock aquifers. Method reporting limits (MRLs) for the PVD samplers were 0.25 μg/L for acetylene (0.0094 μM) and 0.28 μg/L for ethene and ethane (0.0099 and 0.0092 μM, respectively); the MRLs for conventional low‐flow groundwater samples were ~40% higher. For two of the sites, the maximum concentrations of acetylene, ethene, and ethane obtained with the PVD samplers were comparable to the conventional low‐flow samples. The frequency of detection for these gases with the PVD samplers was also comparable to conventional low‐flow groundwater sampling. For one of the sites with higher levels of acetylene (maximum of 13 μg/L), the concentrations from the PVD samplers were approximately twofold higher than those with conventional low‐flow groundwater sampling. Based on robust detection of acetylene, ethene, and/or ethane, it appears likely that degradation of TCE is occurring at the three sites that were sampled. The use of PVD samplers can reduce the possibility of false negative results to provide another line of evidence in support of natural attenuation.

Wastewater contaminants in a fractured bedrock aquifer and their potential use as enteric virus indicators.

Domestic wastewater is a source of persistent organic pollutants and pathogens to the aquatic environment, including groundwater aquifers. Wastewater contaminants include a variety of personal care products, pharmaceuticals, endocrine disrupters, bacteria, and viruses. Groundwater from 22 wells completed in a semi-confined to confined, fractured Silurian dolostone aquifer in southern Wellington County, Ontario, Canada, was analyzed for 14 organic wastewater contaminants (4 artificial sweeteners, 10 pharmaceuticals) as well as E. coli, total coliforms, and 6 human enteric viruses. Enteric viruses were detected in 8.6% of 116 samples, and at least one organic wastewater contaminant was detected in 82% of the wells (in order of decreasing detection frequency: acesulfame, ibuprofen, sulfamethoxazole, triclosan, carbamazepine, and saccharin). Virus indicator metrics [positive and negative predictive values (PPV, NPV), sensitivity, specificity] were calculated at the sample and well level for the organic wastewater compounds, E. coli, and total coliforms. Fecal bacteria were not good predictors of virus presence (PPV = 0%-8%). Of the potential chemical indicators, triclosan performed the best at the sample level (PPV = 50%, NPV = 100%), and ibuprofen performed the best at the well level (PPV = 60%, NPV = 67%); however, no samples had triclosan or ibuprofen concentrations above their practical quantification limits. Therefore, none of the compounds performed sufficiently well to be considered reliable for assessing the potential threat of enteric viruses in wastewater-impacted groundwater in this bedrock aquifer. Future studies need to evaluate the indicator potential of persistent organic wastewater contaminants in different types of aquifers, especially in fractured rock where heterogeneity is strong.IMPORTANCEAssessing the potential risk that human enteric viruses pose in groundwater aquifers used for potable water supply is complicated by several factors, including: (i) labor-intensive methods for the isolation and quantification of viruses in groundwater, (ii) the temporal variability of these viruses in domestic wastewater, and (iii) their potentially rapid transport in the subsurface, especially in fractured rock aquifers. Therefore, aquifer risk assessment would benefit from the identification of suitable proxy indicators of enteric viruses that are easier to analyze and less variable in wastewater sources. Traditional fecal indicators (e.g., E. coli and coliforms) are generally poor indicators of enteric viruses in groundwater. While many studies have examined the use of pharmaceutical and personal care products as tracers of domestic wastewater and fecal pollution in the environment, there is a paucity of data on the potential use of these chemical tracers as enteric virus indicators, especially in groundwater.

Environmental Tracers Tritium 3H and SF6 Used to Improve Knowledge of Groundwater Sustainability of a Crystalline Rock Island Aquifer of Tobago, West Indies

Understanding the residence time, flow velocities, and storage volumes in fractured rock aquifers is essential for assessing the sustainable use of groundwater in island aquifer systems. Environmental tracers such as tritium (3H) and sulfur hexafluoride (SF6) that enter the aquifer systems through natural hydrological processes are effectively used to quantify the “apparent age” of young groundwater (<60 years) to calculate these aquifer characteristics. The island of Tobago is predominantly a fractured rock aquifer system that contains groundwater sub-basins that produce 100–1000% more than their calculated recharge. In this study, we analyzed ten samples from groundwater wells throughout the island for environmental tracers to establish the “apparent age” of the groundwater. We then estimated the active flux and storage volume of the aquifer. The 3H results of ten wells throughout the island range between 0.02 and 0.60 TU. The groundwater samples to the south of the island possess more pre-modern groundwater ages (older ages) than samples in the north of the island. These tritium values reveal age ranges between 18 and 60–plus years. This finding suggests that fractures create flow paths oriented towards the southern parts of the basin that enable the transport of groundwater across significant topographic boundaries and watershed divides. Additionally, the data may support that southern wells are connected to isolated old sources of groundwater. SF6 values range between 0.91 and 7.97 (part per trillion volume) with an interpreted age range between 18 and 39 years. Both tracers support the original hypothesis that longer residence time waters are located south of the island except for three SF6 samples which are believed to be affected by “excess air”.

Experimental study on non-Darcian flow through a single artificial fracture for different fracture apertures and surface roughness

Abstract This study aims to explore the influence of various geometrical and hydraulic parameters on flow behavior and hydraulic conductivity in a single artificial fracture through a series of laboratory experiments. Laboratory experiments were conducted to examine unconfined groundwater flow through an artificially constructed single fracture. The fracture model consisted of varying aperture sizes (3, 9, and 12 mm) and different surface roughness conditions (fine, medium, and coarse sand coatings). Non-Darcian turbulent flow characteristics were observed at different flow rates, and the gradient of Reynolds number versus average flow velocity increased with aperture size. Flow parameters of the Darcian, Izbash, and Forchheimer models were calculated to characterize the flow behavior. Both the Forchheimer and Izbash models were found suitable for describing the non-Darcian flow characteristics under the prevailing conditions. The study revealed that hydraulic conductivity depended on flow length for fractures with different apertures and surface roughnesses, likely due to the presence of 2-D torturous flow within the rough fracture surface. These findings contribute to a better understanding of groundwater flow in fractured rock aquifers and provide valuable insights for modeling and managing such systems.

Benefits and Cautions in Data Assimilation Strategies: An Example of Modeling Groundwater Recharge.

Assimilating recent observations improves model outcomes for real-time assessments of groundwater processes. This is demonstrated in estimating time-varying recharge to a shallow fractured-rock aquifer in response to precipitation. Results from estimating the time-varying water-table altitude (h) and recharge, and their error covariances, are compared for forecasting, filtering, and fixed-lag smoothing (FLS), which are implemented using the Kalman Filter as applied to a data-driven, mechanistic model of recharge. Forecasting uses past observations to predict future states and is the current paradigm in most groundwater modeling investigations; filtering assimilates observations up to the current time to estimate current states; and FLS estimates states following a time lag over which additional observations are collected. Results for forecasting yield a large error covariance relative to the magnitude of the expected recharge. With assimilating recent observations of h, filtering and FLS produce estimates of recharge that better represent time-varying observations of h and reduce uncertainty in comparison to forecasting. Although model outcomes from applying data assimilation through filtering or FLS reduce model uncertainty, they are not necessarily mass conservative, whereas forecasting outcomes are mass conservative. Mass conservative outcomes from forecasting are not necessarily more accurate, because process errors are inherent in any model. Improvements in estimating real-time groundwater conditions that better represent observations need to be weighed for the model application against outcomes with inherent process deficiencies. Results from data assimilation strategies discussed in this investigation are anticipated to be relevant to other groundwater processes models where system states are sensitive to system inputs.

Estimation of specific yield in bedrock near-surface zone of hilly watersheds by examining the relationship between base runoff, storage and groundwater level

A catchment-scale method for estimation of specific yield (Sy) in the zone of groundwater level fluctuation is proposed. It is applicable to hilly watersheds, where deep groundwater discharge – not drained by local streams as baseflow – is small and can be neglected. Therefore, it is mostly employable for bedrock flow systems, dominated by shallow unconfined fractured rock aquifers. Method provides an estimate of specific yield (Sy) by combined analysis of streamflow recession, storage/runoff relationship and groundwater level fluctuation (Q-S-H). For groundwater storage (S) values evaluation, river discharge (Q) master recession curves are constructed and interpreted. The method produces as more reliable results as number of groundwater level observation sites increases. As example, it is demonstrated at the Levočský potok Brook watershed (Western Carpathians, Slovakia), built by fracture porosity dominated Paleogene sediments. Estimated characteristic Sy value is from the interval 0.001–0.002 and 0.002–0.005 for low and medium storage/runoff conditions – or bottom and middle part of GWL fluctuation zone – respectively.

A global inertial permeability for fluid flow in rock fractures: Criterion and significance

Fluid flow in rock fractures is described by Darcy's law and its nonlinear extension, the Forchheimer equation in form of a quadratic polynomial. Since the nonlinearity of the Forchheimer equation, determining its key coefficient of inertial permeability (i.e., the quadratic coefficient representing inertial losses) requires a series of flow tests under different hydraulic gradients. However, how to choose the volumetric flux range for deriving inertial permeability is never scrutinized. Here we report a seemly irrational dependency of inertial permeability on volumetric flux range by massive direct numerical simulations of fluid flow in rock fractures. This variation presents traceability and is closely related to the evolution of microscale eddies and macroscale flow regimes. A concept of global inertial permeability is thus proposed from this traceable variation, which can reproduce the whole flow regime in rock fractures. Two parametrized criterion models are subsequently established for calculating global inertial permeability using the lowest cost and ensuring its accuracy: one in terms of minimum Reynolds number for practical purposes, and the other in terms of minimum eddy volume ratio for mechanistic interpretation. Further discussions indicate that the proposed global inertial permeability is of great importance to non-Darcian flow prediction and flow regime assessment. Moreover, the parameterized criterion models in the form of power-law function for determining global inertial permeability, drawn from 2D fracture simulations, can apply to 3D fractures and large-scale fractured rock aquifers. The findings and results in this study are of great significance for accurately evaluating the hydraulic conductivity of rock fractures, especially at relatively high flow rates, and are useful for many geophysical and industrial applications.

Four-dimensional characterization of a PFOA-contaminated fractured rock aquifer (FRA) in Bennington, Vermont, U.S.A

Between 2016 and 2019, the Vermont Department of Environmental Conservation discovered that groundwater from over 50% (355/628) of wells and springs in the Bennington area had perfluorooctanoic acid (PFOA) concentrations > 20 parts/trillion (ppt) (Vermont action level). The distribution of contaminated wells was complex, with very different PFOA concentrations in closely-spaced wells (<100 m apart). PFOA was introduced into the environment by a factory that produced Teflon™-coated fabrics from about 1978–2002. Manufacturing involved PFOA expulsion from smokestacks followed by atmospheric transport, deposition, and downward leaching to the aquifer system. The fractured rock aquifer was characterized through physical components (geologic mapping, spatial analysis of wells, and geophysical logging) to build a Three-dimensional conceptual model and then map the spatial distribution of chemical groundwater tracers (PFOA, major and trace elements, stable isotopes, and and recharge ages), adding a fourth temporal dimension. The field area consists of four imbricated metamorphic rock slices bounded by thrust faults. The dominant fracture sets that overprint these slices control the N-S and E-W topographic “grain.” Geophysical logging divided wells into those completed internal to or on the margins of thrust slices. Major and trace element analysis shows that groundwater can be spatially discriminated into carbonate (TDS, Ca, Mg, HCO3), shaly (Sr, U, SO4), and siliceous (K, Si) groups. Depleted H2O stable isotope signatures characterize the siliceous group. Recharge-ages determined from CFC11, 12, and 113 and tritium correlate with geochemical groups and were divided into “older” (~1953–73) and “younger” (1963–88), which predate and are synchronous with factory operation, respectively. Our conceptual model shows recharge from highlands to the west and east flowing toward the Bennington valley where it is locally tapped by wells along thrust faults or fracture zones. This water has older recharge ages and anomalously low PFOA concentrations (<6.7 ppt in 12 of 15 samples, mean = 7.9 ppt) compared to the majority of wells in this region, which are completed in the internal portions of thrust slices. These wells contain groundwater with younger recharge ages and high PFOA concentrations (mean = 290 ppt) and occurrence (>20 ppt in 28/36 samples), indicating mixing of younger and older waters.

Hydrogeochemical evolution of the shallow and deep basaltic aquifers in Tamborine Mountain, Queensland (Australia)

Tamborine Mountain, Queensland (Australia), is a prime example of a basalt fractured-rock aquifer. Yet very little is known about the hydrochemistry of this groundwater system. Both analytical (major ions and stable isotopes) and multivariate (hierarchical cluster analysis, principal component analysis and factor analysis) analyses were used in this study to investigate the factors that interact within this aquifer system, in order to determine groundwater hydrogeochemistry. A new approach was applied to the data by classifying hydrographs by water type to clearly identify differing aquifer zones. Three distinct groundwater chemistry types were identified, and they were differentiated by variations in depth. Shallow bores were dominated by Na–Cl waters, deep bores were dominated by Na–HCO3 and Ca–HCO3 waters, and the two deepest bores were dominated by mixed water types. The evaluation of hydrogeochemical data has determined that both mineral weathering processes and groundwater/surface-water interaction had a strong influence on the hydrogeochemistry. Seasonal effects were minimal in the study area based on physicochemical parameters and ion chemistry. However, stable isotopic data show temporal trends. Increased rainfall events during the wet season produced a depletion in δ18O and increased d-excess values. The opposite is found during the dry season as a result of higher evaporation rates that are not hindered by intense rainfall events.

Geometry‐Based Prediction of Solute Transport Process in Single 3D Rock Fractures Under Laminar Flow Regime

AbstractSubsurface substance migration in the fractured rock aquifer is mainly controlled by fractures, and prediction of solute transport in fractures is important to many geophysical processes and engineering activities. This study explores the possibility of predicting transport process in single rock fractures from measurable physical properties. For this purpose, we conducted a large number of pore‐scale simulations of solute transport in single 3D rock fractures with different apertures and various degrees of roughness. The numerically‐obtained breakthrough curves under laminar flow regime were reproduced with high fidelity by the classic analytical solution within the framework of macroscopic advection‐dispersion theory. The fitted transport coefficients (hydrodynamic dispersion coefficient and transport velocity) were normalized by the parallel‐plate model's values, and were further parameterized by aperture b and roughness parameter σ/b in form of ωeb/ψ + λ(mσ/b − 1)/bn and 1 − α(σ/b)/βb, respectively. By incorporating these parameterized transport coefficients and the modified cubic law into the classic analytical solution, we proposed a new transport predictive model. This model was adaptively degenerated into the classic analytical solution of parallel‐plate model, and was proved to be capable of predicting macroscopic solute transport only relying on arithmetic‐mean mechanical aperture and its standard deviation, without solving the velocity and concentration fields. By the established model, we quantitatively revealed the impacts of aperture and roughness on the breakthrough curve under laminar flow regime. This study puts a step toward predicting solute transport process in fractured rock aquifers with measurable geometric properties.

The impact of livestock activities and geochemical processes on groundwater quality of fractured volcanic rock aquifer: Lake Çıldır watershed (NE Turkey)

This paper presents the impact of livestock activities and geochemical processes on the water quality of a fractured volcanic rock aquifer in the Lake Çıldır watershed, located at the northeastern part of Turkey. The existence of a high livestock population and animal grazing activities in meadow and pasturelands of the watershed during the short summer period poses serious stress on both surface and groundwater resources being the only drinking water supply for the local communities. Therefore, understanding the effect of grazing and livestock breeding activities occurring in the recharge areas of the fractured volcanic rock aquifer is vital to take precautions in order to protect limited water supplies at the watershed and vulnerable lake ecosystem as well. The mean nitrate content of the groundwater was measured at 6.4 ± 6.6 (std. dev) mg/L in the wet (before grazing) period and 7.1 ± 5.9mg/L in the dry (after grazing) period. Despite low nitrate concentration levels of groundwater, microbial contamination was observed in the spring waters at alarming levels especially after the animal grazing activities. 56%, 26%, and 11% of the groundwater samples showed bacterial contamination in terms of total coliform, fecal coliform, and fecal streptococci contents, respectively, prior to grazing activity, while in pursuit of intense livestock grazing at highland, these microbial indicators have been increased to 92%, 85%, and 77% in the dry period. A significant increase observed in fecal contamination indicates the negative impact of livestock activities on groundwater quality. Al (200-638µg/L) and Fe (66-218µg/L) enrichments locally observed in groundwater were related to advanced argillic alteration (kaolinization) and hematization zones in pyroclastic rocks.

Legacy pollutants in fractured aquifers: Analytical approximations for back diffusion to predict atrazine concentrations under uncertainty

We present novel analytical approximations for the estimation of travel distance and relative height of solute concentration peaks within a single fracture system for pollutants that have been temporarily applied at a constant rate in the past. These approximations are used to investigate the spatiotemporal evolution of the concentration of atrazine, as an example for many other so-called legacy compounds that are still found in the groundwater of fractured rock aquifers even decades after their application has stopped. This is done in a stochastic framework to account for the uncertainty in relevant parameters, focusing on probabilities of exceeding the given legal concentration limit and the expected length of the recovery period. We specifically consider the properties of the Muschelkalk limestone aquifer in the Ammer river catchment in SW Germany, and the three major types of carbonate rock facies: Shoal, Tempestite, and Basinal limestones. Atrazine sorption parameters have been determined in laboratory experiments. The simulations confirm that diffusion-limited sorption and desorption may cause considerable atrazine levels long after application stop. For the properties of the considered rock facies types, and corresponding parameter ranges, atrazine concentration above the legal limit is supposed to be limited to locations referring to only a few years of travel time. If the concentration exceeds the legal limit by the year 2022, it will take decades to centuries until recovery.

A data-driven modeling approach for the sustainable remediation of persistent arsenic (As) groundwater contamination in a fractured rock aquifer through a groundwater recirculation well (IEG-GCW®)

Persistent arsenic (As) pollution sources from anthropogenic activities pose a serious threat to groundwater quality. This work aims to illustrate the application of an innovative remediation technology to remove As from a heavily contaminated fractured aquifer at a historically polluted industrial site. Groundwater circulation well (GCW) technology was tested to significantly increase and accelerate the mobilization and removal of As in the source area. The GCW extracts and re-injects groundwater at different depths of a vertical circulation well. By pumping out and reinjecting in different screen sections of the well, the resulting vertical hydraulic gradients create recirculation cells and affect and mobilize trapped contaminants that cannot be influenced by traditional pumping systems. The first 45-m deep IEG-GCW® system was installed in 2020, equipped with 4 screen sections at different depths and with an above-ground As removal system by oxidation and filtration on Macrolite (Enki). A geomodeling approach supports both remediation and multi-source data interpretation. The first months of operation demonstrate the hydraulic effectiveness of the IEG-GCW® system in the fractured rock aquifer and the ability to significantly enhance As removal compared to conventional pumping wells currently feeding a centralized treatment system. The recirculation flow rate amounts to about 2 m3/h. Water pumped and treated by the GCW system is reintroduced with As concentrations reduced by an average of 20%–60%. During the pilot test, the recirculating system removed 23 kg As whilst the entire central pump-and-treat (P&T) system removed 129 kg, although it treated 100 times more water volume. The P&T plant removed 259 mg As per m3 of pumped and treated groundwater while the GCW removed 4814 mg As per m3 of the treated groundwater. The results offer the opportunity for a more environmentally sustainable remediation approach by actively attacking the contamination source rather than containing the plume.

Using tritium and radiocarbon activities to constrain regional modern and fossil groundwater mixing in Southern Africa

This study combines historical records of 14C and 3H in the atmosphere and soil with renewal rate and groundwater lumped parameter models to predict the abundance of 14C and 3H in groundwater over time. 624 groundwater samples from numerous studies, over four decades (1978–2019), in South Africa, Namibia, Botswana and Mozambique were collated to compare with predicted groundwater activities of 14C and 3H within the South African Development Community (SADC) region. Spatial datasets of carbonate bearing lithology, C3/C4 vegetation, summer/winter rainfall and coastal proximity were used to apply corrections to 14C and 3H data. Corrected values of 14C and 3H were compared with the theoretical abundance of these tracers, derived from the lumped parameter models, to estimate the general mean residence times and presence of groundwater mixing between modern recharge and older groundwaters. This study found that corrected values produced varying mean residence times derived from 14C ages (∼500–28500 years) and a wide range of potentially mixed waters within each aquifer system (0–100 % of tested wells) across the study area. The largest proportions of mixed groundwater, as well as the youngest mean residence times, were found in alluvial and primary fractured rock aquifers (e.g., western coast of South Africa and southern Mozambique). The smallest proportions of mixed groundwater were predicted in deep confined clay-rich aquifers as well as layered coal bearing carbonate sequences (e.g., Orapa, Malwewe and Serowe, Botswana). Insights into the proportions of mixed groundwater and mean residence times can help assess hydrological resilience on a regional scale. Such information is pertinent in promoting socio-economic development and increased water/food security in the SADC region. By understanding the resilience of groundwater resources, robust and informed strategies for water equality and GDP growth in the SADC region can be envisioned and implemented.