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.

     

A geostatistically based approach to the identification of aquifer transmissivities in Yolo Basin, California

A geostatistically based approach is developed for the identification of aquifer transmissivities in Yolo Basin, California. The approach combines weighted least-squares with universal kriging and cokriging techniques in an overall scheme that (1)considers a prioriknown information on aquifer transmissivity and specific capacities of wells, (2)considers uncertainties in water level and transmissivity data, and (3)estimates the reliability of the generated transmissivity values. Minimization of a global least-squares function that incorporates calibration and plausibility criteria leads to a transmissivity map that shows a good agreement with pumping-test results.     

The use of the marine biotic index AMBI in the assessment of the ecological status of the Obidos lagoon (Portugal)

This study reports the longitudinal distribution of sediment properties, including inorganic and organic contaminants, and the structure of the benthic community in Obidos lagoon, a coastal system permanently connected to the sea and with negligible freshwater sources. Sediments from the upper to central lagoon consist of fine particles (91%) and from the lower lagoon of sands (94%). Chemical composition is strongly correlated to the percentage of fine particles. Contamination is relatively low in those sediments suggesting the effect of diffuse sources. The increase in organic matter content from down- to upstream areas was associated with the dominance of opportunistic species, while sensitive and indifferent species to organic enrichment were mainly associated to the clean sandy downstream area. The marine biotic index (AMBI) was suitable for the discrimination of the biological and environmental gradients in the Obidos lagoon and was highly related with the gradient of organic matter content in this system.     

Syntectonic emplacement of the Middle Jurassic Concón Mafic Dike Swarm, Coastal Range, central Chile (33° S)

The Concón Mafic Dike Swarm (CMDS) consists of basaltic to andesitic dikes emplaced into deformed Late Paleozoic granitoids during the development of the Jurassic arc of central Chile. The dikes are divided into an early group of thick dikes (5–12 m) and a late group of thin dikes (0.5–3 m). Two new amphibole ⁴⁰Ar/³⁹Ar dates obtained from undeformed and deformed dikes, constrain the age of emplacement and deformation of the CMDS between 163 and 157 Ma. Based on radiometric ages, field observations, AMS studies and petrographic data, we conclude that the emplacement of the CMDS was syntectonic with the Jurassic arc extension and associated with sinistral displacements along the NW-trending structures that host the CMDS. The common occurrence of already deformed and rotated xenoliths in the dikes indicates that deformation in the granitoids started previously.The early thick dikes and country rocks appear to have been remagnetized during the exhumation of deep-seated coastal rocks in the Early Cretaceous (around 100 Ma). The remanent magnetization in late thin dikes is mainly retained by small amounts of low-Ti magnetite at high temperature and pyrrhotite at low temperature. The magnetization in these dikes appears to be primary in origin. Paleomagnetic results from the thin dikes also indicate that the whole area was tilted  23° to the NNW during cooling of the CMDS.The NNW–SSE extension vectors deduced from the paleomagnetic data and internal fabric of dikes are different with respect to extension direction deduced for the Middle–Late Jurassic of northern Chile, pointing to major heterogeneities along the margin of the overriding plate during the Mesozoic or differences in the mechanisms driving extension during such period.     

Possible pingo fields in the Utopia basin, Mars: Geological and climatical implications

The presence of pingos on Mars has been hypothesized since the period of the Viking mission. In fact, a diverse range of pingo-like features has been found at various martian sites including Elysium, Chryse and Utopia Planitiae in the northern lowlands. Due to the morphology and the geological setting, some of those features were interpreted in different ways, creating some controversies, as happened in Athabasca Valles. This reflects the complexity of interpreting these features by remote sensing and multiple plausible interpretations of the same feature. With the objective of identifying new possible pingos or rootless cones on Mars, we selected a study area in Utopia Planitia (10-55° N, 210-260° W) where the presence of both features is possible due to its geological history (volcanic and hydrological). We analyzed more than 2100 Mars Orbiter Camera (MOC)-narrow angle images in addition to Viking, Thermal Emission Imaging System (THEMIS), and High Resolution Stereo Camera (HRSC) images, together with Mars Orbiter Laser Altimeter (MOLA)-derived Digital Elevation Models (DEMs) with a Geographic Information System (GIS). We found in 94 MOC-narrow angle images dome, cone, and ring-shaped features. We analyzed them from morphological and morphometrical points of view in order to compare them with relevant features on Mars and Earth. We tested different possible origins for those features following the approach of multiple working hypotheses. We conclude that the dome, cone, and ring-shaped features could be pingos, which is in agreement with their geological settings. Regarding the driving heat source for the formation of the purported pingos, we propose the existence of a heat source, possibly a magma chamber, underneath the surface of the Utopia basin. Together with possible climatic shifts, the past activities of the heat source may have caused melting of ground ice. The pingo growth due to freezing of the water would have occurred during the following cold climatic conditions.     

Trees, lateral shrinkage and building damage

This article presents the study conducted and the actions taken to stabilise the foundations of 25 low-rise (two/three-storey) buildings in Alcázar de San Juan (Ciudad Real, central Spain). Owing to the immediate attention required, it was decided that the study would focus on procuring a synthetic model that would be able to provide a satisfactory explanation for the tendency of the movements, overlooking the detailed model of each building. This analysis pointed out that the mobilisation was probably due to both the shrinkage caused by the existing trees during the dry-growing seasons, and to the loss of lateral confinement caused by the excavation of a service trench adjacent to the footings of the buildings. In keeping with this mechanism, stabilization was carried out by transplanting the trees growing near the damaged buildings and installing a line of piles to provide lateral support to the footings. The evolution of the damage was therefore stopped.     

Upwelling signals in radiocarbon from early 20th-century Peruvian bay scallop (Argopecten purpuratus) shells

We quantified Δ¹⁴C, δ¹⁸O, and δ¹³C cycles along ontogeny within four bay scallop (Argopecten purpuratus) shells collected from Callao Bay, Salaverry, and Sechura Bay, Peru following the 1907–1908 non-El Niño years and the 1925–1926 El Niño. Δ¹⁴C and δ¹³C generally covary; Δ¹⁴C and δ¹⁸O vary inversely. Simultaneous decreases in Δ¹⁴C and increases in δ¹⁸O in non-El Niño shells are followed by constant Δ¹⁴C and gradually decreasing δ¹⁸O, which we interpret as evidence for discrete marine upwelling events followed by warming of the initially cold upwelled water. Upwelling changes from El Niño events are detectable with difficulty in mollusk shell Δ¹⁴C.     

Model prediction capacity of ephemeral gully evolution in conservation tillage systems

Ephemeral gully (EG) erosion has an important impact on agricultural soil losses and increases field surface hydrology connectivity and transport of pollutants to nearby water bodies. Watershed models including an EG component are scarce and not yet properly evaluated. The objective of this study is to evaluate the capacity of one such tool, AnnAGNPS, to simulate the evolution of two EG formed in a conservation tillage system. The dataset for model testing included runoff measurements and EG morphological characteristics during 3 years. Model evaluation focused on EG evolution of volume, width, and length model outputs, and included calibration and testing phases and a global sensitivity analysis (GSA). While the model did not fully reproduce width and length, the model efficiency to simulate EG volume was satisfactory for both calibration and testing phases, supporting the watershed management objectives of the model. GSA revealed that the most sensitive factors were EG depth, critical shear stress, headcut detachment exponent coefficient b, and headcut detachment leading coefficient a. For EG outputs the model was additive, showing low sensitivity to interactions between the inputs. Prediction of EG spatial evolution on conservation tillage systems requires improved development of gully erosion components, since many of the processes were developed originally for traditional tillage practices or larger channel systems. Our results identify the need for future research when EG form within conservation tillage systems, in particular to study gully headcut, soil erodibility, and width functions specific to these practices.     

Damages observed in locations of Oaxaca due to the Tehuantepec Mw8.2 earthquake,Mexico

On September 7, 2017, at 23:49 h (local time), a Mw8.2 intermediate-depth normal-fault earthquake occurred in the Gulf of Tehuantepec, 133 km away from Pijijiapan, Chiapas, and about 700 km away from Mexico City. This event caused 95 fatalities and severe damage to different types of structures located close to the epicenter. The main objective of this work is to present observed damages caused in the state of Oaxaca by this earthquake, which were mainly concentrated in self-built houses and historical and ancient buildings. The locations visited by the reconnaissance team of the Institute of Engineering from UNAM in Oaxaca included Salina Cruz, Tehuantepec, Ixtaltepec, Juchitán, Huatulco and La Ventosa.