Assessment of aquifer vulnerability to industrial waste water using resistivity measurements. A case study, along El-Gharbyia main drain, Nile Delta, Egypt

1D resistivity sounding and 2D resistivity imaging surveys were integrated with geological and hydrochemical data to assess the aquifer vulnerability and saltwater intrusion in the north of Nile Delta, Egypt. In the present study, the El-Gharbyia main drain was considered as a case study to map the sand bodies within the upper silt and clay aquitard. Twenty Schlumberger soundings and six 2D dipole-dipole profiles were executed along one profile close to the western side of the main drain. In addition, 14 groundwater samples and 4 surface water samples from the main drain were chemically analyzed to obtain the major and trace elements concentrations.The results from the resistivity and hydrochemical data were used to assess the protection of the groundwater aquifer and the potential risk of groundwater pollution. The inverted resistivities and thicknesses of the layers above the aquifer layer were used to estimate the integrated electrical conductivity (IEC) that can be used for quantification of aquifer vulnerability. According to the aquifer vulnerability assessment of an underlying sand aquifer, the southern part of the area is characterized by high vulnerability zone with slightly fresh to brackish groundwater and resistivity values of 11-23 Ω.m below the clay cap. The resistivity sections exhibit some sand bodies within the clay cap that lead to increase the recharging of surface waste water (650 mg/l salinity) and flushing the upper part of underlying saltwater aquifer. The region in the north has saltwater with resistivity less than 6 Ω.m and local vulnerable zones within the clay cap. The inverted 2D dipole-dipole profiles in the vulnerable zones, in combination with drilling information have allowed the identification of subsoil structure around the main drain that is highly affected by waste water.     

Geoelectromagnetic measurements across the southern Senegal basin (West Africa)

Magnetotelluric and differential geomagnetic sounding surveys, consisting of nine soundings, were performed in 1984 along a 200-km profile across the southern Senegal basin. They were intended to obtain information concerning the resistivity structure of the crust and upper mantle and the distribution of the induced electric currents. Magnetotelluric data indicate that two-dimensional resistivity models are appropriate for the region. The zone above the basement is highly inhomogeneous in geoelectrical structure. Very conductive sediments (0.6-3 ohm m) appear in the Mesozoic-Cenozoic Senegal basin. These sediments lie at depths of up to 4500 m on the west end of the profile. Below this, a modest resistivity material (10–30 ohm m) extends to a maximum depth of about 3000 m. The material at depth on the cast part of the traverse line is thought to be Palaeozoic sediments of the Bove basin. The depth of the magnetotelluric basement lies between about 250 m (in the east) and 4800 m (in the west). The crust is characterised by a drop in electrical resistivity at a depth of 15 km below the east part of the profile. Considering the total section, we observe a general trend towards lower resistivities at depths in excess of 100 km, the transition from 2000 ohm m to about 2 ohm m occurs in the depth range 100 to 175 km. An analysis of the geomagnetic variation field has identified a concentration of telluric current flow beneath the deep basin. It appears that the additional currents flowing in the striking direction of the Senegal basin are largely controlled by sedimentary rocks of high conductivity lying at depths less than 5 km. Model studies show that the local conductivity distribution is able to explain the currents circulating in the thick well-conducting sediments.     

Characterization of seasonal local recharge using electrical resistivity tomography and magnetic resonance sounding

A groundwater recharge process of heterogeneous hard rock aquifer in the Moole Hole experimental watershed, south India, is being studied to understand the groundwater flow behaviour. Significant seasonal variations in groundwater level are observed in boreholes located at the outlet area indicating that the recharge process is probably taking place below intermittent streams. In order to localize groundwater recharge zones and to optimize implementation of boreholes, a geophysical survey was carried out during and after the 2004 monsoon across the outlet zone. Magnetic resonance soundings (MRS) have been performed to characterize the aquifer and measure groundwater level depletion. The results of MRS are consistent with the observation in boreholes, but it suffers from degraded lateral resolution. A better resolution of the regolith/bedrock interface is achieved using electrical resistivity tomography (ERT). ERT results are confirmed by resistivity logging in the boreholes. ERT surveys have been carried out twice—before and during the monsoon—across the stream area. The major feature of recharge is revealed below the stream with a decrease by 80% of the calculated resistivity. The time‐lapse ERT also shows unexpected variations at a depth of 20 m below the slopes that could have been interpreted as a consequence of a deep seasonal water flow. However, in this area time‐lapse ERT does not match with borehole data. Numerical modelling shows that in the presence of a shallow water infiltration, an inversion artefact may take place thus limiting the reliability of time‐lapse ERT. A combination of ERT with MRS provides valuable information on structure and aquifer properties respectively, giving a clue for a conceptual model of the recharge process: infiltration takes place in the conductive fractured‐fissured part of the bedrock underlying the stream and clayey material present on both sides slows down its lateral dissipation. Copyright © 2007 John Wiley & Sons, Ltd.     

Geophysical surveys for identifying saline groundwater in the semi‐arid region of the central Altiplano,Bolivia

In the central part of the Bolivian Altiplano, the shallow groundwater presents electrical conductivities ranging from 0·1 to 20 mS/cm. In order to study the origin of this salinity pattern, a good knowledge is required of the geometry of the aquifer at depth. In this study, geophysics has been used to complement the sparse data available from drill holes. One hundred time‐domain electromagnetic (TDEM) soundings were carried out over an area of 1750 km². About 20 geological logs were available close to some of the TDEM soundings. Three intermediate results were obtained from the combined data: (i) the relationship between the electrical conductivity of the groundwater and the formation resistivity, (ii) geoelectrical cross‐sections and (iii) geoelectrical maps at various depths. The limited data set shows a relationship between resistivity and the nature of the rock. From the cross‐sections, a conductive substratum with a resistivity of less than 1 Ω·m was identified at most of the sites at depths ranging from 50 to 350 m. This substratum could be a clay‐rich formation containing brines. Using derived relationships, maps of the nature of the formation (sandy, intermediate and clayey sediments) were established at depths of 10 and 50 m. Discrimination between sand and clays was impossible where groundwater conductivity is high (>3 mS/cm). In the central part of the area, where the groundwater conductivity is low, sandy sediments are likely to be present from the surface to a depth of more than 200 m. Clayey sediments are more likely to be present in the south‐east and probably constitute a hydraulic barrier to groundwater flow. In conclusion, the study demonstrates the efficiency of the TDEM sounding method to map conductive zones. Copyright © 2001 John Wiley & Sons, Ltd.     

Electric-hydraulic conductivity correlation in fractured crystalline bedrock: Central Landfill, Rhode Island, USA

Remote sensing and geoelectrical methods were used to find water-bearing fractures in the Scituate granite under the Central Landfill of Rhode Island. These studies were necessary to evaluate the integrity of the sanitary landfill and for planning safe landfill extensions. The most useful results were obtained with fracture trace analysis using Landsat and SLAR imagery in combination with ground-based resistivity measurements using Schlumberger vertical electrical soundings based on the assumption of horizontally layered strata. Test borings and packer tests confirmed, in the presence of a lineament and low bedrock resistivity, the probable existence of high bedrock fracture density and high average hydraulic conductivity. However, not every lineament was found to be associated with high fracture density and high hydraulic conductivity. Lineaments alone are not a reliable basis for characterising a landfill site as being affected by fractured bedrock. Horizontal fractures were found in borings located away from lineaments. High values of hydraulic conductivity were correlated with low bedrock resistivities. Bedrock resistivities between 60 and 700 Ω m were associated with average hydraulic conductivities between 4 and 60 cm/day. In some cases very low resistivities were confined to the upper part of the bedrock where the hydraulic conductivity was very large. These types of fractures apparently become narrower in aperture with depth. Bedrock zones having resistivities greater than 1000 Ω m showed, without exception, no flow to the test wells. Plots of bedrock resistivity versus the average hydraulic conductivity indicate that the resistivity decreases with increasing hydraulic conductivity. This relationship is inverse to that found in most unconsolidated sediments and is useful for estimating the hydraulic conductivity in groundwater surveys in fractured bedrock. In appropriate settings such as the Central Landfill site in New England, this electric-hydraulic correlation relationship, supplemented by lineament trace analysis, can be used effectively to estimate the hydraulic conductivity in bedrock from only a limited number of resistivity depth soundings and test wells.     

Geophysical Assessment of Groundwater Potential: A Case Study from Mian Channu Area,Pakistan

An integrated study using geophysical method in combination with pumping tests and geochemical method was carried out to delineate groundwater potential zones in Mian Channu area of Pakistan. Vertical electrical soundings (VES) using Schlumberger configuration with maximum current electrode spacing (AB/2 = 200 m) were conducted at 50 stations and 10 pumping tests at borehole sites were performed in close proximity to 10 of the VES stations. The aim of this study is to establish a correlation between the hydraulic parameters obtained from geophysical method and pumping tests so that the aquifer potential can be estimated from the geoelectrical surface measurements where no pumping tests exist. The aquifer parameters, namely, transmissivity and hydraulic conductivity were estimated from Dar Zarrouyk parameters by interpreting the layer parameters such as true resistivities and thicknesses. Geoelectrical succession of five‐layer strata (i.e., topsoil, clay, clay sand, sand, and sand gravel) with sand as a dominant lithology was found in the study area. Physicochemical parameters interpreted by World Health Organization and Food and Agriculture Organization were well correlated with the aquifer parameters obtained by geoelectrical method and pumping tests. The aquifer potential zones identified by modeled resistivity, Dar Zarrouk parameters, pumped aquifer parameters, and physicochemical parameters reveal that sand and gravel sand with high values of transmissivity and hydraulic conductivity are highly promising water bearing layers in northwest of the study area. Strong correlation between estimated and pumped aquifer parameters suggest that, in case of sparse well data, geophysical technique is useful to estimate the hydraulic potential of the aquifer with varying lithology.     

Geothermal study of reunion island: Audiomagnetotelluric survey

In 1979 and 1980, 535 magnetotelluric soundings using a frequency range of 1700 Hz-8 Hz were performed on the island of Réunion for geothermal exploration.Few direct hot-water discharges were observed. Consequently, geophysical methods, particularly the audiomagnetotelluric method, were used extensively. Favourable geological conditions for this method were encountered and the results, which were controlled using classical electrical methods on test areas, suggest an unusual distribution of resistivities for lava flows situated above suspected geothermal areas. These layers progressively decrease in resistivity down to a very conductive layer. In areas where these conductive layers were nearest the surface, detailed studies were carried out showing a close correlation between decreasing resistivity and increasing hydrothermal alteration. In addition, gradient wells reveal high geothermal gradients in such areas. The conductive layers revealed by audiomagnetotelluric soundings seem to correlate with thermal effects creating progressive hydrothermal alteration from an inferred hot-water reservoir up toward the surface.     

Joint Audio-Magnetotelluric and Passive Seismic Imaging of the Cerdanya Basin

The structure of Cerdanya Basin (north-east of Iberian Peninsula) is partly known from geological cross sections, geological maps and vintage geophysical data. However, these data do not have the necessary resolution to characterize some parts of Cerdanya Basin such as the thickness of soft soil, geometry of bedrock or geometry of geological units and associated faults. For all these reasons, the main objective of this work is to improve this deficiency carrying out a detailed study in this Neogene basin applying jointly the combination of passive seismic methods (H/V spectral ratio and seismic array) and electromagnetic methods (audio-magnetotelluric and magnetotelluric method). The passive seismic techniques provide valuable information of geometry of basement along the profile. The maximum depth is located near Alp village with a bedrock depth of 500 m. The bedrock is located in surface at both sites of profile. The Neogene sediments present a shear-wave velocity between 400 and 1000 m/s, and the bedrock basement presents a shear-wave velocity values between 1700 and 2200 m/s. These results are used as a priori information to create a 2D resistivity initial model which constraints the inversion process of electromagnetic data. We have obtained a 2D resistivity model which is characterized by (1) a heterogeneous conductivity zone (<40 Ohm m) that corresponds to shallow part of the model up to 500 m depth in the centre of the profile. These values have been associated with Quaternary and Neogene sediments formed by silts, clays, conglomerates, sandstones and gravels, and (2) a deeper resistive zone (1000–3000 Ohm m) interpreted as Palaeozoic basement (sandstones, limestones and slates at NW and conglomerates and microconglomerates at SE). The resistive zone is truncated by a discontinuity at the south-east of the profile which is interpreted as the Alp-La Tet Fault. This discontinuity is represented by a more conductive zone (600 Ohm m approx.) and is explained as a combination of fractured rock and a fluid network. The result highlights that the support between different geophysical methods is essential in producing geophysical meaningful models.     

Electrical resistivity imaging and hydrochemical analysis for groundwater investigation in Kuala Langat,Malaysia

ABSTRACT

Integrated two-dimensional electrical resistivity imaging (ERI) and hydrochemical surveys were used to investigate the groundwater alluvial aquifer in Kuala Langat, Malaysia. The study in the Langat basin considered the thickness of the aquifer, the depth of the bedrock, the regions influenced by seawater intrusion, and the monitoring of water levels. The resistivity imaging results show that the upper layer consists of clay, while the second layer is an aquifer whose thickness varies mostly in the range of 10–30 m, and in some cases extends to 40 m. The bedrock depth varies from 30 to 65 m. The chemical analyses were carried out on groundwater samples from nine boreholes collected between 2008 and 2012. The analyses indicate that the total dissolved solids (TDS) exceed 1000 mg L^-1 near the coastal area and are often less than 500 mg L^-1 further inland. The ERI and hydrochemical analyses reveal that groundwater in the study area, especially towards the coast, is a mixture of brackish and fresh waters.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR M.D. Fidelibus     

An audiomagnetotelluric study of the Meugueur-Meugueur ring structure, Aïr, Niger: ring dyke or cone sheet?

Three radial audiomagnetotelluric (AMT) sounding profiles were carried out across the narrow, 65-km diameter troctolitic Meugueur-Meugueur ring structure, central Aïr, Niger, to study its electrical configuration; one profile extended across the bedrock into the large Ofoud complex situated slightly off geographical centre within the ring. Apparent resistivity data from 27 sites ranged from isotropic to strongly anisotropic. In nearly all soundings, one- and two-dimensional modelling indicated the presence of a major zone of low resistivity (60–600 Ωm), about 200 m thick, dipping steeply inwards at an angle of 65–80° and extending to a depth of at least 2–5 km. This layer, overlain and underlain by rocks of higher resistivities in excess of 5000 Ωm, is taken to be the outer contact. A highly resistive body, about 200 m in width, dipping inwards to a depth of at least 4 km is taken to be the Meugueur-Meugueur intrusion, which is thus interpreted to be a cone sheet.     

Mapping Coastal Aquifers by Joint Inversion of DC and TEM Soundings-Three Case Histories

Electrical and electromagnetic methods are well suited for coastal aquifer studies because of the large contrast in resistivity between fresh water-bearing and salt water-bearing formations. Interpretation models for these aquifers typically contain four layers: a highly resistive unsaturated zone; a surficial fresh water aquifer of intermediate resistivity; an underlying conductive, salt water saturated aquifer; and resistive substratum. Additional layers may be added to allow for variations in lithology within the fresh water and salt water layers. Two methods are evaluated: direct current resistivity and time domain electromagnetic soundings. Use of each method alone produces nonunique solutions for resistivities and/or thicknesses of the different layers. We show that joint inversion of vertical electric and time domain electromagnetic soundings produces a more tightly constrained interpretation model at three test sites than is produced by inversion methods applied to each data set independently.     

Determination of aquifer geometry through geophysical methods: A case study from Quetta Valley, Pakistan

Due to increase in population and agricultural activities, the aquifer of Quetta Valley is under tremendous stress and the water table is declining at an increasing rate. This situation necessitates evaluation of the aquifer system, for which information about geometry of the aquifer is a prerequisite. However, there are no drilling-to-bedrock data available; therefore, electrical resistivity, seismic reflection and gravity methods were employed to determine geometry of the aquifer. Interpretation of vertical electrical soundings provided information about the depth-tobedrock at some specific points, whereas seismic reflection delineated bedrock topography along two lines. The depths to bedrock inferred from electrical resistivity and seismic reflection data were used as constraints in the modeling of gravity data. 2.75D gravity models were constructed along lines with a regular spacing. Map of depth-to-bedrock was prepared by contouring the depth given by the gravity models. Combination of these geophysical methods depicted the geometry of the aquifer. This example shows that in a similar geological setting proper integration of geophysical exploration methods can determine the aquifer geometry with an acceptable reliability and at an appropriate cost.     

Investigation of groundwater resources using controlled-source radio magnetotellurics (CSRMT) in glacial deposits in Heby,Sweden

We have combined tensor radio magnetotelluric- (RMT, 15–250 kHz) and controlled source tensor magnetotelluric (CSTMT, 1–12 kHz) data for the mapping of aquifers in gravel formations lying in between crystalline bedrock and clay rich sediments in the Heby area some 40 km west of Uppsala in Sweden. The estimated transfer functions, the impedance tensor and the tipper vector generally satisfy 1D or 2D necessary conditions except for the lowest CSTMT frequencies where near field effects become more dominant.The data measured from 8 profiles were inverted with the Rebocc code of Siripunvaraporn and Egbert (2000) assuming plane wave conditions. This meant that only 12 frequencies in the range of 4–180 kHz could be used. The four lowest frequencies of CSTMT in the range of 1–2.8 kHz were excluded because of source effects. Data from all profiles were inverted with a starting model of 100 Ω-m and a relative error floor of 0.02 on apparent resistivity, corresponding to less than 1^° on phase. Tipper vectors are generally small except when source effects become dominant in the lowest frequencies of CSTMT and were therefore not used for inversion. Comparing with models derived from vertical electrical soundings, refraction and reflection seismic data as well as ground truth from exploration wells assessed the reliability of the deep part of the models. Furthermore we carried out a non-linear resolution analysis to better quantify the depth extent of the aquifers.The inverted models from the Heby area show well the thickness variations of glacial deposits overlying crystalline bedrock. Generally, the upper 20 m of the models have resistivities below 40 Ω-m, taken to represent clay rich formations. Below the clay layer resistivities increase to about 40–400 Ω-m, interpreted to represent sand/gravel formations with a maximum thickness of about 40 m and a width of several hundred metres. This is a potential aquifer that extends in approximately N–S direction for some kilometres.     

STUDYING THE INTERCONNECTION OF GROUND AND SURFACE WATERS IN REPRESENTATIVE BASINS UNDER THE PROGRAMME OF THE INTERNATIONAL HYDROLOGICAL DECADE

Abstract

Four geoelectrical soundings were measured with a combination of Schlumberger and azimuthal or equatorial dipole electrode arrays on a Carboniferous limestone basin of the Condroz area, Belgium. The measuring technique is briefly outlined as well as the interpretation procedure, which follows a closed-loop scheme with control of calculated model curves. Some general problems of interpretation of geoelectrical sounding curves are tackled, as far as they have a practical bearing on the treatment of Condroz soundings.

The problem of determining the very high resistivity of limestone is approached through ARCHIE's formula, an empirical relation between the bulk rock resistivity, the porosity and the electrolyte resistivity. An evaluation of the latter two parameters, combined with electrical horizontal conductance measurements directly made on resistivity sounding curves, offers a possibility for fast determination of the total water storage in a limestone aquifer. Such storage determinations could be applied whenever an aquifer shows up as a conductive layer interbedded between two highly resistant layers (e.g. nonsaturated limestone and compact, non-fractured limestone).     

A field test of imaging properties of rotational invariants of the magnetotelluric impedance tensor

A part of the Békés Basin (an extensional sub‐basin of the Pannonian Basin, where the basement under thick Pannonian sediments is well known from deep boreholes and from seismic measurements, and where many magnetotelluric (MT) soundings have been carried out for frequencies ranging from 1 to 10^?3 Hz) was selected as a test area to assess the imaging performances of various apparent‐resistivity definitions computed with rotational invariants of either the real part of the complex impedance tensor, or its imaginary part, or both. A comparison (based on earlier 3D numerical studies) has been made between the magnetotelluric images obtained in this way and the depths to the high‐resistivity basement, as known from boreholes and seismic investigations. The correlation coefficient between the series of basement depth values at 39 MT sites and the apparent‐resistivity values was found to be stronger and high correlation appeared at a shorter period when it was computed with apparent resistivities based on the real tensor rather than with apparent resistivities based on the imaginary tensor. In the light of our studies, ρ_{Re Z} and the impedance phase seem to be more informative than any other combination of magnetotelluric interpretation parameters.     

Geoelectrical structure of the central zone of Piton de la Fournaise volcano (Réunion)

 A study of the geoelectrical structure of the central part of Piton de la Fournaise volcano (Réunion, Indian Ocean) was made using direct current electrical (DC) and transient electromagnetic soundings (TEM). Piton de la Fournaise is a highly active oceanic basaltic shield and has been active for more than half a million years. Joint interpretation of the DC and TEM data allows us to obtain reliable 1D models of the resistivity distribution. The depth of investigation is of the order of 1.5 km but varies with the resistivity pattern encountered at each sounding. Two-dimensional resistivity cross sections were constructed by interpolation between the soundings of the 1D interpreted models. Conductors with resistivities less than 100 ohm-m are present at depth beneath all of the soundings and are located high in the volcanic edifice at elevations between 2000 and 1200 m. The deepest conductor has a resistivity less than 20 ohm-m for soundings located inside the Enclos and less than 60–100 ohm-m for soundings outside the Enclos. From the resistivity distributions, two zones are distinguished: (a) the central zone of the Enclos; and (b) the outer zone beyond the Enclos. Beneath the highly active summit area, the conductor rises to within a few hundred meters of the surface. This bulge coincides with a 2000-mV self-potential anomaly. Low-resistivity zones are inferred to show the presence of a hydrothermal system where alteration by steam and hot water has lowered the resistivity of the rocks. Farther from the summit, but inside the Enclos, the depth to the conductive layers increases to approximately 1 km and is inferred to be a deepening of the hydrothermally altered zone. Outside of the Enclos, the nature of the deep, conductive layers is not established. The observed resistivities suggest the presence of hydrated minerals, which could be found in landslide breccias, in hydrothermally altered zones, or in thick pyroclastic layers. Such formations often create perched water tables. The known occurrence of large eastward-moving landslides in the evolution of Piton de la Fournaise strongly suggests that large volumes of breccias should exist in the interior of the volcano; however, extensive breccia deposits are not observed at the bottom of the deep valleys that incise the volcano to elevations lower than those determined for the top of the conductors. The presence of the center of Piton de la Fournaise beneath the Plaine des Sables area during earlier volcanic stages (ca. 0.5 to 0.150 Ma) may have resulted in broad hydrothermal alteration of this zone. However, this interpretation cannot account for the low resistivities in peripheral zones. It is not presently possible to discriminate between these general interpretations. In addition, the nature of the deep conductors may be different in each zone. Whatever the geologic nature of these conductive layers, their presence indicates a major change of lithology at depth, unexpected for a shield volcano such as Piton de la Fournaise. Received: 3 November 1999 / Accepted: 15 September 1999     

A Tracer Test for Detecting Cross Contamination Along a Monitoring Well Column

A tracer test was used to evaluate whether cross contamination exists along a monitoring well completed through a shallow ground water system in fractured clay and screened in a sand and gravel aquifer. The fractured clay is separated from the sand and gravel deposit by a layer of highly plastic unfractured clay. A natural vertical downward hydraulic gradient of approximately 0.5 exists between the shallow system and the sand and gravel aquifer. Ground water contamination was detected in an adjacent monitoring well screened in the fractured clay and in the monitoring well screened in the sand and gravel deposit. No ground water contamination was apparent in an intermediate well screened in the unfractured clay layer. A tracer of sodium bromide was injected into a shallow boring near the monitoring wells. The tracer was detected in the monitoring well in the sand and gravel aquifer after three to seven days. The bromide concentration continued to increase in this well with time while the concentration in the shallow boring declined. This trend of tracer concentration indicates the tracer has in fact migrated downward and possibly traveled along the well column.     

Geomorphology of the Alluvial Sediments and Bedrock in an Intermontane Basin: Application of Variogram Modeling to Electrical Resistivity Soundings

The study describes a methodology used to integrate legacy resistivity data with limited geological data in order to build three-dimensional models of the near subsurface. Variogram analysis and inversion techniques more typically found in the petroleum industry are applied to a set of 1D resistivity data taken from electrical surveys conducted in the 1980s. Through careful integration with limited geological data collected from boreholes and outcrops, the resultant model can be visualized in three dimensions to depict alluvium layers as lithological and structural units within the bedrock. By tuning the variogram parameters to account for directionality, it is possible to visualize the individual lithofacies and geomorphological features in the subsurface. In this study, an electrical resistivity data set collected as part of a groundwater study in an area of the Peshawar basin in Pakistan has been re-examined. Additional lithological logs from boreholes throughout the area have been combined with local outcrop information to calibrate the data. Tectonic activity during the Himalayan orogeny has caused uplift in the area and generated significant faulting in the bedrock resulting in the formation of depressions which are identified by low resistivity values representing clays. Paleo-streams have reworked these clays which have been eroded and replaced by gravel–sand facies along paleo-channels. It is concluded that the sediments have been deposited as prograding fan-shaped bodies and lacustrine deposits with interlayered gravel–sand and clay–silt facies. The Naranji area aquifer system has thus been formed as a result of local tectonic activity with fluvial erosion and deposition and is characterized by coarse sediments with high electrical resistivities.     

Magnetotelluric survey across the active spreading zone in Southwest Iceland

The results of 10 magnetotelluric soundings, performed along a 110-km-long profile crossing the constructive plate boundary in southwest Iceland, are presented. Apparent resistivities are interpreted by a horizontally stratified earth model to yield a pseudo cross-section along the profile. The crust-mantle interface contains a well conductive layer. The depth to the good conductor increases with age of the crust and the distance from the axial zone. This layer is interpreted as partially molten basalt, at a temperature about 1100°C and a volume fraction of the melt phase in the range 10–20%.The high-conductivity layer probably disappears west of the Borgarnes anticlinal axis, which separates the older (to the west) and younger (to the east) flood basalts in western Iceland, indicating that the temperature below the oldest part of the profile lies below the solidus curve of basalt. Recent seismic crustal investigations in the same area indicate a state of partial melting or a magma chamber, which agrees with the results of the magnetotelluric soundings.