A high resolution lithospheric magnetic field model over China

CHAMP satellite data and ground-based magnetic observations are used and combined to map the lithospheric magnetic field over China by means of the revised spherical cap harmonic analysis(R-SHCA)modeling technique.The magnetic field is described to a spatial resolution of 150 km at the mean Earth’s radius,which represents a good compromise between the resolutions afforded by surface and satellite data.We compare the magnetic anomalies modeled at the regional scale with composite regions containing large-scale of tectonic structures.These regions,including the Tarim Basin and the Tibetan Plateau,are correlated with regional magnetic anomalies at satellite altitude but contain a significant number of small-scale and complex magnetic structures at the mean Earth’s radius.These magnetic anomalies are globally consistent with the known geological features in China but also offer a way to delineate the contours of the geological blocks and to discuss the connection between magnetic anomalies and the heat flow distribution in this region.     

汞异常与地震活动关系及预报效能研究进展

汞主要来源于地球深部地幔岩石和岩浆,具有溶解度小、蒸气压低、穿透能力强等特点.地震前,汞在高温和强压力梯度下穿过不同的介质到达地表,导致汞震前异常.汞的前兆异常特征通常被用作地震前预测的重要指标.本文综合国内应用汞的地震观测案例,得出以下结论:(1)与气态汞相比,水溶性汞具有收集方便、不易受传输介质干扰、前兆异常幅度大等特征,易震前观测;(2)水溶性汞观测异常结束到发震时间间隔较长,而气态汞异常一般持续到地震发生;(3)震前汞的异常幅度通常是本底值的数倍至几十倍,这往往与震中距离、震级大小及地质构造有关;(4)对于5级及以下地震,水溶性汞异常范围通常在200 km以内,5~6级地震映震范围在500 km以内.此外,汞震前观测异常与孕震应力场作用下的地下流体活动有关,汞的同位素分馏技术提供了一种新的手段来解决断裂带汞的来源以及汞与地壳深部流体活动的关系.     

Experimental errors of measurements of vertical geomagnetic gradients in the stratosphere

The errors of measurements of vertical geomagnetic gradients at altitudes of 20–40 km, using a balloon magnetic gradiometer with a 6-km-long measuring base oriented along gravity, have been studied in the work. The errors related to the deviation of the measuring base position relative to the vertical have been studied during the real balloon flight with the help of the navigation GPS receivers. The deviations of the measuring base within 5°, which can sometimes reach 15°, have been obtained. This results in a decrease in the magnetic gradient measurement accuracy due to the errors introduced in the specification of the normal magnetic field used to detect magnetic anomalies. To eliminate this error, a GPS receiver was built in each magnetometer in order to observe magnetometers during synchronous measurements and to correct the measurements for the normal magnetic field. It has been indicated that the effect of deviations of the measuring base position on the results is not more than 2% of the measured value at such organization of a gradiometer.     

Geophysical researches (gravity and magnetic) of the Eratosthenes Seamount in the eastern Mediterranean Sea

New free-air gravity and magnetic maps of the Eratosthenes Seamount and its vicinity were regenerated from potential field data. Stages of data processing are power spectrum, upward continuation, filtering on the free-air gravity anomaly data. RTP, pseudo-gravity transformation map, power spectrum, upward continuation, filtering, AS, and HGAS were applied on the magnetic data. A HGAS map shows the images and locations of the Eratosthenes magnetic body. Spectral analysis of the gravity and magnetic anomalies indicates that there is an elliptical elongated structure of the Eratosthenes Seamount in the width of approx. 86 km NW-SE orientation and in the length of 138 km NE-SW orientation, with a strike of N40°E and inclined to NW. It is considered that 22.49 ± 0.08 km obtained from power spectrum of the gravity data may be related to the crust thickness. Also, 15.67 ± 0.02 km obtained from power spectrum of the magnetic data is considered to be related to the magmatic basement of the Eratosthenes Seamount.     

Inversion of magnetic measurements of the CHAMP satellite over the Pannonian Basin

The Pannonian Basin is a deep intra-continental basin that formed as part of the Alpine orogeny. In order to study the nature of the crustal basement we used the long-wavelength magnetic anomalies acquired by the CHAMP satellite. The anomalies were distributed in a spherical shell, some 107,927 data recorded between January 1 and December 31 of 2008. They covered the Pannonian Basin and its vicinity. These anomaly data were interpolated into a spherical grid of 0.5° × 0.5°, at the elevation of 324 km by the Gaussian weight function. The vertical gradient of these total magnetic anomalies was also computed and mapped to the surface of a sphere at 324 km elevation. The former spherical anomaly data at 425 km altitude continued downward to 324 km. To interpret these data at the elevation of 324 km we used an inversion method. A polygonal prism forward model was used for the inversion. The minimum problem was solved numerically by the Simplex and Simulated annealing methods; a L₂ norm in the case of Gaussian distribution parameters and a L₁ norm was used in the case of Laplace distribution parameters. We interpret that the magnetic anomaly was produced by several sources and the effect of the sable magnetization of the exsolution of hemo-ilmenite minerals in the upper crustal metamorphic rocks.     

用最小二乘两步迭代法求解磁性球体几何与磁性参数

通过对常规最小二乘法在求解磁性球体参数过程中产生发散解的原因分析表明:非线性方程组中待定参数过多,特别是角度参数,是影响最小二乘法收敛性的主要因素.为此,提出将磁异常三分量作为观测值,用矢量磁矩作为待定参数,以替代磁化强度磁倾角和磁偏角,从而消除了非线性观测方程中的角度参数影响.根据观测值与磁矩的线性关系以及磁性球体中...     

Topographic heterogeneity of the core-mantle boundary in subduction zones and its effect on kinematics of main geomagnetic field sources

This work considers the kinematics of the source of the main geomagnetic field (MGF) near the core-mantle boundary under the Caribbean region. This source was selected because (i) its trajectory for the reviewed 110 years crosses this boundary, (ii) the region belongs to the so-called cemeteries of the tectonic plates, and (iii) numerous works have studied the structural heterogeneities of the lower mantle in this region with seismic tomography. Our study of the structural heterogeneities of the lower mantle and the trajectory of the MGF source indicates that the relics of the ancient tectonic plates in this region not only reach the coremantle boundary but could penetrate to the liquid core as well to the depth of 300 km. The “cemeteries” of the tectonic plates span significant areas in size. If the topographic heterogeneities of the core-mantle boundary, which are formed by the relics of the ancient tectonic plates, reach several hundreds of kilometers, then they can significantly affect the kinematics of individual structured flows in the liquid core and, consequently, change the spatial structure of secular MGF variations on the Earth’s surface.     

Magnetic studies of basalt fragments recovered by deep drilling in Iceland,and the “magnetic layer” concept

The magnetic susceptibility of 1300 samples of igneous rock drill cuttings obtained from eight deep drill holes in Iceland has been measured, in order to directly provide limits on the thickness of the layer which is the source of the magnetic anomalies over Iceland. The remanent magnetism of some of the material has also been studied, and the variation of magnetic susceptibility in 740 lava flows from eastern Iceland has been analysed as a function of depth of burial.All the results indicate no systematic change of susceptibility with depth up to 2.0 km. The Curie point of all deeply buried basalts in Iceland appears to be close to that of magnetite, so that the magnetic layer may be 5 km or more in thickness when susceptibility contrasts are considered; lateral contrasts in primary remanence may reach to 3 km depth. Derivation of a magnetic layer thickness in Iceland from analyses of magnetic anomalies, using methods which have been conventionally applied to marine magnetic anomalies could, on the other hand, yield much lower apparent thickness values (less than 1 km).We therefore argue that estimates of the magnetic layer thickness in oceanic regions should be based on considerations of magnetite Curie point isotherm behaviour, rather than on anomaly analysis.     

Magnetic fields in a limb flare on July 19, 2012

The magnetic field intensity in the limb solar flare of July 19, 2012, was measured by the Zeeman bisector splitting in the H _α line. An average magnetic field is established to attain 200 G at a measurement error of ±100 G in the top of a radiant flare loop at a height of about 40 Mm. The confinement of a strong magnetic field in the high coronal arcade at a relatively weak external field of the corona (about 1–2 G) is considered using a model of a force-free magnetic flux rope with a fine magnetic structure at a scale of about 300 km.     

Effect of clayey media parameters on the negative response of a coincident loop1

The aim of this work is to consider the phenomenon of negative response and its detectability in a polarizable half-space formed of Clay-water mixture which commonly exists in water aquifers and which can be used as a simple indicator of the presence of Underground water. We make use of an analytical treatment for the transient voltage induced in a coincident loop lying on a half-space as a basis for our computations. A Cole-Cole model is used to represent the Clay with appropriate parameters. As saturated Clay is characterized by high conductivity and moderate chargeability, it is hardly detectable relative to the practical noise level. However, there are optimum Clay parameters that offer a maximum negative response. It is also observed that for a very short time constant τ in the Cole-Cole model, the amplitude of the negative response decreases with the decrease of the time constant due to the fast decay of the polarization current. Finally, it is found that there is an optimum value for the loop radius which gives the largest value of the negative response and this loop radius depends on the model parameters.     

THE TRANSIENT ELECTROMAGNETIC RESPONSE OF A CONDUCTING SPHERE IN AN IMPERFECTLY CONDUCTING HALF-SPACE*

The presence of a conducting environment about a spherical ore body must be considered when calculating the transient electromagnetic response of the ore body due to a step current flowing in a large circular loop at the earth's surface. Failure to do this can easily lead to errors in excess of 10% in numerical calculations. Moreover, there is only a limited time interval in which the response of the spherical conductor is easily seen. In a poorly conducting ground the resonance response of the sphere is the first to be excited. Later, however, the non-resonance or wave-type response is excited. These waves destructively interfere and finally the response of the sphere decays with time as t^?7/2. For a range of times and depths the best loop for detecting the sphere has about the same radius as the sphere.     

利用地电场线性极化特性提取地震前兆异常的方法研究

处理平凉台2008—2017年的地电场观测资料,首先对地电场线性极化和非极化异常进行分析,提出利用地电场线性极化特性提取地震前兆异常的方法。其次,采用相关数据处理和计算方法对周期性质、极化系数、极化方向、垂直极化方向上的投影等进行讨论,结果表明平凉台地电场线性极化正常背景特征为:极化方向变化范围N(66.7±6)°E,极化系数在0.7～1之间,垂直极化方向上的投影值波动范围(-0.53～2.61) mV/km。最后,统计分析平凉台周边中强地震与地电场非线性极化异常的关系,认为平凉台地电场对中强地震有较好的映震能力,地电场非极化异常可以作为判断地震前兆异常的一个重要指标,其异常表现为:从震前非极化到震后恢复极化的一个过程,持续时间一般在4～15天,多出现在震前5天到1个月不等。     

A new least-squares minimization approach to depth and shape determination from magnetic data

We have developed a least‐squares minimization approach to depth determination using numerical second horizontal derivative anomalies obtained from magnetic data with filters of successive window lengths (graticule spacings). The problem of depth determination from second‐derivative magnetic anomalies has been transformed into finding a solution to a non‐linear equation of the form, f(z) = 0. Formulae have been derived for a sphere, a horizontal cylinder, a dike and a geological contact. Procedures are also formulated to estimate the magnetic angle and the amplitude coefficient. We have also developed a simple method to define simultaneously the shape (shape factor) and the depth of a buried structure from magnetic data. The method is based on computing the variance of depths determined from all second‐derivative anomaly profiles using the above method. The variance is considered a criterion for determining the correct shape and depth of the buried structure. When the correct shape factor is used, the variance of depths is less than the variances computed using incorrect shape factors. The method is applied to synthetic data with and without random errors, complicated regionals, and interference from neighbouring magnetic rocks. Finally, the method is tested on a field example from India. In all the cases examined, the depth and the shape parameters are found to be in good agreement with the actual parameters.     

A gradient‐ratio method for the semi‐automatic interpretation of gravity map data sets

The application of semi‐automatic interpretation techniques to potential field data can be of significant assistance to a geophysicist. This paper generalizes the magnetic vertical contact model tilt‐depth method to gravity data using a vertical cylinder and buried sphere models. The method computes the ratio of the vertical to the total horizontal derivative of data and then identifies circular contours within it. Given the radius of the contour and the contour value itself, the depth to the source can be determined. The method is applied both to synthetic and gravity data from South Africa. The Matlab source code can be obtained from the author upon request.     

南海珠江口盆地钻井BHT温度校正及现今地温场特征

温度是盆地现今地温场研究最重要的基础参数之一,针对钻井井底温度(BHT)一般低于地层真实温度的情况,本文对珠江口盆地16口钻井的BHT数据进行了校正,并根据校正后的温度数据采用最小二乘法新增计算了16个地温梯度及大地热流数据.结合前人发表数据,本文绘制了珠江口盆地最新地温梯度等值线图及大地热流等值线图.除此之外,系统收集整理了盆地钻井大地热流、生热率以及热导率等参数,利用一维热传导方程计算了盆地1~5km深度处的温度,并绘制了盆地深部温度等值线图.在此基础上,分析了盆地现今地温场特征及其成因,以期为珠江口盆地地球动力学背景研究及油气资源评价与勘探提供地热学参数.     

A direct method of interpreting gravity and magnetic anomalies: The case of a horizontal cylinder

Summary A new method of interpreting the gravity and magnetic anomalies is introduced with special reference to the magnetic anomalies of a horizontal cylinder. The method consists of calculating the functions of the anomaly and its distance from an arbitrary point. These form a simple linear equation with coefficients related to the parameters defining the body. Since each observation forms a separate linear equation, the required normal equations are formed by the method of least squares and solved for the coefficients and hence for the various parameters defining the target. The discussion here is confined to the vertical magnetic anomalies. The application of the method to horizontal and total field anomalies of two dimensional bodies is also outlined.     

Estimation of the normal geomagnetic field and geological interpretation of the normal field anomalies over middle-northern Croatia

Based on the ground survey of total-field magnetic data at 53 sites, which recorded a time-series over the interval 2003–2005, the core field + long wavelength lithospheric field over the middle-northern Croatia region was calculated. The area of the survey was 18900 km², an average distance between the neighbouring sites being 12 km. The results were reduced to 2004.5 epoch. This “normal” total-field was estimated using the firstorder Taylor polynomial as a function of geographic coordinates, and the polynomial coefficients were calculated with three methods of adjustment: simple and weighted least squares fits and adjustment according to the most frequent value. The stability of the normal field was tested using the Monte Carlo-type test, by decreasing the input data set in each adjustment (up to 90%). All obtained field residuals (measured-“normal”) were mutually compared, as well as fit coefficients of the Taylor polynomials. The residual values indicate the presence of shorter-wavelength anomalies, specifically three major residual anomalies were found (−106 nT, 74 nT and 57 nT). The geostatistical analysis of the ground survey data and the normal field residuals (respectively), using the median absolute deviation method, was further conducted in order to evaluate the calculated anomalies. The geological situation around the anomalies derived by the median absolute deviation method, and around the normal field residual anomalies, is given. The correlation was found between higher value anomalies of the normal field residuals, and shallow volcanic rocks and oil field, respectively.     

The zero-line method of interpreting total field magnetic anomalies of spherical ore-bodies

Summary A method called the zero-line method for interpreting the total field magnetic anomalies of spherical bodies is described. The contour of zero-anomaly is controlled by a linear equation of the typeAx ² +By ² +Cxy +Dx +Ey +F = 0, where the coefficientsA toF are related to the depth of the sphere and the parameters relating to its position and magnetisation. From the coordinates (x, y) of the various points lying on the zero-contour, five normal equations for the above equation are set and solved for the various coefficients. The various parameters relating to the position and magnetisation of the sphere are then determined from these coefficients. It is expected that similar methods can be developed for interpreting magnetic anomalies of other three dimensional bodies also.     

NWC通信台在电离层中激发电磁响应的时变特征

本文利用DEMETER卫星VLF频段电场和磁场频谱数据对DEMETER卫星运行期间2005年至2009年澳大利亚甚低频(Very Low Frequency)通信台NWC发射的通信信号造成的电离层电磁响应的日变化、季节变化及年变化特征进行了统计分析,统计结果表明电磁响应日变化显著,夜间电场强度明显增强可达40dB,磁场变化略小也可为15dB左右,而季节变化不显著,年变化主要受太阳活动的影响,太阳活动越强,电磁响应越小.为解释数据分析结果,对地-电离层电磁波传播过程采用传递矩阵方法进行了模拟计算,模拟结果与数据分析的结果一致.我们认为这种随时间变化的特点可能由250km以下电离层电子密度分布特征导致,因此研究250km以下的电离层电子密度变化可能对寻找地震电离层电磁异常有重要意义.     

Dynamics of escape of high-energy electrons from a spherical plasma volume in the geomagnetic field

Using results of numerical modeling, the dynamics of escape of electrons, produced as a result beta-decay, to the external magnetic field from a spherical plasma volume with an expelled magnetic field is studied. The dependence of the fraction of escaped electrons on time and radius of the plasma volume has been obtained for two kinds of electron sources: a point isotropic source, located at the center of the sphere, and a volume isotropic source. It is shown that for a point source some part of electrons remains in the sphere, whereas for a volume source all electrons, at different values of the magnetic cavern radius, leave it at an identical relative escape rate.