An Assessment of Recently Released High-Degree Global Geopotential Models Based on Heterogeneous Geodetic and Ocean Data

The development of the global geopotential model (GGM) broadens its applications in ocean science, which emphasizes the importance for model assessment. We assess the recently released high-degree GGMs over the South China Sea through heterogeneous geodetic observations and synthetic/ocean reanalysis data. The comparisons with a high resolution (∼3 km) airborne gravimetric survey over the Paracel Islands show that XGM2019e_2159 has relatively high quality, where the standard deviation (SD) of the misfits against the airborne gravity data is ∼3.1 mGal. However, the comparisons with local airborne/shipborne gravity data hardly discriminate the qualities of other GGMs that have or truncated to the same expansion degree. Whereas, the comparisons with the synthetic/ocean reanalysis data demonstrate that the qualities of the values derived from different GGMs are not identical, and the ones derived from XGM2019e_2159 have better performances. The SD of the misfits between the mean dynamic topography (MDT) derived from XGM2019e_2159 and the ocean data is 2.5 cm; and this value changes to 7.1 cm/s (6.8 cm/s) when the associated zonal (meridian) geostrophic velocities are assessed. In contrast, the values derived from the other GGMs show deteriorated qualities compared to those derived from XGM2019e_2159. In particular, the contents computed from the widely used EGM2008 have relatively poor qualities, which is reduced by 3.9 cm when the MDT is assessed, and by 4.0 cm/s (5.5 cm/s) when the zonal (meridian) velocities are assessed, compared to the results derived from XGM2019e_2159. The results suggest that the choice of a GGM in oceanographic study is crucial, especially over coastal zones. Moreover, the synthetic/ocean data sets may be served as additional data sources for global/regional gravity field assessment, which are useful in regions that lack of high-quality geodetic data.

Controlling atmospheric forcing parameters of global ocean models: sequential assimilation of sea surface Mercator-Ocean reanalysis data

In the context of stand alone ocean models, the atmospheric forcing is generally computed using atmospheric parameters that are derived from atmospheric reanalysis data and/or satellite products. With such a forcing, the sea surface temperature that is simulated by the ocean model is usually significantly less accurate than the synoptic maps that can be obtained from the satellite observations. This not only penalizes the realism of the ocean long-term simulations, but also the accuracy of the reanalyses or the usefulness of the short-term operational forecasts (which are key GODAE and MERSEA objectives). In order to improve the situation, partly resulting from inaccuracies in the atmospheric forcing parameters, the purpose of this paper is to investigate a way of further adjusting the state of the atmosphere (within appropriate error bars), so that an explicit ocean model can produce a sea surface temperature that better fits the available observations. This is done by performing idealized assimilation experiments in which Mercator-Ocean reanalysis data are considered as a reference simulation describing the true state of the ocean. Synthetic observation datasets for sea surface temperature and salinity are extracted from the reanalysis to be assimilated in a low resolution global ocean model. The results of these experiments show that it is possible to compute piecewise constant parameter corrections, with predefined amplitude limitations, so that long-term free model simulations become much closer to the reanalysis data, with misfit variance typically divided by a factor 3. These results are obtained by applying a Monte Carlo method to simulate the joint parameter/state prior probability distribution. A truncated Gaussian assumption is used to avoid the most extreme and non-physical parameter corrections. The general lesson of our experiments is indeed that a careful specification of the prior information on the parameters and on their associated uncertainties is a key element in the computation of realistic parameter estimates, especially if the system is affected by other potential sources of model errors.

Controlling atmospheric forcing parameters of global ocean models: sequential assimilation of sea surface Mercator-Ocean reanalysis data

In the context of stand alone ocean models, the atmospheric forcing is generally computed using atmospheric parameters that are derived from atmospheric reanalysis data and/or satellite products. With such a forcing, the sea surface temperature that is simulated by the ocean model is usually significantly less accurate than the synoptic maps that can be obtained from the satellite observations. This not only penalizes the realism of the ocean long-term simulations, but also the accuracy of the reanalyses or the usefulness of the short-term operational forecasts (which are key GODAE and MERSEA objectives). In order to improve the situation, partly resulting from inaccuracies in the atmospheric forcing parameters, the purpose of this paper is to investigate a way of further adjusting the state of the atmosphere (within appropriate error bars), so that an explicit ocean model can produce a sea surface temperature that better fits the available observations. This is done by performing idealized assimilation experiments in which Mercator-Ocean reanalysis data are considered as a reference simulation describing the true state of the ocean. Synthetic observation datasets for sea surface temperature and salinity are extracted from the reanalysis to be assimilated in a low resolution global ocean model. The results of these experiments show that it is possible to compute piecewise constant parameter corrections, with predefined amplitude limitations, so that long-term free model simulations become much closer to the reanalysis data, with misfit variance typically divided by a factor 3. These results are obtained by applying a Monte Carlo method to simulate the joint parameter/state prior probability distribution. A truncated Gaussian assumption is used to avoid the most extreme and non-physical parameter corrections. The general lesson of our experiments is indeed that a careful specification of the prior information on the parameters and on their associated uncertainties is a key element in the computation of realistic parameter estimates, especially if the system is affected by other potential sources of model errors.