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.

     

Utilisation of the image analysis method for the detection of the morphological anisotropy of calcite grains in marble

This paper presents the development and utilisation of an automated image processing algorithm for detection and analysis of grains. Using optical polarising microscopy, a set of colored images are collected from an area on a thin section. A filtering operation, using rotation of a morphological alternating sequence filter (based on a structuring element), is used to remove twinning features within individual grains. Filtering is followed by the watershed segmentation technique to determine grain boundaries. The method is used for the identification of calcite grains in marble and the subsequent analysis of morphological anisotropy.     

Hydrogeological research on intensively exploited deep aquifers in the ‘Loma de Úbeda’ area (Jaén, southern Spain)

The intensive use of groundwater for irrigation in the area of Úbeda (‘Loma de Úbeda’, Jaén, southern Spain) has transformed an area of traditionally rain-fed dry farmland into fields with some of the highest olive oil productivity in the world. Early hydrogeological research studies, initiated just after the beginning of the groundwater exploitation, revealed that the water was collected from three different overlapping aquifers occupying an area of over 1,100 km², with the lower aquifers located at depths from 300 to over 700 m in an area of 440 km². Multidisciplinary research, based on geological characterization, and piezometric, hydrochemical and isotopic data, has led to a conceptual model of functioning in this complex hydrogeological system. The proposed model allows for the identification of the recharge areas, and the discharge, which is at present mainly associated with the groundwater pumping. Areas of mixing of waters from the different aquifers and the main hydrogeochemical processes affecting groundwater quality are described.     

Study of a rockfall in a limy conglomerate canyon (Covarrubias, Burgos, N. Spain)

In this paper, a rockfall of some tenths of cubic metres of rock, occurred in a conglomerate canyon in the northern Castilian plateau (N-Central Spain), is described in detail. This includes the basic data (geology, mapping, meteorological data and rock involved), the mechanical mechanism of the instability, the propagation and reach of the fallen rock boulders and the final location of the pieces. The study is based on detailed geological, aerial photography and geotechnical studies and advance topographical measurements. Geological engineering techniques are used to analyze the instability initiation and rockfall propagation. This study contributes to a better understanding of the geological processes leading to the occurrence of calcareous canyons and cliffs in the drainage network of the upper Castilian plateau. It is also helpful to improve our knowledge on natural rockfall phenomena for risk engineering purposes aiming to quantify, control and mitigate associated hazards.     

Metal-enriched nanoparticles and black carbon: A perspective from the Brazil railway system air pollution

Having a better understanding of air pollutants in railway systems is crucial to ensure a clean public transport. This study measured, for the first time in Brazil, nanoparticles (NPs) and black carbon (BC) on two ground-level platforms and inside trains of the Metropolitan Area of Porto Alegre (MAPA). An intense sampling campaign during thirteen consecutive months was carried out and the chemical composition of NPs was examined by advanced microscopy techniques. The results showed that highest concentrations of the pollutants occur in colder seasons and influenced by variables such as frequency of the trains and passenger densities. Also, internal and external sources of pollution at the stations were identified. The predominance of NPs enriched with metals that increase oxidative stress like Cd, Fe, Pb, Cr, Zn, Ni, V, Hg, Sn, and Ba both on the platforms and inside trains, including Fe-minerals as hematite and magnetite, represents a critical risk to the health of passengers and employees of the system. This interdisciplinary and multi-analytical study aims to provide an improved understanding of reported adverse health effects induced by railway system aerosols.     

Hydrologic models coupled with 2D hydrodynamic model for high-resolution urban flood simulation

Natural Hazards - Floods are the most frequent natural disaster and pose a very challenging threat to many cities worldwide. Understanding the flood dynamic is essential for developing strategies...     

Cohesive‐frictional model for bond and splitting action of prestressing wire

This paper presents a numerical procedure for bond between indented wires and concrete, and the coupled splitting process of the surrounding concrete. The bond model is an interface, non‐associative, plasticity model. It is coupled with a cohesive fracture model for concrete to take into account the splitting of such concrete. Bond between steel and concrete is fundamental for the transmission of stresses between both materials in precast prestressed concrete. Indented wires are used to improve the bond in these structural elements. The radial component of the prestressing force, increased by Poisson's effect, may split the surrounding concrete, decreasing the wire confinement and diminishing the bonding. The combined action of the bond and the splitting is studied with the proposed model. The results of the numerical model are compared with the results of a series of tests, such as those which showed splitting induced by the bond between wire and concrete. Tests with different steel indentation depths were performed. The numerical procedure accurately reproduces the experimental records and improves knowledge of this complex process. Copyright © 2010 John Wiley & Sons, Ltd.     

The Sima de los Huesos (Burgos,northern Spain): palaeoenvironment and habitats of Homo heidelbergensis during the Middle Pleistocene

Interpreting how environmental dynamics respond to global climate change and how this has affected human evolution and dispersal is an on-going topic of debate. During the early Middle Pleistocene (～0.6–0.4 Ma), as compared to earlier, environmental conditions were relatively more stable, with longer climatic cycles alternating between open and forested landscapes. During this interval, humans spread successfully providing an important number of fossil sites where fossils or tools are reported. The Atapuerca-Sima de los Huesos (Burgos, northern Spain) site (Atapuerca-SH) is one of the earliest localities with hominin evidence in the European Middle Pleistocene, with the most important accumulation of Homo heidelbergensis so far. We have analyzed the abundant faunal record from Sima de los Huesos, which is mainly comprised of carnivores, in order to approach an interpretation of the palaeoenvironmental circumstances where these hominids inhabited within the Sierra. Other sites from Sierra de Atapuerca referred to the same Faunal Unit (FU 6), are roughly contemporaneous, and include important ungulates, which are here analyzed with Atapuerca-SH. Additional information provided by isotopic analysis helps elucidate the ancient ecology of taxa present in Sima de los Huesos allowing for an accurate portrayal of the setting in which humans lived. The timing of the spread of Homo heidelbergensis is dominated by a relative climatic and environmental stability and points to a landscape dominated by savannah-like open woodland.     

Structure and stratigraphy of the Teplá–Barrandian Neoproterozoic,Bohemian Massif: A new plate-tectonic reinterpretation

The Teplá–Barrandian unit (TBU) of the Bohemian Massif exposes a section across the once extensive Avalonian–Cadomian belt, which bordered the northern active margin of Gondwana during late Neoproterozoic. This paper synthesizes the state-of-the-art knowledge on the Cadomian basement of the TBU to redefine its principal component units, to revise an outdated stratigraphic scheme, and to interpret this scheme in terms of a recent plate-tectonic model for the Cadomian orogeny in the Bohemian Massif. The main emphasis of this paper is on an area between two newly defined fronts of the Variscan pervasive deformation to the NW and SE of the Barrandian Lower Paleozoic overlap successions. This area has escaped the pervasive Variscan (late Devonian to early Carboniferous) ductile reworking and a section through the Cadomian orogen is here superbly preserved.The NW segment of the TBU consists of three juxtaposed allochthonous belts of unknown stratigraphic relation (the Kralovice–Rakovník, Radnice–Kralupy, and Zbiroh–?árka belts), differing in lithology, complex internal strain patterns, and containing sedimentary and tectonic mélanges with blocks of diverse ocean floor (meta-)basalts. We summarize these three belts under a new term the Blovice complex, which we believe represents a part of an accretionary wedge of the Cadomian orogen.The SE segment of the TBU exposes the narrow Pi?ín belt, which is probably a continuation of the Blovice complex from beneath the Barrandian Lower Paleozoic, and a volcanic arc sequence (the Davle Group). Their stratigraphic relation is unknown. Flysch units (the ?těchovice Group and Svrchnice Formation) overlay the arc volcanics, and both units contain material derived from volcanic arc. The former was also sourced from the NW segment, whereas the latter contains an increased amount of passive margin continental material. In contrast to the Blovice complex, the flysch experienced only weak Cadomian deformation.The new lithotectonic zonation fits the following tectonic scenario for the Cadomian evolution of the TBU well. The S- to SE-directed Cadomian subduction beneath the TBU led to the involvement of turbidites, chaotic deposits, and 605 ± 39 Ma ocean floor in the accretionary wedge represented by the Blovice complex. The accretionary wedge formation mostly overlapped temporally with the growth of the volcanic arc (the Davle Group) at ~ 620–560 Ma. Upon cessation of the arc igneous activity, the rear of the wedge and some elevated portions of the arc were eroded to supply the deep-water flysch sequences of the ?těchovice Group, whereas the comparable Svrchnice Formation (~ 560 to < 544 Ma) was deposited in a southeasterly remnant basin close to the continental margin. The Cadomian orogeny in the TBU was terminated at ~ 550–540 Ma by slab breakoff, by final attachment of the most outboard ~ 540 Ma oceanic crust, and by intrusion of ~ 544–524 Ma boninite dikes marking the transition from the destructive to transform margin during the early/middle Cambrian.     

The Varre-Sai chondrite,a Brazilian fall: petrology and geochemistry

The Varre-Sai meteorite fell along the border of the states of Espirito Santo and Rio de Janeiro, Brazil; on 19 June 2010 at 5:40 pm. Petrography and X-ray powder diffraction (XRD) indicate that the rock is an L5 S4 chondrite, with blastoporphyritic texture that has not been previously described. Geochemical data based on major and rare-earth elements (REEs) show that Varre-Sai is highly similar to the other L chondrites. In Harker diagrams, Varre-Sai, L, and LL chondrites form a single group, suggesting no significant chemical differences between them and contributing to the long-standing debate of whether LL chondrites form a distinct group or whether they are a subset of the L group. Harker diagrams also define a trend from E to H and L/LL chondrites, similar to the cosmochemical trends suggested by other authors. The behaviour of Fe₂O_3t and NiO indicates a relationship with Fe-Ni alloys, and their trend in the diagram suggests some chemical differentiation in the ordinary chondrite parental bodies. The REE content in Varre-Sai, normalized to C chondrites, falls in the field of L chondrites and others, but with slight REE enrichment. The chemical differences in chondrites, mainly in REEs, Fe₂O_3t and NiO could be alternatively interpreted as variations in the inherited agglutinated materials as chondrules, Ca–Al-rich inclusions and Fe–Ni nodules.