Abstract:The gravity and gradient anomalies contain valuable information about the underground geological structures at various depths. Deep and shallow buried source bodies are able to be identified through multi-scale field separation processes, and visual comprehensions of geological structures can be obtained via 3D density inversion techniques. In this study, we propose an improved 3D imaging strategy based on gravitational field separation using the preferential continuation filter. This strategy incorporates the relationship between spectral features and buried depths of source bodies, allowing for a one-step transformation from planar gravity and full-tensor gradient field observations to a 3D density structure in the wave-number domain. Synthetic tests validate the effectiveness and robustness of the gravity and gradient imaging approaches, highlighting their advantages in high vertical resolution and low computational requirements. Nonetheless, it should be noted that the imaging effects of horizontal gradients Γxx and Γyy are unsatisfactory due to their weak noise resistance. Thus, they are not suitable for real data applications. The other imaging approaches are further applied to recover the subsurface 3D density structure beneath the Weishan cone in Wudalianchi Volcanic Field, Northeastern China. Our results provide insights into the possible location and shape of the low-density magma chamber. Also, the potential presence of partial melts is inferred and supported from a gravity perspective. The primary advantage of these approaches is their ability to generate a reasonable geological model in scenarios with limited prior information and physical property constraints. As a result, they have significant practical value in the field of applied geophysics, including mineral exploration and volcanology studies.

What problem does this paper attempt to address?

The paper mainly aims to address the following issues: ### Research Background and Objectives - **Research Background**: Gravity and its gradient anomalies can reflect information about underground geological structures. Through a multi-scale field separation process, it is possible to identify deep and shallow buried sources and obtain a visual understanding of geological structures using 3D density inversion techniques. - **Main Objectives**: Propose an improved 3D imaging strategy based on gravity field separation and optimized continuation filters to achieve a direct transformation from planar gravity and full tensor gradient field observations to 3D density structures in the wavenumber domain. ### Issues Addressed - **Improve Vertical Resolution and Computational Efficiency**: By integrating optimized continuation filtering methods, this strategy can enhance the accurate identification of overlapping field sources and clearly image complex underground density distributions and geological structures. - **Overcome Limitations of Existing Methods**: Traditional field separation methods have two main limitations in the 3D imaging process: difficulty in accurately correlating separation scales with the burial depth of mass layers, and lack of integration between field separation steps and physical property imaging, which may introduce additional errors. - **Application Challenges**: The study also focuses on challenges in practical applications, especially in cases with scarce prior information, by adjusting the field separation effect to enhance the accuracy of inversion results. ### Main Contributions - **Introduction of a New Method**: A new gravity field imaging strategy is proposed, combining optimized continuation filtering methods, which can directly transform planar field data to 3D density structures in the wavenumber domain in a single step. - **Theoretical and Practical Validation**: - Theoretically, the mathematical principles of optimized continuation filtering and the imaging method of the gravity field in the wavenumber domain are detailed. - The effectiveness and robustness of the method are verified through synthetic tests, particularly excelling in high vertical resolution and low computational requirements. - Practical data (gravity and gradient data from the Wudalianchi volcanic area) further validate the practicality of the method, especially in identifying the location and shape of low-density magma chambers. ### Summary The method proposed in this paper aims to improve the accuracy and efficiency of 3D gravity imaging, particularly suitable for geological model construction with limited prior information. By optimizing parameter settings, this method can effectively distinguish geological structures at different depths, offering significant practical value for mineral exploration and volcanology research.

Improved Approaches for 3D Gravity and Gradient Imaging Based on Potential Field Separation: Application to the Magma Chamber in Wudalianchi Volcanic Field, Northeastern China

3D Inversion of Gravity Anomaly – an Application of Technology of Apparent Density Imagery

A Fast Interpretation Method of Gravity Gradiometry Data Based on Magnetic Dipole Localization

The inversion of apparent density based on the separation and continuation of potential field

3D density imaging using gravity and gravity gradient in the wavenumber domain and its application in the Decorah

Two-dimensional Inversion of CSAMT Data and Three-Dimensional Geological Mapping for Groundwater Exploration in Tongkeng Area, Hunan Province, China

3D hybrid imaging based on gravity migration and regularized focusing inversion to predict the Poyang Basin interface

Quantifying the Effects of Near-Surface Viscosity on Seismic Acquisition Geometry: A case study from Chepaizi Exploration Area, Junggar Basin (NW China)

Gravity and Magnetic Data Inversion With Random Projection

Gravity Inversion in Spherical Coordinates With Dynamic Re-Weighting Matrix

3D gravity and magnetic inversion reveal geothermal structure dominated by radioactive heat production of granites in the Weihe basin

A New Body‐wave Amplitude Ratio‐based Method for Imaging Shallow Crustal Structure and Its Application in the Sichuan Basin, Southwestern China

A Density-Weighted and Cross-Gradient Constrained Joint Inversion Method of Gravity and Vertical Gravity Gradient Data in Spherical Coordinates and Its Application to Lunar Data

Joint Inversion of Gravity and Gravity Gradient Data Based on Cross-Gradient Function

Gravity inversion with a binary structure constraint imposed by magnetotelluric data and its application to a Carlin-type gold deposit

3-D probability density imaging of Euler solutions using gravity data: a case study of Mount Milligan, Canada

Supplemental Material: Magma Recharging Beneath the Weishan Volcano of the Intraplate Wudalianchi Volcanic Field, Northeast China, Implied from 3-D Magnetotelluric Imaging

A Deep Learning Gravity Inversion Method Based on a Self-Constrained Network and Its Application

Gravity and magnetic joint imaging based on Gramian constraints

Structure-guided gravity inversion for layered density modeling with an application in the Chezhen Depression, Bohai Bay Basin

Gradient descent fusion for gravity and magnetic data