Interpretation of ground VLF-EM data in terms of inclined sheet-like conductor models

The response of two-dimensional, inclined, sheet-like conductors with, low conductance values to plane wave electromagnetic fields in the very low frequency (VLF) range has been evaluated by using a numerical technique. The conductance values of the conductors considered are appropriate for those produced by water and/or clay-filled fracture and shear zones in the Precambrian crystalline rocks of the Canadian Shield. The surrounding host rock was assumed to be, resistive with resistivities in the 1–10 k.m range to reflect the high resistivities over the shield areas. No overburden was assumed in this analysis.The results of the computations are presented in the form of characteristic interpretation diagrams to interpret ground VLF data in the field, where facilities for direct numerical modelling may not be available. A method for interpreting ground VLF data using such characteristic diagrams has been proposed in this paper which requires a prior knowledge of the host rock resistivity and the inclination of the conductor. These two parameters may be derived from a VLF resistivity survey and from appropriate filtering of the VLF tilt angle response. The interpretation method was applied to a ground VLF anomaly obtained at a research site near Atikokan in NW Ontario, which yielded an interpretation compatible with information from geological mapping.Geological Survey of Canada Contribution No. 51888.     

Interpretation of VLF Resistivity Data for Ground Water Contamination Surveys

Very low frequency (VLF) military communications systems provide a primary field that can be used for shallow geophysical surveys to locate ground water contamination and vertical geologic contacts. Useful properties that can be easily obtained from the interaction of the earth and the primary field are the magnitude of the vertical secondary magnetic field, the surface impedence, and the phase angle between the electrical and magnetic horizontal components. The variations in the secondary magnetic field can be related to vertical geologic contacts, such as the edges of landfill trenches. The surface impedence yields an apparent terrain conductivity, which can be used to locate low-resistivity anomalies often associated with contaminated ground water. The phase angle gives information on vertical variations in resistivity, phase angles less than 45° indicating increasing resistivity with depth. The depth of penetration of the VLF field is about one skin depth. For a frequency of 20 kHz, the skin depth in meters is approximately equal to 3.67 where p is terrain resistivity in ohmmeters.     

A comparison of helicopter-borne electromagnetics in frequency- and time-domain at the Cuxhaven valley in Northern Germany

Two different airborne electromagnetic methods were applied in the same area: the frequency-domain helicopter-borne electromagnetic (HEM) system operated by the Federal Institute for Geosciences and Natural Resources, Germany, and the time-domain SkyTEM system of the HydroGeophysics Group at the University of Aarhus, Denmark. For verification of and comparison with the airborne methods, ground-based transient electromagnetics and 2-D resistivity surveying were carried out. The target of investigation was the Cuxhaven valley in Northern Germany, which is a significant local groundwater reservoir. The course of this buried valley was revealed by drillings and the shape was determined by reflection seismics at several cross sections.We applied electrical and electromagnetic methods to investigate the structure of the valley filling consisting of gravel, sand, silt and clay. The HEM survey clearly outlines a shallow conductor at about 20m depth and a deeper conductor below 40m depth inside the valley. This is confirmed by 2-D resistivity surveying and a drilling. The thickness of the deeper conductor, however, is not revealed due to the limited investigation depth of the HEM system. The SkyTEM survey does not resolve the shallow conductor, but it outlines the thickness of the deeper clay layer inside the valley and reveals a conductive layer at about 180m depth outside the valley. The SkyTEM results are very consistent with ground-based transient electromagnetic soundings.Airborne electromagnetic surveying in general has the advantage of fast resistivity mapping with high lateral resolution. The HEM system is cost-efficient and fast, but the more expensive and slower SkyTEM system provides a higher depth of investigation. Ground-based geophysical surveys are often more accurate, but they are definitively slower than airborne surveys. It depends on targets of interest, time, budget, and manpower available by which a method or combination of methods will be chosen. A combination of different methods is useful to obtain a detailed understanding of the subsurface resistivity distribution.     

Airborne VLF Measurements and Variations of Ground Conductivity: A Tutorial

Airborne very low frequency (VLF) data are routinely collected by national agencies and commercial companies together with other passive geophysical measurements of the static magnetic field and radiometric data. The purpose of this paper is to demonstrate that both standard three-component VLF and tensor VLF (TVLF) data contain a lot of useful quantitative and qualitative information about the electrical conductivity distribution in the upper few hundred meters of the crystalline basement. We first give a new derivation of the fundamental transfer functions (the tipper) used in the TVLF technique. We then show that the tipper can be estimated from simultaneous measurements of the wave magnetic fields from at least two transmitters with somewhat different frequencies, and present a simple model by which the maximum error introduced by the difference in frequencies can be found. Single transmitter scalar VLF maps emphasise those conductive structures that have dominant strikes in the direction of the transmitter. Multiple transmitter transfer functions are dependent only upon the underlying conductivity structure. Two dimensional structures can be quantitatively modelled by modern inversion methods developed originally for deep electromagnetic magnetotelluric (MT) soundings. In such cases three-component VLF measurements can be modelled easily upon appropriate rotation of the co-ordinate system to “strike” co-ordinates. Single frequency transfer functions (tippers) have real and imaginary parts that carry information on not only lateral contrasts in conductivity, as usually stated in text books, but, taken together, they provide a robust tool for determining the background conductivity level away from distinct conductors, and they can also be used to discriminate between deep and shallow conductors. Based upon simulations using multi-frequency data, it can be concluded that such a new development would dramatically increase the resolving power of airborne VLF measurements.     

A PHYSICAL SCALE MODEL STUDY OF THE COMPARATIVE PERFORMANCE OF TWO MODES OF OPERATION FOR FIXED-LOOP TURAM-TYPE EM SYSTEMS1

The results of a physical scale model study of the conventional mode of operation of fixed loop electromagnetic systems and an alternative mode called the tx-parallel mode in which traverses are run parallel to the long axis of the rectangular transmitter loop are presented. The results show that over thick or dipping conductors, the tx-parallel configuration provides coupling with the target which is comparable with that provided by the conventional configuration. In addition, the tx-parallel configuration is shown to provide more consistent indications of the direction and magnitude of conductor dip. Over wide conductors, where separate conventional surveys are needed to define the opposite edges of such conductors, it is shown that only a single tx-parallel survey is needed to locate both edges of the conductor. The tx-parallel results were found to allow better resolution of the individual anomalies due to closely spaced parallel conductors. The tx-parallel response of identical sheet conductors of opposite dip indicated that the response of the separate sheets could be recognized even when the two sheets were placed at zero separation. This was found to be due in part to spatial displacement of the individual current vortices within each conductor owing to their mutual repulsion.     

ELECTROMAGNETIC PROSPECTING FOR GROUNDWATER IN PRECAMBRIAN TERRAINS IN THE REPUBLIC OF UPPER VOLTA*

The “Autorité des amenegements des valées des Voltas (AVV)” is establishing new rural settlements in the Volta valleys. First, a survey of available water supplies is performed. Economic aquifers in Precambrian terrains are deep (15–50 m) and usually occur in fractured zones accompanying faults. Such zones can be identified on aerial photographs, but their precise location on the ground is virtually impossible by visual means. Because of the small size of the aquifers, a location error of 5 m can make the difference between a productive well and a dry hole. Traditionally, resistivity profiling has been used as the means of locating the fractured zones in the field. Our studies suggest that the task can be performed faster, cheaper and more accurately by VLF and EM methods. Because of the limited choice of transmitting stations reccivable in Upper Volta, the VLF method is not sufficiently sensitive to detect conductors with a strike between 45° and 105°. The results obtained with a multifrequency, horizontal-loop EM (HLEM) system were satisfactory in all investigated areas. During the 1980 field season, 35 target areas were surveyed. Of the 24 holes drilled so far, 23 are productive. The weathered layer is a source of distinctive HLEM anomalies, which are characteristic of the underlying rocks. Therefore, different interpretational procedures had to be developed for granitic and volcano-sedimentary areas. Despite the high background level of in-phase and quadrature components, which varied with thickness and conductivity of the weathered layer, aquifers could be detected at a depth greater than 30 m. Attempts were made to interpret the HLEM results quantitatively using two models: a three-layer medium and a valley discontinuity. The latter model is more realistic, but more scale modelling will have to be performed to permit development of viable interpretational procedures. Meanwhile, phasor diagrams based on drilling and resistivity sounding data can be used to estimate the aquifer depth.     

太阳活动低年电离层磁场VLF波的观测特性研究

本文利用搭载在DEMETER卫星上的感应式磁力仪(Instrument Magnetometer Search-Coil,IMSC)探测数据分析了磁场甚低频(Very Low Frequency,VLF)波功率谱的空间分布.在排除地磁扰动影响(Dst≤-30 nT,Kp≥3,AE≥200 nT)的前提下,我们给出200...     

MAGNETIC WAVETILT MEASUREMENTS FOR GEOLOGICAL FRACTURE MAPPING1

The measurement of wavetilt is diagnostic for determining the electrical characteristics of the upper layers of the ground at VLF and LF frequency ranges. Theoretical and field studies have indicated that electric wavetilt using the transverse magnetic (TM) waves detects lateral inhomogeneities virtually instantly as abrupt changes in electrical properties are encountered. Theoretical studies have also indicated that magnetic wavetilt measurements using transverse electric (TE) waves are superior to electric wavetilts for such purposes. An experimental survey was conducted at two locations near Atikokan, Ontario, to verify the theoretical predictions. The survey area, forming a part of a large granitic pluton, was mapped earlier by various geophysical techniques, including the ground VLF-EM method, to detect weak conductors formed either by the presence of fractures in the bedrock filled with water and/or clay, or overburden filling bedrock depressions. A small, multi-turn, horizontal loop was used during the survey as the transmitter to generate TE waves at eleven frequencies from 10.7 to 58.5 kHz. The magnetic wavetilt measurements detected all previously known conductors at the two locations. In addition, the survey detected several weak conductors that were missed by the VLF survey. Thus, the survey indicated the usefulness of magnetic wavetilt results for detection of weak conductors at shallow depths, which may have application in engineering geophysical surveys. The multi-frequency wavetilt data also provided some indications of the depth and depth extent of such conductors.     

STATISTICAL ANALYSIS OF AIRBORNE ELECTROMAGNETIC SIGNATURES OF DIFFERENT CONDUCTOR GROUPS IN TROPICAL TERRAIN OF KARNATAKA,INDIA*

AEM (airborne electromagnetic) anomalies caused by massive sulphide conductors and superficial conductors can be recognized with a statistical method, as shown by an analysis of Input AEM data from Karnataka State. The weathering in the survey area is of tropical type. Parameters, such as various amplitude ratios and time parameters (inverse of decay rate) for exponential and power-law decay were analyzed for sulphide bodies, conducting soil, superficial conductors, and cultural conductors. Time parameters τ₁ (exponential decay) is defined as ratio of time differences between the third and fourth channel to the logarithmic value of the relative amplitude of the two channels. Time parameter K₁ (power-law decay) is defined as ratio of the difference of the logarithmic values of the delay times of the third and fourth channels to the logarithmic value of the relative amplitude of the two channels. The two parameters have been useful in recognizing sulphide conductors. Also the first channel Input amplitude and logarithmic plot of the transients appear to be helpful in conductor identification. Channel ratios seem to be the least effective parameters of conductor identification. In the area studied both power-law and the conventional exponential decay were found equally suitable for approximating Input AEM transients.     

Solar cycle changes in daytime VLF subionospheric attenuation

Observations of field strengths of signals from a number of VLF transmitters, after propagation over long paths in the Earth–ionosphere waveguide, have been used to examine changes in the daytime attenuation rate in the course of a solar cycle. The measurements reported were recorded in the period 1986–1996. The paths studied range in length over about 8–14 Mm; they included NLK (Seattle) and NPM (Hawaii) to Dunedin, NZ, and NAA (Cutler, Maine) and NSS (Annapolis, Maryland) to Faraday, Antarctica. The frequencies monitored were mainly in the range 21–25 kHz but measurements near 10 kHz of Omega Hawaii at Faraday were also used. The daytime VLF attenuation rates at solar minimum were found to be greater by about 0.3 dB/Mm than at solar maximum.     

飞行高度同时反演的固定翼航空瞬变电磁-维反演

航空电磁测量记录中,不仅感生电动势测量数据有观测误差,而且高度计测量数据也有误差,直接进行常规反演往往导致反演结果不可靠,研究飞行高度数据有误差下的反演算法具有实际意义.本文以层状模型的固定翼时间域航空电磁多分量理论响应数据为例,提出了两种针对飞行高度计记录数据有误差时的正则化反演算法,一个是自适应正则化反演方法,另一...     

Resistivity and self-potential changes associated with volcanic activity: The July 8, 2000 Miyake-jima eruption (Japan)

The Miyake-jima volcano abruptly erupted on July 8, 2000 after 17 years of quiet and gave birth to a crater, 1 km in diameter and 250 m deep. This expected unrest was monitored during the years 1995–2000 by electromagnetic methods including DC resistivity measurements and self-potential (SP) surveys. Beneath the 2500 yr old Hatcho-Taira summit caldera audio-magnetotelluric soundings made in 1997–98 identified a conductive medium, 200–500 m thick (within the 50 Ω m isoline) located at a few hundred metres depth. It was associated with the active steady-state hydrothermal system centred close to the 1940 cone and extending southward. A DC resistivity meter set in a Schlumberger array with 600, 1000 and 1400 m long injection lines evidenced strong resistivity changes between September 1999 and July 3, 2000 in the vicinity of the newly formed crater. The apparent resistivity has reached about three times its initial values on the 1400 m long line and has lowered to about 20% on the 600 m line. Just prior to the July 8, 2000 eruption SP mapping made inside the summit Hatcho-Taira caldera revealed negative anomalies where positive ones had occurred during the previous tens of years. The largest negative anomaly, −225 mV in amplitude, mainly took place above the 1940 cone which collapsed in the crater formation. A permanent 1 km long SP line across the caldera suggests accelerating changes during the 3 months preceding the eruption. On a larger scale, the comparison between 1995 and 2000 surveys has shown a global increase of the hydrothermal activity beneath the volcano. Its source could have been 250 m to the south of the crater. These observations suggest that the hydrothermal system was slowly disturbed in the months preceding the eruption while drastic changes have occurred during the 2 weeks before the summit collapse when tectonic and volcanic swarms have appeared.     

A three-dimensional crustal geoelectric model of the Ukrainian Shield

A 3-D crustal geoelectric model of the Ukrainian Shield (USh) is constructed from magnetic variation data. High electrical conductivity anomalies with resistivities of 1–100 Ω m are located at depths of up to 30 km from the basement surface. A high conductivity layer with ρ = 25 Ω m and with its upper boundary at a depth of 70 km is supposed to exist in the upper mantle of the southwestern USh.     

INTERPRETATION OF PULSE EM ANOMALIES OVER GANI CONDUCTORS,KURNOOL, INDIA*

A survey using Crone pulse electromagnetic (PEM) equipment was conducted near the village of Gani in Andhra Pradesh, India. The anomalies obtained had a shape typical for an inclined sheet conductor. From the multichannel TEM response, the multifrequency Slingram response was calculated. A comparison of the Slingram response of the Gani conductors with the free-air model response of a thin sheet conductor shows phase rotation and attenuation of the anomaly vector, thereby confirming the presence of a conductive overburden. Interpretation of the Gani results using free-air phasor diagrams causes an overestimation of depth and conductance of the target conductor unless data obtained at a sufficiently low frequency are used for the interpretation. For areas like Gani this frequency can be as high as 336 Hz.     

Integrated Geophysical Studies in the East-Indian Geothermal Province

Integrated geophysical surveys using vertical electrical sounding (VES), very low frequency (VLF) EM, radiation counting, total magnetic field and self-potential (SP) measurements are carried out to characterize the geothermal area around a hot spring in the Nayagarh district, Orissa, India that lies in the East Indian geothermal province. The study was performed to delineate the fracture pattern, contaminated groundwater movement and possible heating source. VES interpretations suggest a three- to four-layer structure in the area. Resistivity survey near the hot spring suggests that weathered and fractured formations constitute the main aquifer system and extend to 60 m depth. Current flow measured at various electrode separations normalized by the applied voltage suggests that fractures extend to a greater depth. Detailed VLF study shows that fractures extend beyond 70 m depth. VLF anomaly has also very good correlation with the total magnetic field measured along the same profiles. Study results suggest that a gridded pattern of VLF survey could map the underground conductive fracture zones that can identify the movement of contaminated groundwater flow. Therefore, precautionary measures can be taken to check further contamination by delineating subsurface conducting structures. Self potential (SP) measured over the hot spring does not show a large anomaly in favor of the presence of a sulphide mineral body. A small positive (5–15mV) SP anomaly is measured which may be streaming potential due to subsurface fluid flow. A high radiation is measured about four kilometers from the hot spring, suggesting possible radiogenic heating. However, the exact nature of the heating source and its depth is not known in the area. Deep resistivity followed by a magneto-telluric survey could reveal the deeper structures.     

Probing the lower ionospheric perturbations associated with earthquakes by means of subionospheric VLF/LF propagation

There have been reports for many years that the ionosphere is very sensitive to seismic effects, and the detection of ionospheric perturbations associated with earthquakes (EQs) attracts a lot of attention as a very promising candidate for short-term EQ prediction. In this review we present a possible use of VLF/LF (very low frequency (3–30 kHz)/low frequency (30–300 kHz)) radio sounding of seismo-ionospheric perturbations. In order to avoid the overlapping with my own previous reviews, we first show some pioneering results for the Kobe EQ and we try to present the latest results including the statistical evidence on the correlation between the VLF/LF propagation anomalies (ionospheric perturbations) and EQs (especially with large magnitude and with shallow depth), medium-distance (6-8 Mm) propagation anomalies, the fluctuation spectra of subionospheric VLF/LF data (the effect of atmospheric gravity waves, the effect of Earth's tides, etc.), and the mechanism of lithosphere-atmosphere-ionosphere coupling. Finally, we indicate the present situation of this kind of VLF/LF activities going on in different parts of the globe and we suggest the importance of international collaboration in this seismo-electromagnetic study.     

LINEAR FILTERING OF VLF DIP-ANGLE MEASUREMENTS*

The suggested interpretation technique is based on discrete linear filtering of VLF data. The output of the described filtering results is expressed in terms of an equivalent current density at a specified depth that would cause the measured magnetic field. The most practical six-point filter gives an accuracy of 8%. The filter is an extension of the commonly used Fraser filter to process VLF dip-angle data. Filtering the same data set for various depths gives an idea about the change of current density with depth. Areas with high current-density correspond to good conductors. The conductor dip can also be determined. The use of the technique is illustrated on theoretical and field examples. In all cases a good correlation with original models and other types of geophysical measurements was obtained. As shown in the last example, the filtering technique is also applicable in interpretation of other electromagnetic methods.     

VLF RESISTIVITY MAPPING AND VERTICALIZATION OF THE ELECTRIC FIELD1

For over 20 years, powerful VLF transmitters have been used as electromagnetic sources for subsurface investigations in mining exploration. Measurements initially concerned the vertical component of the magnetic field or the inclination of the field and were later extended to measurement of the horizontal electric field in the direction of the transmitter, to determine the resistivity of the terrain. Measurement of the electric field is usually performed with electric lines, grounded or not, with lengths of at least 5 m. This paper presents the concept of a VLF resistivity meter with a very short electric sensor (1 m) and the results obtained with it. This technique improves the measurement of the electric field, which is in principle a point value. It also permits a higher spatial sampling rate and, by closely linking the electric sensor with the magnetic sensor on a lightweight mount, makes it possible for the instrument to be used by a single operator. In addition, transformation of the electric field data, analogous to reduction to the pole in magnetism, is proposed to correct the horizontal deformation of the anomalies created by polarization of the primary field. Comparison with direct current electrical measurements shows highly satisfactory correlations. This transformation, validated for VLF, can be extended to any electrical or electromagnetic method using a uniform primary field, i.e. gradient array in direct current or low-frequency magnetotellurics. We call this verticalization of the electric field. Resistivity measurements and mapping using the VLF frequency range can be applied not only to mining but also to a wide range of shallow geophysical studies (hydrology, civil engineering, etc.) and are not limited to problems concerning the location of conductive targets     

Geophysical techniques to aquifer locating and monitoring for industrial zones in North Hanoi, Vietnam

Geophysical methods were applied for hydrogeological targets in many countries including Vietnam. This paper presents results of using complex geophysical techniques as well as 2D electrical resistivity imaging (ERI), vertical electrical sounding (VES), very low frequency (VLF), and seismic refraction for geological structure investigation for locating the aquifers and assessing the hydrogeological conditions for groundwater potential in industrial zones of North Hanoi, Vietnam. The locations of two aquifers are determined by their depth and thickness on the basis of resistivity and seismic velocity values which were proved by stratifications of three boreholes to 40–60 m of depth on the study area. There are connections from surface water to shallow aquifer by hydraulic windows, as follows from VLF data. The deeper aquifer can be considered as a potential groundwater supply, but the water level is descending in time, as shown by hydrological monitoring. However, with careful use and by reducing sources of pollution, groundwater can continue to be an important natural resource for future.     

Airborne VLF measurements and mapping of ground conductivity in Sweden

Airborne VLF data are routinely collected by The Geological Survey of Sweden (SGU) as part of its bedrock mapping programme. In this paper we demonstrate that the novel Tensor VLF technique developed at Uppsala University and SGU can provide useful qualitative and quantitative information about the electrical conductivity distribution in the upper few hundred meters. Single transmitter scalar VLF maps emphasize those conductive structures that have dominant strikes in the direction of the transmitter. The tensor tipper (essentially the vertical magnetic field from currents along the strike direction) calculated from multiple transmitters is dependent only upon the underlying conductivity structure. Transformation of the tipper into the peaker (the horizontal divergence) has proven to enhance the lateral resolution while the transformation to the apparent resistivity can be used to discriminate different rock types. Two case histories from the application of VLF data are presented in this study. Two dimensional structures can be quantitatively modelled by modern inversion methods developed originally for deep electromagnetic MT soundings. Direct inversion of the real and imaginary parts of the tipper provides more quantitative information about the subsurface resistivity distribution.