Conditioning mean steady state flow on hydraulic head and conductivity through geostatistical inversion

Nonlocal moment equations allow one to render deterministically optimum predictions of flow in randomly heterogeneous media and to assess predictive uncertainty conditional on measured values of medium properties. We present a geostatistical inverse algorithm for steady-state flow that makes it possible to further condition such predictions and assessments on measured values of hydraulic head (and/or flux). Our algorithm is based on recursive finite-element approximations of exact first and second conditional moment equations. Hydraulic conductivity is parameterized via universal kriging based on unknown values at pilot points and (optionally) measured values at other discrete locations. Optimum unbiased inverse estimates of natural log hydraulic conductivity, head and flux are obtained by minimizing a residual criterion using the Levenberg-Marquardt algorithm. We illustrate the method for superimposed mean uniform and convergent flows in a bounded two-dimensional domain. Our examples illustrate how conductivity and head data act separately or jointly to reduce parameter estimation errors and model predictive uncertainty.This work is supported in part by NSF/ITR Grant EAR-0110289. The first author was additionally supported by scholarships from CONACYT and Instituto de Investigaciones Electricas of Mexico. Additional support was provided by the European Commission under Contract EVK1-CT-1999-00041 (W-SAHaRA-Stochastic Analysis of Well Head Protection and Risk Assessment).     

Coupled thermal, hydraulic and geochemical evolution of pyritic tailings in unsaturated column experiments

The evolution of pore-water and the composition of solid phases in the vadose zone of pyritic tailings was studied by means of unsaturated column experiments. Several columns of water-saturated mine tailings were dried during 125 days under controlled laboratory conditions. The columns were dismantled at four successive drying stages and the evolution of pore-water, mineralogy, water content and temperature was characterized.Sulfide and aluminosilicate minerals present in the waste dissolved, releasing sulfate and other solutes (mainly Fe, Zn, Cu, Al, Mg and Ca) to the pore-water. Evaporation caused a crust of efflorescent, water-soluble sulfates to develop over the complete top surface of the columns and into the pores of the underlying waste material. This crust, which has also been identified in the field, changed the hydraulic properties of the tailings and produced a decrease in the evaporation rate of the columns. Moreover, these water-soluble precipitates (mainly rozenite, szomolnokite, halotrichite, hexahydrite, mirabilite and gypsum) acted as temporary sinks for Cd, Pb, Co and Ni, which could be released to the surface run-off or the groundwaters during rainfall events under field conditions.Pore-water evolution was determined not only by geochemical processes (dissolution of sulfides and aluminosilicates, precipitation of secondary phases) but also by thermal and hydraulic processes. Progressive dilution was observed in the lower part of the columns. Dilution was caused by the thermally driven vapor flux from the top of the column to its colder bottom and subsequent condensation therein. This process, which may also occur in tailings under sub-arid climate, played a key role on the evolution of pore-water with increasing drying.     

Dynamics of negative hydraulic barriers to prevent seawater intrusion

Negative hydraulic barriers that intercept inflowing saltwater by pumping near the coast have been proposed as a corrective measure for seawater intrusion in cases where low heads must be maintained. The main disadvantage of these barriers is that they pump a significant proportion of freshwater, leading to contamination with saltwater at the well. To minimize such mixing, a double pumping barrier system with two extraction wells is proposed: an inland well to pump freshwater and a seawards well to pump saltwater. A three-dimensional variable density flow model is used to study the dynamics of the system. The system performs very efficiently as a remediation option in the early stages. Long-term performance requires a well-balanced design. If the pumping rate is high, drawdowns cause saltwater to flow along the aquifer bottom around the seawater well, contaminating the freshwater well. A low pumping rate at the seawards well leads to insufficient desalinization at the freshwater well. A critical pumping rate at the seawater well is defined as that which produces optimal desalinization at the freshwater well. Empirical expressions for the critical pumping rate and salt mass fraction are proposed. Although pumping with partially penetrating wells improves efficiency, the critical pumping rates remain unchanged.     

Inverse problem in hydrogeology

The state of the groundwater inverse problem is synthesized. Emphasis is placed on aquifer characterization, where modelers have to deal with conceptual model uncertainty (notably spatial and temporal variability), scale dependence, many types of unknown parameters (transmissivity, recharge, boundary conditions, etc.), nonlinearity, and often low sensitivity of state variables (typically heads and concentrations) to aquifer properties. Because of these difficulties, calibration cannot be separated from the modeling process, as it is sometimes done in other fields. Instead, it should be viewed as one step in the process of understanding aquifer behavior. In fact, it is shown that actual parameter estimation methods do not differ from each other in the essence, though they may differ in the computational details. It is argued that there is ample room for improvement in groundwater inversion: development of user-friendly codes, accommodation of variability through geostatistics, incorporation of geological information and different types of data (temperature, occurrence and concentration of isotopes, age, etc.), proper accounting of uncertainty, etc. Despite this, even with existing codes, automatic calibration facilitates enormously the task of modeling. Therefore, it is contended that its use should become standard practice.

---

Resumen Se sintetiza el estado del problema inverso en aguas subterráneas. El énfasis se ubica en la caracterización de acuíferos, donde los modeladores tienen que enfrentar la incertidumbre del modelo conceptual (principalmente variabilidad temporal y espacial), dependencia de escala, muchos tipos de parámetros desconocidos (transmisividad, recarga, condiciones limitantes, etc), no linealidad, y frecuentemente baja sensibilidad de variables de estado (típicamente presiones y concentraciones) a las propiedades del acuífero. Debido a estas dificultades, no puede separarse la calibración de los procesos de modelado, como frecuentemente se hace en otros campos. En su lugar, debe de visualizarse como un paso en el proceso de entendimiento del comportamiento del acuífero. En realidad, se muestra que los métodos reales de estimación de parámetros no difieren uno del otro en lo esencial, aunque sí pueden diferir en los detalles computacionales. Se discute que existe amplio espacio para la mejora del problema inverso en aguas subterráneas: desarrollo de códigos amigables al usuario, acomodamiento de variabilidad a través de geoestadística, incorporación de información geológica y diferentes tipos de datos (temperatura, presencia y concentración de isótopos, edad, etc), explicación apropiada de incertidumbre, etc. A pesar de esto, aún con los códigos existentes, la calibración automática facilita enormemente la tarea de modelado. Por lo tanto, se sostiene que su uso debería de convertirse en práctica standard.

---

Résumé Létat du problème inverse des eaux souterraines est synthétisé. Laccent est placé sur la caractérisation de laquifère, où les modélisateurs doivent jouer avec lincertitude des modèles conceptuels (notamment la variabilité spatiale et temporelle), les facteurs déchelle, plusieurs inconnues sur différents paramètres (transmissivité, recharge, conditions aux limites, etc.), la non linéarité, et souvent la sensibilité de plusieurs variables détat (charges hydrauliques, concentrations) des propriétés de laquifère. A cause de ces difficultés, le calibrage ne peut être séparé du processus de modélisation, comme cest le cas dans dautres cas de figure. Par ailleurs, il peut être vu comme une des étapes dans le processus de détermination du comportement de laquifère. Il est montré que les méthodes dévaluation des paramètres actuels ne diffèrent pas si ce nest dans les détails des calculs informatiques. Il est montré quil existe une large panoplie de techniques d ‹inversion : codes de calcul utilisables par tout-un-chacun, accommodation de la variabilité via la géostatistique, incorporation dinformations géologiques et de différents types de données (température, occurrence, concentration en isotopes, âge, etc.), détermination de lincertitude. Vu ces développements, la calibration automatique facilite énormément la modélisation. Par ailleurs, il est souhaitable que son utilisation devienne une pratique standardisée.

     

Introductory review of specific factors influencing urban groundwater, an emerging branch of hydrogeology, with reference to Barcelona, Spain

A number of specific factors must be considered when dealing with groundwater in urban areas. Urbanization significantly affects the natural water cycle, both in terms of quantity and quality. In particular, the main contributors to recharge and discharge clearly differ from those in natural systems. Moreover, water can affect underground structures and infrastructure characteristics of cities such as basements, public transport services (trains, underground railways, etc.), and utility conduits. As a result, urban groundwater is emerging as a distinct branch of hydrogeology. The objective of this paper is to review some of the topics that are specific to urban groundwater. These include (1) fluctuations in groundwater levels caused by changes in land and water uses; (2) pollution problems caused by point or non-point sources in urban areas; (3) characterization and quantification of the components contributing to groundwater recharge and discharge; (4) specific characteristics of groundwater flow and solute transport models in urban areas; and (5) integration of data for sustainable urban water management. Some of these issues are illustrated for the particular case of Barcelona, where a comprehensive hydrogeological study has been carried out during the last few years by both the public and the private sectors.

---

Resumen Se debe considerar una serie de factores específicos al tratar las aguas subterráneas en áreas urbanas. El proceso de urbanización afecta el ciclo natural del agua de manera significativa tanto en lo que se refiere a cantidad como a calidad. En particular, los contribuyentes principales a la recarga y descarga claramente difieren de aquellos aplicables a los sistemas naturales. Es más, el agua puede afectar los estructuras subterráneas de las ciudades tales como sótanos, servicios de transporte público (trenes, trenes subterráneos) y los servicios básicos. Como resultado de esto, la hidrogeologia urbana está emergiendo como una rama específica de la hidrogeología. El objectivo de este artículo es revisar algunos de los temas que son específicos al agua subterránea urbana. Estos incluyen: 1) Fluctuaciones en los niveles de agua subterránea causados por cambios en los usos de terreno y del agua; 2) Problemas de contaminación causados por fuentes puntuales o no puntuales en áreas urbanas; 3) Caracterización y cuantificación de los componentes que contribuyen a la recarga y descarga de agua subterránea; 4) Características especificas de flujo de agua subterránea y modelos de transporte de solutos en áreas urbanas; y 5) Integración de datos para el manejo sostenible de agua urbana. Algunos de estos puntos se han ilustrados para el caso particular de Barcelona en donde se ha llevado a cabo un estudio hidrogeológico amplio durante los últimos años tanto por el sector privado como el público.

---

Résumé Une série de facteurs spécifiques doit être considérée lorsque leau souterraine en région urbaine est étudiée. Lurbanisation affecte significativement le cycle naturel de leau, autant du point de vue quantitatif que qualitatif. Plus particulièrement, les zones de recharge et de décharge diffèrent considérablement des systèmes naturels. Leau souterraine peut affecter les structures souterraines ainsi que les infrastructures municipales telles que les sous-sols, les systèmes de transport publique (trains, métros, ...) et les conduits de service. De ce fait, lhydrogéologie en milieu urbain émerge comme une branche distincte de lhydrogéologie classique. Lobjectif de cet article est de revoir certaines notions qui sont spécifiques à létude de leau souterraine en milieu urbain. Parmi ces notions : 1) les fluctuations de la nappe phréatique causées par des changements dans lutilisation du territoire et de leau, 2) les problèmes de pollution causés par des sources ponctuelles ou diffuses en milieu urbain, 3) la caractérisation et quantification des composants qui contribuent à la recharge et décharge des aquifères, 4) les caractéristiques spécifiques des modèles découlement de leau souterraine et de transport de solutés en milieux urbain, et 5) lintégration des données pour la gestion durable de leau en milieu urbain. Certaines de ces questions sont étudiées pour le cas particulier de la ville de Barcelone où une étude hydrogéologique complète a été menée au cours des dernières années par les secteurs privé et public.

     

Deterioration processes affecting prehistoric rock art engravings in granite in NW Spain

Prehistoric rock art sites are endangered despite conservation efforts. The lack of scientific documentation regarding weathering agents affecting rock art and the absence of specific diagnostic protocols hinder the development of conservation strategies. The aim of this research was to investigate active deterioration processes in a granite petroglyph site located in Mougás (Galicia, NW Spain) by characterizing the granite, conducting a geotechnical study of the outcrop and describing and analysing the main weathering processes. Two main deterioration factors were identified. First, water favours block disjunction at the massif scale and causes pitting and surface erosion at the millimetre scale that affects the readability of the engravings. Second, high temperatures associated with wildfires cause mineral transformations that increase the susceptibility of the rock to weathering. Identifying deterioration factors is a first step in developing appropriate preventive conservation measures, which should aim to reduce rock contact time with water (technically affordable in the short term) and to reduce the probability of wildfire occurrence (technically more complex and possibly with longer‐term results). Copyright © 2018 John Wiley & Sons, Ltd.