Time-domain electromagnetic migration in the solution of inverse problems

Time-domain electromagnetic (TDEM) migration is based on downward extrapolation of the observed field in reverse time. In fact, the migrated EM field is the solution of the boundary-value problem for the adjoint Maxwell's equations. The important question is how this imaging technique can be related to the solution of the geoelectrical inverse problem. In this paper we introduce a new formulation of the inverse problem, based on the minimization of the residual-field energy flow through the surface or profile of observations. We demonstrate that TDEM migration can be interpreted as the first step in the solution of this specially formulated TDEM inverse problem. However, in many practical situations this first step produces a very efficient approximation to the geoelectrical model, which makes electromagnetic migration so attractive for practical applications. We demonstrate the effectiveness of this approach in inverting synthetic and practical TDEM data.     

Natural ring structures on the Baikal ice cover: Analysis of experimental data and mathematical modeling

Hydrophysical studies and mathematical modeling of ring structures during ice cover on Lake Baikal have shown that their existence at the stage of ice cover degradation is due to anticyclonic currents. Such currents can be generated as a result of local upwelling, which we associate with the rise of methane hydrates from the top layer of bottom sediments and their dissociation. Analysis of satellite images shows that the radii of ice rings range from 1300 to 2400 m, which is close to the baroclinic Rossby radius. The measured ice thicknesses in the area of the rings are in agreement with model calculations. Deep water renewal in Lake Baikal can also be associated with the rise of hydrates.     

Variations of the X-ray INTEGRAL spectrum of the active galactic nucleus of NGC 4945

Results of 11-year-long X-ray INTEGRAL observations of the nucleus of Seyfert galaxy NGC 4945 in the 3–500 keV range were processed. A two-component spectrum model, which includes strong radiation absorption in the Compton-thick torus around the AGN “central engine” and secondary radiation reflected from the torus walls, was used in the analysis. The following primary spectrum parameters were determined based on the data accumulated throughout the entire exposure period: photon index Γ = 1.60 ± 0.07, exponential cutoff energy E _c =157 _-22 ⁺²⁹ keV, and column density of the medium that absorbs primary radiation N _H,1 =5.0 _-0.9 ^+1.0 × 10²⁴ cm^–2. The column density of the medium absorbing reflected radiation is two orders of magnitude lower. Both the X-ray flux in the ranges of 20–40, 40–60, and 60–100 keV and the shape of the X-ray spectrum of NGC 4945 vary. The spectrum shape variations may be induced by inhomogeneities of the absorbing medium surrounding the AGN. At the same time, there is some evidence for moderate spectrum variations in the highenergy region, which may be associated with changes in the “central engine.”     

Seasonal variations in the vertical structure of pelagian water stratum in the Southern Baikal

For the first time, T,S-analysis was used to determine the specifics of seasonal variations in the vertical structure of Lake Baikal active layer. In the under-ice period, the active layer includes the under-ice, top winter, and upper intermediate water masses. The under-ice water mass, unlike other masses, shows an increase in mineralization to 100.74 mg/kg, which corresponds to a release of 71.1 g salt under 1 m² of water surface in a layer 0–40 m in the process of salt freezing out during ice cover formation and accretion. In the phases of mixing (homothermy), the water masses of the active layer transform into a surface homogeneous mass. In summer and autumn, surface and upper intermediate water masses, separated by a water mass of summer thermocline can be identified. A specific feature of the summer thermocline water mass is the increased sum of ions because of an increase in HCO ₃ ^- concentration at the decay of organic matter accumulating in the bottom part of the thermocline. The existence of the under-ice water mass and the water mass of summer thermocline was established in Lake Baikal for the first time. In the deep-water zone (>250 m), except for the bottom parts, the lower water masses (the lower intermediate and the deep) are permanent, their characteristics remaining stable during the year. The changes in the bottom water mass are due to the character of the processes of bottom water renewal.     

Restructuring of the Solar Magnetic Fields and Flare Activity Centers in Cycle 24

The development of the solar magnetic activity in cycle 24 has been analyzed. It has been shown that the significant north–south asymmetry of magnetic activity was accompanied by the asynchronous reorganization of solar magnetic fields in the northern and southern hemispheres. The formation of unipolar magnetic regions after the decay of activity centers has been studied. The meridional transport of unipolar magnetic regions leading to changes in the zonal structure of the solar magnetic field has been shown. Long-lived centers of flare activity have been found to exist during the periods of magnetic field restructuring. The spatiotemporal analysis of the flare ensemble making it possible to diagnose non-stationary processes in the solar atmosphere has been shown.     

Incorporating known petrophysical model in the seismic full-waveform inversion using the Gramian constraint

In this paper, we develop a general approach to integrating petrophysical models in three-dimensional seismic full-waveform inversion based on the Gramian constraints. In the framework of this approach, we present an example of the frequency-domain P-wave velocity inversion guided by an electrical conductivity model. In order to introduce a coupling between the two models, we minimize the corresponding Gramian functional, which is included in the Tikhonov parametric functional. We demonstrate that in the case of a single-physics inversion guided by a model of different physical type, the general expressions of the Gramian functional and its gradients become simple and easy to program. We also prove that the Gramian functional has a non-negative quadratic form, so it can be easily incorporated in a standard gradient-based minimization scheme. The developed new approach of seismic inversion guided by the known petrophysical model has been validated by three-dimensional inversion of synthetic seismic data generated for a realistic three-dimensional model of the subsurface.     

电导率各向异性的海洋电磁三维有限单元法正演

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.     

Large-scale three-dimensional inversion of EarthScope MT data using the integral equation method

In this paper we apply 3D inversion to MT data collected in the Northwestern United States as a part of the EarthScope project. By the end of 2009 MT data had been collected from 262 stations located throughout Oregon, Washington, Idaho, and most of Montana and Wyoming. We used data from 139 MT stations in this analysis. We developed fully parallelized rigorous 3D MT inversion software based on the integral equation method with variable background conductivity. We also implemented a receiver footprint approach which considerably reduced the computational resources needed to invert the large volumes of data covering vast areas. The data set used in the inversion was obtained through the Incorporated Research Institutions for Seismology (IRIS). The inversion domain was divided into 2.7 M cells. The inverted electrical conductivity distribution agrees reasonably well with geological features of the region.