Improvement in the Determination of the Marine Geoid by Estimating the Bathymetry from Altimetry and Depth Soundings

• Published in: Marine geodesy 2005-03, Vol.28 (1), p.81-102
• Main Authors: VERGOS, GEORGIOS S., SIDERIS, MICHAEL G.
• Format: Article
• Language: English
• Subjects: bathymetry estimation method; gravity inversion; least-squares collocation; Marine; Marine geoid
• ISSN: 0149-0419; 1521-060X
• DOI: 10.1080/01490410590884593

The contribution of bathymetry to the estimation of gravity field related quantities is investigated in an extended test area in the Mediterranean Sea. The region is located southwest of the island of Crete, Greece, bounded between 33 ˆ ≤ ϕ ≤ 35 ˆ and 15 ˆ ≤ λ ≤ 25 ˆ . Gravity anomalies from the KMS99 gravity field and shipborne depth soundings are used with a priori statistical characteristics of depths in a least-squares collocation procedure to estimate a new bathymetry model. Two different global bathymetry models, namely JGP95E and Sandwell and Smith V8, are used to derive the depth a priori statistical information, while the estimated model is compared against both the global ones and the shipborne depth soundings to assess whether there is an improvement. Various marine geoid models are estimated using ERS1 and GEOSAT Geodetic Mission altimetry and shipborne gravity data. In that process, the effect of the bathymetry is computed using both the estimated and the original depths through a residual terrain modeling reduction. The TOPEX/Poseidon Sea Surface Heights, known for their high accuracy and precision, and the GEOMED solution for the geoid in the Mediterranean are used as control for the validation of the new geoid models and to assess the improvement that the estimated depths offer to geoid modeling. The results show that the newly estimated bathymetry agrees better (by about 30 to 300 m) with the shipborne depth soundings and provides smoother residual geoid heights and gravity anomalies (by about 8-20%) than those from global models. Finally, the achieved accuracy in geoid modeling ranges between 6 and 10 cm (1σ).