New analytical solution and associated software for computing full-tensor gravitational field due to irregularly shaped bodies

Saraswati, Anita Thea; Cattin, Rodolphe; Mazzotti, S.; Cadio, C.

doi:10.1007/s00190-019-01309-y

Cited by 14 publications

(12 citation statements)

References 36 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Three different algorithms, including the GEEC software package developed in Saraswati et al (2019), the vectorized MATLAB code VecWerSch developed by us based on the algorithm in Werner and Scheeres (1997) is assumed for the asteroid. Forward results show that the difference between the two analytical methods is predictably small.…”

Section: Asteroids 433 Eros and 101955 Bennumentioning

confidence: 99%

Modified Parker's Method for Gravitational Forward and Inverse Modeling Using General Polyhedral Models

2021

JGR Solid Earth

View full text Add to dashboard Cite

Parker's method (Parker, 1973), initially introduced by the author for the rapid calculation of gravitational and magnetic anomalies, of 2D and 3D complexly layered models based on the FFT algorithm, is by no doubt one of the most widely applied algorithms in potential field data processing and interpretation (Nabighian, Ander, et al., 2005;. By raising the interface functions of a layered model to the E n th power and transforming into Fourier series using the FFT algorithm, potential anomalies on a horizontal plane above the layered model can be computed much more efficiently than using space-domain numerical integration techniques. The method has been proved particularly useful in solving forward problems such as gravimetric terrain correction and isostatic correction (

show abstract

Section: Asteroids 433 Eros and 101955 Bennumentioning

confidence: 99%

Modified Parker's Method for Gravitational Forward and Inverse Modeling Using General Polyhedral Models

2021

JGR Solid Earth

View full text Add to dashboard Cite

show abstract

“…We model the IBM slab using the GEEC algorithm (Saraswati et al 2019) on Hayes et al (2012) slab geometry. This method uses irregular polyhedrons, which enables to almost perfectly reproduce the slab geometry.…”

Section: Benefit Of the Tensor Componentsmentioning

confidence: 99%

“…As previously, we reproduced the geometry of the IBM slab from the Slab1.0 model (Hayes et al 2012) with the GEEC algorithm (Saraswati et al 2019) and computed the gravity and gravity gradient responses at the satellite height. We computed the gravity signal at the same height (225km) to ease the direct comparison between G and GG inversion results.…”

Section: Comparison At the Satellite Heightmentioning

confidence: 99%

“…We estimate the complementarity between satellite GG data and ground G data on the realistic geometry of the IBM case. The satellite GG data are computed as previously mentioned (Saraswati et al 2019). We build the synthetic ground gravity data set using real data coverage extracted from the BGI (Bureau Gravimetrique International) database.…”

Section: Complementarity Between Satellite Gravity Gradient and Ground Gravity Datamentioning

confidence: 99%

See 1 more Smart Citation

From space to lithosphere: inversion of the GOCE gravity gradients. Supply to the Earth’s interior study

Plasman

Tiberi

Cadio

et al. 2020

Geophysical Journal International

Self Cite

View full text Add to dashboard Cite

SUMMARY The emergence of high resolution satellite measurements of the gravitational field (GOCE mission) offers promising perspectives for the study of the Earth’s interior. These new data call for the development of innovant analysis and interpretation methods. Here we combine a forward prism computation with a Bayesian resolution approach to invert for these gravity gradient data configuration. We apply and test our new method on satellite data configuration, that is 225 km height with a global and homogeneous geographic distribution. We first quantify the resolution of our method according to both data and parametrization characteristics. It appears that for reasonable density contrast values (0.1 g cm−3) crustal structures have to be wider than ∼28 km to be detectable in the GOCE signal. Deeper bodies are distinguishable for greater size (35 km size at 50 km depth, ∼80 km at 300 km depth). We invert the six tensor components, among which five are independent. By carefully testing each of them and their different combinations, we enlighten a trade off between the recovery of data and the sensitivity to inversion parameters. We particularly discussed this characteristic in terms of geometry of the synthetic model tested (structures orientation, 3-D geometry, etc.). In terms of RMS value, each component is always better explained if inverted solely, but the result is strongly affected by the inversion parametrization (smoothing, variances, etc.). On the contrary, the simultaneous inversion of several components displays a significant improvement for the global tensor recovery, more dependent on data than on density variance or on smoothness control. Comparing gravity and gradient inversions, we highlight the superiority of the GG data to better reproduce the structures especially in terms of vertical location. We successfully test our method on a realistic case of a complex subduction case for both gradient and gravity data. While the imaging of small crustal structures requires terrestrial gravity data set, the longest wavelength of the slab is well recovered with both data sets. The precision and homogeneous coverage of GOCE data however, counterbalance the heterogeneous and often quite non-existence coverage of terrestrial gravity data. This is particularly true in large areas which requires a coherent assemblage of heterogeneous data sets, or in high relief, vegetally covered and offshore zones.

show abstract

“…More recent additions to the related literature include new mathematical approaches that eliminate singularities or expand the analytical formulation to include non-homogeneous distributions (D'Urso, 2013(D'Urso, , 2014D'Urso and Trotta, 2017;Ren et al, 2018Ren et al, , 2020 as well as applications covering a wide variety of areas in geodesy, geophysics and astronomy (Schmidt et al, 2011;Bucha et al, 2016;Chanut et al, 2017;Fukushima, 2017;Jiang et al, 2018;Zhang and Chen, 2018;Benedek et al, 2018;Capponi et al, 2018;Saraswati et al, 2019;Vitale and Fedi, 2020).…”

Section: Introductionmentioning

confidence: 99%

A computational review of the line integral analytical formulation of the polyhedral gravity signal

Tsoulis

Gavriilidou

2021

Geophysical Prospecting

View full text Add to dashboard Cite

The generally shaped polyhedron is a widely used model in gravity field modelling and interpretation. Its induced gravity signal -gravitational potential and its derivatives up to second order -has been studied extensively in the geophysical literature. The class of solutions with special interest is the one which leads to closed analytical expressions, as these offer an exact representation of the gravity signal of finite threedimensional distributions while being at the same time linked to a flexible means of geometric modelling. Of the several mathematical algorithms available, the line integral approach involves no approximations and is, however, connected with certain singularities, occurring for specific relative positions of the computation point with respect to the polyhedral source. The present contribution analyses the algorithmic details of the polyhedral line integral approach focusing on its geometric and computational aspects. Following the definitions of the individual coordinate systems and the two-step application of the Gauss divergence theorem for each face and polygonal boundary of the source, the geometric insight of the algorithm is presented, which permits a deeper understanding of the significance of the involved numerical singular terms. An overview of the line integral analytical approach for the polyhedral gravity signal is presented with emphasis on its geometric and computational aspects. Matlab code and results for the two considered case studies, a prismatic source and asteroid Eros, are provided as electronic supplement.

show abstract

New analytical solution and associated software for computing full-tensor gravitational field due to irregularly shaped bodies

Cited by 14 publications

References 36 publications

Modified Parker's Method for Gravitational Forward and Inverse Modeling Using General Polyhedral Models

Modified Parker's Method for Gravitational Forward and Inverse Modeling Using General Polyhedral Models

From space to lithosphere: inversion of the GOCE gravity gradients. Supply to the Earth’s interior study

A computational review of the line integral analytical formulation of the polyhedral gravity signal

Contact Info

Product

Resources

About