Resolution of freshwater and saline water aquifers by composite geophysical data analysis methods

Abstract

The resolution of the freshwater and saline water aquifers in a coastal terrain (Mahanadi Basin, India) is updated. We analysed electrical borehole log data at four sites and compared the water resistivity regime of the freshwater and saline water zones obtained from electrical borehole logging, with the resistivity regime obtained by interpreting vertical electrical sounding (VES) data. The multilayer VES data interpretation is modified to a simple model, containing only the freshwater zone and the saline water zone. The composite geophysical parameters of the freshwater and saline water zones, in particular the resistivity and longitudinal unit conductance regime, are identified. The resolution obtained from the composite geophysical data analyses is very clear and convincing. The composite longitudinal unit conductance regime of the saline water zones is very high compared to that of the freshwater zones. This makes the identification of the two aquifers easy and increases its reliability. A technique which enables analysis of composite geophysical data of freshwater and saline water zones at VES sites in the vicinity of the borehole log sites is proposed. The significance of longitudinal unit conductance in resolving the freshwater and saline water aquifers is illustrated graphically. The proposed technique is validated by correlating the longitudinal unit conductance and resistivity with the total dissolved solids. The efficiency of the technique is validated by carrying out discriminant function analysis.

Citation Hodlur, G. K., Dhakate, R., Sirisha, T. & Panaskar, D. B. (2010) Resolution of freshwater and saline water aquifers by composite geophysical data analysis methods. Hydrol. Sci. J. 55(3), 414–434.     

Groundwater Microbial Communities Along a Generalized Flowpath in Nomhon Area,Qaidam Basin,China

Spatial distribution (horizonal and vertical) of groundwater microbial communities and the hydrogeochemistry in confined aquifers were studied approximately along the groundwater flow path from coteau to plain in the Nomhon area, Qinghai‐Tibet plateau, China. The confined groundwater samples at different depths and locations were collected in three boreholes through a hydrogeological section in this arid and semi‐arid area. The phylogenetic analysis of 16S rRNA genes and multivariate statistical analysis were used to elucidate similarities and differences between groundwater microbial communities and hydrogeochemical properties. The integrated isotopic geochemical measurements were applied to estimate the source and recharge characteristics of groundwater. The results showed that groundwater varied from fresh to saline water, and modern water to ancient water following the flowpath. The recharge characteristics of the saline water was distinct with that of fresh water. Cell abundance did not vary greatly along the hydrogeochemical zonality; however, dissimilarities in habitat‐based microbial community structures were evident, changing from Betaproteobacteria in the apex of alluvial fan to Gammaproteobacteria and then to Epsilonproteobacteria in the core of the basin (alluvial‐lacustrine plain). Rhodoferax, Hydrogenophaga, Pseudomonas, and bacterium isolated from similar habitats unevenly thrived in the spatially distinct fresh water environments, while Sulfurimonas dominanted in the saline water environment. The microbial communities presented likely reflected to the hydrogeochemical similarities and zonalities along groundwater flowpath.     

Modeling Coastal Salinity in Quasi 2D and 3D Using a DUALEM‐421 and Inversion Software

Rising sea levels, owing to climate change, are a threat to fresh water coastal aquifers. This is because saline intrusions are caused by increases and intensification of medium‐large scale influences including sea level rise, wave climate, tidal cycles, and shifts in beach morphology. Methods are therefore required to understand the dynamics of these interactions. While traditional borehole and galvanic contact resistivity (GCR) techniques have been successful they are time‐consuming. Alternatively, frequency‐domain electromagnetic (FEM) induction is potentially useful as physical contact with the ground is not required. A DUALEM‐421 and EM4Soil inversion software package are used to develop a quasi two‐ (2D) and quasi three‐dimensional (3D) electromagnetic conductivity images (EMCI) across Long Reef Beach located north of Sydney Harbour, New South Wales, Australia. The quasi 2D models discern: the dry sand (<10 mS/m) associated with the incipient dune; sand with fresh water (10 to 20 mS/m); mixing of fresh and saline water (20 to 500 mS/m), and; saline sand of varying moisture (more than 500 mS/m). The quasi 3D EMCIs generated for low and high tides suggest that daily tidal cycles do not have a significant effect on local groundwater salinity. Instead, the saline intrusion is most likely influenced by medium‐large scale drivers including local wave climate and morphology along this wave‐dominated beach. Further research is required to elucidate the influence of spring‐neap tidal cycles, contrasting beach morphological states and sea level rise.     

Electromagnetic‐Induction Logging to Monitor Changing Chloride Concentrations

Water from the San Joaquin Delta, having chloride concentrations up to 3590 mg/L, has intruded fresh water aquifers underlying Stockton, California. Changes in chloride concentrations at depth within these aquifers were evaluated using sequential electromagnetic (EM) induction logs collected during 2004 through 2007 at seven multiple‐well sites as deep as 268 m. Sequential EM logging is useful for identifying changes in groundwater quality through polyvinyl chloride‐cased wells in intervals not screened by wells. These unscreened intervals represent more than 90% of the aquifer at the sites studied. Sequential EM logging suggested degrading groundwater quality in numerous thin intervals, typically between 1 and 7 m in thickness, especially in the northern part of the study area. Some of these intervals were unscreened by wells, and would not have been identified by traditional groundwater sample collection. Sequential logging also identified intervals with improving water quality—possibly due to groundwater management practices that have limited pumping and promoted artificial recharge. EM resistivity was correlated with chloride concentrations in sampled wells and in water from core material. Natural gamma log data were used to account for the effect of aquifer lithology on EM resistivity. Results of this study show that a sequential EM logging is useful for identifying and monitoring the movement of high‐chloride water, having lower salinities and chloride concentrations than sea water, in aquifer intervals not screened by wells, and that increases in chloride in water from wells in the area are consistent with high‐chloride water originating from the San Joaquin Delta rather than from the underlying saline aquifer.     

Seawater intrusion and fresh groundwater hydraulics in fjord delta aquifers inferred from ground penetrating radar and resistivity profiles—Sunndalsøra and Esebotn, western Norway

Geophysical surveys have been carried out in two fjord delta aquifers in western Norway. Geophysical methods comprised ground penetrating radar (GPR), shallow reflection seismic and resistivity measurements. The purpose of the investigations was to evaluate hydrogeological properties of the fjord delta aquifers with emphasis on the possibilities of abstracting saline groundwater for use in fish farming. At Sunndalsøra, reflection seismic was helpful in the mapping of the base of the aquifer. Resistivity profiles revealed both the spatial and seasonal variability of seawater intrusion. Penetration depths of ground penetrating radar (GPR-D) revealed the spatial variation of the boundary of the transition zone between saline and fresh groundwater. Maps of GPR-D can be used to indicate direction and magnitude of mean hydraulic gradients and thus to locate optimal sites for both fresh and saline groundwater abstraction.In fjord delta aquifers, both rivers and abandoned river channels constitute important groundwater divides. Optimal locations for saline groundwater abstraction wells are in areas separated from the rivers by abandoned river channels. In areas between the rivers and the abandoned river channels, groundwater abstraction will result in decreasing salinity of the water.     

Geoelectrical and hydrochemical investigations for characterizing the salt water intrusion in the Khanasser valley, northern Syria

An integrated approach of geoelectrical and hydrochemical investigation surveys was proposed for indicating contact regions between saline and fresh groundwater in the Khanasser valley region, northern Syria. The qualitative and quantitative interpretations of 34 vertical electrical soundings (VES) enable to characterize the salt water intrusion laterally and vertically. The established iso-apparent resistivity maps for different AB/2 spacings obviously indicate the presence of a lowresistivity (less than 4 Ohm·m) zone related to the salt water intrusion in the Quaternary and Paleogene deposits. The different hydrochemical and geophysical parameters, such as electrical resistivity, total dissolved solids (TDS) and major ions concentrations used to characterize the salt water intrusion gave almost similar results in locating and mapping the different boundaries of the groundwater salinity. The proposed approach is useful for mapping the interface between different groundwater qualities, and can be therefore used to successfully characterize the salt water intrusion phenomenon in other semi-arid regions. The application of such an approach is a powerful tool and can be used for water resource management in the water scarce areas.     

Distribution of fossil seawater and its role in Neogene sedimentary rock landslides in Niigata,eastern marginal region of the Japan Sea

The Neogene marine sedimentary rock area in the eastern marginal region of the Japan Sea is an area with some of the highest landslide densities in Japan. Some of the landslides in this area have been known to involve saline groundwater, which can be the cause of these landslides. In order to demonstrate the relationships between landslides and saline water, topographic, geological, groundwater, and electromagnetic surveys were performed in the eastern marginal region of the Japan Sea. Many landslides and gravitational slope deformations with linear depressions and small scarps were recognized in the study area. The resistivity profile obtained by an electromagnetic survey suggests that there is a wide zonal distribution of saline water with salt concentrations equivalent to seawater at depths of 50–100 m or more and that the groundwater shallower than 50 m has an electrical conductivity of less than 100 mS/m. The shallow resistive groundwater is inferred to be meteoric water that replaced the saline groundwater, which likely weakened the bedrock, resulting in landslides. A ridge of competent tuff overlying mudstone has many linear depressions from gravitational slope deformation and low‐resistivity water to a depth of 600 m, which suggests that the mudstone was weakened by water replacement and deformed under the tuff caprock. The saline groundwater is inferred to be fossil seawater trapped in pores during sediment deposition, which is brought near the ground surface along with rocks by tectonic movement in the hills. Thus, the saline water and its fresh water replacement are among the important basic causes of the landslides. The oil well data obtained in the eastern marginal region of the Japan Sea suggest that such saline water replacement has occurred widely and that replacement is likely one of the predispositions for the frequent landslides there.     

Interpretation of magnetic resonance soundings in rocks with high electrical conductivity

Magnetic resonance sounding (MRS) is an electromagnetic method designed for groundwater investigations. MRS can be applied not only for studying fresh-water aquifers, but also in areas where intrusion of saline water is rendering the subsurface electrically conductive. In the presence of rocks with a high electrical-conductivity attenuation and a phase shift of the MRS signal may influence the efficiency of the MRS method. We investigated the performance of MRS for allowing us to propose a procedure for interpreting MRS data under these conditions. For numerical modeling, we considered a subsurface with a resistivity between 0.5 and 10 Ω m. The results show that the depth of investigation with MRS depends upon the electrical conductivity of groundwater and surrounding rocks, on the depth of the saline water layer, and on the amount of fresh water above the saline water. For interpreting MRS measurements, the electrical conductivity of the subsurface is routinely measured with an electrical or electromagnetic method. However, due to the equivalence problem, the result obtained with these methods may be not unique. Hence, we investigated the influence of the uncertainty in conductivity distribution provided by transient electromagnetic measurements (TEM) on MRS results. It was found that the uncertainty in TEM results has an insignificant effect on MRS.     

Geoelectrical structure and hydrogeological investigations of the southern Rif Cordillera (Morocco)

An electrical prospecting survey is conducted in the Rharb basin, a semi‐arid region in the southern part of the Rifean Cordillera (Morocco) to delineate characteristics of the aquifer and the groundwater affected by the marine intrusion related to Atlantic Ocean. Analysis and interpretations of electrical soundings, bi‐logarithmic diagrams and the geoelectrical sections highlight a monolayer aquifer in the southern part, a multilayer system in the northern part of the Rharb basin and lenticular semi‐permeable formations. Several electrical layers have been deduced from the analysis of bi‐logarithmic diagrams: resistant superficial level (R₀), conducting superficial level (C₀), resistant level (R), intermediary resistant level (R′), conducting level (Cp) and intermediary layer of resistivity (AT). Spatial distribution of the resistivity deduced from the interpretation of apparent resistivity maps (AB = 400 and 1000 m) and the decreasing of resistivity values (35–10 Ωm), in particular in the coastal zone show that this heterogeneity is related to several anomalies identified in the coastal area, which result from hydraulic and geological processes: (i) heterogeneous hydraulic conductivity in particular in the southern part of the Rharb; (ii) lateral facies and synsedimentary faulting and (iii) the relationship between the electrical conductivity and chloride concentration of groundwater shows that salinity is the most important factor controlling resistivity. The distribution of fresh/salt‐water zones and their variations in space along geoelectrical sections are established through converting subsurface depth‐resistivity models. Copyright © 2010 John Wiley & Sons, Ltd.     

The Interface Configuration of the Fresh‐/Dead Sea Water – Theory and Measurements

The high‐density Dead Sea water (1.235 g/cm³) forms a special interface configuration with the fresh groundwater resources of its surrounding aquifers. The fresh groundwater column beneath its surroundings is around one tenth of its length compared to oceanic water. This fact alone indicates the vulnerability of the fresh groundwater resources to the impacts of changes in the Dead Sea level and to saltwater migration. Ghyben‐Herzberg and Glover equations were used to calculate the volumes of water in coastal aquifers which were replaced by freshwater due to the interface seaward migration as a result of the drop in the level of the Dead Sea. For that purpose, the dynamic equation of Glover approach has been integrated to accommodate that type of interface readjustment. The calculated amounts of freshwater which substituted salt Dead Sea water due to the migration of interface are 3.21 · 10¹¹ m³, from a Dead Sea level of –392 m to τ411 m below sea level. The average porosity of coastal aquifers was calculated to range from 2.8 to 2.94%. Geoelectric sounding measurements showed that areas underlying the coastal aquifers formerly occupied by the Dead Sea water are gradually becoming flushed and occupied by freshwater. The latter is becoming salinized due to the residuals of Dead Sea water in the aquifer matrix, the present salinity of which is lower than that of the Dead Sea water. At the same time salt dissolution from the Lisan Marl formation is causing collapses along the shorelines in the form of sinkholes, tens of meters in diameter and depth.     

Integrated geophysical and chemical study of saline water intrusion

Surface geophysical surveys provide an effective way to image the subsurface and the ground water zone without a large number of observation wells. DC resistivity sounding generally identifies the subsurface formations-the aquifer zone as well as the formations saturated with saline/brackish water. However, the method has serious ambiguities in distinguishing the geological formations of similar resistivities such as saline sand and saline clay, or water quality such as fresh or saline, in a low resistivity formation. In order to minimize the ambiguity and ascertain the efficacy of data integration techniques in ground water and saline contamination studies, a combined geophysical survey and periodic chemical analysis of ground water were carried out employing DC resistivity profiling, resistivity sounding, and shallow seismic refraction methods. By constraining resistivity interpretation with inputs from seismic refraction and chemical analysis, the data integration study proved to be a powerful method for identification of the subsurface formations, ground water zones, the subsurface saline/brackish water zones, and the probable mode and cause of saline water intrusion in an inland aquifer. A case study presented here illustrates these principles. Resistivity sounding alone had earlier failed to identify the different formations in the saline environment. Data integration and resistivity interpretation constrained by water quality analysis led to a new concept of minimum resistivity for ground water-bearing zones, which is the optimum value of resistivity of a subsurface formation in an area below which ground water contained in it is saline/brackish and unsuitable for drinking.     

DETERMINATION OF A SALT-WATER INTERFACE BY ELECTRIC RESISTIVITY DEPTH SOUNDINGS / Détermination d'un interface d'eau salée à l'aide des mesures verticales de résistivité électrique

Abstract

Vertical electrical resistivity soundings were measured near hydrological observation wells in order to ascertain whether geophysical means could be used to map saline water intrusion into a fresh water aquifer in Israel. The soundings showed that the low resistivity layers associated with the salt water are readily discernible. The technique was applied to the entire coastal belt and resulted in a detailed study of the saline water body and its extent. Measurements were repeated six years later and good agreement between the two sets of measurements was noted. This technique is therefore judged to be an accurate tool for the mapping of salt water intrusion in fresh-water aquifers.     

Using chlorofluorocarbons (CFCs) and tritium to improve conceptual model of groundwater flow in the South Coast Aquifers of Laizhou Bay,China

The southern coastal plain of Laizhou Bay, which is the area most seriously affected by salt water intrusion in north China, is a large alluvial depression, which represents one of the most important hydrogeological units in the coastal region of northern China. Chlorofluorocarbons (CFCs, including CFC‐11, CFC‐12 and CFC‐113) and tritium were used together for dating groundwater up to 50 years old in the study area. There are two cones of depression, caused by intensive over‐exploitation of fresh groundwater in the south and brine water in the north. The assigned CFC apparent ages for shallow groundwater range from 8 a to >50 a. A binary mixing model based on CFC‐113 and CFC‐12 concentrations in groundwater was used to estimate fractions of young and pre‐modern water in shallow aquifers and to identify groundwater mixing processes during saltwater intrusion. Discordance between concentrations of different CFC compounds indicate that shallow groundwater around the Changyi cone of depression is vulnerable to contamination. Pumping activities, CFC contamination, mixing and/or a large unsaturated zone thickness (e.g. >20 m) may be reasons for some groundwater containing CFCs without tritium. Saline intrusion mainly occurs because of large head gradients between fresh groundwater in the south and saline water bodies in the north, forming a wedge of saline water below/within fresh aquifer layers. Both CFC and tritium dates indicate that the majority of the saline water is from >50 a, with little or no modern seawater component. Based on the distribution of CFC apparent ages, tritium contents plus chemical and physical data, a conceptual model of groundwater flow along the investigated Changyi‐Xiaying transect has been developed to describe the hydrogeological processes. Three regimes are identified from south to north: (i) fresh groundwater zone, with a mixing fraction of 0.80–0.65 ‘young’ water calculated with the CFC binary mixing model (groundwater ages <34 a) and 1.9–7.8TU of tritium; (ii) mixing zone characterized by a mixing fraction of 0.05–0.65 young groundwater (ages of 23–44 a), accompanied by local vertical recharge and upward leakage of older groundwater; and (iii) salt water zone, mostly comprising waters with ages beyond the dating range of both CFCs and tritium. Some shallow groundwater in the north of the Changyi groundwater depression belongs to the >50a water group (iii), indicating slow velocity of groundwater circulation and possible drawing in of saline or deep groundwater that is tracer‐free. Copyright © 2012 John Wiley & Sons, Ltd.     

Understanding and modelling spatial drain–aquifer interactions in a low‐lying coastal aquifer—the Thurne catchment,Norfolk, UK

Seawater intrusion into fresh groundwater formations generally results inadvertently from human activities, such as over‐abstraction from coastal aquifers. This article describes the data analysis to quantify drain–aquifer interactions in a low‐lying pump‐drained coastal aquifer, which is subject to saline intrusion due to widespread land drainage, and the resulting development and application of a numerical groundwater model to understand the spatial groundwater system behaviour (including groundwater salinity fluxes). Without measured flow data in this pump‐drained catchment, a novel groundwater head‐dependent approach to hydrograph separation is described. Time‐variant and time‐invariant MODFLOW analyses are utilised to examine the flow processes. A new approach to calculate drain coefficients, which represent the extensive network of drainage ditches in the regional model, using field information, is described; the sum of the drainage coefficients are close to the values independently estimated from the head‐dependent hydrograph separation. Results show that (1) the groundwater flows into the drainage systems are well reproduced using the new drain coefficients, (2) particle tracking of fresh and saline water can explain observed spatial salinity distribution within drainage networks and (3) the modelled flow of seawater across the coast is approximately 25% greater than that discharged by the pumps, demonstrating the need for drainage management to be aware of the slow response of groundwater systems to past drainage system changes. The article demonstrates that numerical groundwater modelling can produce the improved understanding needed to inform management decisions in such complex environments. Copyright © 2010 John Wiley & Sons, Ltd.     

Integration of remote sensing and electrical sounding data for hydrogeological exploration—a case study of Bakhar watershed,India / Intégration de données de télédétection et de sondages électriques pour l'exploration hydrogéologique—étude de cas du bassin versant de Bakhar,Inde

Abstract

A geographical information system (GIS) was used for the integration of hydrological data acquired using remote sensing and geoelectrical techniques to understand the groundwater condition of Bakhar watershed, Mirazpur District, UP, India. Indian remote sensing IRS-1D, LISS—III data were used to prepare a geomorphological and lineament map of the Bakhar watershed. Vertical electrical sounding (VES) was carried out in different geomorphic units, and ranges of electrical resistivity values were assigned to the different formations by calibrating electrical resistivity with borehole data. Based on these, a subsurface resistivity map and an aquifer thickness map were prepared. Several layers were superimposed using GIS techniques. Each theme was assigned a weight, depending on its influence on groundwater recharge. Each class or unit in the map was assigned a knowledge-based rank from one to four, depending on its significance in storage and transmittance of groundwater, and these were then multiplied by the layer weighting to produce a score. Based on these scores, the watershed was categorized into different groundwater potential zones. The results indicate that the eastern and northern parts of the study area have very good groundwater potential to meet the demands of water for irrigation and domestic purposes, whereas the southern region has poor groundwater potential zones. Such integrated analysis has not been attempted so far in this region for hydrogeological investigation.     

Correlation of vertical electrical sounding and electrical borehole log data for groundwater exploration

We have correlated the longitudinal unit conductance CL obtained from interpreted vertical electrical sounding data with the formation resistivity R_t and the formation resistivity factor F, obtained by carrying out electrical borehole logging. Interpreted geophysical data of eleven soundings and two electrical borehole log records are used for the analysis. The geophysical data used were acquired in a sedimentary basin. The study area is called Lower Maner Basin located in the province of Andhra Pradesh, India. Vertical electrical soundings were carried out using a Schlumberger configuration with half current electrode separation varying from 600–1000 m. For logging the two boreholes, a Widco logger‐model 3200 PLS was used. True formation resistivity R_t was calculated from a resistivity log. Formation resistivity factor F was also calculated at various depths using R_t values. An appreciable inverse relation exists between the correlated parameters. The borehole resistivity R_t and the formation resistivity factor F decrease with the increase in the longitudinal unit conductance CL. We have shown the use of such a relation in computing borehole resistivity R_t and formation resistivity factor F at sites that posses only vertical electrical sounding data, with a fair degree of accuracy. Validation of the correlation is satisfactory. Scope for updating the correlation is discussed. Significance and applications of the relation for exploration of groundwater, namely to update the vertical electrical sounding data interpretation by translating the vertical electrical sounding data into electrical borehole log parameters, to facilitate correlations studies and to estimate the porosity (φ), permeability (K) and water saturation S_w of water bearing zones are discussed.     

Dynamic river–groundwater exchange in the presence of a saline,semi‐confined aquifer

Understanding groundwater–surface water exchange in river banks is crucial for effective water management and a range of scientific disciplines. While there has been much research on bank storage, many studies assume idealized aquifer systems. This paper presents a field‐based study of the Tambo Catchment (southeast Australia) where the Tambo River interacts with both an unconfined aquifer containing relatively young and fresh groundwater (<500 μS/cm and <100 years old) and a semi‐confined artesian aquifer containing old and saline groundwater (electrical conductivity > 2500 μS/cm and >10 000 years old). Continuous groundwater elevation and electrical conductivity monitoring within the different aquifers and the river suggest that the degree of mixing between the two aquifers and the river varies significantly in response to changing hydrological conditions. Numerical modelling using MODFLOW and the solute transport package MT3DMS indicates that saline water in the river bank moves away from the river during flooding as hydraulic gradients reverse. This water then returns during flood recession as baseflow hydraulic gradients are re‐established. Modelling also indicates that the concentration of a simulated conservative groundwater solute can increase for up to ~34 days at distances of 20 and 40 m from the river in response to flood events approximately 10 m in height. For the same flood event, simulated solute concentrations within 10 m of the river increase for only ~15 days as the infiltrating low‐salinity river water drives groundwater dilution. Average groundwater fluxes to the river stretch estimated using Darcy's law were 7 m³/m/day compared with 26 and 3 m³/m/day for the same periods via mass balance using Radon (²²²Rn) and chloride (Cl), respectively. The study shows that by coupling numerical modelling with continuous groundwater–surface water monitoring, the transient nature of bank storage can be evaluated, leading to a better understanding of the hydrological system and better interpretation of hydrochemical data. Copyright © 2015 John Wiley & Sons, Ltd.     

Electrical earth resistivity surveying for delineating the characteristics of ground water in a semi‐arid region in the Khanasser Valley,northern Syria

A resistivity survey is conducted in Khanasser Valley, a semi‐arid region in northern Syria, to delineate the characteristics of ground water affected by the salt‐water intrusion related to Al‐Jaboul Sabkha. Existing wells were used to measure salinity and conductivity of water samplings. Vertical electrical sounding was carried out near the existing wells. The combination of resistivity and hydrogeological data enables the establishment of empirical relationships between earth resistivity, water resistivity, and the amount of total dissolved solids. These relationships are then used in order to derive salinity maps for electrode spacings of 70, 100, and 150 m. The distribution of fresh, brackish and salt‐water zones and their variations in space along two longitudinal profiles (LP1 and LP2) are established through converting subsurface depth–resistivity models into different ground‐water areas. The constructed ground‐water area maps allow interfaces between different water zones to be determined. Copyright © 2006 John Wiley & Sons, Ltd.     

SPONTANEOUS POTENTIALS AND RESISTIVITY SURVEYS APPLIED TO HYDROGEOLOGY IN A VOLCANIC AREA: CASE HISTORY OF THE CHAÎNE DES PUYS (PUY-DE-DÔME,FRANCE)1

Two geophysical methods [resistivity soundings and spontaneous potentials (SP)] are used to investigate aquifers in the Quaternary volcanic formations of the Chaîne des Puys (France). The interpretation of the resistivity soundings required a specific work of determination of the resistivity of the formations concerned. The basement topography, characterized by axial horst and graben structures and perpendicular palaeovalleys, was revealed. SP was developed experimentally in the field and theoretically. It is shown that the SP anomalies can be considered as double-layer potentials, with their source on the aquifers' interfaces. The most significant of them is the water table, which creates an SP anomaly of similar but inverted form. SP anomalies due to flows in unconfined aquifers, in palaeovalleys, and in volcanic ranges are computed and compared with the observed ones. These two methods have determined the extent and the boundary conditions of the different hydrogeological basins and have determined the major drainage axes and the groundwater flow pattern.     

Calibrating a Salt Water Intrusion Model with Time‐Domain Electromagnetic Data

Salt water intrusion models are commonly used to support groundwater resource management in coastal aquifers. Concentration data used for model calibration are often sparse and limited in spatial extent. With airborne and ground‐based electromagnetic surveys, electrical resistivity models can be obtained to provide high‐resolution three‐dimensional models of subsurface resistivity variations that can be related to geology and salt concentrations on a regional scale. Several previous studies have calibrated salt water intrusion models with geophysical data, but are typically limited to the use of the inverted electrical resistivity models without considering the measured geophysical data directly. This induces a number of errors related to inconsistent scales between the geophysical and hydrologic models and the applied regularization constraints in the geophysical inversion. To overcome these errors, we perform a coupled hydrogeophysical inversion (CHI) in which we use a salt water intrusion model to interpret the geophysical data and guide the geophysical inversion. We refer to this methodology as a Coupled Hydrogeophysical Inversion‐State (CHI‐S), in which simulated salt concentrations are transformed to an electrical resistivity model, after which a geophysical forward response is calculated and compared with the measured geophysical data. This approach was applied for a field site in Santa Cruz County, California, where a time‐domain electromagnetic (TDEM) dataset was collected. For this location, a simple two‐dimensional cross‐sectional salt water intrusion model was developed, for which we estimated five uniform aquifer properties, incorporating the porosity that was also part of the employed petrophysical relationship. In addition, one geophysical parameter was estimated. The six parameters could be resolved well by fitting more than 300 apparent resistivities that were comprised by the TDEM dataset. Except for three sounding locations, all the TDEM data could be fitted close to a root‐mean‐square error of 1. Possible explanations for the poor fit of these soundings are the assumption of spatial uniformity, fixed boundary conditions and the neglecting of 3D effects in the groundwater model and the TDEM forward responses.