ABSTRACT
We have developed the open-source software Tesseroids, a set of command-line programs to perform forward modeling of gravitational fields in spherical coordinates. The software is implemented in the C programming language and uses tesseroids (spherical prisms) for the discretization of the subsurface mass distribution. The gravitational fields of tesseroids are calculated numerically using the Gauss-Legendre quadrature (GLQ). We have improved upon an adaptive discretization algorithm to guarantee the accuracy of the GLQ integration. Our implementation of adaptive discretization uses a “stack-based” algorithm instead of recursion to achieve more control over execution errors and corner cases. The algorithm is controlled by a scalar value called the distance-size ratio () that determines the accuracy of the integration as well as the computation time. We have determined optimal values of for the gravitational potential, gravitational acceleration, and gravity gradient tensor by comparing the computed tesseroids effects with those of a homogenous spherical shell. The values required for a maximum relative error of 0.1% of the shell effects are for the gravitational potential, for the gravitational acceleration, and for the gravity gradients. Contrary to previous assumptions, our results show that the potential and its first and second derivatives require different values of to achieve the same accuracy. These values were incorporated as defaults in the software.
REFERENCES
- , 2012, GOCE satellite derived gravity and gravity gradient corrected for topographic effect in the South Central Andes region: GOCE derivatives in the south central Andes: Geophysical Journal International, 190,
941–959 , doi:10.1111/gji.2012.190.issue-2 .GJINEA 0956-540X - , 2007, Spherical prism gravity effects by Gauss-legendre quadrature integration: Geophysical Journal International, 169,
1–11 , doi:10.1111/gji.2007.169.issue-1 .GJINEA 0956-540X - , 2006, Multiple root finder algorithm for legendre and Chebyshev polynomials via Newton’s method: Annales Mathematicae et Informaticae,
3–13 . - , 2013a, Reference frame transformation of satellite gravity gradients and topographic mass reduction: Journal of Geophysical Research: Solid Earth, 118,
759–774 , doi:10.1029/2012JB009747 . - , 2013b, GOCE gravity gradient data for lithospheric modeling: International Journal of Applied Earth Observation and Geoinformation, 35,
16–30 , doi:10.1016/j.jag.2013.11.001 . - , 2015, Exploration of tectonic structures with GOCE in Africa and across-continents: International Journal of Applied Earth Observation and Geoinformation, 35,
88–95 , doi:10.1016/j.jag.2014.01.013 . - , 2013, The European Alps and nearby orogenic belts sensed by GOCE: Bollettino di Geofisica Teorica e Applicata, 54,
321–334 . - , 2011, The enigmatic Chad lineament revisited with global gravity and gravity-gradient fields: Geological Society, Special Publications, 357,
329–341 , doi:10.1144/SP357.18 . - , 2014, Perturbing effects of sub-lithospheric mass anomalies in GOCE gravity gradient and other gravity data modelling: Application to the atlantic-mediterranean transition zone: International Journal of Applied Earth Observation and Geoinformation, 35,
54–69 , doi:10.1016/j.jag.2014.02.003 . - , 2013, Optimized formulas for the gravitational field of a tesseroid: Journal of Geodesy, 87,
645–660 , doi:10.1007/s00190-013-0636-1 .JOGEF8 0949-7714 - , 2007, A comparison of the tesseroid, prism and point-mass approaches for mass reductions in gravity field modelling: Journal of Geodesy, 81,
121–136 , doi:10.1007/s00190-006-0094-0 .JOGEF8 0949-7714 - , 1987, Introduction to numerical analysis: Dover Publications.
- , 2007, Matplotlib: A 2D graphics environment: Computing in Science & Engineering, 9,
90–95 , doi:10.1109/MCSE.2007.55 . - , 1977, A direct computation of gravity and magnetic anomalies caused by 2-and 3-dimensional bodies of arbitrary shape and arbitrary magnetic polarization by equivalent-point method and a simplified cubic spline: Geophysics, 42,
610–622 , doi:10.1190/1.1440732 .GPYSA7 0016-8033 - , 1991, An exact solution for the gravity curvature (Bullard B) correction: Geophysics, 56,
1179–1184 , doi:10.1190/1.1443138 .GPYSA7 0016-8033 - , 2011, An efficient and adaptive approach for modeling gravity effects in spherical coordinates: Journal of Applied Geophysics, 73,
221–231 , doi:10.1016/j.jappgeo.2011.01.004 .JAGPEA 0926-9851 - , 2013, Explaining the thick crust in Paraná Basin, Brazil, with satellite GOCE gravity observations: Journal of South American Earth Sciences, 45,
209–223 , doi:10.1016/j.jsames.2013.03.008 .JAESE9 0895-9811 - , 2010, Data structures for statistical computing in Python:
Proceedings of the 9th Python in Science Conference ,51–56 . - , 2006, Estimation of distant relief effect in gravimetry: Geophysics, 71, no. 6,
J59–J69 , doi:10.1190/1.2338333 .GPYSA7 0016-8033 - , 2007, IPython: a system for interactive scientific computing: Computing in Science and Engineering, 9,
21–29 , doi:10.1109/MCSE.2007.53 .CSENFA 1521-9615 - , 2011, Mayavi: 3D visualization of scientific data: Computing in Science & Engineering, 13,
40–51 , doi:10.1109/MCSE.2011.35 . - , 2013, Global Moho from the combination of the CRUST2.0 model and GOCE data: Geophysical Journal International, 195,
222–237 , doi:10.1093/gji/ggt247 .GJINEA 0956-540X - , 2015a, A tesserioid (spherical prism) in a geocentric coordinate system with a local-north-oriented coordinate system: figshare, http://dx.doi.org/10.6084/m9.figshare.1495525, accessed 20 May 2016.
UiedaL. , 2015b, Tesseroids: Forward modeling in spherical coordinates. Open-source package, http://tesseroids.leouieda.com/en/latest/ and http://zenodo.org/record/16033, accessed 30 June 2016.- , 2016, Source code repository for “Tesseroids: Forward-modeling gravitational fields in spherical coordinates”, https://github.com/pinga-lab/paper-tesseroids, accessed 20 May 2016.
- , 2008, A comparison of different mass elements for use in gravity gradiometry: Journal of Geodesy, 82,
637–653 , doi:10.1007/s00190-008-0219-8 .JOGEF8 0949-7714
