Low-frequency geoelectrical methods include mainly self-potential, resistivity, and induced polarization techniques, which have potential in many environmental and hydrogeological applications. They provide complementary information to each other and to in-situ measurements. The self-potential method is a passive measurement of the electrical response associated with the in-situ generation of electrical current due to the flow of pore water in porous media, a salinity gradient, and/or the concentration of redoxactive species. Under some conditions, this method can be used to visualize groundwater flow, to determine permeability, and to detect preferential flow paths. Electrical resistivity is dependent on the water content, the temperature, the salinity of the pore water, and the clay content and mineralogy. Time-lapse resistivity can be used to assess the permeability and dispersivity distributions and to monitor contaminant plumes. Induced polarization characterizes the ability of rocks to reversibly store electrical energy. It can be used to image permeability and to monitor chemistry of the pore water-minerals interface. These geophysical methods, reviewed in this paper, should always be used in concert with additional in-situ measurements (e.g. in-situ pumping tests, chemical measurements of the pore water), for instance through joint inversion schemes, which is an area of fertile on-going research.

The complex conductivity of soils remains poorly known despite the growing importance of this method in hydrogeophysics. In order to fill this gap of knowledge, we investigate the complex conductivity of 71 soils samples (including four peat samples) and one clean sand in the frequency range 0.1 Hz to 45 kHz. The soil samples are saturated with six different NaCl brines with conductivities (0.031, 0.53, 1.15, 5.7, 14.7, and 22 S m−1, NaCl, 25°C) in order to determine their intrinsic formation factor and surface conductivity. This data set is used to test the predictions of the dynamic Stern polarization model of porous media in terms of relationship between the quadrature conductivity and the surface conductivity. We also investigate the relationship between the normalized chargeability (the difference of in‐phase conductivity between two frequencies) and the quadrature conductivity at the geometric mean frequency. This data set confirms the relationships between the surface conductivity, the quadrature conductivity, and the normalized chargeability. The normalized chargeability depends linearly on the cation exchange capacity and specific surface area while the chargeability shows no dependence on these parameters. These new data and the dynamic Stern layer polarization model are observed to be mutually consistent. Traditionally, in hydrogeophysics, surface conductivity is neglected in the analysis of resistivity data. The relationships we have developed can be used in field conditions to avoid neglecting surface conductivity in the interpretation of DC resistivity tomograms. We also investigate the effects of temperature and saturation and, here again, the dynamic Stern layer predictions and the experimental observations are mutually consistent.

Seawater intrusion has often resulted in scarce fresh groundwater resources in coastal lowlands. Careful management is essential to avoid the overexploitation of these vulnerable fresh groundwater resources, requiring detailed information on their spatial occurrence. Airborne electromagnetics (EM) has proved a valuable tool for efficient mapping of ground conductivity, as a proxy for fresh groundwater resources. Stakeholders are, however, interested in groundwater salinity, necessitating a translation of ground conductivity to groundwater salinity. This paper presents a methodology to construct a high-resolution (50 × 50 × 0.5 m 3 ) 3D voxel model of groundwater chloride concentration probability, based on a large-scale (1800 km 2 , 9640 flight line kilometres) airborne EM survey in the province of Zeeland, the Netherlands. Groundwater chloride concentration was obtained by combining pedotransfer functions with detailed lithological information. The methodology includes a Monte Carlo based forward uncertainty propagation approach to quantify the inherent uncertainty in the different steps. Validation showed good correspondence both with available groundwater chloride analyses, and with ground-based hydrogeophysical measurements. Our results show the limited occurrence of fresh groundwater in Zeeland, as 75% of the area lacks fresh groundwater within 15 m below ground surface. Fresh groundwater is mainly limited to the dune area and sandy creek ridges. In addition, significant fresh groundwater resources were shown to exist below saline groundwater, where infiltration of seawater during marine transgressions was hindered by the presence of clayey aquitards. The considerable uncertainty in our results highlights the importance of applying uncertainty analysis in airborne EM surveys. Uncertainty in our results mainly originated from the inversion and the 3D interpolation, and was largest at transition zones between fresh and saline groundwater. Reporting groundwater salinity instead of ground conductivity facilitated the rapid uptake of our results by relevant stakeholders, thereby supporting the necessary management of fresh groundwater resources in the region.

International audienceInduced polarization (more precisely the magnitude and phase of impedance of the subsurface) is measured using a network of electrodes located at the ground surface or in boreholes. This method yields important information related to the distribution of permeability and contaminants in the shallow subsurface. We propose a new time-lapse 3-D modelling and inversion algorithm to image the evolution of complex conductivity over time. We discretize the subsurface using hexahedron cells. Each cell is assigned a complex resistivity or conductivity value. Using the finite-element approach, we model the in-phase and out-of-phase (quadrature) electrical potentials on the 3-D grid, which are then transformed into apparent complex resistivity. Inhomogeneous Dirichlet boundary conditions are used at the boundary of the domain. The calculation of the Jacobian matrix is based on the principles of reciprocity. The goal of time-lapse inversion is to determine the change in the complex resistivity of each cell of the spatial grid as a function of time. Each model along the time axis is called a 'reference space model'. This approach can be simplified into an inverse problem looking for the optimum of several reference space models using the approximation that the material properties vary linearly in time between two subsequent reference models. Regularizations in both space domain and time domain reduce inversion artefacts and improve the stability of the inversion problem. In addition, the use of the time-lapse equations allows the simultaneous inversion of data obtained at different times in just one inversion step (4-D inversion). The advantages of this new inversion algorithm are demonstrated on synthetic time-lapse data resulting from the simulation of a salt tracer test in a heterogeneous random material described by an anisotropic semi-variogram

No abstract