Clay fine fissuring monitoring using miniature geo-electrical resistivity arrays

This article describes a miniaturised electrical imaging (resistivity tomography) technique to map the cracking pattern of a clay model. The clay used was taken from a scaled flood embankment built to study the fine fissuring due to desiccation and breaching process in flooding conditions. The potential of using a miniature array of electrodes to follow the evolution of the vertical cracks and number them during the drying process was explored. The imaging technique generated two-dimensional contoured plots of the resistivity distribution within the model before and at different stages of the desiccation process. The change in resistivity associated with the widening of the cracks were monitored as a function of time. Experiments were also carried out using a selected conductive gel to slow down the transport process into the cracks to improve the scanning capabilities of the equipment. The main vertical clay fissuring network was obtained after inversion of the experimental resistivity measurements and validated by direct observations.     

Impact of sandy beach recovery on solute transport in coastal unconfined aquifers

Hydrogeology Journal - Beach recovery describes the processes by which there is a natural restoration of beach material and coastal morphology following storm events, and these processes are common...     

Measurement of vertical environmental-head profiles in unconfined sand aquifers using a multi-channel manometer board

The vertical hydraulic gradient in an unconfined sand aquifer frequently exceeds the horizontal gradient when measured close to a surface discharge boundary, around an abstraction bore or close to a saline interface. The vertical gradient can be efficiently investigated using a multi-channel manometer board connected to a bundled piezometer installed using hollow-stem-auger-drilling techniques. The design and construction of a manometer board is described. Manometer board data from a site in a 30-m-deep fresh-water aquifer adjacent to a pumping bore are described. The bore is screened in the lower part of the aquifer and the impacts of partial penetration are evident on the manometer board results. Data are then described from the centre of a coastal sand spit where fresh water overlies saline water. Fluid electrical conductivity measurements from the individual tubes of the bundled piezometer are used to derive fluid density at each of the bundled piezometer ports. The average fluid density of water above the mini-piezometer port at each depth is then calculated using Gaussian quadrature. The environmental-water head profiles are compared to point-water and fresh-water head profiles and the implications of the incorrect use of point-water or fresh-water head vertical profiles are clearly demonstrated.     

Cation exchange processes and human activities in unconfined aquifers

During the 1999–2002 water years, a hydrogeological research project was carried out on the unconfined aquifer of Trifilia in the Peloponnese. Seawater intrusion due to overpumping, and intensive use of fertilizers caused the groundwater quality degradation that is a typical case for the coastal aquifers in Greece. Isopiezometric maps along with ion distribution balances, ion distribution maps and factor analysis indicate the existence of three zones of groundwater quality. In the first zone of saline water, a cation exchange process between the Ca²⁺ of sediments and the Na⁺ of groundwater contribute to the formation of the water type Na⁺-Ca²⁺-Cl^–. In the second zone, which is considered as a transition zone, dominate the Ca²⁺-Na⁺-HCO₃^–-Cl^– water type. In the third zone of Ca²⁺-HCO₃^–-SO₄^2– water type, relationships among Ca²⁺, SO₄^2–, NO₃^– and NH₄ can be attributed to the dissociation of ammonium nitrate and sulfate fertilizers on one hand, and Ca²⁺ derivation from cation exchange processes between water, rocks and clay minerals, such as smectite and illite, on the other.     

Numerical values of shape factors for field permeability tests in unconfined aquifers

Acta Geotechnica - Field permeability tests are used to evaluate the local hydraulic conductivity. Their interpretation requires knowing the value of a shape factor, c. Regular values for shape...     

Proof of multiplicity of solutions for groundwater seepage in recharged heterogeneous unconfined aquifers

A new analytical proof is presented for steady‐state seepage in recharged heterogeneous unconfined aquifers. The paper also presents a detailed procedure and important rules for performing correctly numerical studies of unsaturated seepage. Once a numerical solution is calibrated with field data, using a set of spatially distributed values for hydraulic conductivity K and effective infiltration EI, any new numerical analysis with a set of αK and αEI values, where α is a constant, yields an equally good calibration. However, if the effective porosities of each layer are unchanged, the groundwater velocities are multiplied by α, whereas the travel times are divided by α, which may help to select α in order to match known travel time data. This is a clear example of multiple solutions to an inverse problem. The paper underlines the role and the need to finely mesh unsaturated zones and also contacts between layers to reach the asymptotic convergence range, as it was carried out to verify the proof and as it should be completed to study any seepage problem. A few consequences of the new analytical proof and the rigorous procedure are shown with examples. Copyright © 2016 John Wiley & Sons, Ltd.     

3D-electrical resistivity tomography monitoring of salt transport in homogeneous and layered soil samples

Monitoring transport of dissolved substances in soil deposits is particularly relevant where safety is concerned, as in the case of geo-environmental barriers. Geophysical methods are very appealing, since they cover a wide domain, localising possible preferential flow paths and providing reliable links between geophysical quantities and hydrological variables. This paper describes a 3D laboratory application of electrical resistivity tomography (ERT) used to monitor solute transport processes. Dissolution and transport tests on both homogeneous and heterogeneous samples were conducted in an instrumented oedometer cell. ERT was used to create maps of electrical conductivity of the monitored domain at different time intervals and to estimate concentration variations within the interstitial fluid. Comparisons with finite element simulations of the transport processes were performed to check the consistency of the results. Tests confirmed that the technique can monitor salt transport, infer the hydro-chemical behaviour of heterogeneous geomaterials and evaluate the performances of clay barriers.     

Comparative analysis between analytical approximations and numerical solutions describing recession flow in unconfined hillslope aquifers

Recession flow of aquifers from a hillslope can be described by the non-linear Boussinesq equation. Under strong assumptions and for specific conceptual formulations, different authors derived analytical approximations or linearized versions to this partial differential equation. A comparative analysis between some analytical approximations of the Boussinesq equation and the numerical solution of the recession flow of an unconfined homogeneous aquifer (horizontal, inclined and concave aquifer floor) was carried out. The objective was to define the range where the analytical solutions approximate the numerical solution. The latter was considered in this study as the reference method, because it requires fewer assumptions. From the considered analytical approximations, exponential decay relationships were found to be mainly valid for fine domain materials when horizontal, mild slopes (less than 2%) and concave aquifer floors were considered, but failed to reproduce coarse aquifer numerical model outflows, in contrast to the quadratic decay relationship, which better reproduce outflows in such domains. On the basis of the comparative analysis, it has been found that recession flows obtained with the considered analytical approximations yield similar values only for certain ranges of aquifer properties and geometries.

---

Résumé L’écoulement de décrue dans les aquifères de versant peut être décrit par l’équation non-linéaire de Boussinesq. En considérant d’importantes hypothèses et des formulations conceptuelles spécifiques, différents auteurs ont dérivé des approximations analytiques ou des versions linéarisées de cette équation partielle différentielle. Une analyse comparative entre des approximations analytiques de l’équation de Boussinesq et une solution numérique de la décrue d’un aquifère homogène libre (limite inférieure de l’aquifère horizontale, inclinée et concave) a été effectuée. L’objectif était de définir la plage pour laquelle les solutions analytiques approchent la solution numérique. Cette dernière a été considérée dans cette étude comme la méthode de référence car elle demande moins d’hypothèses. Parmi les approximations analytiques considérées, les relations de décroissance exponentielle se sont avérées être principalement valides pour des matériaux fins lorsque des couches aquifères horizontales, concaves et de pentes douces (moins de 2%) étaient considérées; en revanche, celles-ci n’ont pas réussi à reproduire les écoulements sortants d’un modèle numérique pour un aquifère à matériaux grossiers, contrairement à la relation de décroissance quadratique qui reproduit mieux les écoulements sortants dans de tels milieux. Sur la base de l’analyse comparative, il s’est avéré que les écoulements en décrue obtenus avec les approximations analytiques considérées présentent des valeurs similaires seulement pour certains critères géométriques et propriétés de l’aquifère.

---

Resumen Puede describirse el flujo de recesión de acuíferos de una ladera mediante la ecuación no-lineal de Boussinesq. Bajo fuertes supuestos y para formulaciones conceptuales específicas, distintos autores han derivado aproximaciones analíticas o versiones lineales para esta ecuación parcial diferencial. Se realizó un análisis comparativo entre algunas aproximaciones analíticas de la ecuación de Boussinesq y la solución numérica del flujo de recesión de un acuífero homogéneo no confinado (piso de acuífero cóncavo, inclinado, y horizontal). El objetivo consistió en definir el rango donde las soluciones analíticas aproximan la solución numérica. En este estudio se consideró la solución numérica como el método de referencia debido a que requiere menores supuestos. De las aproximaciones analíticas consideradas se encontró que las relaciones de desintegración exponencial eran principalmente válidas para materiales de ámbitos finos donde se consideraron pisos de acuíferos horizontal, pendientes medianas (menores a 2%) y cóncavos pero fallaron en reproducir flujos de salida de modelos numéricos de ámbito grueso, en contraste con la relación de desintegración cuadrática la cual reproduce mejor los flujos de salida en tales ámbitos. En base al análisis comparativo se encontró que los flujos de recesión obtenidos con las aproximaciones analíticas consideradas aportaron valores similares solo para ciertos rangos de propiedades y geometrías de acuíferos.

     

Evaluation and mapping of Dead Sea coastal aquifers salinity using Transient Electromagnetic (TEM) resistivity measurements

Evaporite karst has intensively developed recently along the Dead Sea (DS) coastal area in Israel and Jordan. It takes place in very saline groundwater dissolving buried salt layers, causing collapse of the surface. In this paper, groundwater salinity throughout the DS coastal area is investigated using the Transient Electromagnetic (TEM) method. Twenty-eight TEM soundings along the DS coastal area were carried out close to observation boreholes to calibrate resistivity–salinity relationships. Groundwater electrical conductivity was measured in these boreholes, and its salinity was analyzed at the laboratory by the Geological Survey of Israel (GSI). Quantitative relationships between bulk resistivity (ρ_x), water resistivity (ρ_w) and chloride concentration (C_cl) were derived in the resistivity range less than 1.0 Ω·m that enabled to evaluate the salinity of the aquifer in in situ conditions. Average values of the effective porosity of sandy sediments, φ_e = 0.32, and of silty ones, φ_e = 0.44, were used to generate the corresponding Archie equations. The study has shown that a DS aquifer with bulk resistivity in the range of 0.55–1.0 Ω·m contains in pores brine with 50–110 g_chloride/l of (22–50% of that in saturated conditions, respectively), i.e. it keeps the potential to dissolve up to 114–174 g/l of salt.     

Towards a nonperturbation transport theory in hetereogeneous aquifers

A large body of existing theories of flow and contaminant transport in aquifers ignore the presence of recharge, eliminate the boundary conditions, neglect transient conditions in groundwater flow, conceive hydraulic gradients as linear, and require parameter variability to be stationary and Gaussian. The most outstanding and difficult to justify assumption is the subjective small size of the stochastic terms (i.e., small perturbation methods), which usually is forced by considering the logarithm of the hydraulic conductivity. Several problems in flow and contaminant subsurface hydrology, such as the enhanced dispersion parameters with plume size or time after injection, remain to be observed in the light of a stochastic theory that allows a more realistic consideration of physical and hydrologic properties. In this article, an attempt is made to reformulate a contaminant transport equation (the variable dispersion equation, VDE) with transport parameters in terms of regional hydrologic and aquifer hydraulic properties, such as recharge rate, spatially random transmissivity, hydraulic gradient, aquifer thickness, and soil porosity. Subsequently, a general analytic procedure, the method of decomposition, is used to derive a solution to the VDE. This procedure does not require small perturbation, logarithmic transformations, or specific probability law assumptions. Comparison tests with existing theoretical and field results are given. The tests illustrate the enhanced dispersion and shifting concentration effects produced by the variable dispersion equation. Finally a generalization of the method to nonstationary dispersion in three-dimensional domains is proposed.     

Pollutant transport in a shallow unconfined aquifer in Perry,Ohio

The groundwater system in northern Perry Township in Lake County, Ohio, is a shallow, unconfined aquifer consisting of periglacial lake beach deposits and less permeable lacustrine plain deposits. Groundwater flow is generally toward Lake Erie from south to north and is controlled by the top of the Ashtabula Till, but strong, local variations are caused by northward flowing streams During the study period, water levels in most wells exhibited a seasonal fluctuation of less than 0 3 m from their mean values. The areal distributions of chloride and nitrate concentrations indicate that road salt runoff easily infiltrates the aquifer and that nitrate may be sourced from fertilizer application. Ground-water flow and solute transport models indicate that in excess of 27 years are required to obtain chemical steady-state under hydrologic steady-state conditions. The simulations also demonstrate that nitrate loading must occur in more than one cultivated field in order to obtain the observed wide-spread nitrate distribution.     

Estimation of hydraulic characteristics from electrical resistivity data in coastal aquifers of southern India

The applicability of geophysical techniques has been examined for evaluating aquifer properties like transmissivity and hydraulic conductivity of coastal aquifers, Tuticorin, Tamil Nadu. The pumping test data of 10 wells are interpreted by using forward modelling to obtain the aquifer characteristics in the study area. The available vertical electrical soundings (VES) data in the vicinity of the sites of pumping test have been interpreted; and true resistivity and thickness are determined at each site in the study area. Empirical relationships are established for estimating the hydraulic parameters from the electrical data.     

Solute transport characteristics and influencing factors in a coastal unconfined aquifer under tidal action identified by image monitoring in a laboratory experimental setup

Hydrogeology Journal - Studies of solute transport in unconfined coastal aquifers typically neglect oceanic oscillations and assume a static seaward boundary condition defined by the mean annual...     

The transport of U- and Th-series nuclides in sandy confined aquifers

Abundances of ²³⁸U, ²³⁴U, ²³²Th, ²²⁶Ra, ²²⁸Ra, ²²⁴Ra, and ²²²Rn were measured in groundwaters of the Ojo Alamo aquifer in northwest New Mexico. This is an arid area with annual precipitation of ∼22 cm. The purpose was to investigate the transport of U-Th series nuclides and their daughter products in an old, slow-moving groundwater mass as a means of understanding water-rock interactions and to compare the results with a temperate zone aquifer. It was found that ²³²Th is approximately at saturation and supports the view of Tricca et al. (2001) that Th is precipitated irreversibly upon weathering, leaving surface coatings of ²³²Th and ²³⁰Th on aquifer grains. Uranium in the aquifer waters has very high [²³⁴U/²³⁸U] ∼ 9 and low ²³⁸U concentrations. These levels can be explained by low weathering rates in the aquifer (w_238U ∼ 2 × 10⁻¹⁸ to 2 × 10⁻¹⁷s⁻¹) using a continuous flow, water-rock interaction model. The Ra isotopes are roughly in secular equilibrium despite their very different mean lifetimes. The ²²²Rn and ²²⁸Ra isotopes in the aquifer correspond to ∼10% of the net production rate of the bulk rock. This is interpreted to reflect an earlier formed irreversible surface coating of Th that provides Ra and Rn to the aquifer waters. The surface waters that appear to be feeding the aquifer have low [²³⁴U/²³⁸U] and high ²³⁸U concentrations. The flow model shows that it is not possible to obtain the high [²³⁴U/²³⁸U] and low [²³⁸U] values in the aquifer from a source like the present vadose zone input. It follows that the old aquifer waters studied cannot be fed by the present vadose zone input unless they are greatly diluted with waters with very low U concentrations. If the present sampling of vadose zone sources is representative of the present input, then this requires that there was a major change in water input with much larger rainfall some several thousand years ago. This may represent a climatic change in the Southwest.     

Groundwater monitoring for cement kiln dust disposal units in karst aquifers

Well-developed karst aquifers tend to be heterogeneous and consist of variable porosities. Groundwater monitoring and the associated data interpretations in such aquifers are often more complicated than porous medium aquifers. Collection of representative data in karst aquifers often requires monitoring at appropriately located wells and/or springs that are proven to connect to the groundwater system. Water samples are to be collected under different flow conditions, including base flow, high-flow, and low-flow. The sampling frequencies may vary from several months for base flows to minutes in response to recharge events. The groundwater monitoring program presented in this paper is for a cement kiln dust mono-fill site in a karst area of southern Indiana. Following dye tracing and extensive geophysical investigations, one spring was selected as a monitoring location. A second spring should be used as a monitoring location when the last cell of the mono-fill begins receiving the wastes. The paper discusses results from the first spring, at which nine background sampling events were completed to evaluate the natural variations of the water quality. Based on the background data, a statistical evaluation plan was developed for 12 water-quality parameters to determine the integrity of the landfill. The statistical power of the statistical analyses was evaluated by Monte Carlo simulations.     

Effect of tidal fluctuations on transient dispersion of simulated contaminant concentrations in coastal aquifers

Variable-density groundwater models require extensive computational resources, particularly for simulations representing short-term hydrologic variability such as tidal fluctuations. Saltwater-intrusion models usually neglect tidal fluctuations and this may introduce errors in simulated concentrations. The effects of tides on simulated concentrations in a coastal aquifer were assessed. Three analyses are reported: in the first, simulations with and without tides were compared for three different dispersivity values. Tides do not significantly affect the transfer of a hypothetical contaminant into the ocean; however, the concentration difference between tidal and non-tidal simulations could be as much as 15%. In the second analysis, the dispersivity value for the model without tides was increased in a zone near the ocean boundary. By slightly increasing dispersivity in this zone, the maximum concentration difference between the simulations with and without tides was reduced to as low as 7%. In the last analysis, an apparent dispersivity value was calculated for each model cell using the simulated velocity variations from the model with tides. Use of apparent dispersivity values in models with a constant ocean boundary seems to provide a reasonable approach for approximating tidal effects in simulations where explicit representation of tidal fluctuations is not feasible.     

Deep resistivity sounding studies for probing deep fresh aquifers in the coastal area of Orissa, India

Exploration and exploitation of groundwater in sedimentary areas are reasonably simple. However, the problem of salinity in coastal areas makes the job very difficult, especially when the freshwater aquifers are not extensive and are entrapped between saline aquifers. States along the eastern coast of India, particularly Orissa with respect to the Mahanadi basin, have acute problems with groundwater salinity. It has been possible to locate horizons of fresh groundwater entrapped between deep saline aquifers in the southwestern part of Mahanadi delta, with the help of deep resistivity soundings along the Delang-Puri profile. This finding has been validated through boreholes and checked with electrical logs of this region. Three freshwater aquifers have been detected: one at shallow depth between 20 and 60?m, the second in the depth range of 90??60?m, and the third in the fractured/weathered basement. The second freshwater aquifer has the most potential; it has a thickness range of 20??0?m and it could be exploited to overcome problematic salinity issues. In general, the depth to basement is variable and it increases seaward.     

Time-lapse resistivity investigations for imaging saltwater transport in glaciofluvial deposits

Five intersecting resistivity sections have been measured in glaciofluvial deposits hosting an aquifer of regional importance situated along a heavy traffic highway in Sweden. The winter salt spreading has caused a regular salinity increase through the years. For imaging the transport of saltwater in the aquifer, the sections were measured exactly in the same location before and after winter, and interpreted using a time-lapse inverse procedure. Some auger drilling and RCPT data were available for correlation. After winter, the resistivity had generally decreased under the water table and increased above it. The decrease in resistivity in the saturated zone is interpreted as a plume of more saline groundwater created by de-icing salt from the road. The increase in the upper layer can be explained by changes in temperature and soil moisture. The study shows that time-lapse resistivity investigations has potential for imaging hydraulic pathways in complex hydrogeological environments.     

Groundwater response to serial stream stage fluctuations in shallow unconfined alluvial aquifers along a regulated stream (West Virginia,USA)

Groundwater response to stream stage fluctuations was studied in two unconfined alluvial aquifers using a year-long time series of stream stages from two pools along a regulated stream in West Virginia, USA. The purpose was to analyze spatial and temporal variations in groundwater/surface-water interaction and to estimate induced infiltration rate and cumulative bank storage during an annual cycle of stream stage fluctuation. A convolution-integral method was used to simulate aquifer head at different distances from the stream caused by stream stage fluctuations and to estimate fluxes across the stream–aquifer boundary. Aquifer diffusivities were estimated by wiggle-matching time and amplitude of modeled response to multiple observed storm events. The peak lag time between observed stream and aquifer stage peaks ranged between 14 and 95 hour. Transient modeled diffusivity ranged from 1,000 to 7,500 m²/day and deviated from the measured and calculated single-peak stage-ratio diffusivity by 14–82 %. Stream stage fluctuation displayed more primary control over groundwater levels than recharge, especially during high-flow periods. Dam operations locally altered groundwater flow paths and velocity. The aquifer is more prone to surface-water control in the upper reaches of the pools where stream stage fluctuations are more pronounced than in the lower reaches. This method could be a useful tool for quick assessment of induced infiltration rate and bank storage related to contamination investigations or well-field management.