Contribution of geophysical surveys to groundwater modelling of a porous aquifer in semiarid Niger: An overview

Subsurface geophysical surveys were carried out using a large range of methods in an unconfined sandstone aquifer in semiarid south-western Niger for improving both the conceptual model of water flow through the unsaturated zone and the parameterization of numerical a groundwater model of the aquifer. Methods included: electromagnetic mapping, electrical resistivity tomography (ERT), resistivity logging, time domain electromagnetic sounding (TDEM), and magnetic resonance sounding (MRS). Analyses of electrical conductivities, complemented by geochemical measurements, allowed us to identify preferential pathways for infiltration and drainage beneath gullies and alluvial fans. The mean water content estimated by MRS (13%) was used for computing the regional groundwater recharge from long-term change in the water table. The ranges in permeability and water content obtained with MRS allowed a reduction of the degree of freedom of aquifer parameters used in groundwater modelling.     

Contribution of magnetic resonance sounding to aquifer characterization and recharge estimate in semiarid Niger

To improve the knowledge of the regionally important Continental Terminal 3 (CT3) aquifer in south-western Niger, fifteen magnetic resonance soundings (MRS) were carried out in December 2005 in the vicinity of wells and boreholes. The output MRS geophysical parameters, i.e. water content and decay constants versus depth, were compared to hydrogeological characteristics, i.e. water table depth, total porosity, specific yield and transmissivity estimated from direct measurements, pumping tests and transient groundwater modelling. The MRS-determined parameters were then used to estimate the rates of groundwater recharge.Contained in poorly consolidated Tertiary sandstones, the CT3 aquifer's water table has continuously risen by 4 m in total over the past four decades. Additionally, a significant portion of this increase has occurred in the past decade alone, with an annual rise now ranging between 0.1 and 0.3 m depending on the monitored well. Increase in groundwater recharge due to land clearance and deforestation explains this situation. According to previous estimations, the pre-clearing recharge ranged from 1 to 5 mm per year in 1950–60 s, while more recent recharge rates (1990s–2000s) range from 20 to 50 mm per year. These recharge values are directly affected by estimated aquifer specific yield value, while the spatial variation of rates of water table rise can be attributed to large scale hydrodynamic heterogeneities in the aquifer. However, few field measurements were available to confirm these assumptions.The main results of this study are: (1) The water table depth and aquifer transmissivity are estimated from MRS output parameters with an average accuracy of ± 10% and ± 9% respectively. (2) The MRS-determined water content is linked to both the total porosity and the specific yield of the aquifer, but no quantitative formulation can be proposed as yet. (3) Using the average MRS-determined water content over the investigated area, i.e. 13%, the groundwater recharge rates can be estimated to be ~ 2 mm per year in the 1950–1960s (pre-clearing period), and ~ 23 mm per year for the last decade. (4) The variations in specific yield and transmissivity cannot explain by themselves the spatial variability of the rise of the water table. (5) The ranges in transmissivity and water content obtained from MRS are more realistic than the groundwater modelling outputs. Therefore, MRS could be used to better constrain the aquifer parameters in groundwater modelling with a dense site network.Finally, this work illustrates how MRS can successfully improve characterisation and transient multi-year groundwater balance of commonly found sedimentary aquifers, particularly when integrated with well observations and pumping tests.     

Characterising groundwater and surface-water interconnections using hydrogeology,hydrochemistry and stable isotopes in the Ouémé Delta,southern Benin

The Mio-Pliocene aquifer of the coastal sedimentary basin of Benin is the most exploited aquifer for water supply to the urbanised region in the southern part of the country. The population explosion is putting increasing pressure on quantitative and qualitative aspects of the groundwater resources. Preventing groundwater contamination caused by surface waters requires a thorough understanding of surface-water/groundwater interactions, especially the interactions between the Mio-Pliocene aquifer and surface waters. This study aimed to investigate the interactions between groundwater and surface waters along the major rivers (Sô River and Ouémé Stream) and brooks in the Ouémé Delta. Field campaigns identified 75 springs located in the valleys which feed the rivers, and thus maintain their base flow. The piezometric results indicated, through flow direction assessment, that the Mio-Pliocene aquifer feeds Ouémé Stream and Sô River. Chemical analyses of groundwater and surface waters show similar chemical facies, and changes in the chemical composition in groundwater are also observed in the surface waters. Moreover, the isotopic signatures of surface waters are similar to those of the groundwater and springs, which led to the identification of potential groundwater discharge areas. As a result of groundwater discharge into surface waters, the fraction of groundwater in the surface water is more than 66% in the brooks, regardless of the season. In the Ouémé Stream and Sô River, the fraction of groundwater is 0–21% between June and September, while from October to March it is 47–100%.

     

Multi-scale imaging of a slow active fault zone: contribution for improved seismic hazard assessment in the Swiss Alpine foreland

Seismic hazard assessment of slow active fault zones is challenging as usually only a few decades of sparse instrumental seismic monitoring is available to characterize seismic activity. Tectonic features linked to the observed seismicity can be mapped by seismic imaging techniques and/or geomorphological and structural evidences. In this study, we investigate a seismic lineament located in the Swiss Alpine foreland, which was discussed in previous work as being related to crustal structures carrying in size the potential of a magnitude M 6 earthquake. New, low-magnitude (?2.0 ≤ M_L ≤ 2.5) earthquake data are used to image the spatial and temporal distribution of seismogenic features in the target area. Quantitative and qualitative analyses are applied to the waveform dataset to better constrain earthquakes distribution and source processes. Potential tectonic features responsible for the observed seismicity are modelled based on new reinterpretations of oil industry seismic profiles and recent field data in the study area. The earthquake and tectonic datasets are then integrated in a 3D model. Spatially, the seismicity correlates over 10–15 km with a N–S oriented sub-vertical fault zone imaged in seismic profiles in the Mesozoic cover units above a major decollement on top of the mechanically more rigid basement and seen in outcrops of Tertiary series east of the city of Fribourg. Observed earthquakes cluster at shallow depth (<4 km) in the sedimentary cover. Given the spatial extend of the observed seismicity, we infer the potential of a moderate size earthquake to be generated on the lineament. However, since the existence of along strike structures in the basement cannot be excluded, a maximum M 6 earthquake cannot be ruled out. Thus, the Fribourg Lineament constitutes a non-negligible source of seismic hazard in the Swiss Alpine foreland.     

Magnetic resonance sounding applied to aquifer characterization

Magnetic resonance sounding (MRS) is distinguished from other geophysical tools used for ground water investigation by the fact that it measures a magnetic resonance signal generated directly from subsurface water molecules. An alternating current pulse energizes a wire loop on the ground surface and the MRS signal is generated; subsurface water is indicated, with a high degree of reliability, by nonzero amplitude readings. Measurements with varied pulse magnitudes then reveal the depth and thickness of water saturated layers. The hydraulic conductivity of aquifers can also be estimated using boreholes for calibration. MRS can be used for both predicting the yield of water supply wells and for interpolation between boreholes, thereby reducing the number of holes required for hydrogeological modeling. An example of the practical application of MRS combined with two-dimensional electrical imaging, in the Kerbernez and Kerien catchments area of France, demonstrates the efficiency of the technique.     

Estimation of anisotropic hydraulic conductivity using geophysical data in a coastal aquifer of Karnataka,India

Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models.     

Estimating aquifer thickness using multiple pumping tests

A method to estimate aquifer thickness and hydraulic conductivity has been developed, consisting of multiple pumping tests. The method requires short-duration pumping cycles on an unconfined aquifer with significant seasonal water-table fluctuations. The interpretation of several pumping tests at a site in India under various initial conditions provides information on the change in hydrodynamic parameters in relation to the initial water-table level. The transmissivity linearly decreases compared with the initial water level, suggesting a homogeneous distribution of hydraulic conductivity with depth. The hydraulic conductivity is estimated from the slope of this linear relationship. The extrapolation of the relationship between transmissivity and water level provides an estimate of the aquifer thickness that is in good agreement with geophysical investigations. The hydraulically active part of the aquifer is located in both the shallow weathered and the underlying densely fractured zones of the crystalline basement. However, no significant relationship is found between the aquifer storage coefficient and initial water level. This new method contributes to filling the methodological gap between single pumping tests and hydraulic tomography, in providing information on the variation of the global transmissivity according to depth. It can be applied to any unconfined aquifer experiencing large seasonal water-table fluctuations and short pumping cycles.     

Estimating specific yield and transmissivity with magnetic resonance sounding in an unconfined sandstone aquifer (Niger)

The unconfined aquifer of the Continental Terminal in Niger was investigated by magnetic resonance sounding (MRS) and by 14 pumping tests in order to improve calibration of MRS outputs at field scale. The reliability of the standard relationship used for estimating aquifer transmissivity by MRS was checked; it was found that the parametric factor can be estimated with an uncertainty ≤150% by a single point of calibration. The MRS water content (θ _MRS) was shown to be positively correlated with the specific yield (Sy), and θ _MRS always displayed higher values than Sy. A conceptual model was subsequently developed, based on estimated changes of the total porosity, Sy, and the specific retention Sr as a function of the median grain size. The resulting relationship between θ _MRS and Sy showed a reasonably good fit with the experimental dataset, considering the inherent heterogeneity of the aquifer matrix (residual error is ～60%). Interpreted in terms of aquifer parameters, MRS data suggest a log-normal distribution of the permeability and a one-sided Gaussian distribution of Sy. These results demonstrate the efficiency of the MRS method for fast and low-cost prospection of hydraulic parameters for large unconfined aquifers.