The instantaneous apparent resistivity tensor: a visualization scheme for LOTEM electric field measurements

Long-offset transient electromagnetic (LOTEM) data have traditionally been represented as early- and late-time apparent resistivities. Time-varying electric field data recorded in a LOTEM survey made with multiple sources can be represented by an 'instantaneous apparent resistivity tensor'. Three independent, coordinate-invariant, time-varying apparent resistivities can be derived from this tensor. For dipolar sources, the invariants are also independent of source orientation. In a uniform-resistivity half-space, the invariant given by the square root of the tensor determinant remains almost constant with time, deviating from the half-space resistivity by a maximum of 6 per cent. For a layered half-space, a distance–time pseudo-section of the determinant apparent resistivity produces an image of the layering beneath the measurement profile. As time increases, the instantaneous apparent resistivity tensor approaches the direct current apparent resistivity tensor. An approximate time-to-depth conversion can be achieved by integrating the diffusion depth formula with time, using the determinant apparent resistivity at each instant to represent the resistivity of the conductive medium. Localized near-surface inhomogeneities produce shifts in the time-domain apparent resistivity sounding curves that preserve the gradient, analogous to static shifts seen in magnetotelluric soundings. Instantaneous apparent resistivity tensors calculated for 3-D resistivity models suggest that profiles of LOTEM measurements across a simple 3-D structure can be used to create an image that reproduces the main features of the subsurface resistivity. Where measurements are distributed over an area, maps of the tensor invariants can be made into a sequence of images, which provides a way of 'time slicing' down through the target structure.     

Magnetotelluric Soundings In the Canadian Rocky Mountains

Summary. Magnetotelluric soundings have been made at 25 stations in the Rocky Mountain Trench (RMT) and Main Ranges near 53° N, close to the centre of a major conductivity anomaly which had been mapped in a magnetovariation array study. Most stations covered the frequency range 0.01–500 Hz and three stations 0.0002–500 Hz. the resistivity tensor shows low to moderate anisotropy in the RMT, but is strongly 2-D or 3-D in the Rocky Mountains. Apparent resistivities as a function of frequency are displayed in pseudosections along the Trench and along a transverse profile across the RMT and into the Main Ranges. In preparation for 2-D modelling, 1-D inversions have been used to construct resistivity-depth sections satisfying both magnitudes and phases of the MT responses. These show very low resistivities, in the range 1—10Ωm, in the upper crust under the RMT and even lower values under the Main Ranges. the latter values give strong confirmation of the Northern Rockies conductor reported by Bingham, Cough & Ingham and are in agreement with models of the conductors fitted to long-period magnetovariation fields by Ingham, Gough & Parkinson. the MT results here reported add some essential depth and resistivity information. It is suggested that the conductors beneath the Rocky Mountains Main Ranges and Trench constitute a thickening at the edge of the Canadian Cordilleran Regional (CCR) conductor. Gough has argued that a wide variety of geophysical and geological parameters indicate high temperatures and partial melting in the mantle under the CCR conductor. At the upper crustal depths penetrated in this magnetotelluric study, it is considered more probable that the high conductivity is caused by hot, saline water of mantle origin rather than silicate melt. the CCR in general may have two layers of fluid producing its high conductivity, silicate melt below and saline hot water above.     

Three-dimensional topography corrections of magnetotelluric data

Topographic effects due to irregular surface terrain may prevent accurate interpretation of magnetotelluric (MT) data. Three-dimensional (3-D) topographic effects have been investigated for a trapezoidal hill model using an edge finite-element method. The 3-D topography generates significant MT anomalies, and has both galvanic and inductive effects in any polarization. This paper presents two different correction algorithms, which are applied to the impedance tensor and to both electric and magnetic fields, respectively, to reduce topographic effects on MT data. The correction procedures using a homogeneous background resistivity derived from a simple averaging method effectively decrease distortions caused by surface topography, and improve the quality of subsurface interpretation. Nonlinear least-squares inversion of topography-corrected data successfully recovers most of structures including a conductive or resistive dyke.     

Results of magnetotelluric and gravimetric measurements in western Nicaragua, Central America

Magnetotelluric and gravity data have been collected within a ca. 170 km long traverse running from the Pacific coast of Nicaragua in the west to the Nicaraguan Highland in the east. This part of Nicaragua is characterized by sedimentary rocks of the Pacific Coastal Plain, separated from the Tertiary volcanic rocks of the Highland by the NW-SE-trending Nicaraguan Depression. 2-D interpretation of the magnetotelluric (MT) data, collected at 13 stations, indicates four regions of high electrical conductivity in addition to the conductive coastal region. Two of these are associated with conducting sediments and pyroclastics in the upper part of the crust. Two other conductive structures have been defined at depth around 20 km and the one best defined is located below the depression. From the distribution of seismic events, volcanic activity in the depression and the similarity in geophysical characteristics with areas such as the Rio Grande Rift, this conductor is interpreted as a melt layer or a complex of magma chambers. Models of the upper lithosphere, constrained by the MT model, vertical electrical sounding (VES) data, seismic data and densities, have been tested using gravity data. A model that passes this test shows a gradual thickening of the crust eastwards from the Pacific coast. An anomaly centred over the depression is interpreted to have its origin in a thinning of the crust. In this model the melt layer is situated on top of the bulge of the lower lithosphere. A change in the composition of the crust, from the Pacific Coastal Plain to the Highland, is indicated from the change in character of the MT response and from the density distribution in the gravity model. This may support the hypothesis that the Pacific region is an accreted terrane. MT and gravity data indicate a depth to a resistive and high-density basement in the depression of ca. 2 km. On the basis of this, the vertical setting in the depression is estimated to be of the order of 2.5 km.     

Interpretation of long-offset transient electromagnetic data from the Odenwald area, Germany, using two-dimensional modelling

The standard 1-D inversion approach for the interpretation of transient electromagnetic (TEM) data usually fails in the presence of near-surface conductivity anomalies. Since multidimensional inversion codes are not routinely available, the only alternative to discarding the data may be trial-and-error forward modelling. We interpret data from a long-offset transient electromagnetic (LOTEM) survey which was carried out in 1995 in the Odenwald area, using 2-D finite-difference modelling. We focus on a subsegment of the LOTEM profile, which was shot with two different electric dipole transmitters. A model is found which consistently explains the electric and magnetic field data at eight locations for both transmitters. First, we introduce a conductive dyke under the receiver spread to explain sign reversals in the magnetic field transients. A conductive slab under one of the transmitters is required to obtain a reasonable quantitative fit for that transmitter. Consideration of the electric field data then requires a modification of the layered earth background. Finally, we study the response of a crustal conductor, which was the original target of the survey. The data are sensitive to the conductor, and for the investigated subset of the data the fits are slightly better without the conductive layer.     

Deep Schlumberger sounding and the crustal resistivity structure of central Australia

Summary. Three 200 km Schlumberger resistivity soundings have been conducted over the central Australian shield, using telephone lines to obtain the large electrode spacings. These represent the first crustal scale controlled source electrical study to be carried out in this continent. A computer controlled data acquisition system was used which allowed precise measurements to be made with only modest emission currents (0.1–0.5 A).
The three soundings, centred on the towns of Renner Springs, Wauchope and Aileron, showed the southern part of the study area (the Arunta Block) to be an order of magnitude more resistive than the more northerly section (the Tennant Creek Block). This difference correlates with the higher heat flow of the Tennant Creek Block. A lowering of apparent resistivity at large electrode spacings for one sounding (Wauchope) is taken to indicate the presence of a low resistivity layer in the middle crust, at a depth less than 20 km. However, the effect of the highly conductive overburden characteristic of inland Australia, combined with the large transverse resistance of the crust, prevented the other two soundings from detecting such a layer. Without support from these two soundings, it is impossible to be sure that the lowered resistivity at Wauchope is not caused merely by lateral variations in near-surface resistivity.
The data also show that crustal resistivities are much lower than the expected values for dry rock, whether or not a low resistivity layer is included in the model. This implies a widespread occurrence of free water in the crust, with greater amounts occurring at depth if the low resistivity zone exists.     

3-D magnetotelluric inversion and model validation with gravity data for the investigation of flood basalts and associated volcanic rifted margins

20 magnetotelluric (MT) soundings were collected on the Isle of Skye, Scotland to provide a high-resolution three-dimensional (3-D) electrical resistivity model of a volcanic province within the framework of a project jointly interpreting gravity, seismic, geological and MT data. The full 3-D inversion of the MT data jointly interpreted with gravity data reveals upper crustal structure. The main features of the model are interpreted in conjunction with previous geological mapping and borehole data. Our model extends to 13 km depth, several kilometres below the top of the Lewisian basement. The top of the Lewisian basement is at approximately 7–8 km depth and the topography of its surface was controlled by Precambrian rifting, during which a 4.5 km thick sequence of Torridonian sediments was deposited. The Mesozoic sediments above, which can reach up to 2.2 km thick, have small-scale depocentres and are covered by up to 600 m of Tertiary lava flows. The interpretation of the resistivity model shows that 3-D MT inversion is an appropriate tool to image sedimentary structures beneath extrusive basalt units, where conventional seismic reflection methods may fail.     

Magnetotelluric image of the crust and upper mantle in the backarc of the northwestern Argentinean Andes

Magnetotelluric data from the backarc of the Central Andes in NW Argentinawere re-examined by employing impedance tensor decomposition and 2-D inversion and modelling techniques. The data in the period range of 50–15 000 s were collected on a profile of 220 km length reaching from the Eastern Cordillera across the Santa Barbara System to the Andean foreland of the Argentinean Chaco.
After a dimensionality analysis, data from most sites were treated as regional 2-D. The exception was the eastern section of the profile, where the magnetotelluric transfer functions for periods ≤ 1000 s reflect a 3-D earth. Application of two tensor decomposition schemes yielded a regional strike direction of N–S, which is the azimuth of the Central Andean mountain chains. Several 2-D models were obtained by pseudo- and full 2-D Occam inversion schemes. Special emphasis was placed on the inversion of phase data to reduce the influence of static shifts in the apparent resistivity data. The smooth inversion models all show a good conductor at depth. A final model was then calculated using a finite element forward algorithm.
The most prominent feature of the resulting model is a conductor which rises from depths of 180 km below the Chaco region to 80 km beneath the Santa Barbara System and the Eastern Cordillera. Its interpretation as a rise of the electrical asthenosphere is supported by seismic attenuation studies. Magnetotelluric results, surface heat-flow distribution in the area, and the electrical properties of crustal and mantle rocks suggest that the upper mantle is predominantly ductile beneath the Eastern Cordillera and the western Santa Barbara System. This generally agrees with anelastic seismic attenuation models of the area and is useful in discriminating between models of Q quality factor distribution.     

Investigation of a Tertiary maar structure using three-dimensional resistivity imaging

This paper presents the application of the electrical resistivity tomography (ERT) method to the investigation of the Tertiary maar structure of Baruth (Germany) known from previous gravimetric surveys. ERT was applied to support the optimum location for a palaeoclimatological drill hole.
  Special modifications of data acquisition, signal processing and inversion are introduced to adapt the method of ERT to the special requirements for the 3-D investigation of structures with horizontal extensions of 1  km or more. More than 5000 dipole–dipole combinations were recorded at three concentric circular electrode arrangements using stand-alone transient data acquisition systems (RefTek).
  We present a fast approximate imaging technique based on the simultaneous iterative reconstruction technique (SIRT). As the complete calculation of the inverse Frechét matrix is avoided, the algorithm is especially suitable for large data and model spaces, where complete inversion is beyond the limits of available computing hardware. The single-step method is applicable to arbitrary irregular electrode layouts. Synthetic tests show that the imaging procedure reconstructs the main features of the subsurface.
  A low-resistivity body could be interpreted as limnic sediments filling the interior of the Tertiary maar crater. Considering the horizontal resistivity gradient, estimates for the lateral and depth extents of the structure were made. An optimum position for a palaeoclimatological borehole was found, and was in good agreement with the gravimetric minimum.     

A magnetotelluric study of the Alpine Fault, New Zealand

Magnetotelluric soundings have been made at seven locations on a 4  km profile crossing the Alpine Fault in the South Island of New Zealand. The 'distortion' techniques of Groom & Bailey (1989 ) and Lilley (1998a , b ) have been used to derive regional apparent resistivity and phase curves that correspond to electromagnetic induction in orientations parallel and perpendicular to the fault. 2-D inversion of the regional responses reveals that a narrow (<1  km wide) conductive zone is associated with the Alpine Fault. This conductor is most probably related to the heating of deep circulating meteoric water in a region in which enhanced temperatures occur at shallow depth due to the tectonic uplift of the Southern Alps.     

A comparison of the stability of stationary and non-stationary magnetotelluric analysis methods

A Magnetotelluric (MT) sounding was carried out at a site in south-east Queensland, in the Clarence-Moreton Basin. The synoptic recordings were taken over a period of four months at sampling frequencies from 500 Hz to 5 × 10-5 Hz. The resulting data was analysed by the stationary cross-frequency and the Cone kernel time-frequency distribution (TFD) methods of MT analysis. The results were compared as apparent resistivities on a daily basis for frequencies above 1 Hz, as well as over all the available data. The TFD MT apparent-resistivity results were more stable and less noisy on an daily basis than the cross-frequency results. Similarly the TFD analysis gave less noisy results than the cross-frequency analysis when all available data was processed. Application of these new non-stationary analysis techniques to MT processing should decrease the bias error problem of the MT methods and so increase reliability and repeatability of MT soundings.     

The time-domain magnetotelluric response over three-layer earth models

Magnetotelluric data observed in frequency and time domains are expressed as apparent resistivity. the apparent resistivity is a weighted spatial average of the subsurface resistivity distribution. In this paper, we develop analytical expressions to compute the apparent resistivity in the time domain for various three-layer earth models. the present approach to computing the magnetotelluric response in the time domain is found to overcome the problems encountered by the method of images. the magnetotelluric response in the time domain for various three-level models have been computed and shown graphically. the time-domain responses show a characteristic behaviour with a small change in layered parameters (resistivity and thickness of the layers), whereas frequency-domain responses do not show such behaviour. This characteristic behaviour of time-domain magnetotelluric sounding curves will be useful in the qualitative interpretation of field data.     

A geomagnetic scattering theory for evaluation of earth structure

Summary. Structural features of the Earth's lower crust and upper mantle can be mapped by the analysis of temporal geomagnetic fluctuations using the electromagnetic scattering theory developed in this paper. Decomposing geomagnetic field fluctuations at the Earth's surface into an excitation part and a scattered part forms the basis of a power series development. The vertical field component is interpreted as a scattering of the excitation field. The horizontal gradient and geomagnetic depth sounding methods are special cases of the theory developed. The horizontal gradient sounding method has a tensorial aspect which has not been recognized before; it should be included to obtain correct penetration depth parameter evaluations from field data.     

Interpreting magnetic data by integral moments

The use of the integral moments for interpreting magnetic data is based on a very elegant property of potential fields, but in the past it has not been completely exploited due to problems concerning real data. We describe a new 3-D development of previous 2-D results aimed at determining the magnetization direction, extending the calculation to second-order moments to recover the centre of mass of the magnetization distribution. The method is enhanced to reduce the effects of the regional field that often alters the first-order solutions. Moreover, we introduce an iterative correction to properly assess the errors coming from finite-size surveys or interaction with neighbouring anomalies, which are the most important causes of the failing of the method for real data. We test the method on some synthetic examples, and finally, we show the results obtained by analysing the aeromagnetic anomaly of the Monte Vulture volcano in Southern Italy.     

基于OLI遥感影像的错戳龙错盐湖水深反演研究

传统的水深测量方法是利用测量船上安装的测深设备和定位设备测出水域各点的水深。用遥感方法进行水深反演,基于OLI遥感影像对错戳龙错盐湖进行水深反演研究,确定OLI4波段是研究错戳龙错盐湖遥感水深反演的最佳波段,建立了此盐湖的水深遥感模型。经过对遥感数据和实测点数据的回归分析,建立了拟合方程。根据拟合方程对水深进行反演,反演结果与实测结果基本一致,从而证明了遥感测水深方法的快速有效性,弥补了传统实测水深方法的不足。     

Magnetotelluric interpretation of crustal and mantle structure in the Grenville Province

Summary Accurate determinations of depths and conductivities of electrical structures in shield regions are often difficult because of the inhomogeneity of the uppermost crust. A magnetotelluric (MT) station (BAT) in the Grenville Province of the Precambrian Shield in eastern Canada has been in operation since 1975 for time-dependency studies of electrical resistivity changes related to earthquakes. The MT response of the station displays low skew with small to moderate anisotropy. One-dimensional inversion of the apparent resistivity and phase reveals two well-defined conductors in the crust, one at 10 km and the second at the base of the crust. The latter has a resistivity less than 50 Ω m. These results are substantiated by three additional MT stations located up to 40 km distant.
Data from other new MT stations and from stations previously published in the literature are compared with two-dimensional computer model results and with the three-dimensional analogue scale model results of Dosso et al. While additional data for periods less than 100 s would be desirable the results from a number of the MT stations are not inconsistent with a widespread occurrence of a conducting zone at the base of the crust in the Grenville. The inversion analysis also indicates the existence of a conductor at some depth greater than 100 km with a resistivity less than 30 Ω m. This may coincide with a seismic low-velocity zone observed in the mantle under the Canadian Shield.     

Application of Magnetotelluric (MT) Resistivity to Imaging of Regional Three-Dimensional Geologic Structures and Groundwater Systems

Widespread definition of a groundwater system in three dimensions is necessary for the management and maintenance of groundwater resources. A magnetotelluric (MT) survey can be an effective geophysical prospecting method for imaging regional geological structures by measuring both shallow and deep resistivity. To demonstrate the capability of an MT survey to characterize a groundwater system, the Kumamoto area of central Kyushu in southwestern Japan was selected as a case study site because of its rich groundwater resources. Three-dimensional (3D) MT resistivity structure to a depth of 5?km was modeled by 1D inversion analysis of raw MT data and 3D interpolation of the resultant resistivity column data by the optimization principle method. Consequently, both deep and shallow aquifers were detected. A high-resistivity zone appears at depths between 500 and 2,000?m between the Futagawa?CHinagu faults and the Usuki?CYatsushiro tectonic line, which supports the existence of an aquiclude under the aquifer. The most important characteristic inferred from the 3D resistivity model is that the deep groundwater system below a depth of 1,000?m has two main flow paths. One path is likely to be through porous rocks because the low resistivity zone is regarded as tuff with sand and gravel, and the other flow path is interpreted to be through fractured zones along the Hinagu faults. Hence, the path and direction of the groundwater flows are probably controlled by geologic structures and the configuration of the active faults. These findings support the effectiveness of the MT method for investigating groundwater systems.     

Forward modelling of direct current and low-frequency electromagnetic fields using integral equations

We present a semi-analytical, unifying approach for modelling the electromagnetic response of 3-D bodies excited by low-frequency electric and magnetic sources. We write the electric and magnetic fields in terms of power series of angular frequency, and show that to obey Maxwell's equations, the fields must be real when the exponent is even, and imaginary when it is odd. This leads to the result that the scattering equations for direct current fields and for fields proportional to frequency can both be explicitly formulated using a single, real dyadic Green's function. Although the underground current flow in each case is due to different physical phenomena, the interaction of the scattering currents is of the same type in both cases. This implies that direct current resistivity, magnetometric resistivity and electric and magnetic measurements at low induction numbers can all be modelled in parallel using basically the same algorithm. We make a systematic derivation of the quantities required and show that for these cases they can all be expressed analytically. The problem is finally formulated as the solution of a system of linear equations. The matrix of the system is real and does not depend on the type of source or receiver. We present modelling results for different arrays and apply the algorithm to the interpretation of field data. We assume the standard dipoledipole resistivity array for the direct current case, and vertical and horizontal magnetic dipoles for induction measurements. In the case of magnetometric resistivity we introduce a moving array composed of an electric dipole and a directional magnetometer. The array has multiple separations for depth discrimination and can operate in two modes. The mode where the predominant current flow runs along the profile is called MMR-TM. This mode is more sensitive to lateral variations in resistivity than its counterpart, MMR-TE, where the mode of conduction is predominantly perpendicular to the profile.     

Effects of galvanic distortion on magnetotelluric data over a three-dimensional regional structure

Although the galvanic distortion due to local, near-surface inhomogeneities is frequency-independent, its effect on the magnetotelluric data becomes, in a 3-D structure, frequency-dependent. Therefore, both the apparent resistivity and the phase responses are disturbed, and a correction should be carried out prior to the 3-D interpretation in order to retrieve the 3-D regional impedance tensor. In many cases, the structure is 2-D for depths corresponding to a first range of periods and 3-D for longer periods (called 2-D/3-D). For these cases, a simple method which allows us to retrieve the 3-D regional impedance tensor (except the static shift) is presented. The method proposed uses the Groom & Bailey decomposition of the distortion matrix for the short periods. Three examples are presented: two using synthetic data and one employing real data. These examples show the effect of the galvanic distortion over a regional 2-D/3-D model and the retrieval of the regional transfer functions from the distorted ones.     

冰雷达探测数据处理方法研究——以南极冰盖Dome A地区的数据处理为例

提出了冰雷达数据处理的流程,详细论述了冰雷达数据处理中常规修正的关键技术：静态校正、增益控制、带通滤波以及偏移处理等。以中国第24次南极科学考察（CHINARE 24）所获取的Dome A 地区30 km×30 km范围冰雷达原始数据处理作为实例,提取得到了该区域冰盖内部等时层埋深以及冰盖的冰厚数据。通过插值展现了冰盖等时层以及冰岩界面的三维空间形态。结合该区域ICESat冰盖表面高程数据,构建出了Dome A地区涵盖冰盖表面、冰盖等时层以及冰岩界面形态特征的三维模型。