Truncated plurigaussian simulations to characterize aquifer heterogeneity

Integrating geological concepts, such as relative positions and proportions of the different lithofacies, is of highest importance in order to render realistic geological patterns. The truncated plurigaussian simulation method provides a way of using both local and conceptual geological information to infer the distributions of the facies and then those of hydraulic parameters. The method ( Le Loc'h and Galli 1994 ) is based on the idea of truncating at least two underlying multi-Gaussian simulations in order to create maps of categorical variable. In this article, we show how this technique can be used to assess contaminant migration in highly heterogeneous media. We illustrate its application on the biggest contaminated site of Switzerland. It consists of a contaminant plume located in the lower fresh water Molasse on the western Swiss Plateau. The highly heterogeneous character of this formation calls for efficient stochastic methods in order to characterize transport processes.     

Analytical approximation to characterize the performance of in situ aquifer bioremediation

The performance of in situ bioremediation to remove organic contaminants from contaminated aquifers depends on the physical and biochemical parameters. We characterize the performance by the contaminant removal rate and the region where biodegradation occurs, the biologically active zone (BAZ). The numerical fronts obtained by one-dimensional in situ bioremediation modeling reveal a traveling wave behavior: fronts of microbial mass, organic contaminant and electron acceptor move with a constant velocity and constant front shape through the domain. Hence, only one front shape and a linear relation between the front position and time is found for each of the three compounds. We derive analytical approximations for the traveling wave front shape and front position that agree perfectly with the traveling wave behavior resulting from the bioremediation model. Using these analytical approximations, we determine the contaminant removal rate and the BAZ. Furthermore, we assess the influence of the physical and biochemical parameters on the performance of the in situ bioremediation technique.     

The use of kelp sieve tube sap metal composition to characterize urban runoff in southern California coastal waters

This study introduces an innovative method for biomonitoring using giant kelp (Macrocystis pyrifera) sieve tube sap (STS) metal concentrations as an indication of pollution influence. STS was sampled from fronds collected from 10 southern California locations, including two reference sites on Santa Catalina Island. Using ICP-MS methodology, STS concentrations of 17 different metals were measured (n = 495). Several metals associated with pollution showed the highest STS concentrations and most seasonal variation from populations inside the Port of Los Angeles/Long Beach. Lowest concentrations were measured at less-urbanized areas: Santa Catalina Island and Malibu. Some metals showed a spatial gradient in STS metal concentration with increasing distance from point sources (i.e. Los Angeles River). Cluster analyses indicate that polluted seawater may affect kelp uptake of metals essential for cellular function. Results show that this method can be useful in describing bioavailable metal pollution with implications for accumulation within an important ecosystem.     

Direct and indirect methods for groundwater investigations: A case-study of MRS and VES in the southern part of Sweden

Magnetic Resonance Sounding (MRS) has been successfully tested for detecting groundwater in two areas in southern Sweden. Measurements of Schlumberger VES have been conducted in the same place as the MRS and the results are generally consistent. Low resistivity layers interpreted as clay are sometimes identified close to the surface. The MRS result in site 2 is a good example of signals penetrating through the clay and deeper aquifer still being detected. The MRS data suggest aquifers that are not only hosted in soft sediment materials (moraine, sand, and mixed materials), but hosted in basement rocks. Based on the MRS and borehole pumping test data, the results agree with yield, average water content and subsurface geological data.     

Optimization of groundwater abstraction from a coastal aquifer

ABSTRACT

An analytical mathematical model, based on Jacobian elliptic functions, has been used to identify feasible wellfield locations and pumping rates for large-scale abstraction from an unconfined coastal aquifer. The choice of optimum feasible wellfield strategy has been made using a simple economic model which calculates the cost of the pipelines required to transport the abstracted groundwater to a large coastal city which forms the demand centre. Results indicate that the cheapest wellfield design would be a single large wellfield. However, because of the need to maintain at least a minimum supply in the city until a new surface water source is developed, a better solution may well be to develop two smaller wellfields pumping a greater total abstraction.     

Mathematical modelling and parameter identification for a coastal aquifer

The identification of parameters is considered for an unconfined coastal underground aquifer which is subject to the action of the sea tide on its boundary. First, by introducing a suitable mathematical model, it is possible to separate the problem of determining the permeability from that of determining the specific yield. Then the identification procedures employed are described in detail and some computational results are presented.     

AMT 、TEM 、VES地层响应特征模拟分析及其联合反演探讨

通过正演模拟方式分析音频大地电磁法(AMT)、中心回线瞬变电磁法(TEM)、直流电测深法(VES)对地层的响应特征,结果表明:TEM、AMT对地层的响应特征类似,即对低阻层反应灵敏,而对高阻层的分辩能力较差;VES对高阻层及低阻层的分辩能力相当,相对于TEM、AMT,VES对低阻层的分辩能力差,而对高阻层的分辨能力较高.基于AMT或TEM与VES对地层的响应具有互补的特点,探讨了对AMT或TEM与VES数据进行的联合反演.结果表明,AMT或TEM与VES进行联合反演,相比单一方法进行反演,结果能得到明显的改善.     

Modification of the GALDIT framework using statistical and entropy models to assess coastal aquifer vulnerability

Recently, groundwater vulnerability assessment of coastal aquifers using the GALDIT framework has been widely used to investigate the process of groundwater contamination. This study proposes multi-attribute decision-making (MADM) entropy and Wilcoxon non-parametric statistical test methods to improve the vulnerability index of coastal aquifers. The rates and weights of this framework were modified using Wilcoxon non-parametric and entropy methods, respectively, and a combined framework of GALDIT-entropy, Wilcoxon-GALDIT, and Wilcoxon-entropy was obtained. Pearson correlation coefficients between the mentioned vulnerability indices and total-dissolved solids (TDS) of 0.51, 0.66 and 0.75, respectively, were obtained. According to the results, the Wilcoxon-entropy index had the highest correlation with TDS. Generally, it can be concluded that the proposed frameworks provide a more accurate estimation of vulnerability distribution in coastal aquifers.     

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.     

Groundwater response to tidal fluctuations in a leaky confined coastal aquifer with a finite length

The study on the hydraulic properties of coastal aquifers has significant implications both in hydrological sciences and environmental engineering. Although many analytical solutions are available, most of them are based on the same basic assumption that assumes aquifers extend landward semi‐infinitely, which does not necessarily reflect the reality. In this study, the general solutions for a leaky confined coastal aquifer have been developed that consider both finitely landward constant‐head and no‐flow boundaries. The newly developed solutions were then used to examine theoretically the joint effects of leakage and aquifer length on hydraulic head fluctuations within the leaky confined aquifer, and the validity of using the simplified solution, which assumes the aquifer is semi‐infinite. The results illustrated that the use of the simplified solution may cause significant errors, depending on joint effects of leakage and aquifer length. A dimensionless characteristic parameter was then proposed as an index for judging the applicability of the simplified solution. In addition, practical application of the general solution for the constant‐head inland boundary was used to characterize the hydraulic properties of a leaky confined aquifer using the data collected from a field site at the Seine River estuary, France, and the versatility of the general solution was further justified.     

The use of satellite data to characterize phytoplankton in Sea of Okhotsk water

For the identification of natural processes taking place in surface water layer in the Sea of Okhotsk and nearby water areas, satellite data were used to determine seasonal variations of chlorophyll “a” concentration and the problem of water typification by satellite data in the visible light range was solved. Analysis of the annual cycle of phytoplankton concentration variations confirmed the presence of a short explosion in phytoplankton population in the southern Sea of Okhotsk during ice cover formation, which had been predicted before by model calculations. Classes of water reflecting different geoenvironmental situations were identified for characteristic regions of the Sea of Okhotsk. The dependences of chlorophyll “a” concentration growth on temperature were calculated, and corrections were introduced in the standard relationships used in biooptic algorithms for different parts of the Sea of Okhotsk.     

Hydraulic diffusivity in a coastal aquifer: spectral analysis of groundwater level in responses to marine system

The hydraulic diffusivity gives a measure of diffusion speed of pressure disturbances in groundwater system; large values of hydraulic diffusivity lead to fast propagation of signals in aquifer. This research provides a novel design and derives spectral representation to determine hydraulic diffusivity using spectral analysis of groundwater levels coupled with time-dependent boundary adjacent to marine system and no flow boundary in aquifer system. To validate the proposed method, water levels of fluctuated boundary and groundwater well in a sandy confined aquifer were collected. The hydraulic diffusivity is then obtained by an inverse process in the non-linear complex form of spectral relationship. The method essentially is constructed on the conceptual design of natural forcing transmitted in large aquifer. It is unlike the conventional field pumping test which is only used to determine hydraulic properties of groundwater in small range around the well. Hydraulic diffusivity of the confined aquifer is determined using real observation and then checked by comparing to the published range. It suggests that without local aquifer test to estimate hydraulic diffusivity in a coastal aquifer using spectral representation with its relevant flow system and boundary has become feasible.     

Closed-form analytical solutions incorporating pumping and tidal effects in various coastal aquifer systems

Pumping wells are common in coastal aquifers affected by tides. Here we present analytical solutions of groundwater table or head variations during a constant rate pumping from a single, fully-penetrating well in coastal aquifer systems comprising an unconfined aquifer, a confined aquifer and semi-permeable layer between them. The unconfined aquifer terminates at the coastline (or river bank) and the other two layers extend under tidal water (sea or tidal river) for a certain distance L. Analytical solutions are derived for 11 reasonable combinations of different situations of the L-value (zero, finite, and infinite), of the middle layer’s permeability (semi-permeable and impermeable), of the boundary condition at the aquifer’s submarine terminal (Dirichlet describing direct connection with seawater and no-flow describing the existence of an impermeable capping), and of the tidal water body (sea and tidal river). Solutions are discussed with application examples in fitting field observations and parameter estimations.     

Determination of aquifer parameters based on measurements of tidal effects on a coastal aquifer near Beihai,China

Measurements of the tide and groundwater levels in coastal zones are of importance in determining the properties of coastal aquifers. The solution to a one‐dimensional unsteady groundwater flow model in a coastal confined aquifer with sinusoidal fluctuation of the tide shows that the tidal efficiency decreases exponentially with distance and the time lag increases linearly with distance from the coast. The aquifer property described by the ratio of storage coefficient to transmissivity is determined if the damping constant of the tidal efficiency or the slope of the time lag with distance are obtained on the basis of tidal measurements. Hourly observations of the tide and groundwater levels at 10 wells on the northern coast near Beihai, China show that with distance from the coast, tidal efficiency decreases roughly exponentially and the time lag increases roughly linearly. The estimated ratio of storage coefficient to transmissivity of the confined aquifer ranges from 1·169 × 10^?6 d m^?2 to 1·83 × 10^?7 d m^?2. For a given transmissivity of 750 m² d^?1, the storage coefficient of the aquifer is 8·7675 × 10^?4 with the tidal efficiency method and 1·3725 × 10^?4 with the time lag method. The damping constant of the tidal efficiency with distance can be defined as the tidal propagation coefficient. The value of the confined aquifer is determined as 0·0018892 m^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

Effects of sea-level rise on ground water flow in a coastal aquifer system

The effects of sea-level rise on the depth to the fresh water/salt water interface were simulated by using a density-dependent, three-dimensional numerical ground water flow model for a simplified hypothetical fresh water lens that is similar to shallow, coastal aquifers found along the Atlantic coast of the United States. Simulations of sea-level rise of 2.65 mm/year from 1929 to 2050 resulted in an increase in water levels relative to a fixed datum, yet a net decrease in water levels relative to the increased sea-level position. The net decrease in water levels was much greater near a gaining stream than farther from the stream. The difference in the change in water levels is attributed to the dampening effect of the stream on water level changes in response to sea-level rise. In response to the decreased water level altitudes relative to local sea level, the depth to the fresh water/salt water interface decreased. This reduction in the thickness of the fresh water lens varied throughout the aquifer and was greatly affected by proximity to a ground water fed stream and whether the stream was tidally influenced. Away from the stream, the thickness of the fresh water lens decreased by about 2% from 1929 to 2050, whereas the fresh water lens thickness decreased by about 22% to 31% for the same period near the stream, depending on whether the stream was tidally influenced. The difference in the change in the fresh water/salt water interface position is controlled by the difference in the net decline in water levels relative to local sea level.     

Modelling groundwater levels in an urban coastal aquifer using artificial neural networks

The prediction of groundwater levels in a basin is of immense importance for the management of groundwater resources, especially in coastal regions where the water table fluctuations are to be limited to avoid seawater intrusion. In this paper, an Artificial Neural Network (ANN) methodology is presented to predict groundwater levels in individual wells with one month lead. Groundwater levels were also predicted in neighboring wells using model parameters from the best network of a well. This methodology is applied to an urban coastal aquifer in Andhra Pradesh state, India. The results suggest that the feed forward neural network with Levenberg Marquardt (LM) algorithm is a good choice for predicting groundwater levels in individual wells. Bayesian Regularization (BR) model parameters of Balaji Nagar well are also used successfully to predict groundwater levels in the study area. It was observed that the ANN‐based algorithms were a better choice for the prediction of groundwater levels with limited hydrological parameters. Copyright © 2007 John Wiley & Sons, Ltd.     

The impact of highly permeable layer on hydraulic system in a coastal aquifer

This study is aimed to understand the hydraulic mechanism of coastal aquifer systems that include highly permeable layers (HPLs). These hydrologic conditions can be found in many volcanic islands that are composed of a series of lava flows discharged into sea or other standing body of water. In the first part, we developed a numerical model based on the geologic and hydrologic data obtained from the eastern Jeju Island, Korea, of which the aquifer contains clinker and hyaloclastite layers. The simulation results reproduced spatial location of fresh‐saline water interface, especially the abrupt decline of interface at the inland part and the thickness variation of transition zone along the cross‐section observed at the eastern Jeju coastal aquifer. We were able to find out that these phenomena are strongly related to the presence of the HPL. In the second part, quantitative analyses were conducted with the use of hypothetical models in order to understand the dynamic characteristics of coastal system that includes HPLs. A series of sensitivity studies were conducted to assess the effect of the horizontal length and vertical depth of HPL on the spatial location of the interface toe and the configuration of transition zone. Various case studies have shown that the seawater intruded into the inland more as the horizontal length of HPL was increased and its vertical depth was decreased. In other simulations including two HPLs, the vertical distance between these two HPLs primarily controlled the flow regime, flux variations, and the configuration of the transition zone. Finally, we performed simulations to evaluate the effect of a rising sea‐level. This study provides more understanding of how the presence of HPL controls the seawater intrusion processes, and the spatial configurations of fresh‐saline water interface at coastal aquifers. Copyright © 2012 John Wiley & Sons, Ltd.     

The effects of the 2004 tsunami on a coastal aquifer in Sri Lanka

On December 26, 2004, the earthquake off the southern coast of Sumatra in the Indian Ocean generated far-reaching tsunami waves, resulting in severe disruption of the coastal aquifers in many countries of the region. The objective of this study was to examine the impact of the tsunami on groundwater in coastal areas. Field investigations on the east coast of Sri Lanka were carried out along a transect located perpendicular to the coastline on a 2.4 km wide sand stretch bounded by the sea and a lagoon. Measurements of groundwater table elevation and electrical conductivity (EC) of the groundwater were carried out monthly from October 2005 to August 2007. The aquifer system and tsunami saltwater intrusion were modeled using the variable-density flow and solute transport code HST3D to understand the tsunami plume behavior and estimate the aquifer recovery time. EC values reduced as a result of the monsoonal rainfall following the tsunami with a decline in reduction rate during the dry season. The upper part of the saturated zone (down to 2.5 m) returned to freshwater conditions (EC < 1000 μS/cm) 1 to 1.5 years after the tsunami, according to field observations. On the basis of model simulations, it may take more than 15 years for the entire aquifer (down to 28 m) to recover completely, although the top 6 m of the aquifer may become fresh in about 5 years.