Bootstrapped discrete scale invariance analysis of geomagnetic dipole intensity

The technique of bootstrapped discrete scale invariance allows multiple time-series of different observables to be normalized in terms of observed and predicted characteristic timescales. A case study is presented using the SINT2000 time-series of virtual axial dipole moment, which spans the past 2 Myr. It is shown that this sequence not only bears a clear signature of a preferred timescale of about 55.6 Ka, but additionally predicts similar features (of shorter and longer duration) that are actually observed on the timescales of historical secular variation and dipole reversals, respectively. In turn, the latter two empirical sources both predict the characteristic timescale found in the dipole intensity sequence. These communal scaling characteristics suggest that a single underlying process could be driving dynamo fluctuations across all three observed timescales, from years to millions of years.     

Application of the normalized surface magnetic charge model to UXO discrimination in cases with overlapping signals

This paper presents an application of the normalized surface magnetic charge (NSMC) model to discriminate objects of interest, such as unexploded ordnance (UXO), from innocuous items in cases when UXO electromagnetic induction (EMI) responses are contaminated by signals from other objects. Over the entire EMI spectrum considered here (tens of Hertz up to several hundreds of kHz), the scattered magnetic field outside the object can be produced mathematically by equivalent magnetic charges. The amplitudes of these charges are determined from measurement data and normalized by the excitation field. The model takes into account the scatterer's heterogeneity and near- and far-field effects. For classification algorithms, the frequency spectrum of the total NSMC is proposed and investigated as a discriminant. The NSMC is combined with the differential evolution (DE) algorithm in a two-step inversion procedure. To illustrate the applicability of the DE–NSMC algorithm, blind test data are processed and analyzed for cases in which signals from nearby objects frequently overlap. The method was highly successful in distinguishing UXO from accompanying clutter.     

2D magnetic resonance tomography applied to karstic conduit imaging

Karstic conduits play a crucial role for water supply in many parts of the world. However, the imaging of such targets is generally a difficult task for most geophysical methods. Magnetic Resonance Sounding (MRS) is a geophysical method designed for imaging of water bearing structures. Initially, MRS was developed for characterizing horizontally stratified aquifers. However, when applying a 1D MRS measuring setup to the imaging of 2D–3D targets, the size of which may be much smaller than the loop, the accuracy and the lateral resolution may not be sufficient. We have studied the possibility of simultaneously processing several MRS aligned along a profile to perform a Magnetic Resonance Tomography (MRT). This work emphasizes the gain of resolution for 2D–3D imagery of MRT versus the interpolation of 1D inversion results of MRS along the same profile. Numerical modelling results show that the MRT response is sensitive to the size and location of the 2D target in the subsurface. Sensitivity studies reveal that by using the coincident transmitting/receiving (TX/RX) setup and shifting the loop around the anomaly area, the depth, section and position of a single karstic conduit with a size smaller than the MRS loop size can be resolved. The accuracy of the results depends on the noise level and signal level, the latter parameter being linked to the depth and volume of the karstic conduit and the water content in the limestone matrix. It was shown that when applying MRT to the localization of 2D anomalies such as karstic conduits, the inclination of the geomagnetic field, the orientation of the MRT profile and the angle of crossover of the conduit by the MRT profile must be taken into account. Otherwise additional errors in interpretation should be expected. A 2D inversion scheme was developed and tested. Both numerical and experimental results confirm the efficiency of the developed approach.     

王艳，廖明生，李德仁等.利用长时间序列相干目标获取地面沉降场.地球物理学报，2007

方法，相干目标分析是利用时间序列多景影像获取地表稳定反射体——高相干点目标的新方法.实验区选择上海市主城区约100 km2 区域，InSAR的实验数据采用覆盖1992～2000年间的25景ERS-1/2的单视复影像.在无先验沉降场模型的情况下，实验结果显示该实验区的地面沉降场在时间范围上跨越8年，平均观测时间采样率约为4个月，垂直形变精度优于5 mm，满足我国对城市地面沉降形变观测的要求.研究结果表明该方法可以用于今后大时间/空间尺度上的高精度地表形变观测.     

Automatic Interpretation of Magnetic Data Using Euler Deconvolution with Nonlinear Background

The voluminous gravity and magnetic data sets demand automatic interpretation techniques like Naudy, Euler and Werner deconvolution. Of these techniques, the Euler deconvolution has become a popular choice because the method assumes no particular geological model. However, the conventional approach to solving Euler equation requires tentative values of the structural index preventing it from being fully automatic and assumes a constant background that can be easily violated if the singular points are close to each other. We propose a possible solution to these problems by simultaneously estimating the source location, depth and structural index assuming nonlinear background. The Euler equation is solved in a nonlinear fashion using the optimization technique like conjugate gradient. This technique is applied to a published synthetic data set where the magnetic anomalies were modeled for a complex assemblage of simple magnetic bodies. The results for close by singular points are superior to those obtained by assuming linear background. We also applied the technique to a magnetic data set collected along the western continental margin of India. The results are in agreement with the regional magnetic interpretation and the bathymetric expressions.     

Effect of surface stress on magneto-elastic surface waves in finitely conducting media

This paper deals with the investigation of the effect of surface stress and conductivity on the propagation of surface wave in isotropic, homogeneous, elastic media under the action of a primary magnetic field. Formulation of the general surface wave propagation problem has been made, and the corresponding frequency equation has been derived. Frequency equations for Rayleigh wave, surface shear wave and Stoneley wave have been deduced from that of general surface wave as special cases. The effects of surface stress and magnetic field on the wave velocities and attenuation factors of Rayleigh wave and surface shear wave are shown by numerical calculations and graphs. Some important wave velocity equations, as obtained by other authors, have been deduced as special cases from the wave velocity equation for Stoneley wave. It is found that the combined effect of surface stress and magnetic field modulates the wave velocity ratios and attenuation factors of Rayleigh wave and surface shear wave to a considerable extent.     

Use of variograms for field magnetometry analysis in Upper Silesia Industrial Region

Measurements of topsoil magnetic susceptibility are often used for quick assessment of soil contamination of anthropogenic origin, with heavy metals or other pollutants. However, because of complicated correlations between low-field magnetic susceptibility (shortened to magnetic susceptibility) of topsoil and soil pollution, the outcome of a field magnetometry survey can not be related directly to soil pollution. For each case study, the results should be interpreted on their own taking into account not only the type of pollution but also pedogenic, biogenic and environmental factors. In practice, it is very difficult to measure and consider all these factors. Here we illustrate the merit of geostatistical methods, which are focused on the spatial variability of a phenomenon, in the interpretation of soil magnetometry results. This article presents the analysis of spatial variability of top soil layers magnetic susceptibility-within the Upper Silesia Industrial Region (USIR)-using semivariance analysis. It also explains how to adjust the sampling density of field magnetometry measurements to spatial variability of the soil pollution as well as to the spatial scale of the investigated area. For this purpose, the values of magnetic susceptibility have been measured by using various sampling densities at areas of different size located within USIR. This enabled to determine the main scales of magnetic susceptibility spatial variability of soils within USIR using semivariance. A few distinct scales of variability were found from the site scale to a more regional scale. Variability ranges of 30 km, 12 km, and 5 km refer to the large regional scale, whereas smaller ranges of few hundreds down to a few tens of meters, can be attributed to the local (site) scale. In addition, the precision of the measuring campaigns, performed within USIR with different sampling densities, was compared through the analysis of the spatial variability of the soil magnetic susceptibility signal by using ordinary kriging. jarek97@yahoo.com, piotr.fabijanczyk@is.pw.edu.pl     

Two-dimensional power series solution for non-axisymmetrical transport in a radially convergent tracer test with scale-dependent dispersion

It has been known for many years that dispersivity increases with solute travel distance in a subsurface environment. The increase of dispersivity with solute travel distance results from the significant variation of hydraulic properties of heterogeneous media and was identified in the literature as scale-dependent dispersion. This study presents an analytical solution for describing two-dimensional non-axisymmetrical solute transport in a radially convergent flow tracer test with scale-dependent dispersion. The power series technique coupling with the Laplace and finite Fourier cosine transform has been applied to yield the analytical solution to the two-dimensional, scale-dependent advection–dispersion equation in cylindrical coordinates with variable-dependent coefficients. Comparison between the breakthrough curves of the power series solution and the numerical solutions shows excellent agreement at different observation points and for various ranges of scale-related transport parameters of interest. The developed power series solution facilitates fast prediction of the breakthrough curves at any observation point.