Summary. We regard the gridding of aeromagnetic data as a geophysical estimation problem and propose to use the inverse interpolation method for gridding aerom

We propose the 3‐D correlation imaging for gravity anomalies and gravity gradiometry data. We also propose to use the 3‐D correlation imaging based on anomaly separation to improve the resolution of the imaging. The test on the gravity anomalies and the gravity gradiometry data of the synthesized Y‐type dykes model and the synthesized multiple rectangular prism model show that the method can image the spatial extension states and the distribution of equivalent excess masses of geological anomalous bodies with high resolution both vertically and horizontally.

The technique of teleseismic receiver function H–κ stacking is popular for estimating the crustal thickness and Vp/Vs ratio. However, it has large uncertainty or ambiguity when the Moho multiples in receiver function are not easy to be identified. We present an improved technique to estimate the crustal thickness and Vp/Vs ratio by joint constraints of receiver function and gravity data. The complete Bouguer gravity anomalies, composed of the anomalies due to the relief of the Moho interface and the heterogeneous density distribution within the crust, are associated with the crustal thickness, density, and Vp/Vs ratio. According to their relationship formulas presented by Lowry and Pérez-Gussinyé (2011), we invert the complete Bouguer gravity anomalies by using a common algorithm of likelihood estimation to obtain the crustal thickness and Vp/Vs ratio, and then utilize them to constrain the receiver function H–κ stacking result. We verified the improved technique on three synthetic crustal models and evaluated the influence of selected parameters, the results of which demonstrated that the novel technique could reduce the ambiguity and enhance the accuracy of estimation. Real data test at two given stations in the NE margin of Tibetan Plateau illustrated that the improved technique provided reliable estimations of crustal thickness and Vp/Vs ratio.

Three-dimensional magnetic inversion, based on the least-square and regularization algorithm in the space domain, is an important tool for quantitative interpretation of magnetic data. However, the common 3D inversion approaches usually require great numbers of forward and inversion calculations and cause low efficiency for inverting large-scale data. Three-dimensional imaging is an alternate rapid tool for qualitative and quantitative interpretation of magnetic data. In this paper, we present a wavenumber-domain iterative approach for 3D imaging of magnetic anomalies and gradients, which could increase imaging efficiency and is suitable for rapidly imaging large-scale data. The wavenumber-domain formulas for forward modeling and imaging of total magnetic anomaly, three magnetic components, magnetic gradients and magnetic full-tensor gradients are deduced and provided. A depth-scale factor and the constraints of magnetic interface are included into the imaging formulas to enhance depth resolution. An iterative algorithm is adopted for the imaging to reduce the fitting error and improve the imaging accuracy. Tests on synthetic and real data from the Sichuan basin, China, verified the feasibility of the presented approaches.