High-resolution monitoring of biogeochemical gradients in a tar oil-contaminated aquifer

The detailed understanding of in situ biodegradation of petroleum hydrocarbons in porous aquifers requires knowledge on biogeochemical gradients, the distribution of individual redox species and microorganisms. The generally limited spatial resolution of conventional monitoring wells, however, hampers appropriate characterization of small-scale gradients and thus localization of the relevant processes. Groundwater sampling across a BTEX plume in a sandy aquifer by means of a novel high-resolution multi-level well (HR-MLW) is presented here. The presence of distinct and steep biogeochemical gradients is demonstrated in the centimeter and decimeter scale, which could not be resolved with a conventional multi-level well. The thin BTEX plume with a vertical extension of only 80 cm exhibited a decline of contaminant concentrations by two orders of magnitude within a few centimeters in the upper and lower fringe zone. The small-scale distribution of sulfate, sulfide and Fe(II) in relation to the contaminants and elevated δ³⁴S and δ¹⁸O values of groundwater sulfate strongly indicated sulfate and iron reduction to be the dominant redox processes involved in biodegradation. High microbial activities and biomass especially at the plume fringes and the slope of chemical gradients supported the concept that the latter are regulated by microbial processes and transverse dispersion, i.e. vertical mixing of electron donors and acceptors. Transverse dispersion therefore was suggested to be a driving factor controlling biodegradation in porous aquifers, but not exclusively limiting natural attenuation processes at this site. Broad overlapping zones of electron donors and electron acceptors point towards additional factors limiting anaerobic biodegradation in situ. The identification of small-scale gradients substantially contributed to a better understanding of biodegradation processes and hence is a prerequisite for the development of reliable predictive mathematical models and future remediation strategies.     

Geochemical processes and solute transport at the seawater/freshwater interface of a sandy aquifer

Geochemical processes occurring at a seawater/freshwater interface were studied in a shallow coastal siliclastic aquifer containing minor amounts of calcite. Data were collected from 106 piezometers in a 120-m transect from the coastline and landward. In the first 40 m from the coastline, a wedge of saltwater is intruding below the freshwater aquifer. The aquifer is strongly reduced with mineralization of organic matter by methanogenesis in the freshwater aquifer, and sulfate reduction dominating in the most seaward part of the saline aquifer. The spatial separation of cations in the aquifer indicated a slow freshening process where Ca²⁺ from freshwater displaced the marine cations Na⁺ and Mg²⁺ from the exchanger complex. The resulting loss of Ca²⁺ from solution decreases the saturation state for calcite and possibly causes calcite dissolution. A storm-flooding event was recorded where pulses of dense seawater sank through the fresh aquifer. As a result, the terminal electron accepting process switched from methanogenesis to sulfate reduction. The pulses of sinking seawater also triggered cation exchange reactions where Ca²⁺ was expelled from the exchanger by seawater Na⁺ and Mg²⁺. The released Ca²⁺ is being flushed from the aquifer by groundwater flow, and this export of Ca²⁺ will, in the long term, cause decalcification of the sediment. The water composition in the aquifer is in a transient state as the result of various processes that operate on different timescales. Oxidation of organic matter occurs continuously but at a rate decreasing on a geological time scale. The freshening of the aquifer operates on the timescale of a few years. The episodic flooding and sinking of seawater through the aquifer proceeds in the course of days to weeks, but occurs irregularly with years in between.     

利用观测孔中地下水位潮汐效应计算天津滨海新区含水层参数

在沿海地区,尤其是围海造陆工程形成的陆域地区地下水水位受潮汐影响较大,使传统水文地质试验求取含水层参数存在较大误差。因此通过合理概化地下水在潮汐作用下运动规律,建立数学模型,推导解析公式求取沿海含水层参数具有重要意义。分析天津滨海新区两处观测孔地下水位及潮汐波动特征,在滞后时间不明显的情况下,利用观测孔水位变幅数据计算了含水层水头扩散系数,并根据承压含水层储水系数经验值进一步获得含水层渗透系数。通过两个观测孔分别计算,对比计算结果互相验证发现,该方法取得了令人满意的结果。利用地下水潮汐效应计算含水层参数可以广泛应用于沿海地区水文地质工作中。     

Groundwater flow and residence time in a karst aquifer using ion and isotope characterization

In order to identify the origin of the main processes that affect the composition of groundwater in a karstic aquifer, a hydrogeochemical and isotopic study was carried out of water from numerous observation wells located in Sierra de Gador, a semiarid region in SE Spain. Several natural and anthropogenic tracers were used to calculate groundwater residence time within this complex aquifer system. Analysis of major ions enabled the principal geochemical processes occurring in the aquifer to be established, and the samples were classified into four distinctive solute groups according to this criterion. Dissolution of carbonate rocks determines the chemical composition of less mineralized water. In another group, the concurrent dissolution of dolomite and precipitation of calcite in gypsum-bearing carbonate aquifer, where the dissolution of relatively soluble gypsum controls the reaction, are the dominant processes. Marine intrusion results in highly mineralized waters and leads to base exchange reactions. The groundwater enrichment of minor and trace elements allowed classification of the samples into two classes that are linked to different flow patterns. One of these classes is influenced by a slow and/or deep regional flow, where the temperature is generally elevated. The influence of sulphate reduces by up to 40 % the barium concentration due to the barite precipitation. Isotope data (T, ¹⁴C) confirm the existence of recent local flows, and regional flow system, and ages of ground water may reach 8000 years. The importance of gypsum dissolution in this aquifer is proved by the δ³⁴S content.     

Hydrogeological characterization of aquifer in palla flood plain of Delhi using integrated approach

The Yamuna flood plains spread across the northern part of Indian subcontinent is home to millions of people. The ever-growing population in these plains make it difficult to sustain the demand of freshwater resources. However, the productive aquifers of flood plains could provide solution for these issues. In this context, it is necessary to understand the aquifer characteristics. Thus, the paper attempts to characterize the aquifer in Palla area of the flood plain using integrated approach. Besides, grain size analysis and site-litholog study, the nature of aquifer material was also ascertained from bulk mineralogy of the sediments using X-Ray Diffraction. The aquifer parameters were estimated with help of long duration pumping test data. Moreover, the effect of pumping on salinity variation and hydrochemical facies evolution was also examined. The sand dominant, unconfined aquifer was estimated to have horizontal hydraulic conductivity in the range of 25 m/day and vertical hydraulic conductivity of 6–7 m /day. While the specific yield of the aquifer was estimated in the range of 0.07–0.08. It is observed that under conducive active flood plain environment, the given sand mineralogy at the site does not allow salinity increase in groundwater even after more than a decade of groundwater pumping. In fact, over years, hydrochemical facies have evolved towards Bicarbonate type. These things put together make the active flood plain aquifer a sustainable groundwater resource.     

Use of sinkhole and specific capacity distributions to assess vertical gradients in a karst aquifer

The carbonate-rock aquifer in the Great Valley, West Virginia, USA, was evaluated using a database of 687 sinkholes and 350 specific capacity tests to assess structural, lithologic, and topographic influences on the groundwater flow system. The enhanced permeability of the aquifer is characterized in part by the many sinkholes, springs, and solutionally enlarged fractures throughout the valley. Yet, vertical components of subsurface flow in this highly heterogeneous aquifer are currently not well understood. To address this problem, this study examines the apparent relation between geologic features of the aquifer and two spatial indices of enhanced permeability attributed to aquifer karstification: (1) the distribution of sinkholes and (2) the occurrence of wells with relatively high specific capacity. Statistical results indicate that sinkholes (funnel and collapse) occur primarily along cleavage and bedding planes parallel to subparallel to strike where lateral or downward vertical gradients are highest. Conversely, high specific capacity values are common along prominent joints perpendicular or oblique to strike. The similarity of the latter distribution to that of springs suggests these fractures are areas of upward-convergent flow. These differences between sinkhole and high specific capacity distributions suggest vertical flow components are primarily controlled by the orientation of geologic structure and associated subsurface fracturing.     

Stability assessment and characterization of colloids in coastal groundwater aquifer system at Kalpakkam

Natural colloids have a potential role in facilitating the transport of radionuclides in groundwater. To assess the role of mobile colloidal phases in radionuclide transport, characterization of colloids for size, size distribution, zeta potential, surface charge and elemental composition is required. Groundwater samples were collected from 12 borewells in the study area and were characterized with respect to physicochemical parameters. Water quality parameters such as pH, temperature, specific conductance, TDS and dissolved O₂ were measured in situ. Based on salinity, two groundwater types were identified in the study area: (1) freshwater type and (2) brackish type. Laboratory and field-scale turbidity measurements in addition to quantitative analysis of major ions were carried out. It was observed that the colloid concentrations are reduced at higher salinity. Zetasizer, particle size analyzer, scanning electron microscopy, energy dispersive X-ray analysis and X-ray diffraction techniques were used to characterize the colloidal particles in groundwater. Colloids were present in all the groundwater samples, the concentration varying between 0.05 and 6 mg/L. Colloid concentration was greater in wells with low ionic strength and the number of colloidal particles varied between 3 × 10⁹ and 4 × 10¹¹ particles/L. The average colloid size varied between 200 and 350 nm for various samples. The zeta potential of the colloidal particles varied between −25.5 and –34.0 mV. SEM analysis of colloidal particles revealed the presence of clays particularly kaolinite and that the mineral composition of colloids reflects the mineralogy of the aquifer. For proper risk assessment and remediation studies, the role of colloids in radionuclide transport assumes greater importance. This study highlights the need for, and relevance of, detailed colloid characterization to assess its role in the migration of radionuclides from near surface disposal facilities.     

Field characterization and data integration to define the hydraulic heterogeneity of a shallow granular aquifer at a sub-watershed scale

Providing a sound basis for aquifer management or remediation requires that hydrogeological investigations carried out to understand groundwater flow and contaminant transport be based on representative data that capture the heterogeneous spatial distribution of aquifer hydraulic properties. This paper describes a general workflow allowing the characterization of the heterogeneity of the hydraulic properties of granular aquifers at an intermediate scale of a few km². The workflow involves characterization and data integration steps that were applied on a 12-km² study area encompassing a decommissioned landfill emitting a leachate plume and its main surface water receptors. The sediments composing the aquifer were deposited in a littoral–sublittoral environment and show evidence of small-scale transitional heterogeneities. Cone penetrometer tests (CPT) combined with soil moisture and electrical resistivity (SMR) measurements were thus used to identify and characterize spatial heterogeneities in hydraulic properties over the study area. Site-specific statistical relationships were needed to infer hydrofacies units and to estimate hydraulic properties from high-resolution CPT/SMR soundings distributed all over the study area. A learning machine approach was used due to the complex statistical relationships between colocated hydraulic and CPT/SMR data covering the full range of aquifer materials. Application of this workflow allowed the identification of hydrofacies units and the estimation of horizontal hydraulic conductivity, vertical hydraulic conductivity and porosity over the study area. The paper describes and discusses data acquisition and integration methodologies that can be adapted to different field situations, while making the aquifer characterization process more time-efficient and less labor-intensive.     

Computational mechanics of the steel–concrete interface

Concrete cracking in reinforced concrete structures is governed by two mechanisms: the activation of bond forces at the steel–concrete interface and the bridge effects of the reinforcement crossing a macro‐crack. The computational modelling of these two mechanisms, acting at different scales, is the main objective of this paper. The starting point is the analysis of the micro‐mechanisms, leading to an appropriate choice of (measurable) state variables describing the energy state in the surface systems: on the one side the relative displacement between the steel and the concrete, modelling the bond activation; on the other hand, the crack opening governing the bridge effects. These displacement jumps are implemented in the constitutive model using thermodynamics of surfaces of discontinuity. On the computational side, the constitutive model is implemented in a discrete crack approach. A truss element with slip degrees of freedom is developed. This degree of freedom represents the relative displacement due to bond activation. In turn, the bridge effect is numerically taken into account by modifying the post‐cracking behaviour of the contact elements representing discrete concrete cracks crossed by a rebar. First simulation results obtained with this model show a good agreement in crack pattern and steel stress distribution with micro‐mechanical results and experimental results. Copyright © 2001 John Wiley & Sons, Ltd.     

Multi-scale characterization of a complex karst and alluvial aquifer system in southern Germany using a combination of different tracer methods

Water suppliers face major challenges such as climate change and population growth. To prepare for the future, detailed knowledge of water resources is needed. In southern Germany, the state water supplier Zweckverband Landeswasserversorgung provides 3 million people with drinking water obtained from a complex karst and alluvial aquifer system and the river Danube. In this study, a combination of different tracing techniques was used with the goal of a multi-scale characterization of the aquifer system and to gain additional knowledge about groundwater flow toward the extraction wells in the Danube Valley. For the small-scale characterization, selected groundwater monitoring wells were examined using single-borehole dilution tests. With these tests, a wide range of flow behavior could be documented, including fast outflow within just a few hours in wells with good connection to the aquifer, but also durations of many weeks in low-permeability formations. Vertical flow, caused by multiple flow horizons or uprising groundwater, was detected in 40% of the tested wells. A regional multi-tracer test with three injections was used to investigate the aquifer on a large scale. For the highly karstified connection between a swallow hole and a spring group, high flow velocities of around 80 m/h could be documented. Exceptionally delayed arrivals, 250 and 307 days after the injection, respectively showing maximum velocities of 0.44 and 0.39 m/h, were observed in an area where low-permeability sediments overlay the karst conduits. With the chosen methods, a distinct heterogeneity caused by the geological setting could be documented on both scales.

     

Interpreting sustainable yield of an aquifer using a fuzzy framework

Reliable estimates for how much water can be safely withdrawn from aquifers without harming the environment is crucial for identifying new water supply sources and fostering sustainable growth. Methodologies to estimate groundwater availability that are rooted in science and yet accomplishable with minimal data are particularly useful for effectual aquifer management. Also, as groundwater management is increasingly becoming a participatory process, these methodologies must be transparent and easily understood by a wide range of audiences. In addition, proposed approaches must also reconcile imprecision and uncertainties arising from lack of data, differences in stakeholders’ perceptions and limitations associated with incomplete aquifer characterization. In this study, the fundamental concept of water balance is coupled with fuzzy regression to develop a scheme for assessing regional-scale groundwater availability. Using the mass-balance approach, anthropogenic water demands (municipal, industrial and agricultural) and ecological demands (baseflows to rivers) can be incorporated into the availability estimation process. The use of fuzzy regression enables the specification of decision makers’ preferences to the adopted procedure and renders the parameter estimation to be more robust in the presence of extreme values. The methodology is illustrated by using it to estimate groundwater availability in the Gulf coast aquifer, underlying Refugio County, TX, USA.     

Simulation of arch dam – foundation interaction with a new friction interface element

A new model for three-dimensional non-linear contact problems with irreversible friction is presented here for the interaction between the rock foundation and an arch dam structure. Based on the finite element method and load incremental theory, a constraint contact element with displacements and contact stresses as node parameters is developed. In this approach, four contact conditions are considered, i.e. fixed, slip, free and mixed. This model can simulate frictional slippage, decoupling and re-bonding of two bodies initially mating at a common interface or with any initial gaps. Furthermore boundary conditions for this element are discussed and treatment measures proposed. This method is shown to be effective and to have the advantage of being easily implemented into standard finite element programs. Solutions are obtained for a centrally loaded, simply supported composite beam and for an end-loaded elastica with initial gaps in regional contact with a rigid surface. The results obtained for these examples are compared to the plane stress solutions by contact friction analysis. As an application example, Quanshui arch dam located in Ruyuan County of Guangdong Province in southern China is simulated with the new element.     

Spatial characterization of environmental gradients in a coastal lagoon, Chincoteague Bay

Spatial patterns of environmental processes are intrinsic yet complex components of estuaries. Spatial characterization of environmental gradients is a necessary step to better understand and classify estuarine environments. A geographic information system is developed to analyze the major abiotic environmental processes, to evaluate accuracy and spatial uncertainty, and to analyze potential zonation within the choked coastal lagoon of Chincoteague Bay in Maryland and Virginia, USA. Spatially extensive grid-based models of environmental gradients are constructed from existing geospatial and environmental databases, including tidal prism, bathymetry, salinity, wave exposure, and Secchi disk depth. Integration of wetland boundaries and bathymetric data provide for full basin analysis of flushing and tidal prism. Multivariate Principal Components Analysis demonstrates the covariation among gradients and provides an empirical approach to mapping multidimensional zones within the lagoon. The project documents the development of an estuarine geographic information system that can be used to analyze and compare estuarine environments and provide data for environmental decision making.     

Separation of groundwater-flow components in a karstified aquifer using environmental tracers

Groundwater discharges from the intensively karstified Taurus Mountains to the Mediterranean Sea, either along the contact zone between the mountains and the Travertine Plateau (the Kirkgozler Springs, 15 m³/s), or through the travertine (e.g. the Dudenbasi Spring, 18 m³/s) and underneath it (unnamed submarine springs, unknown discharges). In an attempt to identify the hydraulic connections between the various outlet points, groundwater was analyzed for stable and radioactive isotopes, CFCs and He. The upgradient springs, belonging to the Kirkgozler–Dudenbasi system, were proven to be a mixture of recent and older water on the basis of their low ¹⁴C values (12–22.4 pmc), their exceptionally high He content (429–991 μcc/kg) and ³He:⁴He (R:Ra) ratios (1.471–2.602) and their measurable ³H and CFC contents (1.9–5.9 TU and 0.84 to 3.27 pmoles/kg, respectively). The older component probably contains an even lower amount of modern C. However, the undersaturation of the mixture with respect to calcite, its high CO₂ content (up to 83 mg/L) and its enriched ¹³C values (−2.2 to −4.1‰) suggest intensive water/rock interactions, which would contribute ¹⁴C-devoid bicarbonates to the solution. Downgradient springs discharging along the Mediterranean coast contain groundwater contributions from higher altitudes, as evidenced by their depleted δ¹⁸O and δD composition with respect to the local precipitation; however, a larger portion of the recent water component could be contributed from direct precipitation on the travertine. This larger component is reflected in the increased ³H (3.4 to 8.4 TU) and ¹⁴C (32.7–63.6 pmc) contents, atmospheric He (43–82 pmoles/kg), R:Ra values (1.006–1.198) and CFC contamination of the water.     

Development of a sharp interface model that simulates coastal aquifer flow with the coupled use of GIS

A simulation program, which works seamlessly with GIS and simulates flows in coastal aquifers, is presented in the present paper. The model is based on the Galerkin finite element discretization scheme and it simulates both steady and transient freshwater and saltwater flow, assuming that the two fluids are separated by a sharp interface. The model has been verified in simple cases where analytical solutions exist. The simulation program works as a tool of the GIS program, which is the main database that stores and manages all the necessary data. The combined use of the simulation and the GIS program forms an integrated management tool offering a simpler way of simulating and studying saline intrusion in coastal aquifers. Application of the model to the Yermasogia aquifer illustrates the coupled use of modeling and GIS techniques for the examination of regional coastal aquifer systems.

---

Resumen Se ha desarrollado un modelo casi tridimensional de elementos finitos para simular el flujo de agua dulce y salada, tanto en régimen estacionario como en transitorio, en sistemas acuíferos costeros, bajo la hipótesis de separación por medio de una interfaz abrupta. Las ecuaciones del modelo han sido discretizadas mediante un esquema de Galerkin de discretización en elementos finitos. El modelo ha sido verificado en casos sencillos para los que existe solución analítica. Todos los datos necesarios se introducen y gestionan con un Sistema de Información Geográfica [SIG] por ordenador. El programa de simulación forma parte del programa de SIG, constituyendo una herramienta integrada de gestión para estudiar la intrusión salina en acuíferos costeros. La aplicación del modelo al acuífero de Yermasogia ilustra el uso acoplado de las técnicas de modelación y de SIG con el fin de examinar sistemas acuíferos costeros a escala regional.

---

Résumé Pour étudier un système aquifère côtier, nous avons développé un modèle aux éléments finis en quasi 3-D qui simule les écoulements deau douce et deau salée en régime aussi bien permanent que transitoire. Les équations qui les régissent sont discrétisées par un schéma de discrétisation de Garlekin aux éléments finis. Le modèle a été vérifié dans des cas simples où il existe des solutions analytiques. Toutes les données nécessaires sont introduites et gérées grâce à un logiciel de gestion de SIG. Le programme de simulation est utilisé comme un outil du logiciel de SIG, constituant ainsi un outil de gestion intégrée dont le but est de simuler et détudier lintrusion saline dans les aquifères côtiers. Lapplication du modèle à laquifère de Yermasogia illustre lutilisation couplée de la modélisation et des techniques de SIG pour létude des systèmes aquifères côtiers régionaux.

     

A study of the pollution of the Aries River (Romania) using capillary electrophoresis as analytical technique

 This paper examines two issues, the extensive pollution occurring in the Aries River, NW Romania, as a result of unchecked discharge of mining effluents into the river system, and the suitability of capillary electrophoresis (CE) as an analytical method for investigations into water chemistry. The results confirm the first objective by providing details on the pollution of the Aries River and its geochemical system and demonstrate the usefulness of CE. In its upper reaches, the river system is characterized by high contents of SO₄ ^2– as a direct result of acid mine effluents and the oxidation of sulphide minerals on mine dumps as well as inflows from settling ponds. Although continuous dilution by natural branch waters and natural water-rock interaction reduces the pollution to some extend, the total level of SO₄ ^2– remains above European averages. The waters of the Aries River, by comparison, contain contents of Cu²⁺ and Zn²⁺ up to 100 times higher than those of unpolluted river water. Received: 1 November 1999 · Accepted: 3 April 2000     

Emplacement of Semail–Emirates ophiolite at ridge–trench collision

The Oman‐Emirates is the largest and best‐exposed ophiolite; consequently, it has attracted significant interest among scientists, together with serious conflicts. Most geologists regard this ophiolite as having formed in an intra‐oceanic subduction zone before being accreted to the Arabian continent. Here, we propose an alternative scenario, supported by detailed field observations and integrated geophysics. The smaller Emirates part of the ophiolite was forced into a nearby continent, in the pre‐collision stage of Tethyan closure. The contraction led to the exhumation of the mantle floor of segmented basins accreted in a rifted system similar to the present‐day Gulf of California. The implied high temperature–high pressure metamorphism and the range of geochemical signatures were introduced during the process of rifting, whereas the larger Oman ophiolite was emplaced by obduction onto and along the subducting continental shore. This Ridge–Trench–Transform system might call for a new process to obduct over continents in particular Tethyan ophiolites.     

Evidence of ancient life at 207 m depth in a granitic aquifer

The results of electron-microscopy investigations of calcite precipitated in a water-conducting fracture in a ca. 1800 Ma granitic rock from 207 m below sea level at the island of Aspo on the southeastern (Baltic) coast of Sweden are compared with measurements of carbon, oxygen, and sulfur isotope composition of the calcite and embedded pyrite. Parts of the calcite had extremely low delta 13C values, indicative of biological activity, and contained bacteria-like microfossils occurring in colonies and as typical biofllms. X-ray microanalysis demonstrated these fossils to be enriched in carbon. Our results provide evidence for ancient life in deep granitic rock aquifers and suggest that the modern microbial life found there is intrinsic. Modeling historical and present geochemical processes in deep granitic aquifers should, therefore, preferably include biologically catalyzed reactions. The results also suggest that the search for life on other planets, e.g., Mars, should include subsurface material.