...
Easy Heading Macro | ||
---|---|---|
|
Preamble
History of modifications
Expand | ||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| ||||||||||||||||||||||||
|
List of datasets covered by this document
The present document applies to the C3S altimeter sea level Climate Data Record (CDR) and the following temporal extensions (Interim CDR). The current version corresponds to the reprocessed DUACS Delayed-Time vDT2021 products.
Related documents
Expand | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
| ||||||||||||
|
Acronyms
Expand | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
General definitions
Expand | ||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| ||||||||||||||
|
Scope of the document
This document is the Product User Guide and Specification (PUGS) document for Version DT2021 sea level products, distributed through the Copernicus Climate Change Service (C3S). This document has been developed within the frame of the 2018/C3S_312b_Lot3_CLS/SC2 contract. It provides the end user with practical information regarding the use of these products.
Executive summary
The current version of the C3S sea level product is the DUACS vDT2021 reprocessed delayed-time altimeter sea level products. It benefits from several advances over the previous version. They include:
- New L2P altimeter standards following expert recommendations (Lievin et al., 2020).
- Improved editing for L3 product for mapping
- More precise definition of the error budgets associated with the different altimeter measurements for the Optimal Interpolation process
This Product User Guide and Specification explains the basic altimetry principles that allow the computation of the altimeter sea level product, and provides a brief description of the associated production system. The details of the input data are provided, including their origin. The technical characteristics of each altimeter mission used in the production system are described, as well as the level 2 altimeter algorithms (geophysical standards and orbit solutions). An empirical correction of the TOPEX-A instrumental drift observed during 1993-1998 is included in the data files. The characteristics of the satellite constellation are described, and the principle of the sea level mapping procedure is provided. Finally, the product characteristics are described (format, nomenclature and data handling variables) and a description of the file content is provided in the Annex.
Product description
The sea level products distributed through the Copernicus Climate Change Service (C3S) comprise:
...
The present document refers to the C3S altimeter sea level Climate Data Record (CDR) and the following temporal extensions (Interim CDR). The current version corresponds to the reprocessed DUACS Delayed-Time vDT2021 product.
Common variables in Altimetry
Altimetry gives access to the Sea Surface Height (SSH) above the reference ellipsoid (see Figure 1) as in Eq (1):
LaTeX Formatting |
---|
Eq (1) $SSH = Orbit - Altimetric Range$ |
...
The reference period N considered can be changed as described in Pujol et al (2016).
Anchor | ||||
---|---|---|---|---|
|
Figure 1: Different concepts of sea surface height used in altimetry.
Variables in the daily sea level product
The variables disseminated to users as part of the C3S sea level product, are the Sea Level Anomalies (variable name in netcdf: 'sla'), Absolute Dynamic Topography (variable name in netcdf: 'adt'), formal mapping error from sla (variable name in netcdf: 'err_sla'), the geostrophic velocities anomalies (variable names in netcdf: 'ugosa' and 'vgosa'), the formal mapping error on zonal/meridional velocity anomalies (variable names in netcdf: 'err_ugosa' and 'err_vgosa') and the absolute geostrophic velocities (variable names in netcdf: 'ugos' and 'vgos').
Geostrophic current1 datasets are also disseminated to users as part of the C3S sea level product, and are generated from the SLAs and the ADTs. It is computed using a nine-point stencil width methodology (Arbic et al., 2012) for latitudes outside the 5°S/5°N band. In the equatorial band, the Lagerloef methodology (Lagerloef et al., 1999) is used. A specific variable is also available (variable name in netcdf: 'tpa_correction') and can be added to the SLA to correct for the observed instrumental drift during the lifetime of the TOPEX-A mission (the correction is null after this period). See Section 4.2.1 for more details.
A variable (variable name in netcdf: 'flag_ice') has been added to flag data using OSI SAF CDR sea ice concentration products (OSI-450) until 2016 and ICDR sea ice concentration (OSI-430-b) from 2016 (also distributed in the Climate Data Store, more info in Lavergne et al., 2019). The flag corresponds to the limit of 15% in sea ice concentration.
Info | ||
---|---|---|
| ||
1The geostrophic velocities are directly derived from the first derivative of the altimeter heights and correspond to the balance between the forces of pressure and those related to the Earth's rotation (the ageostrophic part of the velocities linked to acceleration is not included). |
Variables in the monthly sea level product
The variables disseminated to users are the Sea Level Anomalies (variable name in netcdf: ‘sla’) and the Eddy Kinetic Energy issued from the geostrophic velocities anomalies (variable name in netcdf: ‘eke’). It is computed using Eq (6):
...
LaTeX Formatting |
---|
Eq (6) $ EKE= (U*U+V*V)*1/2 $ |
Processing
The Delayed-Time DUACS component maintains a consistent and user-friendly altimeter database using state-of-the-art recommendations from the altimetry community.
...
- Data acquisition
- Input data quality control
- Intercalibrate and unify
- Along-track products generation
- Gridded merged products generation
- Final quality control
Input data and corrections
The altimeter measurements used to compute the C3S sea level product consist of Level-2 products from different missions called Delayed-Time Geophysical Data Records (GDR) or Non Time Critical (NTC) products. Details of the different L2 altimeter products sources and delay of availability are given in Table 1.
Anchor | ||||
---|---|---|---|---|
|
...
The auxiliary products (altimeter standards, geophysical corrections) used in the production are described in Table 2. They are the most up-to-date standards (whose timeliness is compatible with the C3S production planning) and most of them follow the recommendations of the ESA Sea Level CCI project (Quartly et al. 2017; Legeais et al., 2018). More details on the description of these standards can be found in Lievin et al., 2020.
Anchor | ||||
---|---|---|---|---|
|
Poseidon Topex | Jason-1 | OSTM/Jason-2 | Jason-3 | ERS-1 | ERS-2 | Envisat | Cryosat-2 | SARAL AltiKa | Sentinel-3A | Sentinel-6 MF | |
Orbit | GSFC STD18 | POE-E | POE-F | POE-F | Reaper | POE-E | POE-F | POE-F | POE-F | POE | |
Ionospheric Correction | Filtered dual-frequency altimeter range measurements [Guibbaud et al. 2015]; DORIS on Poseidon | Filtered dual-frequency altimeter range [Guibbaud et al. 2015] (from SSB C-band) | Filtered dual-frequency altimeter from [Guibbaud et al. 2015] & c> 170 from L2 GDRF | Reaper NIC09 model [Scharroo and Smith, 2010] | GIM [Ijima et al., 1999] | Filtered from L2; c>65: GIM [Ijima et al., 1999] corrected for 8mm bias | GIM [Ijima et al., 1999] | Filtered from L2 | L2 Filtered dual frequency | ||
Sea State Bias | Non parametric [Tran et al. 2010] ; BM4 on Poseidon | Non parametric [Tran 2015] | Non parametric [ Tran 2012] | Non parametric from J2 [ Tran 2012] & c>170 from [Tran 2020 report] J3 GDRF | BM3 [Gaspar and Ogor, 1994] | Non parametric [Mertz et al., 2005] | Non parametric [ Tran 2017] | Non parametric [ Tran 2018] Baseline C | Non parametric [ Tran 2018] | Non parametric [ Tran 2012] | Non parametric SSB [Tran 2021] from J3 GDR |
Wet Troposphere | GPD+ [Fernandes and Lazaro, 2015] | JMR (GDRE) radiometer | AMR radiometer | AMR radiometer (c>170 from L2 GDRF) | GPD+ [Fernandes and Lazaro, 2015] | MWR radiometer reprocessed | GPD+ [Fernandes and Lazaro, 2015] | Neural Network (5 entries) V4 | MWR 3 radiometer | MWR radiometer | |
Dry Troposphere | ERA5 (1-hour) model based | ||||||||||
Dynamical Atmospheric Correction | TUGO High frequencies forced with analysed ERA5 pressure and wind field + inverse barometer Low frequencies | TUGO HF forced with analysed ERA 5 pressure an d wind field; and after 02/2016 MOG2D HF forced with analysed ECMWF pressure and wind field + inverse barometer Low Frequencies | MOG2D HF forced with analysed ECM WF pressure and wind [Carrere and Lyard, 2003; operational version 3.2.0] + inverse barometer Low Frequencies | TUGO High frequencies forced with analysed ERA5 pressure and wind field + inverse barometer Low frequencies | TUGO High frequencies forced with analysed ERA5 pressure and wind field; and after 02/2016 MOG2D High frequencies forced with analysed ECMWF pressure and wind field + inverse barometer Low frequencies | TUGO HF forced with analysed ERA5 pressure an d wind field; and after 02/2016 MOG2D HF forced with analysed ECMWF pressure and wind field + inverse barometer Low Frequencies | MOG2D High frequencies forced with analysed ECMWF pressure and wind field [Carrere and Lyard, 2003; operational version 3.2.0] + inverse barometer Low frequencies | ||||
Ocean Tide | FES 2014 B [Carrère et al. 2016] | ||||||||||
Internal Tide | Zaron 2019 (HRETv8.1 tidal frequencies: M2, K1, S2, O1) | ||||||||||
Pole Tide | Desai et al., 2015 ; Mean Pole Location 2017 | ||||||||||
Solid Tide | Elastic response to tidal potential [Cartwright and Tayler, 1971; Cartwright and Edden, 1973] | ||||||||||
Mean Sea Surface | Composite (SIO,CNES/CLS15,DTU15) [Sandwell et al.,2017 ; Ole et al.; Pujol et al.,2018] | ||||||||||
Mean Dynamic Topography | CNES_CLS18 (Mulet et al, 2021) combined with CMEMS_2020 | ||||||||||
Glacial Isostatic Adjustment (GIA) | The DUACS L4 products are not corrected from GIA effects |
...
Warning | ||
---|---|---|
| ||
Between 1993 and 1998, the retrievals of global mean sea level (MSL) have been known to be affected by an instrumental drift in the TOPEX-A measurements, which has been quantified by several studies as discussed in the C3S Product Quality Assessment Report ([C3S_PQAR], section 3.2) and in Legeais et al. (2020). The altimeter sea level community agrees that it is necessary to correct the TOPEX-A record for the instrumental drift to improve the accuracy and the uncertainty of the total sea level record. An empirical correction of this drift based on a global comparison between altimetry and _in situ_ tide gauge measurements (WCRP sea level budget group, 2018) has been proposed in the data files. The correction value included in the dedicated variable can be added to the gridded SLA, to correct for the observed instrumental drift during the lifetime of the TOPEX-A mission (the correction is null after this period). This is a global correction to be added a posteriori (and not before) to the global mean sea level estimate derived from the gridded sea level data. It can be applied at regional or local scales as a best estimate (better than no correction, since the regional variation of the instrumental drift is unknown). However, even if the corrections proposed by the different studies available lead to similar global MSL trends and accelerations (in agreement with climate models), there is not yet a consensus on the best approach to estimate the drift correction at global and regional scales. The recommendation of the Ocean Surface Topography Science Team (OSTST) is to wait for the future release of a reprocessed TOPEX dataset. Therefore, the TOPEX-A correction has been proposed as a separate variable within the C3S sea level data files vDT2021 (and not directly included in the SLA estimate). See the sea level Product Quality Assessment Report [C3S_PQAR] for further details. |
Altimetry constellation
The complete altimetry satellite constellation used in the C3S sea level product is illustrated in Figure 2.
Anchor | ||||
---|---|---|---|---|
|
Figure 2: Overview of the L2P products (input for DUACS system) availability period for each altimetric mission.
The C3S sea level altimeter product is based on a satellite constellation with a stable number of altimeters in order to ensure the long-term stability of the ocean observation system. The different altimeter satellites included in the product are the reference missions and the complementary missions as well as missions of opportunity, as illustrated in Figure 3 and described below:
- the reference missions are the TOPEX/Poseidon, Jason-1, Jason-2, Jason-3 and Sentinel-6MF, which have been successively introduced into the production system. These missions are essential for the computation of the long- term trend of the MSL since they are used to wedge complementary missions in terms of sea level drift. Sentinel-6MF is the current reference mission used in the system and it has replaced Jason-3 in February 2022.
- the complementary missions provide additional information for the estimation of mesoscale signal variabilities (>200-300 km) and also increase the observing capacity at high latitudes, which is of great interest for climate. The missions that have successively been included in the C3S product are ERS-1, ERS-2, Envisat, SARAL/Altika and presently Sentinel-3A. Note that the ERS-1 mission was operated in an ice phase (phase D) from 21/12/1993 to 10/04/1994; no ERS-1 altimeter measurements have been used as input to the sea level production system during this period. As no other altimeter data are available, this means that the C3S product is based on TOPEX/Poseidon data only during this 3.5-month period. During the following two successive geodetic phases (phase E, 10/04/1994 – 28/09/1994 and phase F, 28/09/1994 – 21/03/1995), the changes to the ERS-1 mission operations (declared as a new mission: ERS-1 geodetic) have been taken into account in sea level data production.
- In addition, after the loss of the Envisat mission in April 2012, only the opportunity CryoSat-2 mission has been available. Thus, this opportunity mission was included in the C3S product until SARAL/AltiKa delayed-time measurements become available in March 2013.
...
Note |
---|
Note that the information about the satellites used to compute each map is given in the global attribute "platform" of each file. The use of such a constant number in the satellite constellation contributes to ensuring the long- term Mean Sea Level (MSL) stability, which is not the case when using all satellites available throughout the altimeter period (see section 3.1 of [C3S_PQAR]). |
Gridded merged product generation
The gridded merged product is based on the along-track altimeter measurements, which have undergone several processing steps, (as described in detail in [C3S_ATBD]). First of all, global and regional inter-mission biases are removed. Then, the along-track measurements are cross-calibrated following Le Traon and Ogor (1998), which allows for the reduction of the long wavelength errors (LWE) and also considers geographically-correlated errors. Along-track high frequency aliased signals are also removed. In addition, the data are filtered (Dufau et al., 2016) with 65km cut-off length low-pass filtering. The along-track measurements are also subsampled for the mapping procedure, keeping one along-track point out of two. All the details are described in section 3.6 of [C3S_ATBD], in Taburet et al. (2019) and Pujol et al (2016). These procedures ensure the long-term stability of the sea level record. An optimal interpolation method is used for the mapping procedure following Ducet et al. (2000) and Le Traon et al. (2003). This ensures mesoscale signal reconstruction. The parameters used for the mapping procedure are a compromise between the characteristics of the physical field to be focused on, and the sampling capabilities associated with the altimeter constellation.
Mean and reference period
The along-track and gridded sea surface heights (sea level anomalies and absolute dynamic topography) are computed with respect to a 20-year reference period (1993-2012). In addition to the reference period, a mean reference convention has been adopted in the DUACS products: the sea level time series has been arbitrarily referenced so that the mean sea level averaged during the year 1993 is set to zero (see Figure 4). This convention explains why the DUACS global mean SLA during the reference period (1993-2012) is different from zero. The obtained value (about 2.5cm without a Glacial Isostatic Adjustment (GIA) correction) is directly related to global sea level rise (see Figure 4, right). The most recently calculated Global Mean Sea Level (1993 to end of 2021) is displayed in Figure 5.
Note |
---|
Note that the proposed correction of the TOPEX-A instrumental drift has been chosen so that the correction is null after the end of the lifetime of the TOPEX mission in 1999 (dashed line in Figure 4 and in Figure 5). With this approach, the corrected Global MSL does not equal to zero in 1993. This approach is the preferred approach to ensure the continuity of the initial and corrected GMSL after 1999 (e.g. for ocean modellers). |
...
Figure 5: Global mean sea level progression for the period 1993-2021 (without GIA correction) deduced from DUACS L4 gridded products. The dashed line represents the sea level time-series corrected for the TOPEX-A instrumental drift (between 1993-1998).
Specifications and target requirements
Spatial and temporal coverage
The daily time series begins on 01/01/1993. The time series benefits from regular temporal extensions approximately 3 times per year (ICDR production), and the timeliness of the product is of 5 months at the minimum. Such a delay depends on:
...
The characteristics of the different missions used in the C3S sea level product are described in Table 3.
Anchor | ||||
---|---|---|---|---|
|
...
Info | ||
---|---|---|
| ||
* ERS-1: No ERS-1 data between 23 December 23,1993 and April 10, 1994 (ERS-1 phase D - 2nd ice phase). |
Validation and uncertainty estimates
Validation activities are carried out to assess the quality of the product. The validation method is described in the Product Quality Assurance Document [C3S_PQAD] and details of the validation results are provided in the Product Quality Assessment Report [C3S_PQAR].
The description of the altimeter errors and characterization of the uncertainties are available in [C3S_PQAR].
Data usage information
Grid characteristics
The product is delivered in a Cartesian grid with the coverage definition detailed in the table below:
...
Note |
---|
Note that the values taken into account to generate a map are ocean values. The mapping process (see section 4.2.3) computes some slight extrapolation to the coasts. This avoids the production of gaps in the data that can occur near the coast, and it also allows for a more precise computation of velocities. |
Format
The product is stored and delivered to users using the NetCDF (Network Common Data Form) using CF (Climate and Forecast) Metadata convention.
File nomenclature
Daily product
The nomenclature of the file is the following:
...
Code Block |
---|
<DateMap>=the date of the map in the form YYYYMMDD |
Monthly product
The nomenclature of the file is the following:
...
Code Block |
---|
<MonthMap>=the month of the map in the form YYYYM |
Data Handling Variables
Daily sea level product
4 dimensions are defined:
...
The variables are listed in Table 5:
Anchor | ||||
---|---|---|---|---|
|
Type | Name | Content | Unit | Scale Factor |
float | time(time) | Time of measurement | days since 1950-01-01 | none |
float | latitude(latitude) | Latitude of measurement | degrees_north | none |
float | longitude(longitude) | Longitude of measurement | degrees_east | none |
float | lat_bnds (latitude,nv) | latitude values at the north and south bounds of each pixel. | degrees_north | none |
float | lon_bnds(longitude,nv) | longitude values at the west and east bounds of each pixel. | degrees_east | none |
int | nv(nv) | Useful for grid definition | none | none |
int | crs | Describes the grid_mapping used by the data in this file. This variable does not contain any data; only information about the geographic coordinates system. | none | none |
int | sla(time,latitude,longitude) | Sea level Anomaly | meters | 10-4 |
int | err_sla(time,latitude,longitude) | Formal mapping error | meters | 10-4 |
int | ugosa(time,latitude,longitude) | Geostrophic velocity anomalies: eastward zonal component | m/s | 10-4 |
int | vgosa(time,latitude,longitude) | Geostrophic velocity anomalies: northward meridian component | m/s | 10-4 |
int | err_ugosa(time,latitude,longitude) | Formal mapping error on zonal velocity anomalies | m/s | 10-4 |
int | err_vgosa(time,latitude,longitude) | Formal mapping error on meridional velocity anomalies | m/s | 10-4 |
int | adt(time,latitude,longitude) | Absolute dynamic topography | meters | 10-4 |
int | ugos(time,latitude,longitude) | Absolute geostrophic velocity: eastward zonal component | m/s | 10-4 |
int | vgos(time,latitude,longitude) | Absolute geostrophic velocity: northward meridian component | m/s | 10-4 |
int | tpa_correction | TOPEX-A instrumental drift correction derived from altimetry and tide gauges global comparisons | m | 10-4 |
int | flag_ice | Ice Flag based on CDR OSI-SAF products until 2016 (OSI-450), Interim products from 2016 (OSI- 430-b) | - | 10-4 |
Monthly sea level product
4 dimensions are defined:
...
The variables are listed in Table 6:
Anchor | ||||
---|---|---|---|---|
|
Type | Name | Content | Unit | Scale Factor |
float | time(time) | Time of measurement | days since 1950-01-01 | none |
float | latitude(latitude) | Latitude of measurement | degrees_north | none |
float | longitude(longitude) | Longitude of measurement | degrees_east | none |
float | lat_bnds (latitude,nv) | latitude values at the north and south bounds of each pixel. | degrees_north | none |
float | lon_bnds(longitude,nv) | longitude values at the west and east bounds of each pixel. | degrees_east | none |
float | climatology_bnds(time,nv) | Useful for grid definition | meters | |
int | nv(nv) | Useful for grid definition | none | none |
int | crs | Describes the grid_mapping used by the data in this file. This variable does not contain any data; only information about the geographic coordinates system. | none | none |
int | sla(time,latitude,longitude) | Sea level Anomaly | meters | 10-4 |
int | eke(time,latitude,longitude) | Eddy Kinetic Energy | cm2/s2 | 10-4 |
Appendix A - Specifications of the daily sea level product
Expand | ||
---|---|---|
| ||
|
Appendix B - Specifications of the monthly sea level product
Expand | ||
---|---|---|
| ||
|
Anchor references references
References
references | |
references |
Arbic, B. K., Scott, R. B., Chelton, D. B., Richman, J.G., and Shriver, J. F.: Effects on stencil width on surface ocean geostrophic velocity and vorticity estimation from gridded satellite altimeter data, J. Geophys. Res., 117, C03029, doi:10.1029/2011JC007367, 2012.
...
Info | ||
---|---|---|
| ||
This document has been produced in the context of the Copernicus Climate Change Service (C3S). The activities leading to these results have been contracted by the European Centre for Medium-Range Weather Forecasts, operator of C3S on behalf of the European Union (Delegation Agreement signed on 11/11/2014 and Contribution Agreement signed on 22/07/2021). All information in this document is provided "as is" and no guarantee or warranty is given that the information is fit for any particular purpose. The users thereof use the information at their sole risk and liability. For the avoidance of all doubt , the European Commission and the European Centre for Medium - Range Weather Forecasts have no liability in respect of this document, which is merely representing the author's view. |
icus Climate Change Service
Related articles
Content by Label | ||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
|