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.

     

Origins of high pH mineral waters from ultramafic rocks, Central Portugal

This paper reviews the geochemical, isotopic (²H, ¹⁸O, ¹³C, ³H and ¹⁴C) and numerical modelling approaches to evaluate possible geological sources of the high pH (11.5)/Na–Cl/Ca–OH mineral waters from the Cabeço de Vide region (Central-Portugal). Water–rock interaction studies have greatly contributed to a conceptual hydrogeological circulation model of the Cabeço de Vide mineral waters, which was corroborated by numerical modelling approaches. The local shallow groundwaters belong to the Mg–HCO₃ type, and are derived by interaction with the local serpentinized rocks. At depth, these type waters evolve into the high pH/Na–Cl/Ca–OH mineral waters of Cabeço de Vide spas, issuing from the intrusive contact between mafic/ultramafic rocks and an older carbonate sequence. The Cabeço de Vide mineral waters are supersaturated with respect to serpentine indicating that they may cause serpentinization. Magnesium silicate phases (brucite and serpentine) seem to control Mg concentrations in Cabeço de Vide mineral waters. Similar δ²H and δ¹⁸O suggest a common meteoric origin and that the Mg–HCO₃ type waters have evolved towards Cabeço de Vide mineral waters. The reaction path simulations show that the progressive evolution of the Ca–HCO₃ to Mg–HCO₃ waters can be attributed to the interaction of meteoric waters with serpentinites. The sequential dissolution at CO₂ (g) closed system conditions leads to the precipitation of calcite, magnesite, amorphous silica, chrysotile and brucite, indicating that the waters would be responsible for the serpentinization of fresh ultramafic rocks (dunites) present at depth. The apparent age of Cabeço de Vide mineral waters was determined as 2790 ± 40 a BP, on the basis of ¹⁴C and ¹³C values, which is in agreement with the ³H concentrations being below the detection limit.     

Magmatic history and geophysical signature of a post-collisional intrusive center emplaced near a crustal-scale shear zone: the Plechý granite pluton (Moldanubian batholith,Bohemian Massif)

The Plechy pluton, southwestern Bohemian Massif, represents a late-Variscan, complexly zoned intrusive center emplaced near the crustal-scale Pfahl shear zone; the pluton thus provides an opportunity to examine the interplay among successive emplacement of large magma batches, magmatic fabric acquisition, and the late-Variscan stress field associated with strike-slip shearing. The magmatic history of the pluton started with the emplacement of the porphyritic Plechy and Haidmühler granites. Based on gravity and structural data, we interpret that the Plechy and Haidmühler granites were emplaced as a deeply rooted, ∼NE–SW elongated body; its gross shape and internal fabric (steep ∼NE–SW magmatic foliation) may have been controlled by the late-Variscan stress field. The steep magmatic foliation changes into flat-lying foliation (particularly recorded by AMS) presumably as a result of divergent flow. Magnetic lineations correspond to a sub-horizontal ∼NE–SW finite stretch associated with the divergent flow. Subsequently, the Třístoličník granite, characterized by steep margin-parallel magmatic foliation, was emplaced as a crescent-shaped body in the central part of the pluton. The otherwise inward-younging intrusive sequence was completed by the emplacement of the outermost and the most evolved garnet-bearing granite (the Marginal granite) along the southeastern margin of the pluton. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

CPLat: first operational experimental processing center for SIRGAS in Argentina

The SIRGAS permanent GPS network which is in fact the IGS network densification for the American continent, consists today of more than 200 stations covering the continent and islands. It is currently processed by the IGS RNAAC SIR centre at Deutsches Geodätisches Forschungsinstitut producing weekly free solutions relying on IGS final orbits and EOP that contribute to the ITRF through IGS. By August 2006, the SIRGAS Working Group I had accepted five proposals for experimental processing centers within the region that would collaborate with IGS RNAAC SIR. One of them, Centro de Procesamiento La Plata (CPLat) in Argentina, began processing 60 stations on October 2006. By January 2007 CPLat reached operational capability, delivering weekly free solution SINEX files, with an internal consistency of 1.5 mm average for the horizontal components, and 3 mm in the vertical. Comparisons with IGS global and IGS RNAAC SIR weekly solutions were taken as external consistency indications, showing average RMS residuals of 1.8, 2.4 and 5 mm for the north, east, and vertical component, respectively. Analysis and comparison of adjusted solution time series from CPLat and other processing centers has proved to be highly valuable for solution QC, namely detection and identification of station anomalous behavior or modelling problems. These procedures will ensure the maintenance of the performance specifications for CPLat solutions. Action is being taken in order to guarantee the continuity of this effort beyond the experimental phase.     

The impact of agriculture on solute and suspended sediment load on a Mediterranean watershed after intense rainstorms

The solute and suspended‐sediment load following five rainstorms (2005–2007) with varied intensities were studied at the Vernegà experimental watershed, north‐western Spain. Two land‐use areas are located within this watershed, the upstream one (forest) with 160 ha a 100% forested area, and the downstream one (agricultural) with 97 ha being 9 ha conventional agricultural field and 88 ha forest. This study investigates the capacity of each land‐use to yield water, suspended sediment concentration (SSC) and dissolved solid concentration (DSC). The hypothesis is that DSC and SSC from the agricultural area are greater than DSC and SSC of the forest area. Results showed that the agriculture area produced significantly greater mean DSC than in the forest area, the main contribution was the Ca²⁺ (24·68 ± 46·52 mg l^?1) ion at the agricultural area. A long‐term sediment production rate at the agricultural outlet was calculated (69·1 tonnes per 100 years) based on the total sediment discharge (TSD) and the recurrence interval of the largest event of the five rainstorms (October 2005). Geographic information system (GIS) spatial data layers of the watershed were produced to determine the relation of tracks, landforms, slopes and forest management to SSC yield in the forest outlet (133·89 ± 308·14 mg l^?1) during the five rainstorms. Agriculture practices are the main cause of soil erosion at the study area. Copyright © 2010 John Wiley & Sons, Ltd.     

Effect of sampling and linkage on fault length and length–displacement relationship

The power-law exponent (n) in the equation: D=cL ⁿ , with D = maximum displacement and L = fault length, would be affected by deviations of fault trace length. (1) Assuming n=1, numerical simulations on the effect of sampling and linkage on fault length and length–displacement relationship are done in this paper. The results show that: (a) uniform relative deviations, which means all faults within a dataset have the same relative deviation, do not affect the value of n; (b) deviations of the fault length due to unresolved fault tip decrease the values of n and the deviations of n increase with the increasing length deviations; (c) fault linkage and observed dimensions either increase or decrease the value of n depending on the distribution of deviations within a dataset; (d) mixed deviations of the fault lengths are either negative or positive and cause the values of n to either decrease or increase; (e) a dataset combined from two or more datasets with different values of c and orders of magnitude also cause the values of n to deviate. (2) Data including 19 datasets and spanning more than eight orders of fault length magnitudes (10⁻²–10⁵ m) collected from the published literature indicate that the values of n range from 0.55 to 1.5, the average value being 1.0813, and the peak value of n _d (double regression) is 1.0–1.1. Based on above results from the simulations and published data, we propose that the relationship between the maximum displacement and fault length in a single tectonic environment with uniform mechanical properties is linear, and the value of n deviated from 1 is mainly caused by the sampling and linkage effects.     

Denudation and retreat of the Serra do Mar escarpment in southern Brazil derived from in situ‐produced 10Be concentration in river sediment

The Serra do Mar escarpment, located along the southeastern coast of Brazil, is a high‐elevation passive margin escarpment. This escarpment evolved from the denudation of granites, migmatites and gneisses. The granites outcrop in the form of a ridge along the escarpment crest, due to its differential erosion (‘sugarloaf’ hills) from the surrounding lithologies. Several studies suggest that the passive margin escarpments are actively retreating toward the interior of the continent. However, no prior study has calculated the long‐term denudation rates of Serra do Mar to test this hypothesis. In this study, we measured the in situ‐produced ¹⁰Be concentration in fluvial sediments to quantify the catchment‐wide long‐term denudation rates of the Serra do Mar escarpment in southern Brazil. We sampled the fluvial sediments from ten watersheds that drain both sides of the escarpment. The average long‐term denudation rate of the oceanic side is between 2.1‐ and 2.6‐fold higher than the rate of the continental side: 26.04 ± 1.88 mm ka^‐1 (integrating over between 15.8 ka^‐1 and 46.6 ka^‐1) and 11.10 ± 0.37 mm ka^‐1 (integrating over between 52.9 ka^‐1 and 85.4 ka^‐1), respectively. These rates indicate that the coastal base level is controlling the escarpment retreat toward the continental high lands, which is consistent with observations made at other high‐elevation passive margins around the globe. The results also demonstrate the differential erosion along the Serra do Mar escarpment in southern Brazil during the Quaternary, where drainages over granites had lower average denudation rates in comparison with those over migmatites and gneisses. Moreover, the results demonstrate that the ocean‐facing catchments have been eroded more intensely than those facing the continent. The results also reveal that drainage over the granites decreases the average denudation rates of the ocean‐facing catchments and the ‘sugarloaf’ hills therefore are natural barriers that slowly retreat once they are exhumed. Copyright © 2013 John Wiley & Sons, Ltd.     

Testing the consistency of the Starburst scenario for active galactic nuclei

In the frame of the Starburst Model, we show that the evolution of a massive stellar cluster in a high metallicity environment can reproduce the observed emission-line spectrum and the UV-optical Spectral Energy Distribution of the Seyfert 2 galaxies and LINERs. We apply the results of our models to three objects: NGC 5506, NGC 5643 and Mk 348.     

Detailed <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar dating of geologic events associated with the Mantos Blancos copper deposit,northern Chile

The ⁴⁰Ar/³⁹Ar geochronological method was applied to date magmatic and hydrothermal alteration events in the Mantos Blancos mining district in the Coastal Cordillera of northern Chile, allowing the distinction of two separate mineralization events. The Late Jurassic Mantos Blancos orebody, hosted in Jurassic volcanic rocks, is a magmatic-hydrothermal breccia-style Cu deposit. Two superimposed mineralization events have been recently proposed. The first event is accompanied by a phyllic hydrothermal alteration affecting a rhyolitic dome. The second mineralization event is related to the intrusion of bimodal stocks and sills inside the deposit. Because of the superposition of several magmatic and hydrothermal events, the obtained ⁴⁰Ar/³⁹Ar age data are complex; however, with a careful interpretation of the age spectra, it is possible to detect complex histories of successive emplacement, alteration, mineralization, and thermal resetting. The extrusion of Jurassic basic to intermediate volcanic rocks of the La Negra Formation is dated at 156.3 ± 1.4 Ma (2σ) using plagioclase from an andesitic lava flow. The first mineralization event and associated phyllic alteration affecting the rhyolitic dome occurred around 155–156 Ma. A younger bimodal intrusive event, supposed to be equivalent to the bimodal stock and sill system inside the deposit, is probably responsible for the second mineralization event dated at ca. 142 Ma. Other low-temperature alteration events have been dated on sericitized plagioclase at ca. 145–146, 125, and 101 Ma. This is the first time that two distinct mineralization events have been documented from radiometric data for a copper deposit in the metallogenic belt of the Coastal Cordillera of northern Chile. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Putative Meteoritic Craters in Río Cuarto (Central Argentina) Interpreted as Eolian Structures

Ten oblong aligned depressions in the Río Cuarto area (provincia de Córdoba, Argentina) were supposed to be the result of very-low-angle Holocene meteoroid impacts. However, we consider that authors that studied the structures did not demonstrated their extraterrestrial origin. We suggest that an eolian origin for the structures of Río Cuarto is more likely. Actually, these landforms integrate large systems of similar deflation/accumulation geoforms aligned according to predominant winds during different periods.