Fast 3D joint inversion of gravity and magnetic data based on cross gradient constraint

The gravity and magnetic data can be adopted to interpret the internal structure of the Earth.To improve the calculation efficiency during the inversion process and the accuracy and reliability of the reconstructed physical property models,the triple strategy is adopted in this paper to develop a fast cross-gradient joint inversion for gravity and magnetic data.The cross-gradient constraint contains solving the gradients of the physical property models and performing the cross-product calculation of their gradients.The sparse matrices are first obtained by calculating the gradients of the physical property models derived from the first-order finite difference.Then,the triple method is applied to optimize the storages and the calculations related to the gradients of the physical property models.Therefore,the storage compression amount of the calculations related to the gradients of the physical property models and the cross-gradient constraint are reduced to one-fold of the number of grid cells at least,and the compression ratio increases with the increase of the number of grid cells.The test results from the synthetic data and field data prove that the structural coupling is achieved by using the fast cross-gradient joint inversion method to effectively reduce the multiplicity of solutions and improve the computing efficiency.

Joint inversion of gravity and vertical gradient data based on modified structural similarity index for the structural and petrophysical consistency constraint

Joint inversion is one of the most effective methods for reducing non-uniqueness for geophysical inversion.The current joint inversion methods can be divided into the structural consistency constraint and petrophysical consistency constraint methods,which are mutually independent.Currently,there is a need for joint inversion methods that can comprehensively consider the structural consistency constraints and petrophysical consistency constraints.This paper develops the structural similarity index(SSIM)as a new structural and petrophysical consistency constraint for the joint inversion of gravity and vertical gradient data.The SSIM constraint is in the form of a fraction,which may have analytical singularities.Therefore,converting the fractional form to the subtractive form can solve the problem of analytic singularity and finally form a modified structural consistency index of the joint inversion,which enhances the stability of the SSIM constraint applied to the joint inversion.Compared to the reconstructed results from the cross-gradient inversion,the proposed method presents good performance and stability.The SSIM algorithm is a new joint inversion method for petrophysical and structural constraints.It can promote the consistency of the recovered models from the distribution and the structure of the physical property values.Then,applications to synthetic data illustrate that the algorithm proposed in this paper can well process the synthetic data and acquire good reconstructed results.

Characteristics of subsurface density variations before the 4.20 Lushan M_s7.0 earthquake in the Longmenshan area: inversion results

The 4.20 Lushan Ms7.0 earthquake occurred on the southwest segment of the Longmenshan fault on 20 April 2013. Some meaningful information on the prepa- ration and occurrence of this earthquake was found based on the dynamic variation of gravity （DVG）. To examine the great progress of the Lushan earthquake, we obtained the density variation （DENV） derived from the DVG using the compact gravity inversion method in this article. The inversion results reveal three main findings： （1） the DENV in the crust in the Jinshajiang fault area changed from positive in 2010-2011 to negative in 2011-2012. （2） The DENV in the Xianshuihe fault area decreased continuously from 2010 to 2012. （3） The DENV of the uppermost mantle of South China decreased in 2010-2011 and increased in 2011-2012. We propose that the flow/expansion of the middle-lower crust beneath the Bayan Har block and Moho subsidence on the southwest margin of the Chuan-Dian block may have been the major causes of the Lushan earthquake.

Isostatic anomaly characteristics and tectonism of the New Britain Trench and neighboring Papua New Guinea

In this paper, we calculated Bouguer gravity anomalies and Airy-Heiskanen isostatic anomalies in the New Britain Trench and surrounding areas of Papua New Guinea. The calculations are based on a topographic model and a gravity anomaly model from the Scripps Institute of Oceanography. Characteristics of the isostatic anomaly and the earthquake dynamic environment of this region are analyzed,and the results show that there are obvious differences in the isostatic state between each block in the region. Tectonic movements are very intense in regions with high positive or negative isostatic gravity anomalies, and a number of sub-plates in this area are driven by external tectonic action, such as plate subduction and thrusting, of the Pacific, Indo-Australian and Eurasian plates. From the distribution of the isostatic gravity anomaly, the tectonic action of anti-isostatic movement obviously is the main source of power. Based on isostatic gravity and spatial distribution of earthquakes in the region, it is obvious that further contraction of the Indo-Australian Plate will result in the southwestern part of the Solomon Plate becoming part of the Owen Stanley thrust belt, and the northern part will enter the lower part of the Bismarck Plate. The eastern part of the Solomon Plate will enter the front of the Pacific Plate, resulting in northward and eastward migration of significant earthquakes along the Solomon Plate.

Crustal structure of the Qiangtang and Songpan-Ganzi terranes(eastern Tibet) from the 2-D normalized full gradient of gravity anomaly

Numerous geophysical studies have revealed the lithospheric structure of the Qiangtang and the Songpan-Ganzi terranes in the eastern Tibetan Plateau.However,crust-mantle evolution and crustal response to the Indian lithospheric subduction are still controversial.Answering these questions requires additional information regarding crustal structure.In this study,the 2-D normalized full gradient(NFG)of the Bouguer gravity anomaly was used to investigate anomalous sources and interpret the crustal structure underneath the Qiangtang and Songpan-Ganzi terranes.The NFG-derived structures with loworder harmonic numbers(N=33 and N=43)showed that an anomalous source beneath the southern Qiangtang terrane had a characteristic northeastward-dipping shape,suggesting the northeastward motion of the crustal material induced by underthrusting Indian lithospheric mantle.The NFG images with harmonic number N=53 showed a large-scale anomalous source in the lower crust of the transformational zone from the Qiangtang terrane to the Songpan-Ganzi terrane,consistent with thickening crust and resistance of lower crustal flow.The anomalous source demonstrated by the NFG results with harmonic number N=71,located in the upper crust underneath the Ganzi-Yushu fault,suggested a seismogenic body of the 2010 M6.9 Yushu event.