The Combined Mean Dynamic Topography MDT CNES-CLS18

The point of estimating a global Mean Dynamic Topography (MDT) is to reference the altimeter Sea Level Anomalies (SLA), computed relative to a 20 year (1993-2012) mean profile, in order to obtain the full signal that results from the dynamic of the ocean : the Absolute  Dynamic Topography (ADT).

Note that the Ssalto/Duacs (M)SLA products are computed relative to 1993-2012 period since April 2014 (Duacs 2014 v15.0 version and posterior version). Thus it is mandatory to be consistent and use (M)SLA and MDT referenced to the same time period. 

The MDT CNES-CLS calculation combines several steps:

  1. We first compute a first guess: a medium scale of Mean Dynamic Topography using only satellite data (gravimetry and altimetry). The raw difference of Mean Sea Surface with geoid model is filtered using an optimal filter to remove noise and residual geodetic signal.
  2. Beside we compute synthetic estimate of MDT and mean geostrophic current from oceanic in-situ data (T/S profiles and drifters).
    1. First, the data are processed to have the same physical content as altimetry:
      1. the drifters (SVP and Argo drifting at the surface) are processed to only extract the geostrophic component
      2. also deep baroclinic and barotropic components are added to dynamic height computed from T/S profiles.
    2. Then sea level anomalies are interpolated and subtracted from in-situ measurements of the ocean state to estimate a mean: MDT and the corresponding mean geostrophic surface velocities.
  3. Then, using a multivariate objective analysis, the synthetic estimates of mean heights and mean velocities are combined to improve the first guess of the MDT. 
  4. Finally, a validation is done by comparing this version to other MDTs and independent in-situ data.

The Table below lists the input data and improvement done for the new CNES-CLS18 compared with the previous CNES-CLS13.

MDT CNES-CLS18 MDT CNES-CLS13
First guess MSS used CNES-CLS15 (Pujol et al, 2018) CNES-CLS11 (Schaeffer et al, 2012)
Geoïd GOCO05S (Mayer-Gürr,et al. 2015)
  Complete GOCE mission (Nov 2009-October 2013)
+ 10.5 years of GRACE data
EGM-DIR-R4 (Bruinsma et al, 2012)
2 years of reprocessed GOCE data
+7 years of GRACE data
First Guess filtering Optimal filter (Rio et al, 2011) with updated parameters Optimal filter (Rio et al, 2011) 
Synthetic estimate from in-situ data Hydrological data CTD and ARGO Pref variable 200/400/900/1200/1900  from CORA4.2 (1993-2013), CORA5.0 (2014-2015) and CORA5.1 (2016)
Period 1993-2016
CTD (Cora3.4), ARGO Pref variable 200/400/900/1200/1900
Period 1993-2012
Drifter Data SD-DAC drifter, both drogued and undrogued: 1993-2016
Argo floats surface velocities: 1997-2016
SD-DAC drifter, both drogued and undrogued: 1993-2012
Argo floats surface velocities: 1997-2013
Ekman model Parameters fitted over the period 1993-2016 by latitude and Mixed Layer Depth (from ARMOR3D)
Two levels: 0m and 15m
Parameters fitted over the period 1993-2012, by longitude, latitude and month (Rio et al, 2014)
Two levels:  0m and 15m 
Wind Slippage correction Update of Rio et al, 2012 in order not to discard the trajectories beginning/end Rio et al, 2012
Drifter filtering Max (24 hours, Inertial Period) 3 days
Altimeter data Delayed-Time CMEMS-DUACS 2018 (Taburet et al, under review) Delayed-Time DUACS-2010 (Dibarboure et al, 2011)
Inversion Resolution Global 1/8° Global 1/4°