Laurent Longuevergne

Extracting coherent regional information from local measurements with Karhunen‐Loève transform: Case study of an alluvial aquifer (Rhine valley, France and Germany)

We investigate the ability of combining the Karhunen-Loeve transform (KLT) with the kriging method to extract regional information from a set of point measurements. This method was applied to a set of 195 piezometric head time series over a period of 17 years from observation wells distributed within the French and German area of the Rhine valley alluvial groundwater body. Piezometric head time series are analyzed with KLT in order to highlight characteristic temporal signals, classified from the most energetic (global) to the least energetic (local) signals. The first five signals amount to 80% of the global variance of the system and are inferred to represent different hydrological contributions (exchanges with rivers and rainfall), but they also represent a significant anthropogenic component. Kriging is then used to regionalize the signals and to build a reconstruction model of the behavior of the whole aquifer containing only filtered information coming from identified source signals.

2D characterization of near-surface V P/V S: surface-wave dispersion inversion versus refraction tomography

The joint study of pressure (P-) and shear (S-) wave velocities (Vp and Vs ), as well as their ratio (Vp /Vs), has been used for many years at large scales but remains marginal in near-surface applications. For these applications, and are generally retrieved with seismic refraction tomography combining P and SH (shear-horizontal) waves, thus requiring two separate acquisitions. Surface-wave prospecting methods are proposed here as an alternative to SH-wave tomography in order to retrieve pseudo-2D Vs sections from typical P-wave shot gathers and assess the applicability of combined P-wave refraction tomography and surface-wave dispersion analysis to estimate Vp/Vs ratio. We carried out a simultaneous P- and surface-wave survey on a well-characterized granite-micaschists contact at Ploemeur hydrological observatory (France), supplemented with an SH-wave acquisition along the same line in order to compare Vs results obtained from SH-wave refraction tomography and surface-wave profiling. Travel-time tomography was performed with P- and SH- wave first arrivals observed along the line to retrieve Vtomo p and Vtomo s models. Windowing and stacking techniques were then used to extract evenly spaced dispersion data from P-wave shot gathers along the line. Successive 1D Monte Carlo inversions of these dispersion data were performed using fixed Vp values extracted from Vtomo p the model and no lateral constraints between two adjacent 1D inversions. The resulting 1D Vsw s models were then assembled to create a pseudo-2D Vsw s section, which appears to be correctly matching the general features observed on the section. If the pseudo-section is characterized by strong velocity incertainties in the deepest layers, it provides a more detailed description of the lateral variations in the shallow layers. Theoretical dispersion curves were also computed along the line with both and models. While the dispersion curves computed from models provide results consistent with the coherent maxima observed on dispersion images, dispersion curves computed from models are generally not fitting the observed propagation modes at low frequency. Surface-wave analysis could therefore improve models both in terms of reliability and ability to describe lateral variations. Finally, we were able to compute / sections from both and models. The two sections present similar features, but the section obtained from shows a higher lateral resolution and is consistent with the features observed on electrical resistivity tomography, thus validating our approach for retrieving Vp/Vs ratio from combined P-wave tomography and surface-wave profiling.

Physical Modelling To Remove Hydrological Effects At Local And Regional Scale: Application To The 100-M Hydrostatic Inclinometer In Sainte-Croix-Aux-Mines (France)

New inclinometers devoted to hydrological studies were set up in the Vosges Mountains (France). Two orthogonal 100-m base hydrostatic inclinometers were installed in December 2004 as well as a hydrometeorological monitoring system for the 100-km2 hydrological unit around the inclinometer. As inclinometers are very sensitive to environmental influences, this observatory is a test site to confront hydrological modelling and geodetic observations. Physical modelling to remove hydrological effects without calibrating on geodetic data is tested on these instruments. Specifically, two deformation processes are most important: fluid pressure variations in nearby hydraulically active fractures and surface loading at regional scale

Temporal Variations of Near-surface Seismic Data at the Ploemeur (France) Hydrogeological Observatory

Near-surface seismic methods are mainly used to determine the geometrical characteristics of hydrosystems (and to provide elements that are interesting for hydrogeologists such as separating aquifer layers, setting up systems boundaries, highlighting fractures etc.). Recent methodological advances suggest the high potential of seismic methods to investigate the mechanical properties of the Critical Zone (CZ), by exploiting the full wealth of seismic records. Indeed, the behavior of Shear (S) and Pressure (P) waves in the presence of water is partially decoupled, so that the ratio of their propagation velocities VP/VS is strongly linked to water saturation. We propose here a time-lapse application of this approach. Two seismic acquisitions were carried out under distinct hydrogeological conditions along the same line at the Ploemeur hydrogeological observatory (South Brittany, France). Vertical component seismic data were recorded to extract: (i) P-wave first arrival times and (ii) Rayleigh-wave phase velocities. The significant variations with time and space, of both datasets, indicate marked changes in mechanical properties of the CZ that have to be compared to soil moisture variations in the unsaturated zone and groundwater level variations.

Contribution of Seismic Methods to Hydrogeophysics

The characterisation and monitoring of aquifer systems mainly rely on piezometric and log data. Delineating spatial variations of lithology between piezometers is a delicate task, which inevitably generates errors possibly propagating into hydrogeological models. Seismic methods have been proposed to: (i) improve the low spatial resolution of borehole data, (ii) provide a characterisation of the subsurface geometry, and (iii) estimate the physical parameters of the medium influenced by the presence of water and the associated flow and transport processes. The joint study of pressure (P-) and shear (S-) wave seismic velocities (VP and VS, respectively), whose evolution is strongly decoupled in the presence of fluid, has been proposed through the estimation of the VP/VS ratio and Poissons ratio. A specific methodology has been developed for the combined exploitation of P- and surface waves present on single seismic records. The use of this methodology in several geological and hydrogeological contexts allowed for estimating VP/VS ratio lateral and temporal variations in good agreement with a priori geological information and existing geophysical and piezometric data. Laser-based ultrasonic techniques were also proposed to put these processing techniques in practice on perfectly controlled physical models and study elastic wave propagation in partially saturated porous media.

Surface-wave Dispersion Stacking on a Granite- micaschists Contact at Ploemeur Hydrological Observatory, France

In the context of a geophysical survey at the Ploemeur hydrological observatory (France), we performed surface-wave profiling for the characterisation of shallow subsurface Shear-wave velocities. Since we anticipated lateral variations but needed great investigation depth, we deployed multifold acquisition geometries and used roll-along dispersion stacking to enable efficient measurements of multi-modal dispersion data. Several offset moving windows have been tested. Represented as pseudo-sections, the phase velocities extracted using a 12-trace window clearly showed three areas coherent with field observation and interestingly consistent with electrical conductivities and P-wave first arrival times. This cross-quality control has been of great help in the choice of the moving window size and revealed itself to be a rewarding step prior to the inversion process.