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

Abstract

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σ).

Keywords:

• Marine geoid
• bathymetry estimation method
• gravity inversion
• least-squares collocation

Acknowledgments

The financial assistance provided to the authors by NSERC and the GEOIDE Network of Centers of Excellence is gratefully acknowledged. The bathymetry estimation was based on a program developed and provided by Per Knudsen; for this, he is gratefully acknowledged. We extensively used the Generic Mapping Tools Version 3.1 (Wessel and Smith 1998) in displaying our results and the GRAVSOFT software (Tscherning et al. 1992) package for the RTM reductions. BGI and the IfE are gratefully acknowledged for providing the shipborne gravity data.