Combining Electrical Resistivity Tomography and Ground Penetrating Radar to study geological structuring of karst Unsaturated Zone

This paper highlights the efficiency and complementarity of a light package of geophysical techniques to study the structure of karst Unsaturated Zone (UZ) in typical Mediterranean environment where soil cover is thin or absent. Both selected techniques, 2D Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT), are widely used in environmental studies and their application is accessible for a lot of scientists/engineers. However, GPR or ERT alone is not able to provide an enhanced characterization of geological features in karst media. In the present study, GPR results supply a near surface high resolution imaging and thus can provide relevant geological information such as stratifications and fractures. Despite the quality of the results GPR's investigation depth remains limited to around 12 m. Apparent and inverted resistivity provided by ERT surveys shows strong lateral and vertical variations. These variations can inform about general geological structuring and feature orientation. ERT is able to prospect down to 40 m but it's a low resolution integrative technique. In the study area the investigated limestone is a commonly electrical resistive formation (more than 2000 Ω.m). However deeper than 5–7 m, the ERT profiles reveal several zones of moderate resistivity (around 900 Ω.m). In these zones a stratification change corresponding to slanted bedding is clearly identified by GPR results. The combination of both GPR and ERT results can allow a well-established geological interpretation. These moderate resistivity zones with slanted beddings can explain the presence of a perennial water flow point 35 m below the surface of the studied site within the underground gallery of the Low-Noise Underground Laboratory (LSBB).     

Structural controls on the 1998 volcanic unrest at Iwate volcano: Relationship between a shallow, electrically resistive body and the possible ascent route of magmatic fluid

Magnetotelluric (MT) measurements were conducted at Iwate volcano, across the entirety of the mountain, in 1997, 1999, 2003, 2006, and 2007. The survey line was 18 km in length and oriented E–W, comprising 38 measurements sites. Following 2D inversion, we obtained the resistivity structure to a depth of 4 km. The surface resistive layer (~ several hundreds of meters thick) is underlain by extensive highly conductive zones. Based on drilling data, the bottom of the highly conductive zone is interpreted to represent the 200 °C isotherm, below which (i.e., at higher temperatures) conductive clay minerals (smectite) are rare. The high conductivity is therefore mainly attributed to the presence of hydrothermally altered clay. The focus of this study is a resistive body beneath the Onigajo (West-Iwate) caldera at depths of 0.5–3 km. This body appears to have impeded magmatic fluid ascent during the 1998 volcanic unrest, as inferred from geodetic data. Both tectonic and low-frequency earthquakes are sparsely distributed throughout this resistive body. We interpret this resistive body as a zone of old, solidified intrusive magma with temperatures in excess of 200 °C. Given that a similar relationship between a resistive body and subsurface volcanic activity has been suggested for Asama volcano, structural controls on subsurface magmatic fluid movement may be a common phenomenon at shallow levels beneath volcanoes.     

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.     

Shallow karst exploration using MT-VLF and DC resistivity methods1

A geophysical test was carried out over a well-located and fairly embedded karstic conduit. The MT-VLF method was selected because of its high resolution and its ability to provide a resistivity parameter sensitive to water and Clay. This method was used together with DC resistivity methods which allow vertical adjustment of the VLF data and show consistency between the investigation and target depths. After correcting the deformations due to the polarization of the primary field, the MT-VLF data show clearly, in the Central part of the site, that the conduit does not coincide with an anomaly axis but coincides with the boundary between a conductive area and a resistive area. 2D modelling confirms that direct detection of the conduit is not feasible and that the conduit is located close to a conductive zone corresponding to a completely Clay-filled fractured zone. This situation was observed on the whole site and the conduit seemed systematically to thread its way between the conductive zones to join the outlet. The distribution of the conductive fractured zones and the direction of the hydraulic gradient were two important elements in predicting the location of the conduits. A 3D approach would increase the probability of finding the conduits in such a case.     

Potential of Electrical Resistivity Tomography to Detect Fault Zones in Limestone and Argillaceous Formations in the Experimental Platform of Tournemire, France

The Experimental platform of Tournemire (Aveyron, France) developed by IRSN (French Institute for Radiological Protection and Nuclear Safety) is located in a tunnel excavated in a clay–rock formation interbedded between two limestone formations. A well-identified regional fault crosscuts this subhorizontal sedimentary succession, and a subvertical secondary fault zone is intercepted in the clay–rock by drifts and boreholes in the tunnel at a depth of about 250 m. A 2D electrical resistivity survey was carried out along a 2.5 km baseline, and a takeout of 40 m was used to assess the potential of this method to detect faults from the ground surface. In the 300 m-thick zone investigated by the survey, electrical resistivity images reveal several subvertical low-resistivity discontinuities. One of these discontinuities corresponds to the position of the Cernon fault, a major regional fault. One of the subvertical conductive discontinuities crossing the upper limestone formation is consistent with the prolongation towards the ground surface of the secondary fault zone identified in the clay–rock formation from the tunnel. Moreover, this secondary fault zone corresponds to the upward prolongation of a subvertical fault identified in the lower limestone using a 3D high-resolution seismic reflection survey. This type of large-scale electrical resistivity survey is therefore a useful tool for identifying faults in superficial layers from the ground surface and is complementary to 3D seismic reflection surveys.     

Resistivity mapping using the VLF-MT method around surface fault ruptures of the 1995 Hyogo-ken Nanbu earthquake, Japan

Abstract Distinctive fault ruptures, the Nojima Fault and Ogura Fault, appeared along the northwestern coast of Awaji Island at the time of the 1995 Hyogo-ken Nanbu earthquake (Kobe earthquake). In order to delineate the shallow resistivity structures around the faults just after they formed, Very Low Frequency Magnetotelluric (VLF-MT) surveys were made at five sites along the Nojima Fault and at one site along the Ogura Fault. Fourteen transects were made at the one site on the Ogura Fault, and another transect covers the area between the two faults. Changes in apparent resistivity or phase, or both, commonly occur when crossing the surface location of one of the faults, except for the northern transects at OGR-0 on the Ogura Fault. Apparent resistivity values of less than 100 Ωm were observed for Tertiary and Quaternary sediments and values larger than 200 Ωm for granitic rocks. The resistivity structures are related to the morphological characteristics of the fault ruptures. Remarkably conductive zones (less than 10 Ωm in apparent resistivity and 30–40 m in width) were found where the surface displacement is distinct and prominent along a single fault plane. If remarkably conductive zones were formed at the time of the 1995 Hyogo-ken Nanbu earthquake, the results provide a good constraint on the dimensions of a conductive zone near the surface that was made by one earthquake. Alternatively, if characteristic resistivity structures existed prior to the earthquake, the conductive zone was probably formed by some tens of earthquakes in relatively modern times. In this case, this phenomenon is inferred to be a concentration of fracturing in a narrow zone and is associated with the formation of clay minerals, which enhance rock conductivity.     

A review on developments in the electrical structure of craton lithosphere

Cratons have a long history of evolution. In this paper, applications of the magnetotelluric method used in the study of craton lithosphere over the past 30 years were reviewed, examining case studies of cratons in North America, South America, Asia, Australia, and Africa. The nuclei of the Archean cratons, for example the Kalahari Craton and Rae Craton, are usually characterized by thick and highly resistive lithospheric roots. During or after the formation of the cratons, tectonothermal events, such as collision, mantle plume, and asthenosphere upwelling led to the formation of high-conductivity zones in the craton lithosphere, which could be attributed to the increased hydrogen content (of nominally anhydrous minerals), higher iron content, and formation of graphite films or sulfides along the grain boundary of minerals. These conductive zones are characterized by resistivity discontinuities in craton lithosphere. In particular, the conductive zones include (1) large-scale lithospheric mantle conductors beneath the Slave Craton, Gawler Craton, and central part of North China Craton(Trans-North China Orogen); (2) near-vertical high-conductivity zone associated with the fossil subduction zone beneath the Dharwar Craton and Slave Craton; and (3) regional lateral electrical discontinuities, such as a conductive anomaly under the Bushveld Complex of the Kaapvaal Craton. The eMoho refers to the electrical discontinuity in the crust-mantle boundary. In existing research, this has been detected under the condition of extremely high lithospheric resistivity with only a slight decrease in the lower crust, and in the case of a very thin conductive lower crust or the lack thereof. In the resistivity model, the unique “mushroom-like” lower crust-lithosphere mantle conductor and very thin lower crust layer of the North China Craton may represent lithosphere destruction and/or thinning. We also find that some of the cratons are still not well understood. Therefore, extensive three-dimensional inversion and joint interpretation of geochemical, geophysical, and geologic data are necessary to understand the tectonic evolutionary history of craton lithosphere.

     

Detection of alteration at the Millennium uranium deposit in the Athabasca Basin: a comparison of data from two airborne electromagnetic systems with ground resistivity data

The Millennium uranium deposit is located within the Athabasca Basin of northern Saskatchewan. The basement rocks, comprised primarily of paleo‐Proterozoic gneisses, are electrically resistive. However, the deposit is associated with highly conductive graphitic metasediments that are intercalated with the gneisses. An unconformity separates the basement rocks from the overlying, horizontally stratified, Proterozoic sandstones of the Athabasca Group (which are also highly resistive). The strike extents of the graphitic metasedimenary packages are extensive and therefore electromagnetic (EM) survey techniques are successful at identifying these zones but do not identify the specific locations where they are enriched in uranium. Through drilling it has been noted that hydrothermal processes associated with mineralization has altered the rocks in the vicinity of the deposits, which should in theory result in a resistivity low. A significant resistivity low has been mapped coincident with the Millennium deposit using ground resistivity survey techniques. However, a comparison of the airborne EM and ground resistivity results reveals that the two data sets have imaged different features. The resistive‐limit (on‐time) windows of the MEGATEM data show conductive features corresponding to lakes located to the west and south of the deposit. The late‐time windows show a feature to the east of the deposit, interpreted as being associated with the east‐dipping graphitic basement conductors (similar to that observed in historical ground EM data collected in this area). The early‐time TEMPEST windows (delay times less than 0.2 ms) show a broad resistivity low located at approximately the same location as where the alteration has been identified through drilling. Modelling the data is not easy but a response that decays prior to 0.3 ms is consistent with 500 Ωm material in the sandstone, a resistivity value close to the lower limit with respect to the hydrothermally altered Athabasca group sediments in this area. The MEGATEM system does not see a conductive zone over the alteration as clearly but the high signal‐to‐noise ratio in the late‐time MEGATEM data means that the conductive material at a greater depth is more coherently imaged.     

跨安宁河断裂带浅孔综合地球物理测井成果分析

利用综合地球物理测井技术测得跨川西安宁河断裂带冕宁—西昌段两个钻孔的电阻率、声波速度、自然电位、自然伽马等4种综合测井曲线,并结合钻孔取芯资料得到了钻孔周围所钻遇岩土层的综合物性参数．本文得到的综合测井解释成果客观真实地反映了安宁河断裂带东、西两盘的介质物性和岩性特征,揭示了断层东、西两盘地层的纵波速度变化、视电阻率值等物性变化特征,为研究安宁河断裂带的物质组成、介质物性、断层结构及其活动状态提供了基础资料．综合地球物理测井技术有望成为监测活动断裂带附近介质物性参数动态变化的一种新手段．      

Characterization of recharge through complex vadose zone of a granitic aquifer by time-lapse electrical resistivity tomography

The vadose zone is the main region controlling water movement from the land surface to the aquifer and has a very complex structure. The use of non-invasive or minimally invasive geophysical methods especially electrical resistivity imaging is a cost-effective approach adapted for long-term monitoring of the vadose zone. The main aim of this work is to know the fractures in the vadose zone, of granitic terrene, through which the recharge or preferred path recharge to the aquifer takes place and thus to relate moisture and electrical resistivity. Time lapse electrical resistivity tomography (TLERT) experiment was carried out in the vadose zone of granitic terrene at the Indian Geophysical Research Institute, Hyderabad along two profiles to a depth of 18 m and 13 m each. The profiles are 300 m apart. Piezometric, rainfall and soil moisture data were recorded to correlate with changes in the rainfall recharge. These TLERT difference images showed that the conductivity distribution was consistent with the recharge occurring along the minor fractures. We mapped the fractures in hard rock or granites to see the effect of the recharge on resistivity variation and estimation of moisture content. These fractures act as the preferred pathways for the recharge to take place. A good correlation between the soil moisture and resistivity is established in the vadose zone of granitic aquifer. Since the vadose zone exhibits extremely high variability, both in space and time, the surface geophysical investigations such as TLERT have been a simple and useful method to characterize the vadose zone, which would not have been possible with the point measurements alone. The analyses of the pseudosection with time indicate clearly that the assumption of the piston flow of the moisture front is not valid in hard rocks. The outcome of this study may provide some indirect parameters to the well known Richard's equation in studying the unsaturated zone.     

Geophysical Discovery of a New LNAPL Plume at the Former Wurtsmith AFB, Oscoda, Michigan

Alight nonaqueous phase liquid (LNAPL) ground water contaminant plume has been discovered by purely geophysical means at the former Wurtsmith Air Force Base (AFB) near Oscoda, Michigan. It is located near another plume called FT-02, which is a well-studied area undergoing natural bioremediation. The plume was discovered by ground penetrating radar (GPR) profiling while extending a long line from FT-02 to establish background variability around that plume. The new plume was apparent because of a high-conductivity "shadow' or GPR reflection attenuation observed below the conductive zone at the top of the aquifer, identical to the pattern observed at the FT-02 plume. Further GPR surveys were conducted by students of a Western Michigan University geophysics field course to outline the proximal part of the plume. The GPR survey was supplemented by an electromagnetic induction (EM) survey which showed a group of four cables crossing the area. Finally, a magnetometer survey was conducted to search for any buried steel objects which might have been missed by the EM survey. The results of the three geophysical surveys were then used by students of a University of Michigan field course to guide subsurface soil and fluid sampling, which verified the presence of residual LNAPL product and ground water with conductivities 2.5 to 3.3 times above background. The plume source is in the vicinity of a vaulted underground storage tank (UST) formerly used for the collection of waste solvents and fuels for subsequent use in the fire training exercises at FT-02. This newly discovered LNAPL plume, along with other "mature' plumes, fits the electrical model which predicts conductive ground water below the decomposing but electrically resistive LNAPLs. Finally, this is a fine example of the cooperative use of a dedicated research site for training by students of two different universities.     

Study of the Santa Catarina aquifer system (Mexico Basin) using magnetotelluric soundings

A tensor magnetotelluric test survey was carried out in the region of Santa Catarina, located in the Chalco sub-basin of the Mexico Basin. The objective was to define the stratification at depth with an emphasis on the geometry of the main aquifer of that region which is partially known from DC resistivity soundings and drilling. High-quality magnetotelluric soundings could be recorded in the immediate vicinity of large urban zones because the sub-surface is very conductive. Interpretation shows that the solid bedrock is located at a depth of at least 800 m to the south and 1300 m to the north; it could, however, be much deeper. Using complementary DC resistivity sounding and well-logging data, three main layers have been defined overlying the bedrock. These layers are, from surface to bottom, an unsaturated zone of sand, volcanic ash and clay about 10 m thick, followed by a very conductive (1.5 ohm·m) 200 m thick layer of sand and ash with intercalated clay, saturated with highly mineralized water, and finally a zone with resistivity increasing gradually to 60 ohm·m. The investigated deep aquifer constitutes most of this third layer. It consists of a sequence of sand, gravel, pyroclastites and mainly fractured basalts. MT resistivity soundings and magnetic transfer functions also indicate that a shallow resistive structure is dipping, from the northwest, into the lacustrine deposits of the basin. This geologic feature is likely to be highly permeable fractured basaltic flows, which provide a channel by which water contaminated by the Santa Catarina landfill may leak into the basin.     

Exploration of geothermal structure in Puga geothermal field, Ladakh Himalayas, India by magnetotelluric studies

To understand the crustal electric structure of the Puga geothermal field located in the Ladakh Himalayas, wide band (1000 Hz–0.001 Hz) magnetotelluric (MT) study have been carried out in the Puga area. Thirty-five MT sites were occupied with site spacing varying from 0.4 to 1 km. The measurements were carried out along three profiles oriented in east–west direction. After the preliminary analysis, the MT data were subjected to decomposition techniques. The one-dimensional inversion of the effective impedance data and the two-dimensional inversion of the TE (transverse electric) and TM (transverse magnetic) data confirm the presence of low resistive (5–25 Ω m) near surface region of 200–300 m thick in the anomalous geothermal part of the area related to the shallow geothermal reservoir. Additionally, the present study delineated an anomalous conductive zone (resistivity less than 10 Ω m) at a depth of about 2 km which is possibly related to the geothermal source in the area. A highly resistive basement layer separates the surface low resistive region and anomalous conductive part. The estimated minimum temperature at the top of conductive part is about 250 °C. The significance of the deeper conductive zone and its relation to the geothermal anomaly in the area is discussed.     

Fiber Optic Measurements of Soil Moisture in a Waste Rock Pile

The design and construction of a waste rock pile influences water infiltration and may promote the production of contaminated mine drainage. The objective of this project is to evaluate the use of an active fiber optic distributed temperature sensing (aFO-DTS) protocol to measure infiltration and soil moisture within a flow control layer capping an experimental waste rock pile. Five hundred meters of fiber optic cable were installed in a waste rock pile that is 70 m long, 10 m wide, and was covered with 0.60 m of fine compacted sand and 0.25 m of non-reactive crushed waste rock. Volumetric water content was assessed by heating the fiber optic cable with 15-min heat pulses at 15 W/m every 30 min. To test the aFO-DTS system 14 mm of recharge was applied to the top surface of the waste rock pile over 4 h, simulating a major rain event. The average volumetric water content in the FCL increased from 0.10 to 0.24 over the duration of the test. The volumetric water content measured with aFO-DTS in the FCL and waste rock was within ±0.06 and ±0.03, respectively, compared with values measured using 96 dielectric soil moisture probes over the same time period. Additional results illustrate how water can be confined within the FCL and monitored through an aFO-DTS protocol serving as a practical means to measure soil moisture at an industrial capacity.     

Architecture of the deep critical zone in the Río Icacos watershed (Luquillo Critical Zone Observatory,Puerto Rico) inferred from drilling and ground penetrating radar (GPR)

In the critical zone, surficial bedrock interactions result in the formation of a mantle of chemically‐ and physically‐altered material defined here as regolith. In the watershed of the Río Icacos, an upland river draining the Luquillo Mountains in tropical Puerto Rico, we explored the influence of lithology (quartz diorite versus hornfels‐facies volcaniclastic rock) on weathering. Regolith profiles were studied by drilling boreholes and imaging the subsurface using ground penetrating radar (GPR). Overall, the regolith structure is not laterally continuous but rather is punctuated by zones of deep fractures that host in situ weathering, corestones, and colluvial material. GPR images of these vertical zones show reflectors at 15–20 m depth. Thus, the architecture of the critical zone in the upper Luquillo Mountains is highly dependent on lithology and its influence on fracture development. At the highest elevations where hornfels overlies quartz diorite, positive feedbacks occur when the water table drops so that oxidative weathering of biotite in the more felsic rock creates microfractures and allows deeper infiltration of meteoric waters. Such exposure results in some of the fastest weathering rocks in the world and may contribute to formation of the knickpoint in the Río Icacos watershed. This work represents the first study combining GPR and drilling to look at the structure of the deep critical zone and demonstrates: (1) the importance of combining direct methods (such as drilling) with indirect methods (such as GPR) to understand the architecture of the critical zone in tropical systems; (2) the interplay of the surficial stress regime, lithology and climate in dictating the architecture of weathering. Copyright © 2016 John Wiley & Sons, Ltd.     

The main structure and tectonic features of the Chernorud-Mukhor site on the western shore of Lake Baikal from TEM and SP measurements

Transient electromagnetic (TEM), self-potential (SP) and geoelectrical mapping measurements were carried out at the Chernorud-Mukhor site in the Priolkhonje area on the western shore of the Lake Baikal. All measurements were made along several profiles across the main strike of the regional Primorsky fault. TEM measurements were carried out in a time range from a few tens of microseconds to several tens of milliseconds. The most important result of the 1D modelling of TEM soundings is the discovery of nearly horizontal boundaries that divide high resistive overlying and well conducting underlying rocks. The resistivity of the former is in the range from 100 Ωm to 1000 Ωm, while the resistivity of the latter varies from less than 1 Ωm to several tens of Ωm. This good conductive zone could also be verified by geoelectrical mapping using Schlumberger array (AB/2=100 m). Due to high conductivity of the underlying rocks only the upper boundary of the conductive layer could be determined by TEM soundings. A regional SP anomaly with amplitude of about −450 mV has also been observed above the low resistivity zone indicating the electron nature of its conductance. Geologically, the conductive zone is represented by a graphite-bearing layer within the region of archean rocks. Since that layer extends over a large area, it may be used as a key in studying structures and tectonics of the Priolkhonje area. A 1D TEM geoelectric section shows a wide, gently sloping syncline as a probable base structure of the Chernorud-Mukhor site. Neotectonic faults divide the syncline into vertically displaced blocks that form a wide complicated graben with a total amplitude of about 250 m.     

Assessing the contribution of electrical resistivity tomography (ERT) and self-potential (SP) methods for a water well drilling program in fractured/karstified limestones

ERT and SP investigations were conducted in carbonate rocks of the Dinant Synclinorium (Walloon Region of Belgium) to find suitable locations for new water wells in zones with little hydrogeological data. Since boreholes information needed to be representative of the area, large fractured zones were searched for the drillings. Large ERT profiles (320 to 640 m) allowed us to image the resistivity distribution of the first 60 m of the subsurface and to detect and characterize (in terms of direction, width and depth) fractured zones expected to be less resistive. Data errors, depth of investigation (DOI) indexes and sensitivity models were analyzed in order to avoid a misinterpretation of the resulting images. Self-potential measurements were performed along electrical profiles to complement our electrical results. Some negative anomalies possibly related to preferential flow pathways were detected. A drilling campaign was conducted according to geophysical results. ‘Ground truth’ geological data as well as pumping tests information gave us a way to assess the contribution of geophysics to a drilling program. We noticed that all the wells placed in low resistivity zones associated with SP anomalies provide very high yields and inversely, wells drilled in resistive zones or outside SP anomalies are limited in terms of capacity. An apparent coupling coefficient between SP signals and differences in hydraulic heads was also estimated in order to image the water table.     

Conductivity structure of crust in the Tangshan seismic area and the possibility of exploring potential seismic sources by magnetotelluric method

It is shown by the result of digital magnetotelluric soundings in the Tangshan seismic area and its surrounding regions that the crust under the surface conductive sediments is divided into two layers,i. e., the resistive upper crust and the conductive lower crust. The upper crust wherein the Tangshan main shock and most of the aftershocks occurred is a convex lens-like body which is cut by faults at the east, south and west sides. The focus of the mainshock was located at the position of maximum thickness of the resistive upper crust while the spatial variation of Curie point isothermal surface and the deepest limit of the depths of aftershocks coincide with the downward depression of the bottom of the resistive upper crust. Thus, the Tangshan main shock and most of its aftershocks were related closely to the resistive upper crust from the view points of either vertical layering or lateral variations. And there were only a very few aftershocks in the conductive lower crust. The mechanical property of the rocks transforms from being brittle in the upper crust into ductile in the lower crust mainly due to the combination of different factors, e.g., increase of confining pressure, change in minerals, rise in temperature as well as stabilization of slips by high pore—pressure. A small amount of water and a rise in temperature may lead to a decrease of the electric resistivity within the rock while a change in the static pressure and mineral content within the rock causes very little change in the electric resistivity. Thus it is deduced that a resistive upper crust and a more conductive lower crust from the view points of either vertical or lateral variations are related to the brittle and ductile properties respectively. Hence it is possible that there is a relationship between the electric structure of the crust and zones of potential seismic hazard. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 354–363, 1991. The project is supported by the Chinese Joint Seismological Science Foundation, and the Swedish International Development Authority.     

Imaging the Mariánské Lázně Fault (Czech Republic) by 3-D ground-penetrating radar and electric resistivity tomography

Geodynamic activity in the area of West Bohemia is typified by the occurrence of earthquake swarms, Quaternary volcanism and high flux of mantle-derived CO₂. The highest swarm activity occurs beneath the eastern edge of the Cheb basin, which is delineated by the NW-SE trending morphologically pronounced Mariánské Lázn?? Fault (MLF) controlling the formation of the basin. The previous trenching survey across the MLF zone has identified several fault strands with possible Quaternary activity. In this paper we present the results of the geophysical survey focused to trace the faults signatures in geophysical sections and to build an image of near surface tectonics. The method of electric resistivity tomography (ERT) along two profiles parallel to the trench identified a strong resistivity contrast between the bodies of sandy gravels in the middle and conductive clayey sands to the west and weathered crystalline basement to the east. The 2-D ground penetration radar (GPR) sections show direct correlation of reflections with lithological boundaries identified in the trench. As expected, the GPR signal amplitudes increase with the resistivities found in the ERT sections. Two of the four faults identified in the trench are indicated in the resistivity and GPR sections. A 3-D GPR measurement has identified a spot of high amplitudes elongated parallel to the MLF trend, which coincides with the high resistivity body. To improve the signal-to-noise ratio of the time slices we stacked the GPR time slices within vertically homogeneous blocks. This provided a contrast image of the sand-gravel body including its boundaries in three dimensions. The detailed analysis of the 3-D GPR cube revealed additional fault that limits the highly reflective sands and appears to be offset by another younger fault. Our results suggest a complex fault pattern in the studied area, which deserves a further study.     

Resistivity imaging of Pleistocene alluvial aquifers in a contractional tectonic setting: A case history from the Po plain (Northern Italy)

In this work we present the hydrogeophysical imaging of a key sector of the Quaternary Po foreland basin (northern Italy), focussing on the reconstruction of clastic aquifers and aquitards in a complex tectono-sedimentary subsurface architecture. The study area includes the relic reliefs of Casalpusterlengo and Zorlesco, two smooth morphological features involving uplifted and gently folded Pleistocene marine to alluvial sediments, plausibly linked to the buried Northern Apennines thrust and fold belt. The geophysical data include 35 Direct Current Vertical Electrical Soundings collected over a 37 km² wide area, acquired with Schlumberger array and maximum half-spacing of 500 m. 1-D resistivity-depth profiles were computed for each VES. An integrated hydrostratigraphic approach was applied, to constrain the interpretation of the geophysical data along several cross-sections, including the comparison of resistivity soundings to stratigraphic logs, borehole electric logs and the pore-water properties.The resistivity interfaces, traceable with the same laterally continuous vertical polarity, were used to develop an electrostratigraphic model in order to portray the stacking of electrostratigraphic units down to 200 m below ground surface. Their vertical associations show a general upward increase of electrical resistivity. This assemblage mimics the regional coarsening upwards depositional trend, from the conductive units of the Plio-Pleistocene marine-to-transitional depositional systems to the resistive units of the Middle–Late Pleistocene fluvial and alluvial plain depositional systems. Middle Pleistocene depositional systems host an alternation of North-dipping, high-to-intermediate permeability aquifer systems (70–180 Ωm, thickness of 5–70 m) separated by low permeability aquitards (20–50 Ωm, thickness up to 40 m). These units pinch out against the Casalpusterlengo and Zorlesco relic reliefs, where they cover the uplifted and folded regional aquitard (20–50 Ωm) formed by Pliocene-Lower Pleistocene clays to sandy silts with gravel lenses in agreement with borehole data. In the deepest part of the local stratigraphy, a broad low-resistivity anomaly (< 10 Ωm) was clearly mapped through the study area. By comparison with electrical borehole logs in deep oil-wells, it could be interpreted as the fresh–saltwater interface due to the presence of connate waters and brines hosted by the marine-to-transitional shales.