Clogging process by suspended solids during groundwater artificial recharge: Evidence from lab simulations and numerical modelling

Permeability reduction of infiltration media due to suspended solid (SS) clogging is the bane of groundwater artificial recharge. To overcome the clogging problem and advance the understanding of the process‐based spatial‐temporal evolution of SS clogging, a 1D laboratory column simulation was carried out, followed by numerical modelling of the experimental data in this study. It was found that clogging caused a reduction in the hydraulic conductivity (K) in the upper layer at the beginning and extended deeper to approximately 50 cm, and no reduction in K was detected below 52 cm throughout the experimental period of 129 hr. The most clogged layer spanned from the surface to a depth of 11 cm, and the middle 11–52 cm was characterized by a slight decrease in K. The clogging rates of the different layers decreased with the depth, which was based on data analysis, with the largest value of 0.038 hr^?1 in the upper 1 cm. The overall K began to decrease from the surface layer and was increasingly affected by clogging with time. A mathematical model was established to simulate the SS clogging process evolution based on considerations of the attachments and detachments of particles. Then the model was applied to perform several scenario analyses after calibration and validation using the data obtained in the experiment. The simulation results indicated that the SS concentration was much more sensitive than the groundwater depth below the land surface, and 10 days of constant recharge is recommended as the disposal cycle of the clogged layer under the given conditions.     

Characterizing Heterogeneity in Infiltration Rates During Managed Aquifer Recharge

Infiltration rate is the key parameter that describes how water moves from the surface into a groundwater aquifer during managed aquifer recharge (MAR). Characterization of infiltration rate heterogeneity in space and time is valuable information for MAR system operation. In this study, we utilized fiber optic distributed temperature sensing (FO‐DTS) observations and the phase shift of the diurnal temperature signal between two vertically co‐located fiber optic cables to characterize infiltration rate spatially and temporally in a MAR basin. The FO‐DTS measurements revealed spatial heterogeneity of infiltration rate: approximately 78% of the recharge water infiltrated through 50% of the pond bottom on average. We also introduced a metric for quantifying how the infiltration rate in a recharge pond changes over time, which enables FO‐DTS to be used as a method for monitoring MAR and informing maintenance decisions. By monitoring this metric, we found high‐spatial variability in how rapidly infiltration rate changed during the test period. We attributed this variability to biological pore clogging and found a relationship between high initial infiltration rate and the most rapid pore clogging. We found a strong relationship (R² = 0.8) between observed maximum infiltration rates and electrical resistivity measurements from electrical resistivity tomography data taken in the same basin when dry. This result shows that the combined acquisition of DTS and ERT data can improve the design and operation of a MAR pond significantly by providing the critical information needed about spatial variability in parameters controlling infiltration rates.     

Clogging of saturated porous media by silt‐sized suspended solids under varying physical conditions during managed aquifer recharge

Managed aquifer recharge is an effective method for utilizing excess flood flows, but clogging of porous media is a limiting factor in the implementation of this water storage technique. In recent years, much research on the physical clogging of porous media during artificial recharge has been conducted. However, the understanding of clogging due to silt‐sized suspended solids (SS) is still inadequate, especially under varying physical conditions. Here, we subjected sand columns to controlled rates of flow and SS suspensions to investigate the influence of media size, SS size, SS concentration, and flow velocity on the clogging of porous media by silt‐sized SS. The results show that the diameter ratio of SS particles to sand grains is the dominant factor influencing the position of physical clogging. As pore velocity increased, the mobility of silt‐sized SS was enhanced and retention in the porous media decreased noticeably. The spatial retention profiles in the porous media were found to vary greatly at different flow velocities. The SS concentration of the infiltrating suspension also dramatically influenced the mobility and deposition of silt‐sized SS particles, such that high concentrations accelerated the clogging process. As the different physical factors changed, the breakthrough curves and retention profiles of silt‐sized SS particles changed obviously and the mechanisms of retention differed. On the whole, clogging position is mainly determined by particle size ratio, but clogging rate is dominated by a variety of factors including particle size ratio, SS concentration, and flow velocity.     

Trench infiltration for managed aquifer recharge to permeable bedrock

Managed aquifer recharge to permeable bedrock is increasingly being utilized to enhance resources and maintain sustainable groundwater development practices. One such target is the Navajo Sandstone, an extensive regional aquifer located throughout the Colorado Plateau of the western United States. Spreading‐basin and bank‐filtration projects along the sandstone outcrop's western edge in southwestern Utah have recently been implemented to meet growth‐related water demands. This paper reports on a new cost‐effective surface‐infiltration technique utilizing trenches for enhancing managed aquifer recharge to permeable bedrock. A 48‐day infiltration trench experiment on outcropping Navajo Sandstone was conducted to evaluate this alternative surface‐spreading artificial recharge method. Final infiltration rates through the bottom of the trench were about 0·5 m/day. These infiltration rates were an order of magnitude higher than rates from a previous surface‐spreading experiment at the same site. The higher rates were likely caused by a combination of factors including the removal of lower permeability soil and surficial caliche deposits, access to open vertical sandstone fractures, a reduction in physical clogging associated with silt and biofilm layers, minimizing viscosity effects by maintaining isothermal conditions, minimizing chemical clogging caused by carbonate mineral precipitation associated with algal photosynthesis, and diminished gas clogging associated with trapped air and biogenic gases. This pilot study illustrates the viability of trench infiltration for enhancing surface spreading of managed aquifer recharge to permeable bedrock. Published in 2010 by John Wiley & Sons, Ltd.     

Groundwater composition and recharge origin in the shallow aquifer of the Djerid oases,southern Tunisia: implications of return flow

Abstract

Major ions and stable isotopes in groundwaters of the Plio-Quaternary shallow aquifer of the Djerid oases, southern Tunisia, were investigated to elucidate the origin of groundwater recharge and the mineralization processes. It has been demonstrated that the groundwater composition is mainly controlled by the water–rock interaction, the encroachment of brines from the Chotts and the return flow of irrigation waters. The isotopically depleted groundwater samples suggest that the recharge waters derive from an old palaeoclimatic origin. However, the enriched groundwater samples reflect the presence of evaporated recharge water. Furthermore, the large negative deuterium-excess values indicate the effect of secondary evaporation processes, probably related to the return flow of irrigation waters pumped from the underlying aquifer.

Editor D. Koutsoyiannis; Associate editor E. Custodio

Citation Tarki, M., Dassi, L. and Jedoui, Y., 2012. Groundwater composition and recharge origin in the shallow aquifer of the Djerid oases, southern Tunisia: implications of return flow. Hydrological Sciences Journal, 57 (4), 790–804.     

Soil profile evolution following land‐use change: implications for groundwater quantity and quality

Soil and vadose zone profiles are used as an archive of changes in groundwater recharge and water quality following changes in land use in an area of the Loess Plateau of China. A typical rain‐fed loess‐terrace agriculture region in Hequan, Guyuan, is taken as an example, and multiple tracers (chloride mass balance, stable isotopes, tritium and water chemistry) are used to examine groundwater recharge mechanisms and to evaluate soil water chloride as an archive for recharge rate and water quality. Results show that groundwater recharge beneath natural uncultivated grassland, used as a baseline, is about 94–100 mm year^?1 and that the time it takes for annual precipitation to reach water table through the thick unsaturated zone is from decades to hundreds of years (tritium free). This recharge rate is 2–3 orders of magnitude more than in the other semiarid areas with similar annual rainfall but with deep‐rooted vegetation and relatively high temperature. Most of the water that eventually becomes recharge originally infiltrated in the summer months. The conversion from native grassland to winter wheat has reduced groundwater recharge by 42–50% (50–55 mm year^?1 for recharge), and the conversion from winter wheat to alfalfa resulted in a significant chloride accumulation in the upper soil zone, which terminated deep drainage. The paper also evaluates the time lag between potential recharge and actual recharge to aquifer and between increase in solute concentration in soil moisture and that in the aquifer following land‐use change due to the deep unsaturated zone. Copyright © 2012 John Wiley & Sons, Ltd.     

Applicability of spectral analysis to determine hydraulic diffusivity

This study is to evaluate the applicability of estimating the one-dimensional horizontal hydraulic diffusivity of an unconfined aquifer with time-dependent fluctuation of lateral head and vertical recharge boundaries using observed water level spectra. Different models of boundary condition are imposed to evaluate the statistical significance between the calculated hydraulic diffusivity (ξ^′) with the given hydraulic diffusivity (ξ). The auto-spectra of the water level in observation wells tapping the same aquifer are closely related to those at the disturbed boundaries. For an aquifer with a constant hydraulic diffusivity, the water level fluctuation in the monitoring wells is linearly related to the water level spectra observed at the boundaries. The spectral density function of aquifer hydraulic head varies inversely with specific yield (S _y) and directly with recharge. Given small variation in water level spectra at the disturbed boundaries, the water level fluctuation in the aquifer is affected by the recharge condition and the aquifer spectral density function is sensitive to S _y. Using an iterative technique to estimate ξ from 1400 sets of given parameters, 99% of the ξ^′/ξ values deviated within only one order of magnitude with the model length (L) being equal to 1 km and 10 km. For L equal to 100 m, approximately 82% of the ξ^′/ξ population falls within two orders of magnitude. Therefore, spectral analysis of aquifer hydraulic head response can be used to estimate the hydraulic diffusivity of an unconfined aquifer which is affected by periodic variations in recharge and head at boundaries.     

Radially convergent groundwater flow in sloping terrain

Abstract

The groundwater flow equation governing the elevation (h) of the steady-state phreatic surface in a sloping aquifer fed by constant recharge over a bi-circular sector is rhh′ ? r ² Bh′ + Pr ² ? PR ² = 0, where r is the radial coordinate, P is a constant involving recharge and aquifer properties, and B is the slope of the aquifer—bedrock boundary. The derived flow equation describes radially convergent flow through a sloping aquifer that discharges to a water body of fixed head. One important simplification is that in which the width of the bi-circular sector is constant, and the draining land becomes a rectangular aquifer. The bi-circular sector and rectangular-strip groundwater flow problems are solved in terms of implicit equations. The solutions for the steady-state phreatic surfaces depend on the ratio of recharge to hydraulic conductivity, the slope of the aquifer-bedrock, and the downstream constant-head boundary. Computational examples illustrate the application of the solutions.     

Effects of variations of river stage and hydraulic conductivity on temporal scaling of groundwater levels: numerical simulations

It was found in previous studies that groundwater levels may fluctuate as a temporal fractal. In this study numerical simulations of groundwater level fluctuations in an unconfined aquifer near a river were conducted to investigate the effects of aquifer heterogeneity and river stage variations on the fractal behavior of the water levels, h(t). Groundwater recharge was taken to be a white-noise process. The aquifer heterogeneity was simulated with a second-order stationary field of hydraulic conductivity (K) with an exponential variogram model. The results showed that groundwater levels fluctuate as a temporal fractal in both homogeneous and heterogeneous aquifers as long as K is less than 10 m/d. Most aquifers may indeed act as a fractal filter which takes a random non-fractal recharge inputs and produces a fractal responses of groundwater level fluctuations. A crossover in temporal scaling of h(t) may appear in more permeable aquifers. Fluctuations of the groundwater level in a homogeneous aquifer are dominated by the recharge process when the river stage is constant or by the river stage variations when the river stage varies in highly permeable aquifers. Heterogeneity plays an important role in the temporal scaling of h(t) in more permeable aquifers: the stronger the heterogeneity, the stronger the temporal scaling of h(t).     

The Sustainable Pumping Rate Concept: Lessons from a Case Study in Central Italy

The term sustainable pumping rate (SPR) is defined as the maximum pumping rate that can be maintained indefinitely without mining an aquifer, and is different from the concept of safe yield (SY), which takes into account also aspects related to the much wider concept of sustainability. The assessment of the SPR for the case study of Petrignano d'Assisi, an alluvial aquifer located in Central Italy, shows the need for a reliable estimate of the global water budget of the aquifer, particularly of the recharge under undisturbed conditions; however, the latter is not sufficient, because the SPR is affected also by the geometry of the aquifer, the hydraulic conductivity pattern, the variation of recharge/discharge ratio induced by the abstractions, and so on. All these aspects are analyzed by means of a numerical flow model calibrated both under undisturbed conditions (1974) and under exploitation conditions (1998 to 2004). The steady-state modeling results show that the relation between recharge and abstractions both at local and global scale is a key point in order to estimate a long-term SPR. Moreover, as it could be necessary to overexploit the aquifer for short periods, e.g., during drought episodes, the estimate of SPR must be performed also in transient conditions, in order to take into account the characteristic time of depletion and the successive recovery.     

Estimation of the percentage of annual groundwater recharge with bomb tritium using a cumulative mass balance method

The bomb tritium (³H) distribution patterns in the aquifer beneath an abandoned landfill at the Canadian Forces Base (CFB) Borden, Ontario, and in a sandy aquifer at Whiteshell Nuclear Research Establishments (WNRE) Pinawa, Manitoba, all in Canada, were delineated in great detail. A sampling and monitoring network of multilevel samplers and bundle piezometers were used. The directions of groundwater flow were established, and the boundary between the tritiated and non-tritiated zones of the two aquifers were closely demarcated. Using a cumulative mass balance method, the³H input mass into the aquifers was compared with the³H mass in groundwater storage to estimate the percentages of annual groundwater recharge from 1953 to 1978. Two recharge calculations for theeffective recharge zone and thetotal recharge area of the aquifers as established from the flownet analysis, and the distributions of dissolved geochemical constitutents showed that theeffective recharge zone calculations gave higher values of 30.6 cm/yr for CFB Borden and 20.1 cm/yr for WNRE while thetotal recharge areas gave lower values of 19.1 and 10.1 cm/yr for the Borden and WNRE aquifers respectively. The two recharge values provide possible minimum and maximum recharge estimates for the two study areas.     

Multi-method groundwater recharge estimation at Eshito micro-watershed,Rift Valley Basin in Ethiopia

ABSTRACT

Understanding recharge processes is fundamental to improve sustainable groundwater resource management. Shallow groundwater (SGW) is being developed for multiple purposes in Ethiopia without consideration of monitoring. We established a citizen science-based hydro-meteorological monitoring network, with a focus on SGW recharge estimation, in Eshito micro-watershed, Ethiopia. Citizen scientists collected rainfall, groundwater-level and stream water-level data. We characterized the shallow aquifer using pumping tests. The data were used to estimate SGW recharge using three methods: chloride mass balance, water-level fluctuation (WLF) and baseflow separation. Approximately 20% and 35% of annual rainfall amount contributes to recharge based on the chloride mass balance and WLF results, respectively. Baseflow separation showed recharge values for the watershed vary from 38% to 28% of annual rainfall at the upstream and downstream gauging stations, respectively. This study shows that the recharge in previously unmonitored micro-watersheds can be studied if citizens are involved in data generation.     

Sensitivity analysis in groundwater vulnerability assessment based on GIS in the Mahdia-Ksour Essaf aquifer,Tunisia: a validation study

Abstract

The assessment of groundwater vulnerability to pollution has proved to be an effective tool for water resource management, especially in arid and semi-arid regions like Mahdia and Ksour Essaf. The main objective of this study is to assess the aquifer vulnerability by applying the DRASTIC method as well as using sensitivity analysis to evaluate the effect of each DRASTIC parameter on the final vulnerability map. An additional objective is to demonstrate the role of the GIS techniques in the vulnerability assessment. The DRASTIC method assigns a high vulnerability to the coast of the Mahdia-Ksour Essaf. The lowest values are observed in the southern part of the study area. A sensitivity analysis applied in this study suggests that net recharge, aquifer media and depth of groundwater are the key factors determining vulnerability. The model is validated with groundwater quality data and the results have shown strong relationships between modified DRASTIC Vulnerability Index and nitrate and chloride concentrations.

Citation Saidi, S., Bouri, S. & Ben Dhia, H. (2011) Sensitivity analysis in groundwater vulnerability assessment based on GIS in the Mahdia-Ksour Essaf aquifer, Tunisia: a validation study. Hydrol. Sci. J. 56(2), 288–304.     

Impact of Turbidity on TCE and Degradation Products in Ground Water

Elevated particulate concentrations in ground water samples can bias contaminant concentration data. This has been particularly problematic for metal analyses where artificially increased turbidity levels can affect metals concentrations and confound interpretation of the data. However, few studies have been conducted to determine the impact of particulates on trichloroethylene (TCE), cis-dichloroethylene (c-DCE), and vinyl chloride concentrations.
Laboratory batch studies and field investigations were conducted to evaluate the effects of suspended solids on VOC concentrations in ground water samples analyzed by purge-and-trap gas chromatography. Three different solids were used to assess the effects of suspended particulates. The solids were aquifer material from a field site in North Carolina and two reference clay minerals (kaolinite and Namontimorillonite). During the laboratory portion of this study, the solids were used to determine effects on TCE concentrations under controlled laboratory conditions.
The same solids were used in a field study to compare the laboratory results with field results. Solids were added to the sample vials prior it) sample collection to intentionally increase turbidity levels in the water samples. Results of the study indicate essentially no decrease in TCE, c-DCH, or vinyl chloride concentrations due to increased turbidity levels.     

Transient Recharge Estimability Through Field‐Scale Groundwater Model Calibration

The estimation of recharge through groundwater model calibration is hampered by the nonuniqueness of recharge and aquifer parameter values. It has been shown recently that the estimability of spatially distributed recharge through calibration of steady‐state models for practical situations (i.e., real‐world, field‐scale aquifer settings) is limited by the need for excessive amounts of hydraulic‐parameter and groundwater‐level data. However, the extent to which temporal recharge variability can be informed through transient model calibration, which involves larger water‐level datasets, but requires the additional consideration of storage parameters, is presently unknown for practical situations. In this study, time‐varying recharge estimates, inferred through calibration of a field‐scale highly parameterized groundwater model, are systematically investigated subject to changes in (1) the degree to which hydraulic parameters including hydraulic conductivity (K) and specific yield (S_y) are constrained, (2) the number of water‐level calibration targets, and (3) the temporal resolution (up to monthly time steps) at which recharge is estimated. The analysis involves the use of a synthetic reality (a reference model) based on a groundwater model of Uley South Basin, South Australia. Identifiability statistics are used to evaluate the ability of recharge and hydraulic parameters to be estimated uniquely. Results show that reasonable estimates of monthly recharge (<30% recharge root‐mean‐squared error) require a considerable amount of transient water‐level data, and that the spatial distribution of K is known. Joint estimation of recharge, S_y and K, however, precludes reasonable inference of recharge and hydraulic parameter values. We conclude that the estimation of temporal recharge variability through calibration may be impractical for real‐world settings.     

Estimating groundwater residence time and recharge patterns in a saline coastal aquifer

A detailed study using environmental tracers such as chloride (Cl^?) and tritium (³H), deuterium (²H) and oxygen (¹⁸O) isotopes was performed in an alluvial coastal aquifer in two contrasting environments (urban and agricultural). These environmental tracers combined with a high‐resolution multi‐level sampling approach were used to estimate groundwater residence time and recharge patterns and to validate the hydrogeochemical conceptual model already proposed in previous studies. δ¹⁸O and δ²H combined with Cl^? data proved that the hypersaline groundwater present in the deepest part of the aquifer was sourced from the underlying hypersaline aquitard via an upward flux. Both chemical and isotopic data were employed to calibrate a density‐dependent numerical model based on SEAWAT 4.0, where ³H and Cl^? were helped quantifying solutes transport within the modelled aquifer. Model results highlighted the differences on estimated recharge in the two contrasting environments, with the urban one exhibiting concentrated recharge because of preferential infiltration associated to the storm water drains network, while scarce local recharge characterized the agriculture setting. In the urban field site, is still possible to recognize at 9 m b.g.l. the input of the atmospheric anthropogenic ³H generated by testing of thermonuclear weapons, while in the agricultural field site, the ³H peak has been washed out at 6 m b.g.l. because the groundwater circulation is restricted only to the upper fresh part of the aquifer, drained by the reclamation system. The presented approach that combined high‐resolution field monitoring, environmental tracers and numerical modelling, resulted effective in validating the conceptual model of the aquifer salinization. Copyright © 2016 John Wiley & Sons, Ltd.     

Comparative performance of soil water balance models in computing semi-arid aquifer recharge

Abstract

Estimating groundwater recharge is essential to ensure the sustainable use of groundwater resources, particularly in arid and semi-arid regions. Soil water balances have been frequently advocated as valuable tools to estimate groundwater recharge. This article compares the performance of three soil water balance models (Hydrobal, Visual Balan v2.0 and Thornthwaite) in the Ventós-Castellar aquifer, Spain. The models were used to simulate wet and dry years. Recharge estimates were transformed into water table fluctuations by means of a lumped groundwater model. These, in turn, were calibrated against piezometric data. Overall, the Hydrobal model shows the best fit between observed and calculated levels (r² = 0.84), highlighting the role of soil moisture and vegetation in recharge processes.

Editor D. Koutsoyiannis; Associate editor X. Chen

Citation Touhami, I., et al., 2014. Comparative performance of soil water balance models in computing semi-arid aquifer recharge. Hydrological Sciences Journal, 59 (1), 193–203.     

Groundwater recharge amidst focused stormwater infiltration

Distributed, infiltration‐based approaches to stormwater management are being implemented to mitigate effects of urban development on water resources. One of the goals of this type of storm water management, sometimes called low impact development or green infrastructure, is to maintain groundwater recharge and stream base flow at predevelopment levels. However, the connection between infiltration‐based stormwater management and groundwater recharge is not straightforward. Water infiltrated through stormwater facilities may be stored in soil moisture, taken up by evapotranspiration or contribute to recharge and eventually base flow. This study focused on a 1.1 km² suburban, low impact development watershed in Clarksburg, Maryland, USA, that was urbanized and contained 73 infiltration‐based stormwater facilities. Continuous water table measurements were used to quantify the movement of infiltrated stormwater. Time series analyses were performed on hydrographs of 7 wells, and the episodic master recession method was used. Persistence in water levels, as measured by autocorrelation function, was found to be positively related to depth to water. Storm properties (precipitation rate and duration) and well location (proximity to the nearest stream) were significant in driving episodic recharge to precipitation ratios. The well that had the highest recharge to precipitation ratios and water table rises of up to 1.5 m in response to storm events was located furthest from the stream and down gradient of stormwater infiltration locations. This work may be considered in evaluating the effects of planned watershed‐scale infiltration‐based stormwater management on groundwater flow systems.     

Groundwater recharge and flow in a small mountain catchment in northern Ethiopia

Abstract

The hydrodynamic behaviour of a sloped phreatic aquifer in the Tigray Highlands in northern Ethiopia is described. The aquifer is situated in the soils of a plateau on top of a basalt sequence and lies on steep slopes; the latter lead to hydraulic gradients that can cause high discharge fluxes. Distinct wet and dry seasons characterize the climate of the Tigray Highlands and recharge is absent during the dry season. Because of the fertile vertisols that have developed, the plateau is heavily cultivated and thus has great local economic, and hence social, importance. Water for land irrigation is almost exclusively delivered by rainfall, which is largely restricted to the period June—September. During the dry season, the water table drops dramatically and the aquifer drains nearly completely, under the strong gravity-driven, sustained discharges. This study strives to give insights into recharge and discharge mechanisms of the aquifer, in order to improve the effectiveness of the implemented water conservation measures.     

Variation of hydrogeochemical characteristics of water in surface flows,shallow wells,and boreholes in the coastal city of Douala (Cameroon)

ABSTRACT

Groundwater is used by 3?million inhabitants in the coastal urban city of Douala, Cameroon, but comprehensive data are too sparse for it to be managed in a sustainable manner. Hence this study aimed to (1) assess the potability of the groundwater; (2) evaluate the spatial variation of groundwater composition; and (3) assess the interaction and recharge mechanisms of different water bodies. Hydrogeochemical tools and methods revealed the following results in the Wouri and Nkappa formations of the Douala basin, which is beneath Douala city: 30% of water samples from hand-dug wells in the shallow Pleistocene alluvium aquifer were saline and highly mineralized. However, water from boreholes in the deeper (49–92 m depth) Palaeocene aquifer was saline-free, less mineralized and potable. Water in the shallow aquifer (0.5–22 m depth) was of Na⁺-K⁺-Cl^?-NO₃^? type and not potable due to point source pollution, whereas Ca⁺-HCO₃^? unpolluted water dominates in the deeper aquifer. Water in the deep and shallow aquifers indicates the results of preferential flow pass and evaporative recharge, respectively. Possible hydrogeochemical processes include point source pollution, reverse ion exchange, remote recharge areas and mixing of waters with different chemical signatures.

EDITOR D. Koutsoyiannis ASSOCIATE EDITOR M.D. Fidelibus