Discovering Crustal Deformation Bands by Processing Regional Gravity Field

Objectives： This article presents a new computational procedure to discover scratches buried in the earth＇s crust. We also validate this new interdisciplinary analysis method with regional gravity data located in a well-known Dabie orogenic zone for test. Methods： Based on the scratch analysis method evolved with mathematical morphology of surfaces, we present a procedure that extracts information of the crustal scratches from regional gravity data. Because the crustal scratches are positively and highly correlated to crustal deformation bands, it can be used for delineation of the crustal deformation belts. The scratches can be quantitatively characterized by calculation of the ridge coefficient function, whose high value traces delineate the deformation bands hidden in the regional gravity field. In addition, because the degree of crustal deformation is an important indicator of tectonic unit divisions, so the crust can be further divided according to the degree of crustal deformation into some tectonic units by using the ridge coefficient data, providing an objective base map for earth scientists to build tectonic models with quantitative evidence. Results： After the ridge coefficients are calculated, we can further enhance the boundary of high ridge-coefficient blocks, resulting in the so-called ridge-edge coefficient function. The high-value ridge-edge coefficients are well correlated with the edge faults of tectonic units underlay, providing accurate positioning of the base map for compilation of regional tectonic maps. In order to validate this new interdisciplinary analysis method, we select the Dabie orogenic zone as a pilot area for test, where rock outcrops are well exposed on the surface and detailed geological and geophysical surveys have been carried out. Tests show that the deformation bands and the tectonic units, which are conformed by tectonic scientists based on surface observations, are clearly displayed on the ridge and ridge-edge coefficient images obtained in this article. Moreover, these computer-generated images provide more accurate locations and geometric details. Conclusions： This work demonstrates that application of modern mathematical tools can promote the quantitative degree in research of modern geosciences, helping to open a door to develop a new branch of mathematical tectonics.

Investigation of the Tikhonov Regularization Method in Regional Gravity Field Modeling by Poisson Wavelets Radial Basis Functions

The application of Tikhonov regularization method dealing with the ill-conditioned problems in the regional gravity field modeling by Poisson wavelets is studied. In particular, the choices of the regularization matrices as well as the approaches for estimating the regularization parameters are investigated in details. The numerical results show that the regularized solutions derived from the first-order regularization are better than the ones obtained from zero-order regularization. For cross validation, the optimal regularization parameters are estimated from L-curve, variance component estimation（VCE） and minimum standard deviation（MSTD） approach, respectively, and the results show that the derived regularization parameters from different methods are consistent with each other. Together with the firstorder Tikhonov regularization and VCE method, the optimal network of Poisson wavelets is derived, based on which the local gravimetric geoid is computed. The accuracy of the corresponding gravimetric geoid reaches 1.1 cm in Netherlands, which validates the reliability of using Tikhonov regularization method in tackling the ill-conditioned problem for regional gravity field modeling.

Determination of the degree 120 time-variable gravity field in the Sichuan-Yunnan region using Slepian functions and terrestrial measurements

The terrestrial time-variable gravity measurements are characterized by a high signal-to-noise ratio and sensitivity to the sources of mass change in the Earth's crust.These gravity data have many applications,such as surface deformation,groundwater storage changes,and mass migration before and after earthquakes.Based on repeated terrestrial gravity measurements at 198 gravity stations in the Sichuan-Yunnan region(SYR)from 2015 to 2017,we determine a time series of degree 120 gravity fields using the localized spherical harmonic(Slepian)basis functions.Our results show that adopting the first 6 Slepian basis functions is sufficient for effective localized Slepian modeling in the SYR.The differences between two gravity campaigns at the same time of year show an obvious correlation with tectonic features.The degree 120 timevariable gravity models presented in this paper will benefit the study of the regional mass migration inside the crust of the SYR and supplement the existing geophysical models for the China Seismic Experimental Site.