An Empirical Method for the Optimal Setting of the Potential Fields Inverse Problem

The use of potential field methods for geophysical exploration purposes is nowadays quite common: these techniques consent to retrieve geological knowledge over extended regions and can give complementary information where other invasive or expensive techniques, such as seismic acquisitions, fail (e.g., in the recovery of geometries of geological horizons beneath a thick salt layer). Recent dedicated satellite gravity and magnetic missions, such as GRACE, GOCE and SWARM together with the exploitation of offshore satellite altimetry and airborne/shipborne surveys, have paved the way to the realization of a variety of global models, characterized by spatial resolutions of about 4 km (both for gravity anomaly and lithosphere magnetic anomalies) and high accuracy (about 3-5$3\text{-}5$ mGal and 20 nT). These models are a valuable source of information to study the geological evolution and characterization of the lithosphere structure, especially at a regional scale. In the present work, some preliminary technical aspects related to the use of these models to perform three-dimensional inversion are discussed, thus defining an empirical but rigorous procedure to set up gravity and magnetic inversion. In particular, we address the questions whether the classical planar approximation is acceptable for regional inversions or if a spherical one is required. We also provide guidance for choosing the best gravity functional (e.g., gravity anomalies or second radial derivative of the anomalous potential) and the optimal sizing of the three-dimensional volume area to be modelled depending on the specific target investigated. The application of the proposed methods to the Mediterranean case study is also presented.