Japan Geoscience Union Meeting 2023

Presentation information

[J] Oral

S (Solid Earth Sciences ) » S-GD Geodesy

[S-GD01] Geodesy and Global Geodetic Observing System

Wed. May 24, 2023 9:00 AM - 10:30 AM 304 (International Conference Hall, Makuhari Messe)

convener:Yusuke Yokota(Institute of Industrial Science, The University of Tokyo), Yuta Mitsui(Department of Geosciences, Shizuoka University), Koji Matsuo(Geospatial Information Authority of Japan), Chairperson:Toshimichi Otsubo(Hitotsubashi University), Takahito Kazama(Graduate School of Science, Kyoto University)

9:30 AM - 9:45 AM

[SGD01-13] Second-order smoothness prior over the Delaunay Tessellation and its application to gravity Bayesian inversion.

*Yuanyuan Niu1, Jiancang Zhuang2,1 (1.The Graduate University for Advanced Studies, 2.The Institute of Statistical Mathematics)


Prior information is always used to form up additional restrictions in geophysical inversions to solve the non-uniqueness problem of the solution. The smoothness (second-order derivative) of the model is one of such important restrictions. Smoothness is usually calculated through interpolation over the regular grids for the reason of easy implementation in numerical calculation. When observed data are irregularly distributed such as in geodetic inversions, Delaunay Tessellation (DT) based interpolation is popularly used to avoid additional interpolations. However, the numerical calculation of the second-order derivatives (smoothness) of a function based on the DT interpolation is more difficult than that of the first-order derivative (flatness). We propose a new method for calculating the smoothness with DT-based interpolation: the quadratic interpolators. The new method is tested through numerical experiments in the framework of full Bayesian inversion and applied to a gravity Bayesian inversion problem.