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.

     

Summer boundary-layer height at the plateau site of Dome’C,antarctica

Measurements of the mean and turbulent structure of the planetary boundary layer using a sodar and a sonic anemometer, and radiative measurements using a radiometer, were carried out in the summer of 1999–2000 at the Antarctic plateau station of Dome C during a two-month period. At Dome C strong ground-based inversions dominate for most of the year. However, in spite of the low surface temperatures (between −50 and −20 °C), and the surface always covered by snow and ice, a regular daytime boundary-layer evolution, similar to that observed at mid-latitudes, was observed during summertime. The mixed-layer height generally reaches 200–300 m at 1300–1400 LST in high summer (late December, early January); late in the summer (end of January to February), as the solar elevation decreases, it reduces to 100–200 m. A comparison between the mixed-layer height estimated from sodar measurements and that calculated using a mixed-layer growth model shows a rather satisfactory agreement if we assign a value of 0.01–0.02 m s⁻¹ to the subsidence velocity at the top of the mixed layer, and a value of 0.003–0.004 K m⁻¹ to the potential temperature gradient above the mixed layer.     

Impact of copper smelting on lakes in the southern Ural Mountains,Russia, inferred from chironomids

Karabash (52°2 N, 60°10 E) is a copper smelting town in the southern Ural Mountains of Russia. The town is affected by sulphur dioxide emissions and deposition of metal-rich particulates from the smelter, acid drainage from old mine workings, and leachates from disused waste dumps and tailings dams. The close proximity of houses to these sources of pollution is of concern to human health and has devastated terrestrial vegetation in the environs. The environmental impact of the smelter on lakes in the area has been assessed using chironomids. Short sediment cores were taken from 16 lakes within a 50 km radius of the smelter and the composition of the chironomid fauna from the bottom of each core, representing conditions prior to the commissioning of the smelter in 1910, was compared with the present chironomid fauna in the surface sediments. Redundancy analysis (RDA) showed that changes in the chironomid fauna of most lakes were driven by trophic change, independent of the industrial activity. Lakes and ponds adjacent to the smelter and waste dumps, which directly receive contaminated waters were devoid of macro- and mesofauna and flora, but there is no evidence that other lakes have been severely impacted by smelter emissions. Local geology ensures that the lakes are well-buffered to the effects of acid deposition which will limit the bioavailability of metals in the water column and sediment.     

Perturbation equations of celestial mechanics

Perturbation equations of celestial mechanics in terms of orbital elements are completely derived in application to the motion of interplanetary dust particle in the gravational field of the Sun and under the action of disturbing forces. Consideration of change of mass of interplanetary dust particle is the most important feature of this derivation. The results obtained are completely general in the case of constant masses.     

The halo population G and K subdwarfs

The photometry of G and K subdwarfs in the Vilnius photometric system permits us to recognize them photometrically and to determine their metallicity, despite the presence of interstellar reddening. The weakening of the strong metallic lines in their spectra can also be used to estimate their metallicity.     

Détermination indirecte des limites d'Atterberg par gravimétrie dynamique

This study presents a new method for Atterberg limits measurements using dynamic drying. Indeed, for some polluted clays, or clay mixtures, classical methods cannot be used. The principle of liquid $(W_{\mathrm{L}})$ and plasticity $(W_{\mathrm{P}})$ limits determination by this method is presented. Experimental values for classical clays are close to those obtained with Casagrande's methods, allowing for larger physical mud characterisation. To cite this article: C. Proust et al., C. R. Geoscience 336 (2004).     

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.     

Harmful microalgae blooms (HAB); problematic and conditions that induce them

HAB occurrence is becoming more frequent and problematic in marine recreational waters. However, the exploitation of the coastal area for recreational use is promoting the necessary conditions for the HAB increase. In terms of the harmful effects, we can consider two types of causative organism: the toxic producers and the high-biomass producers. Toxic events can be produced by a very low concentration of the causative organism. This characteristic implies a difficulty for the monitoring programs in relation to human health. It is important to point out in the context of human health and HAB events, that in some coastal regions (e.g. the Mediterranean basin) HABs are an emerging problem. In these regions, the local population and visitors may face a health risk that is difficult to measure. The monitoring of toxic species has mainly been associated -with shellfish farming. However, the risk of intoxication could become even greater in areas not subject to legislation of aquaculture activities.     

Spatial distribution and ecotoxicity of petroleum hydrocarbons in sediments from the Galicia continental shelf (NW Spain) after the Prestige oil spill

The distribution of aliphatic and aromatic hydrocarbons was determined in surface sediments collected at 36 stations along the Galicia continental shelf (NW Spain), following the Prestige oil spill. Sampling was performed in December 2002, just after the accident, and in February and September 2003. Concentrations of PAHs (summation operator 13 parent components) were in the range of 0.9-422 microg/kgdw, the highest values being close to coastal urban areas (e.g. Pontevedra and A Coru?a), whereas in the stations of the area most heavily impacted by the spill (off Costa da Morte) concentrations were in the range of 14.8-89.6 microg/kgdw, with a certain predominance of alkylated compounds, which may suggest a mixture of petrogenic and pyrolytic sources. The detailed study of petrogenic molecular markers (e.g. steranes and triterpanes) showed the occurrence of an old (weathered) petrogenic chronic pollution in the shelf sediments but not of the Prestige oil, with the possible exception of few stations in the area of Costa da Morte. This was attributed to the heavy nature of the spilled oil that was barely dispersed in the water column and mainly stranded on the coast or sedimented in the form of oil patches. The addition of increasing amounts of fuel oil to a representative sediment sample showed that the molecular indices were indicative of the presence of the Prestige oil when the amount was above 1g/kg of sediment. The toxicity of selected samples (showing the higher PAH concentrations) was tested using the bivalve embryogenesis bioassay. Embryogenesis success reached high values in all cases (80-88%, with 86% in the control), indicating a lack of toxicity in the sediments and supporting the conclusion that the patchiness of the fuel eventually reaching the seafloor reduced its impact on the benthic communities of the Galician shelf.     

Conduit-vent structures and related proximal deposits in the Las Cañadas caldera, Tenerife, Canary Islands

The Las Cañadas caldera wall and the outer slopes of the caldera provide three-dimensional exposures of numerous proximal-welded fallout deposits and have been mapped in detail. As a result, some parts of the Ucanca and Guajara Formations of the stratigraphy of Martí et al. (1994) have been divided into members that correspond to individual eruptions. Mapping has also revealed the occurrence of conduit-vent structures associated with proximal-welded fallout deposits. Conduit-vent structures consist of an upper flaring area and a lower narrow conduit. Conduit-vent geometry and dimensions include cylindrical plugs and eruptive fissures steeply dipping towards the caldera depression and elongated vents. The flaring area can be rather asymmetric and is usually filled by down-vent rheomorphic flow of the proximal fallout deposit. The lower conduits are filled by lava plug, agglutination of juveniles onto conduit walls and dyke intrusion with eventual dome extrusion. The eruption dynamics of welded fallout deposits and magma fragmentation within the conduit are consistent with an evolution from explosive to effusive. In this context conduit flow regimes evolve from turbulent to annular flow in which the conduit is progressively choked, and laminar flow leading to the final conduit closure.