The changes in the electric conductivity of the lithosphere in the source region of the strongest Olyutora earthquake in the Koryak highlands

The results of the magnetotelluric (MT) soundings before and after an earthquake are analyzed. The interpretation is based on the longitudinal and transverse MT curves (along and across the strike of the main tectonic elements, respectively). The MT curves are distorted by the ρ- and coast effect. The distortions due to the coast effect are estimated by the testing three-dimensional (3D) model. It is established that the coast effect distortion at the periods up to 1000 s is small and can be disregarded. The divergence of the longitudinal and transverse MT curves, which points to the presence of the deep faults, is thoroughly studied. The inversion of the MTsounding curves is carried out by the REBOCC program of the numerical two-dimensional modeling. This program implements the procedures of elimination of the ρ-effect and the joint inversion of the longitudinal and transverse MT sounding curves. The obtained geoelectrical cross sections provide an insight into the structure of electrical conductivity of the lithosphere before and after the earthquake. The more intense variations in the electric conductivity are observed in the zone of the deep faults. These variations are related to the changes in the porosity and saturation of the rocks by the highly mineralized fluids.     

The geoelectrical properties of the earth around the Karymshina multidisciplinary observation site in southern Kamchatka

This paper is concerned with the Earth’s electromagnetic field, a deep geoelectrical section, and the dynamics of earth conductivity. This analysis is based on MTS curves along directions that to a first approximation go along and across southern Kamchatka. It is shown that the longitudinal and transverse MTS curves are subject to the influence of local and regional geoelectrical inhomogeneities. The coast effect was studied by 3-D numerical modeling. The patterns we have found were used for interpretation of generalized MTS curves. The resulting parameters of the geoelectrical section were refined by reducing the longitudinal curve to the standard curve of apparent resistivity. The results from this interpretation are to be refined as more geoelectrical information is forthcoming. The dynamics of lithospheric conductivity were studied from data acquired during the 2005–2008 monitoring of magnetotelluric impedance in the range of periods between several hundred and several thousand seconds. The basis for the analysis was assumed to be the transverse impedance and its phase. The latter was found to exhibit anomalous bay fluctuations that might have been related to earthquakes. A possible origin of these anomalies is discussed.     

The deep structure of the southern kamchatka volcanic zone from geophysical data

We discuss the interpretation method and results of magnetotelluric soundings in combination with other geological and geophysical data. The interpretation method was developed by studying possible distortions in MTS curves using 3D-numerical modeling of the magnetotelluric field. Deep conductivity was studied by using longitudinal MTS curves below a period of 400 s, which are nearly unaffected by the induction effect due to marine electrical currents. Transverse curves were used to obtain more detail for the geoelectric model. Inversion of average longitudinal MTS curves resulted in a geoelectric section of the lithosphere down to a depth of 60 km. Anomalies of high conductivity in the lithosphere were detected and were found to produce certain effects in gravity and seismic velocities. MTS and seismic tomography data were used to determine the possible origin of the high conductivity anomaly and to estimate rock porosity and the concentration of magma melts.     

Magnetotelluric sounding of Kamchatka

The MTS data acquired in Kamchatka during the last 30 years have been analyzed and summarized. Our interpretation is based on curves oriented along and across Kamchatka. Longitudinal and transverse curves can be affected by local geoelectric inhomogeneities. These were suppressed by conformal averaging. A bimodal interpretation of average longitudinal and transverse curves yielded a deep geoelectric model, which can be adopted as a starting point to be subsequently refined by 3D numerical modeling. The model involves a crustal conductive layer extending along central Kamchatka. In the east of the peninsula this layer is connected with crustal transverse conductive zones as wide as 50 km. Those zones have extensions toward the Pacific Ocean. Major centers of present-day volcanism occur in the transverse zones. The upper mantle contains an asthenospheric conductive layer forming an uplift beneath the present-day volcanic belt of Kamchatka.     

A deep geoelectric model of the Bol’she-Bannyi hydrothermal system,Kamchatka

This paper reports the results of magnetotelluric sounding (MTS) in the area where the Bol’she-Bannyi hydrothermal springs are discharged. The MTS curves were inverted on the assumption of a twodimensional inhomogeneous model using longitudinal and transverse curves of apparent resistivity. It was found that the geoelectric section contains a nearly vertical anomaly of high electrical conductivity at depths of 5.5–8 km, which is the signature of a deep-seated fault. The resulting geoelectric section for the upper crust and data from regional magnetotelluric soundings were used as a basis for developing a conceptual deep model of the Bol’she-Bannyi hydrothermal system. We quote an approximate estimate of rock porosity. According to the model, deep fluids come from a crustal layer into the subvertical deep-seated fault then penetrate via fissures into the sedimentary–volcanogenic cover, and finally arrive at the ground surface in zones of high rock permeability. We provide a recommendation for drilling a deep well in order to determine the potential of the Bol’she-Bannyi hydrothermal field.     

3D magnetotelluric modelling including surface topography

An edge finite‐element method has been applied to compute magnetotelluric (MT) responses to three‐dimensional (3D) earth topography. The finite‐element algorithm uses a single edge shape function at each edge of hexahedral elements, guaranteeing the continuity of the tangential electric field while conserving the continuity of magnetic flux at boundaries. We solve the resulting system of equations using the biconjugate gradient method with a Jacobian preconditioner. The solution gives electric fields parallel to the slope of a surface relief that is often encountered in MT surveys. The algorithm is successfully verified by comparison with other numerical solutions for a 3D‐2 model for comparison of modelling methods for EM induction and a ridge model. We use a 3D trapezoidal‐hill model to investigate 3D topographic effects, which are caused mainly by galvanic effects, not only in the Zxy mode but also in the Zyx mode. If a 3D topography were approximated by a two‐dimensional topography therefore errors occurring in the transverse electric mode would be more serious than those in the transverse magnetic mode.     

海岸效应对近海地区大地电磁测深数据畸变作用研究

在近海地区采集的大地电磁测深数据通常受到海岸效应的影响,使得大地电磁测深数据发生畸变,因而很难利用大地电磁测深资料较为可靠地获得地下深部的电性结构.本文通过正演模拟方法,分析和总结海水深度变化和海底地形变化对近海地区大地电磁测深数据的畸变影响.当测区与海岸线的距离小于目标频率的大地电磁场趋肤深度时,高导海洋的存在会严重影响测区内电磁场的分布.由于海岸效应的影响,大地电磁测深视电阻率曲线和相位曲线均会发生不同程度的畸变,在低频部分,这种畸变作用尤为明显.大地电磁测深一维Occam反演方法和二维非线性共轭梯度反演方法,对近海地区浅部地层具有较好的反演效果.随着海水深度的增加和海底地形的复杂变化,两种反演方法均会出现不同程度的假异常,为地质解释工作造成了影响.近渤海地区的实测大地电磁测深数据在低频部分可能受到海岸效应的影响而导致视电阻率曲线的严重畸变.     

Analysis of observed and model mt fields based on three-dimensional mathematical modeling: Case study of the Verkhnee Penzhino-Korf profile

We propose a system for the analysis of magnetotelluric (MT) data, which makes use of the invariant characteristics of the impedance tensor such as the maximum and minimum induction curves and the phase tensor. We examine the coefficients of the appearance and normalization of principal values of the impedance tensor. By the case study for Koryakiya, it is shown that the three-dimensional (3D) mathematical modeling and the Wiese-Parkinson vectors allow one to correct the results of one-dimensional (1D) and two-dimensional (2D) inversion of MT curves. Comparison between model and observed data based on the 1D inversion of MTS curves provides a pictorial view of the distortions of MT curves and their sensitivity to the parameters of a geological cross section.     

Coast effect distortion of marine magnetotelluric data: Insights from a pilot study offshore northeastern Japan

We report on strong coast effect distortions observed for broadband marine magnetotelluric (MT) data collected on the forearc offshore northeastern Japan. Eight days of horizontal electric and magnetic fields recorded at eight seafloor stations and the horizontal magnetic fields from a land remote station were processed with a robust multiple-station algorithm, yielding good MT responses and inter-station transfer functions at periods of 7–10,000 s. Transverse electric (TE) mode responses have cusps in apparent resistivity and negative phases at periods around 1000 s, while the transverse magnetic (TM) mode responses are galvanically depressed below the TE responses. An analysis of inter-station transfer functions confirms that the apparent resistivity cusps are a magnetic field, rather than electric field, phenomenon, consisting of an amplitude minimum and rapid phase change around a characteristic frequency. Poynting vectors for a TE coast effect model study illustrate that the anomalous phases are associated with energy diffusing back up to the seafloor from below, after being turned around from its usual downward propagating trajectory by inductive coupling between the conductive ocean and the resistive seafloor along the continental margin. We show that the characteristic frequency and position of the TE mode apparent resistivity cusps are determined by a relatively simple combination of the electrical resistivity of the seafloor, the depth of the ocean, and the distance from the coastline. By including coastlines and bathymetry in 2D inversion, we recover the seafloor conductivity structure along the forearc, demonstrating that broadband data can constrain the thickness of conductive forearc sediments and the underlying high resistivity associated with the mantle wedge and subducting oceanic lithosphere.     

The 3D profile inversion of the MTS–MVS data by the example of the Shiveluch volcanic zone in the Kamchatka Peninsula

The joint use of the magnetovariational (MV) and magnetotelluric (MT) data and the a priori information about the areal distribution of the deep electrical conductivity within the studied area and the integral conductance of the sedimentary cover enabled us to conduct a three-dimensional (3D) inversion of these data and construct the 3D model of the northeastern part of Kamchatka. The results show that within the Shiveluch volcano, north and south of it, there are northwesterly striking conductive faults penetrating up to the mantle which ascends from a depth of 41 to 35 km. The electrical conductivity of these crustal structures which accommodate high seismic and volcanic activity is probably accounted for by fluid saturation.     

基于非结构化网格的任意复杂2D RMT有限元模拟

射频大地电磁法(RMT)是以潜艇天线发射的射线源等作为场源的一种地球物理勘探方法,目前被广泛应用于近地表工程和环境地球物理勘探.RMT数据解释常采用基于静态假设的大地电磁法(MT)程序,往往会反演出不真实的浅层目标体.为解决这一长期困扰RMT资料解释的问题,本文实现了考虑位移电流的RMT有限元数值模拟.为了处理任意复杂模型,如起伏地形,非结构的三角形网格被用于离散RMT模型.首先通过算例对比,验证了程序的正确性和可靠性.通过Dike模型讨论了空气层厚度对RMT数值解的影响,结果表明当空气层厚度大于1/4波长即可满足精度要求.以山脊模型为例计算了位移电流对RMT响应的影响,表明位移电流的影响会随着山脊高程的增加而增大.最后通过舒家店实际模型进一步验证了位移电流在RMT中的重要性.     

3-D Global Induction in the Oceans and Solid Earth: Recent Progress in Modeling Magnetic and Electric Fields from Sources of Magnetospheric,Ionospheric and Oceanic Origin

Electromagnetic induction in the Earth’s interior is an important contributor to the near-Earth magnetic and electric fields. The oceans play a special role in this induction due to their relatively high conductivity which leads to large lateral variability in surface conductance. Electric currents that generate secondary fields are induced in the oceans by two different processes: (a) by time varying external magnetic fields, and (b) by the motion of the conducting ocean water through the Earth’s main magnetic field. Significant progress in accurate and detailed predictions of the electric and magnetic fields induced by these sources has been achieved during the last few years, via realistic three-dimensional (3-D) conductivity models of the oceans, crust and mantle along with realistic source models. In this review a summary is given of the results of recent 3-D modeling studies in which estimates are obtained for the magnetic and electric signals at both the ground and satellite altitudes induced by a variety of natural current sources. 3-D induction effects due to magnetospheric currents (magnetic storms), ionospheric currents (Sq, polar and equatorial electrojets), ocean tides, global ocean circulation and tsunami are considered. These modeling studies demonstrate that the 3-D induction (ocean) effect and motionally-induced signals from the oceans contribute significantly (in the range from a few to tens nanotesla) to the near-Earth magnetic field. A 3-D numerical solution based on an integral equation approach is shown to predict these induction effects with the accuracy and spatial detail required to explain observations both on the ground and at satellite altitudes. On leave from Institute of Terrestrial Magnetism, Ionosphere and Radiowave Propagation, Russian Academy of Sciences, 142190 Troitsk, Moscow region, Russia.     

3D geoelectric model of the Mutnov steam hydrothermal deposit

The article discusses the results of the interpretation of magnetotelluric sounding (MTS) data generated in the vicinity of the Mutnov steam hydrothermal deposit during the period from 2004 to 2006. Approaches to an interpretation of magnetotelluric (MT) data in a situation characterized by significant 3D distortions and the coastal effect phenomenon are discussed. The resulting 3D geoelectric model is presented. Based on analysis of the created geoelectric model, the principal geoelectric horizons are characterized, which includes the parameters of the integral conductivity and the morphology of the geoelectric boundaries. Zones of the highest values of the electrical conductivity of the watered tuffaceous sediment horizon have been identified, as well as zones of the maximum lateral electrical conductivity contrasts that are interpreted as discharge zones of hydrotherms. The results of the drilling carried out based on MTS data are presented.     

大地电磁局部畸变异常特征三维数值模拟研究

局部畸变问题曾经困扰大地电磁资料反演解释几十年,大地电磁三维数值模拟技术的发展为剖析局部畸变特点和得到可靠的反演成像结果提供了技术基础。本文采用三维数值模拟成像方法对典型三维局部畸变模型进行模拟分析。三维数值模拟结果显示:电场分量垂直电性分界面的极化模式视电阻率曲线(对应二维情况下TM模式)在穿越低阻异常体界面时,曲线会先上移后下移,而在穿越高阻异常体界面时,曲线会先下移后上移,这与电性分界面处积累面电荷产生的二次电场有关。三维模型中XY模式、YX模式视电阻率和相位在三维异常体附近的水平变化是呈现近似垂向对称的,该现象与电场垂直跨越电性界面时视电阻率的变化规律是吻合的,当测线分别沿X方向和Y方向展布时,三维情况下的XY和YX模式分别对应二维情况下的TM模式。低阻小异常体对区域构造响应的畸变影响比高阻小异常体要严重。低阻小异常体对二维区域响应的两种极化模式视电阻率和相位都有非常明显的畸变影响,相比较而言对TE模式的畸变要大于TM模式,因此我们在做二维反演解释时,可优先考虑拟合TM模式数据。位于小异常体中心上方测点的三维畸变响应虽然与对应真实二维区域响应的差异比较大,但可以等效于某种二维模型响应,这种由局部畸变造成的假二维响应在实际野外数据的解释中是需要注意的。      

The comparison of the experimental MT data with the results of numerical modelling for the Kamchatka Peninsula

The experimental MT data and theoretical data calculated by the thin sheet technique for a low-frequency telluric field are analysed and compared for the Kamchatka Peninsula. The comparison shows a good agreement between the experimental and theoretical data and results in a conclusion that the DC-type distortions of the telluric field play the most important role in the structure of the MT field. The component of the telluric field which is parallel to the general direction of the Kamchatka Peninsula is less distorted and should be used for the estimation of deep geoelectric structures.     

Electrical conductivity anomalies in the earth

Anomalies of short-period geomagnetic variations have been found in various regions over the world. It is known that such anomalies arise from electromagnetic induction within an electrical conductivity anomaly or from local perturbation of induced electric currents by a conductivity anomaly. In order to investigate a regional electric state in the Earth, conductivity anomaly (CA) studies based on anomalous behaviors of geomagnetic variations have been extensively undertaken, as well as studies based on magnetotelluries in which induced currents are directly used.Some of the geomagnetic variation anomalies, however, turned out to be caused by surface conductors, such as sea water and sediments. Anomalies of this sort have been intensively studied and classified into coast, island, peninsula, and strait effects in the case of sea effects. Three-dimensional conduction or channelling of induced electric currents is sometimes observed in the cases of sediments and some crustal conductivity anomalies. However, anomalies of such surface origins often provide some information of the underground conductivity structure.Electrical conductivity anomalies can be classified into two types: anomalies originating in the crust and in the upper mantle. Many of crustal anomalies are well correlated with metamorphic belts, fracture zones, and hydrated layers, and magnetic and gravity anomalies are also often found over the conductivity anomalies. Most of mantle anomalies have been interpreted mainly in terms of high temperature and partial melting, since conductivity anomalies coincide well with anomalies in heat flow and seismic wave velocities.     

大地电导率横向突变处磁暴感应地电场的邻近效应

大地电性结构的横向变化会对磁暴时的感应地电流和地面电磁场产生影响.本文假设扰动地磁场变化的源为地面以上一定高度的面电流,以某一典型层状大地电导率结构为基础,构造含有电导率横向突变的地电模型.针对感应电流的方向与横向分界面平行的情形,采用伽辽金有限元法对电导率横向突变处的感应地电场进行了分析,揭示了电导率横向差异产生的趋肤效应和邻近效应的机理,针对与电性结构分界面平行的输电线路,从评估地磁感应电流的角度讨论了影响的严重程度和范围.     

Transient electromagnetic surveys with unimodal transverse magnetic field: ideas and results

The theory behind transient electromagnetic surveys can be well described in terms of transverse magnetic and transverse electric modes. Soundings using transverse magnetic and transverse electric modes require different source configurations. In this study, we consider an alternating transverse magnetic field excitation by a circular electric dipole. The circular electric dipole transmitter is a horizontal analogue of the vertical electric dipole. Offshore surveys using circular electric dipole might represent an alternative to the conventional marine controlled‐source electromagnetic method at shallow sea and/or for exploring relatively small targets. Field acquisition is carried out by recording either electric or magnetic responses. Electric responses bear information on the 1D structure of a layered earth and successfully resolve high‐resistivity targets in marine surveys. Land‐based circular electric dipole soundings are affected by induced polarisation. On the contrary, magnetic responses are absent on the surface of a 1D earth, and as a result, they are very sensitive to any and even very small 3D conductivity perturbations. In addition, they are sensitive to induced polarisation or some other polarisation effects in the subsurface. At present, circular electric dipole transmitters and magnetic receivers are successfully used in on‐land mineral and petroleum exploration.     

郯庐断裂带中段(沂沭断裂带)电性结构研究与孕震环境

穿过郯庐断裂带中段（沂沭断裂带,36°N）所做的大地电磁测深（MT）剖面长约150km.使用Robust技术和远参考道大地电磁方法处理观测数据.通过分析视电阻率、阻抗相位、Swift二维偏离度和区域走向,定性确定测区的电性结构.二维反演解释中选择非线性共轭梯度（NLCG）方法,使用TE、TM两种模式资料联合反演,沿剖面的二维电性结构显示：自西向东,鲁西隆起、郯庐断裂带、胶莱坳陷及鲁东隆起4个电性区块分别对应,鲁东和鲁西隆起区为高阻,郯庐断裂带电性结构复杂,高、低阻相间,胶莱坳陷为低阻（高导）区.沿MT剖面附近曾发生3个地震,其震源区处在电性变化剧烈部位,并在震源区附近存在高导体.