An integrative geological and geophysical approach to characterize a superficial deltaic aquifer in the Camargue plain,France

Deltaic aquifers are complex due to the important heterogeneity of their structure and their hydrogeological functioning. Auger drilling provides localized, but very robust geological and hydrogeological information, while geophysical surveys can provide integrated subsurface information. An effective, easy-to-use and low-cost methodology combining geological/hydrogeological information from Auger drillings and the results from three geophysical techniques (Electromagnetic mapping, Electrical Sounding and Electrical Resistivity Tomography) is being developed to characterize the structure of a typical Mediterranean deltaic aquifer. A first hypothesis about hydrodynamic properties of the aquifer is also obtained. The study area is located in the Rhone delta (Middle Camargue/southern France). Integration of geophysical and geological techniques allowed identifying the presence and lateral extension of the Saint-Ferréol paleochannel, the vertical contact between lagoonal–fluvial deposits and the marine clayed silt that separate the superficial aquifer from the deeper aquifer. Likewise, high north–south heterogeneity and east–west homogeneity were highlighted in the study area. Presence of clay in sandy deposits in the low areas implies changes in lateral hydraulic permeability. This fact, jointed to the low hydraulic gradient, suggests a slow groundwater flow in the local system. The Rhone delta presents a typical configuration of a Mediterranean deltaic aquifer, thus this methodological approach can be used for similar deltaic Mediterranean systems.     

Geophysical characterisation of karstic networks – Application to the Ouysse system (Poumeyssen,France)

In the framework of the management of karstic aquifers, geophysical reconnaissance can be used to locate conduits and caves, and to characterise the surrounding limestone matrix. Suitable characterisation of heterogeneities in the karstic environment is, however, challenging for ground-based geophysical methods. The present article describes the results, and evaluates the response and accuracy of combined geophysical measurements carried out at the Poumeyssen test site in France, involving electrical resistivity imaging (ERI), magnetic resonance sounding (MRS), “mise-à-la-masse” electrical mapping, and seismic tomography. This site provides the opportunity to study a relatively wide, shallow, water-filled conduit whose location and shape are known from topographic work carried out by cave divers. Seismic and MRS provided the exact location and width of the conduit, to within a few meters. The seismic and electrical data suggest that the limestone medium surrounding the conduit is not homogeneous.     

Joint use of TEM and MRS methods in a complex geological setting

Transient Electromagnetic (TEM), known also as Time Domain Electromagnetic (TDEM) and Magnetic Resonance Sounding (MRS) methods were applied jointly to investigate variations in lithology and groundwater salinity in the Nahal Hever South area (Dead Sea coast of Israel). The subsurface in this area is highly heterogeneous and composed of intercalated sand and clay layers over a salt formation, which is partly karstified. Groundwater is very saline, with a chloride concentration of 100–225 g/l. TEM is known as an efficient tool for investigating electrically conductive targets like saline water, but it is sensitive to both the salinity of groundwater and the porosity of rocks. MRS, however, is sensitive primarily to groundwater volume, but it also allows delineating of lithological variations in water-saturated formations. MRS is much less sensitive to variations in groundwater salinity in comparison with TEM. We show that MRS enables us to resolve the fundamental uncertainty in TEM interpretation caused by the equivalence between groundwater resistivity and lithology. Combining TEM and MRS, we found that the sandy Dead Sea aquifer filled with Dead Sea brine is characterized by a bulk resistivity of ρ_x > 0.4 Ωm, whereas zones with silt and clay in the subsurface are characterized by a bulk resistivity of ρ_x < 0.4 Ωm. These observations are confirmed by calibration of the TEM method performed near 18 boreholes.     

Chemistry of metal–humic complexes contained in Megalopolis lignite and potential application in modern organomineral fertilization

Lignite samples from two deposits located in the Megalopolis Basin, Southern Greece, were evaluated for their potential applicability as raw materials for the production of organomineral fertilizers. Fundamental chemical analyses were carried out to demonstrate high humic substances and metal contents. To determine their relative distribution in the Megalopolis lignite extract, eight elements, namely Na, K, Cd, Mn, Mg, Pb, Zn, and Cu, were studied both in H₂O and in Na₄P₂O₇/NaOH solutions. The behavior of these metals showed significant variations; Zn, Pb, Cd, and Cu associate mostly to the humic substances and proved scarce in the water extract. Contrarily, K and Mg gave a significantly low total yield in the Na₄P₂O₇/NaOH solution, while Mn was classified among the least extracted elements. Further enrichment of Megalopolis humic substances in these metals was achieved; Pb and Mg proved the most and least retained metal, respectively. Decomplexation titration curves of humic matter saturated with these metal ions demonstrated that novel organomineral fertilizing materials may develop based on optimized metal ion and humate contents, which can retain metals in a soluble form within a wide pH range. Formation of complexes between humic substances and Zn, Cd, and Mg was clearly indicated.     

Multi-model multi-analysis ensemble weather forecasting on the grid for the South Eastern Mediterranean Region

Weather forecasting is based on the use of numerical weather prediction (NWP) models that are able to perform the necessary calculations that describe/predict the major atmospheric processes. One common problem in weather forecasting derives from the uncertainty related to the chaotic behaviour of the atmosphere. A solution to that problem is to perform in addition to “deterministic” forecasts, “stochastic” forecasts that provide an estimate of the prediction skill. A computationally feasible approach towards this aim is to perform “ensemble forecasts”. Indeed, in the frame of SEE-GRID-SCI EU funded project a Regional scale Multi-model, Multi-analysis ensemble forecasting system (REFS) was built and ported on the Grid infrastructure. REFS is based on the use of four limited area models (namely BOLAM, MM5, ETA, and NMM) that are run using a multitude of initial and boundary conditions over the Mediterranean. This paper presents the tools and procedures followed for developing this application at a production level.     

A BME solution of the stochastic three-dimensional Laplace equation representing a geothermal field subject to site-specific information

This work develops a model of the geothermal field in the Nea Kessani region (Greece) by means of the Bayesian maximum entropy (BME) method, which describes the temperature variations across space in the underground geological formations. The geothermal field is formed by a thermal reservoir consisting of arcosic sandstones. The temperature distribution vs depth was first investigated by the Greek Institute of Geology and Mineral Exploration (IGME) using measurements in a set of vertical drill holes. These measurements showed that hot fluids rising from the deep enter the reservoir in a restricted area of the field and flow towards local thermal springs. The field modelling, which was based on the powerful BME concept, involves the solution of a stochastic partial differential equation that assimilates important site-specific information. The stochastic three-dimensional steady-state Laplace equation was considered as general knowledge and the drilling exploration data were used to construct the specificatory knowledge base in the BME terminology. The produced map is more informative and, in general, it gives higher temperature estimates compared to previous studies of the same region. This is also in agreement with the quartz geothermometry analysis carried out by IGME.     

Assimilating remotely sensed snow observations into a macroscale hydrology model

Accurate forecasting of snow properties is important for effective water resources management, especially in mountainous areas like the western United States. Current model-based forecasting approaches are limited by model biases and input data uncertainties. Remote sensing offers an opportunity for observation of snow properties, like areal extent and water equivalent, over larger areas. Data assimilation provides a framework for optimally merging information from remotely sensed observations and hydrologic model predictions. An ensemble Kalman filter (EnKF) was used to assimilate remotely sensed snow observations into the variable infiltration capacity (VIC) macroscale hydrologic model over the Snake River basin. The snow cover extent (SCE) product from the moderate resolution imaging spectroradiometer (MODIS) flown on the NASA Terra satellite was used to update VIC snow water equivalent (SWE), for a period of four consecutive winters (1999–2003). A simple snow depletion curve model was used for the necessary SWE–SCE inversion. The results showed that the EnKF is an effective and operationally feasible solution; the filter successfully updated model SCE predictions to better agree with the MODIS observations and ground surface measurements. Comparisons of the VIC SWE estimates following updating with surface SWE observations (from the NRCS SNOTEL network) indicated that the filter performance was a modest improvement over the open-loop (un-updated) simulations. This improvement was more evident for lower to middle elevations, and during snowmelt, while during accumulation the filter and open-loop estimates were very close on average. Subsequently, a preliminary assessment of the potential for assimilating the SWE product from the advanced microwave scanning radiometer (AMSR-E, flown on board the NASA Aqua satellite) was conducted. The results were not encouraging, and appeared to reflect large errors in the AMSR-E SWE product, which were also apparent in comparisons with SNOTEL data.     

Seismicity characteristics of a potentially active Quaternary volcano: The Tatun Volcano Group, northern Taiwan

The Tatun Volcano Group (TVG) is located at the northern tip of Taiwan, near the capital Taipei and close to two nuclear power plants. Because of lack of any activity in historical times it has been classified as an extinct volcano, even though more recent studies suggest that TVG might have been active during the last 20 ka. In May 2003 a seismic monitoring project at the TVG area was initiated by deploying eight three-component seismic stations some of them equipped with both short-period and broadband sensors. During the 18 months observation period local seismicity mainly consisted of high frequency earthquakes either occurring as isolated events, or as a continuous sequence in the form of spasmodic bursts. Mixed and low frequency events were also present during the same period, even though they occurred only rarely. Arrival times from events with clear P-/S-wave phases were inverted in order to obtain a minimum 1D velocity model with station corrections. Probabilistic nonlinear earthquake locations were calculated for all these events using the newly derived velocity model. Most high frequency seismicity appeared to be concentrated near the areas of hydrothermal activity, forming tight clusters at depths shallower than 4 km. Relative locations, calculated using the double-difference method and utilising catalogue and cross-correlation differential traveltimes, showed insignificant differences when compared to the nonlinear probabilistic locations. In general, seismicity in the TVG area seems to be primarily driven by circulation of hydrothermal fluids as indicated by the occurrence of spasmodic bursts, mixed/low frequency events and a b-value (1.17 ± 0.1) higher than in any other part of Taiwan. These observations, that are similar to those reported in other dormant Quaternary volcanoes, indicate that a magma chamber may still exist beneath TVG and that a future eruption or period of unrest should not be considered unlikely.     

Normal faults associated with volcanic activity arc

Volcanic centers (volcanoes, fumaroles or solfatara fields), epicenters of strong shallow earthquakes (with focal depths up to 20 km) and epicenters of intermediate depth strong earthquakes (with focal depths between 120 and 160 km) in the southern Aegean volcanic arc can be grouped into five, well defined, linear clusters trending about N60°E. This lineation of shallow earthquakes and volcanic activity is attributed to five corresponding normal faults which are named after the five corresponding volcanic centers (Sousaki, Methana, Milos, Santorini and Nisyros). This is supported by a similar trend of the geomorphological features (grabens and islands) and of geophysical features (Bouguer anomalies), as well as by other seismological data (fault plane solutions and the origins of tsunamis) and geological information on the Santorini caldera. The greater volcanic activity in the eastern volcanic centers (Santorini and Nysiros) compared to the western volcanic centers (Sousaki, Methana and Milos) is attributed to the higher rate of extensional crustal deformation. In the eastern part of the volcanic arc it is 26 mm/yr: in the west it is 2 mm/yr. The delineation of the epicenters of the intermediate depth earthquakes along the same five lines indicates the existence of five corresponding rupture zones in the lower (leading) part of the descending lithospheric slab (at depths of 120–180 km). These deep zones are probably the sources of hot material which is ascending vertically upwards and intrudes into the crust along its fracture zones. The orientation of these zones explains the focusing of the macro-seismic results of these deep shocks at narrow regions of the sedimentary arc (Peloponnesus, Crete, etc).