English

Computational algorithm for downward continuation of gravity anomalies

Numerical Analysis 2025-07-14 v1 Numerical Analysis

Abstract

The downward continuation of potential fields from the Earth's surface into the subsurface is a critical task in gravity exploration, as it helps to identify the sources of gravity anomalies. This problem is often addressed by solving a first-kind integral equation using regularization techniques to stabilize an inherently unstable process. A similar approach is used in our work, where the continued field is represented as the potential of a simple layer or its vertical derivative. The constancy of the density sign of this equivalent simple layer preserves the sign of anomalies, provided that the layer's surface encloses all anomalous sources. This constraint is a key feature of our algorithm for the downward continuation of potential fields. To enforce, for instance, non-negativity in the simple layer density, we employ the NNLS (Non-Negative Least Squares) method. The efficiency of the proposed method is demonstrated on model examples.

Keywords

Cite

@article{arxiv.2507.08506,
  title  = {Computational algorithm for downward continuation of gravity anomalies},
  author = {D. K. Ivanov and L. N. Temirbekova and P. N. Vabishchevich},
  journal= {arXiv preprint arXiv:2507.08506},
  year   = {2025}
}

Comments

17 pages, 11 figures

R2 v1 2026-07-01T03:56:26.492Z