Magnetotelluric strike rules

Geophysics ◽  
1987 ◽  
Vol 52 (3) ◽  
pp. 267-278 ◽  
Author(s):  
P. Zhang ◽  
R. G. Roberts ◽  
L. B. Pedersen

The distortion of the magnetotelluric impedance tensor by complex “near‐surface” structure leads to leakage between the elements of the tensor. The magnetotelluric impedance tensor for our principal model, which has both a local and a regional strike, can be written in the long‐period limit as a sum of the regional, undistorted impedance and a perturbed impedance. The latter can be written as a product of a local distortion (which can be regarded as thin‐sheet distortion in the long‐period range) tensor and the regional impedance. Local and regional strikes are found by rotating the impedance tensor into directions in which diagonal elements are proportional and column elements are proportional, respectively. The regional impedance tensor is calculated assuming that the strikes are uniquely defined. An example from a crystalline area with well conducting fracture zones illustrates the model concepts. A weighted least‐squares procedure is used for the estimation of distortion parameters.

Geophysics ◽  
1990 ◽  
Vol 55 (3) ◽  
pp. 372-375 ◽  
Author(s):  
Philip E. Wannamaker

Modeling of three‐dimensional (3-D) thin‐sheet telluric anomalies has been a popular means of estimating the distortions of magnetotelluric (MT) response functions in the vicinity of an upper crustal resistivity inhomogeneity (Kaikkonen, 1986) because structures with an arbitrary 3-D shape in plan view can be simulated with modest computational effort. In the class of thin‐layer problems of this note, the sheet structure overlies an infinitely resistive basement; and anomalous fields are computed in the long‐period or zero‐frequency (dc) limit (e.g., Hermance, 1982). Under the assumption of no vertical current flow into the basement, this 3-D problem reduces to a two‐dimensional (2-D) dc solution across the sheet. Since only the long‐period limit is considered, however, the conclusions that can be drawn about the effects of general 3-D structures are restricted. The purpose of this note is to clarify restrictions on the use of 2-D MT interpretations in 3-D areas based on just thin‐sheet telluric modeling.


Geophysics ◽  
2012 ◽  
Vol 77 (4) ◽  
pp. E301-E315 ◽  
Author(s):  
Thomas Kalscheuer ◽  
Juliane Hübert ◽  
Alexey Kuvshinov ◽  
Tobias Lochbühler ◽  
Laust B. Pedersen

Magnetotelluric (MT), radiomagnetotelluric (RMT), and, in particular, controlled-source audiomagnetotelluric (CSAMT) data are often heavily distorted by near-surface inhomogeneities. We developed a novel scheme to invert MT, RMT, and CSAMT data in the form of scalar or tensorial impedances and vertical magnetic transfer functions simultaneously for layer resistivities and electric and magnetic galvanic distortion parameters. The inversion scheme uses smoothness constraints to regularize layer resistivities and either Marquardt-Levenberg damping or the minimum-solution length criterion to regularize distortion parameters. A depth of investigation range is estimated by comparing layered model sections derived from first- and second-order smoothness constraints. Synthetic examples demonstrate that earth models are reconstructed properly for distorted and undistorted tensorial CSAMT data. In the inversion of scalar CSAMT data, such as the determinant impedance or individual tensor elements, the reduced number of transfer functions inevitably leads to increased ambiguity for distortion parameters. As a consequence of this ambiguity for scalar data, distortion parameters often grow over the iterations to unrealistic absolute values when regularized with the Marquardt-Levenberg scheme. Essentially, compensating relationships between terms containing electric and/or magnetic distortion are used in this growth. In a regularization with the minimum solution length criterion, the distortion parameters converge into a stable configuration after several iterations and attain reasonable values. The inversion algorithm was applied to a CSAMT field data set collected along a profile over a tunnel construction site at Hallandsåsen, Sweden. To avoid erroneous inverse models from strong anthropogenic effects on the data, two scalar transfer functions (one scalar impedance and one scalar vertical magnetic transfer function) were selected for inversion. Compared with a regularization of distortion parameters with the Marquardt-Levenberg method, the minimum-solution length criterion yielded smaller absolute values of distortion parameters and a horizontally more homogeneous distribution of electrical conductivity.


2020 ◽  
Vol 13 (7) ◽  
pp. 3521-3542
Author(s):  
Ian Ashpole ◽  
Aldona Wiacek

Abstract. We compare MOPITT Version 7 (V7) Level 2 (L2) and Level 3 (L3) carbon monoxide (CO) products for the 1∘×1∘ L3 grid box containing the coastal city of Halifax, Canada (longitude −63.58∘, latitude 44.65∘), with a focus on the seasons DJF and JJA, and highlight a limitation in the L3 products that has significant consequences for the temporal trends in near-surface CO identified using those data. Because this grid box straddles the coastline, the MOPITT L3 products are created from the finer spatial resolution L2 products that are retrieved over both land and water, with a greater contribution from retrievals over water because more of the grid box lies over water than land. We create alternative L3 products for this grid box by separately averaging the bounded L2 retrievals over land (L3L) and water (L3W) and demonstrate that profile and total column CO (TCO) concentrations, retrieved at the same time, differ depending on whether the retrieval took place over land or water. These differences (ΔRET) are greatest in the lower troposphere (LT), where mean retrieved volume mixing ratios (VMRs) are greater in L3W than L3L, with maximum mean differences of 11.6 % (14.3 ppbv, p=0.001) at the surface level in JJA. Retrieved CO concentrations are more similar, on average, in the middle and upper troposphere (MT and UT), although large differences (in excess of 40 %) do infrequently occur. TCO is also greater in L3W than L3L in both seasons. By analysing L3L and L3W retrieval averaging kernels and simulations of these retrievals, we demonstrate that, in JJA, ΔRET is strongly influenced by differences in retrieval sensitivity over land and water, especially close to the surface where L3L has significantly greater information content than L3W. In DJF, land–water differences in retrieval sensitivity are much less pronounced and appear to have less of an impact on ΔRET, which analysis of wind directions suggests is more likely to reflect differences in true profile concentrations (i.e. real differences). The original L3 time series for the 1∘×1∘ grid box containing Halifax (L3O) corresponds much more closely to L3W than L3L, owing to the greater contribution from L2 retrievals over water than land. Thus, in JJA, variability in retrieved CO concentrations close to the surface in L3O is suppressed compared to L3L, and a declining trend detected using weighted least squares (WLS) regression analysis is significantly slower in L3O (strongest surface level trend identifiable is −1.35 (±0.35) ppbv yr−1) than L3L (−2.85 (±0.60) ppbv yr−1). This is because contributing L2 retrievals over water are closely tied to a priori CO concentrations used in the retrieval, owing to their lack of near-surface sensitivity in JJA, and these are based on monthly climatological CO profiles from a chemical transport model and therefore have no yearly change (surface level trend in L3W is −0.60 (±0.33) ppbv yr−1). Although our analysis focuses on DJF and JJA, we demonstrate that the findings also apply to MAM and SON. The results that we report here suggest that similar analyses be performed for other coastal cities before using MOPITT surface CO.


Geophysics ◽  
1992 ◽  
Vol 57 (9) ◽  
pp. 1127-1137 ◽  
Author(s):  
Andreas Hördt ◽  
Vladimir L. Druskin ◽  
Leonid A. Knizhnerman ◽  
Kurt‐Martin Strack

The interpretation of long‐offset transient electromagnetic (LOTEM) data is usually based on layered earth models. Effects of lateral conductivity variations are commonly explained qualitatively, because three‐dimensional (3-D) numerical modeling is not readily available for complex geology. One of the first quantitative 3-D interpretations of LOTEM data is carried out using measurements from the Münsterland basin in northern Germany. In this survey area, four data sets show effects of lateral variations including a sign reversal in the measured voltage curve at one site. This sign reversal is a clear indicator of two‐dimensional (2-D) or 3-D conductivity structure, and can be caused by current channeling in a near‐surface conductive body. Our interpretation strategy involves three different 3-D forward modeling programs. A thin‐sheet integral equation modeling routine used with inversion gives a first guess about the location and strike of the anomaly. A volume integral equation program allows models that may be considered possible geological explanations for the conductivity anomaly. A new finite‐difference algorithm permits modeling of much more complex conductivity structures for simulating a realistic geological situation. The final model has the zone of anomalous conductivity aligned below a creek system at the surface. Since the creeks flow along weak zones in this area, the interpretation seems geologically reasonable. The interpreted model also yields a good fit to the data.


Geophysics ◽  
1978 ◽  
Vol 43 (3) ◽  
pp. 563-575 ◽  
Author(s):  
H. O. Seigel ◽  
D. H. Pitcher

The Tridem vertical coplanar airborne electromagnetic system provides simultaneous in‐phase and quadrature information at frequencies of 500, 2000 and 8000 Hz. The system can map a broad range of earth conductors of simple geometry and provide quantitative estimates of their conductivities and dimensions. Computer programs have been developed to automatically interpret the six channels of Tridem data, plus the output of an accurate radar altimeter, to determine the depth of burial, conductivity and thickness of a near‐surface, flat‐lying conducting horizon. In limiting cases, the interpretation provides the conductance (conductivity‐thickness product) of a thin sheet (ranging from 100 mmhos to 100 mhos) or the conductivity of a homogeneous earth (ranging from 1 mmhos/m to 10 mhos/m). Two actual field examples are presented from Ontario, Canada; one relating to the mapping of overburden conditions (sand, clay and rock, etc) and the other to the mapping of the distribution of a buried lignite deposit. Other areas of potential application of the system to surficial materials would include groundwater mapping, permafrost investigations, and civil engineering studies for roads and pipelines.


2019 ◽  
Vol 218 (2) ◽  
pp. 1008-1031 ◽  
Author(s):  
Shunguo Wang ◽  
Mehrdad Bastani ◽  
Steven Constable ◽  
Thomas Kalscheuer ◽  
Alireza Malehmir

SUMMARY Boat-towed radio-magnetotelluric (RMT) measurements using signals between 14 and 250 kHz have attracted increasing attention in the near-surface applications for shallow water and archipelago areas. A few large-scale underground infrastructure projects, such as the Stockholm bypass in Sweden, are planned to pass underneath such water zones. However, in cases with high water salinity, RMT signals have a penetration depth of a few metres and do not reach the geological structures of interest in the underlying sediments and bedrock. To overcome this problem, controlled source signals at lower frequencies of 1.25 to 12.5 kHz can be utilized to improve the penetration depth and to enhance the resolution for modelling deeper underwater structures. Joint utilization of boat-towed RMT and controlled source audio-magnetotellurics (CSAMT) was tested for the first time at the Äspö Hard Rock Laboratory (HRL) site in south-eastern Sweden to demonstrate acquisition efficiency and improved resolution to model fracture zones along a 600-m long profile. Pronounced galvanic distortion effects observed in 1-D inversion models of the CSAMT data as well as the predominantly 2-D geological structures at this site motivated usage of 2-D inversion. Two standard academic inversion codes, EMILIA and MARE2DEM, were used to invert the RMT and CSAMT data. EMILIA, an object-oriented Gauss–Newton inversion code with modules for 2-D finite difference and 1-D semi-analytical solutions, was used to invert the RMT and CSAMT data separately and jointly under the plane-wave approximation for 2-D models. MARE2DEM, a Gauss–Newton inversion code for controlled source electromagnetic 2.5-D finite element solution, was modified to allow for inversions of RMT and CSAMT data accounting for source effects. Results of EMILIA and MARE2DEM reveal the previously known fracture zones in the models. The 2-D joint inversions of RMT and CSAMT data carried out with EMILIA and MARE2DEM show clear improvement compared with 2-D single inversions, especially in imaging uncertain fracture zones analysed in a previous study. Our results show that boat-towed RMT and CSAMT data acquisition systems can be utilized for detailed 2-D or 3-D surveys to characterize near-surface structures underneath shallow water areas. Potential future applications may include geo-engineering, geohazard investigations and mineral exploration.


Geophysics ◽  
1975 ◽  
Vol 40 (5) ◽  
pp. 790-797 ◽  
Author(s):  
Larry R. Grillot

A method of analyzing the tensor nature of magnetotelluric fields from short independent sections of filtered MT data is presented in which: (1) horizontal field components are filtered at chosen periods; (2) short sections which show a well defined signal over 3–4 cycles are selected from the records; (3) the polarizations of these short record sections are plotted; and (4) Fourier coefficients calculated for section pairs showing differing magnetic field polarizations are then used to calculate the elements of the impedance tensor. Results of this analysis performed on MT field data show that the use of MT signal pairs yields uniform principal resistivities for calculations using several different section pairs of the same period. These results also illustrate internal consistency of the resistivity estimates over the complete period range of the data (30–7500 sec). A brief comparison of the MT pairs calculations and results derived from a least‐squares approach at a selected period shows good agreement between the two methods. The overall results suggest that the use of independent signal pairs in the determination of the tensor impedance could complement more traditional methods in some cases.


1980 ◽  
Vol 88 ◽  
pp. 383-387 ◽  
Author(s):  
Douglas S. Hall ◽  
Jerzy M. Kreiner ◽  
Steven N. Shore

We have collected all available times of minimum, over 1650, for 23 regular, 7 short-period, and 4 long-period eclipsing RS CVn binaries using the card catalogues at the Sternberg Astronomical Institute and the Krakow Astronomical Observatory, as well as unpublished data. We examined them critically, discarded those grossly in error, and assigned weights. For AD Cap, RV Lib, and ε UMi there are virtually no data. For the remaining 31 we determined quadratic ephemerides by weighted least squares. Values of d1nP/dt are given in the Table with their errors. Roughly 2/3 are variable at greater than the 2 σ level. Significant decreases outnumber significant increases by almost 2:1. We show that light curve asymmetry produced by the distortion wave probably accounts for the relatively small, rapid fluctuations in the O-C curve but not the long-term period changes.


Geophysics ◽  
2006 ◽  
Vol 71 (2) ◽  
pp. G27-G34 ◽  
Author(s):  
Yusen Ley-Cooper ◽  
James Macnae ◽  
Terry Robb ◽  
Julian Vrbancich

We investigate the properties of EM signals in several different response-parameter domains to identify calibration errors in helicopter electromagnetic (HEM) data. In particular, we define a dimensionless response parameter α, derived from frequency-domain data, that is numerically identical to the historic wire-loop response parameter, and is closely related to the thin-sheet and half-space response parameters. The arctangent of α is the phase ϕ of the secondary field. We further define a dimensionless amplitude response parameter β, calculated as the ratio between inductive limits estimated from the data and from system geometry. The inductive limit calculated from geometry provides an initial altitude correction to the data amplitudes. Additional data corrections further correct phase effects and altimeter variations. The amplitude and phase errors in calibration become independent differences between the data and the fitted model in the ϕβ domain. This investigation was undertaken in the response-parameter domain rather than in the data domain, allowing the analysis to be independent of absolute values of conductivity and system frequencies. Statistical analysis in the ϕβ domain of recently acquired data suggests that amplitude calibration errors in HEM data may cause fitted models to have systematic depth errors of 1 to 2 m for near-surface conductors; variable altitude measurement errors are about 1.5 m over seawater.


Sign in / Sign up

Export Citation Format

Share Document