Contributors: Lin Gilbert (University of Leeds), Sebastian Simonsen (Technical University of Denmark), Jan Wuite (ENVEO IT GmbH)
Issued by: University of Leeds / Lin Gilbert
Date: 04/05/2021
Ref: C3S_312b_Lot4.D3.IS.9-v3.0_SEC_Product_User_Guide_Specification_i1.0.docx
Official reference number service contract: 2018/C3S_312b_Lot4_EODC/SC2
History of modifications
List of datasets covered by this document
Related documents
Acronyms
Scope of the document
This document is the Product User Guide and Specification for ice sheet surface elevation change (SEC), part of the Copernicus Ice Sheets and Ice Shelves service. It describes the datasets produced, and their specifications.
Executive summary
The service addresses three essential climate variables (ECVs) by providing four separate products.
- Ice velocity is given for Greenland in product D3.IS.4
- Gravimetric mass balance is given for Greenland and Antarctica in product D3.IS.5
- Surface elevation change is given for:
- Antarctica in product D3.IS.6.1
- Greenland in product D3.IS.6.2
We document here the description of each of the CDR v3.0 for the two Polar region SEC products with a guide to their usage and specifications.
Product Change Log
The same product changes are applicable to both Antarctic and Greenland products.
Version | Product Changes |
---|---|
V1 | Initial product, using data from ERS1, ERS2, EnviSat and CryoSat-2. |
V2 | Data from Sentinel-3A added. |
V3 | Data from Sentinel-3B added. |
1. Surface elevation change, Antarctica – D3.IS.6.1
1.1. Product description
The product contains surface elevation change rates and their uncertainties from the Antarctic ice sheet, ice shelves, ice rises and islands on a regular geographic grid at regular time intervals. The change rate is calculated over a 5-year period. An example plot from one data grid is shown in Figure 1 below.
Figure 1: Example data from Antarctic surface elevation change product D3.IS.6.1
The initial CDR input data is radar altimetry surface elevation measurements from ERS1, ERS2, Envisat, CryoSat-2, Sentinel-3A and Sentinel-3B.
The product is updated monthly and will always contain the whole-time range to date. There is a time lag of 2 months between new measurements being taken and their being added to the product.
Validation is performed with respect to observation campaigns by the Airborne Topographic Mapper, a scanning laser altimeter flown on board aircraft by Operation IceBridge (Studinger 2014). Each campaign lasts for one Antarctic spring season (October to December), and so validation is performed annually.
Specification summary:
- Sensors: ERS1 RA, ERS2 RA, Envisat RA-2, CryoSat-2 SIRAL, and Sentinel-3A/B SRAL
- Grid: 25km by 25km polar stereographic projection based on 0°E, 71°S on the WGS84 ellipsoid (EPSG: 3031)
- Time range: 1992 to present
- Time intervals: 5yr window moving in monthly steps
- Validation frequency: annual
- Quality flagging: flags provided for steep terrain and missing data
Details of the product file format and contents are below. Details of methodology may be found in the related document, the Algorithm Theoretical Basis Document.
1.2. Target requirements
Targets are set by two separate bodies. The Global Climate Observing System (GCOS) maintains definitions of Essential Climate Variables and their requirements at https://gcos.wmo.int/en/essential-climate-variables/ice-sheets-ice-shelves/ecv-requirements. The C3S project itself provides key performance indicator targets.
Details of the CDR v3 performances against these targets may be found in the Product Quality Assessment Report, a related document [D2].
The requirements summarized for the SEC products, based on the combined GCOS and C3S performance indicator targets are shown in Table 1.
Table 1: SEC products targets and performance indicators
Statistic | Target | Target source |
---|---|---|
Stability at pixel-level | 0.1 m/y | GCOS |
Accuracy at basin-level | 0.1 m/y | GCOS |
Proportion of area within 0.1 m/y accuracy at pixel-level | No target stated | C3S project |
Surface coverage, aggregated over one year | 65% ERS1, ERS2, Envisat, Sentinel-3A/B | C3S project |
1.3. Data usage information
1.3.1. Product data format and content
The product is provided as a NetCDF file containing stacked grids of the surface elevation change rate and associated uncertainty, and validity flags. The grids cover the Antarctic ice sheets, shelves, rises and islands at 25km x 25km resolution in a polar stereographic projection with centre longitude 0E and true scale latitude 71S. There is one grid per month, each cell containing the rate of surface elevation change derived from a 5 year period centred on that gird's timestamp. The change rate and its uncertainty are given in m/year. Missing data is indicated by a floating point NaN (not a number) value. Single-layer flag grids are provided for surface type (eg ocean, land ice, ice shelf, ice rise or island) and high slope (< 2⁰, 2⁰ <= slope <= 5⁰, > 5⁰), on the same projection as the data. The main variables are listed in Table 2 below.
Table 2: Antarctic surface elevation change main data variables.
Variable name | Variable description | Type |
---|---|---|
x | Centre of grid cell on X axis, in m | 32-bit float |
y | Centre of grid cell on Y axis, in m | 32-bit float |
longitude | Longitude of grid cell centre, in degrees east | 32-bit float |
latitude | Latitude of grid cell centre, in degrees north | 32-bit float |
time | Central time of surface elevation change rate derivation, in hours since 1990.0 | 32-bit float |
sec | Surface elevation change rate, in m/year | 32-bit float |
sec_uncert | Uncertainty on surface elevation change rate, in m/yr | 32-bit float |
sec_ok | Validity flag for surface elevation change rate | Byte |
surface_type | Flag for geographical surface type in cell | Byte |
high_slope | Flag for geographical slope class (ie low/medium/high) in cell | Byte |
1.3.2. Product known limitations
Please note: the gridded data and its uncertainties are well-defined at pixel level. If they are to be combined to drainage basin-level the user should be sure to consider the varying terrain and ice dynamics within the basin - a simple mean value will not be representative.
1.3.3. Product nomenclature
The dataset filename as delivered to EODC is of the format
C3S_AntIS_RA_SEC_versx_yyyy-mm-dd.nc, where:
- C3S is the overall project
- AntIS indicates the Antarctic ice sheet region
- RA indicates that the data source is radar altimetry
- SEC is the ECV addressed, ie surface elevation change
- versx is the version number of the processing code
- yyyy-mm-dd is the creation date of the dataset
The latest file created should always be used, as the files are accumulative – each one contains all previous data as well as its monthly updates.
1.3.4. Structure of netCDF files
The header data for an example netCDF data file is given here.
Table 3: Sample of the structure of the provided NetCDF file for the Antarctic surface elevation change and exemplified by the November 2020 test dataset
netcdf C3S_AntIS_RA_SEC_vers3_2020-11-30 { :institution = "Copernicus Climate Change Service" ; :contact = "copernicus-support@ecmwf.int" ; :project = "C3S_312b_Lot4_ice_sheets_and_shelves" ; :creation_date = "2020-10-27T:14:59:21Z" ; :comment = "Data is geophysically corrected, instruments 60 month power corrected (ERS1/2, EnviSat, CryoSat-2) or 36 month (Sentinel-3A) or 12 month (Sentinel-3B). SEC uses crossover method, cross-calibration by elevation regression. Time coverage 1994.83 to 2018.17 using data from 2.5 years further at each end of the range. Longitude 0 to 360 degrees, latitude -90 to -57.664 degrees. Grid_projection ESPG: 3031, ie PS 0E 71S WGS84, grid bottom left at -2.6e6m in x and -2.2e6m in y, grid cell width in x and y 25km" ; :history = "Product version 3.0" ; :summary = "Surface elevation change rate derived for Antarctica in 25km by 25km grid cells over a 5 year window moving at a monthly cadence." ; :keywords = "EARTH SCIENCE CLIMATE INDICATORS CRYOSPHERIC INDICATORS GLACIAL MEASUREMENTS GLACIER ELEVATION/ICE SHEET ELEVATION, EARTH SCIENCE\tCRYOSPHERE GLACIERS/ICE SHEETS\tGLACIER ELEVATION/ICE SHEET ELEVATION" ; :license = "C3S general license" ; } |
2. Surface elevation change, Greenland – D3.IS.6.2
2.1. Product description
The product contains surface elevation change rates and their associated uncertainties for the Greenland ice sheet and is provided on a regular grid at monthly resolution. The basis for the elevation change estimate for the older satellites (ERS-1, ERS-2 and Envisat) is a running 5-year mean, whereas for the ongoing satellite missions (Cryosat-2, Sentinel-3 A and B) the elevation change estimate is based on the monthly evaluation of a 3-year baseline. An example plot from one data grid is shown in Figure 2.
Figure 2: Example accumulated surface elevation change map produced from product D3.IS.6.2
The current version 3 of the CDR product builds on the legacy of its preceding versions:
CDRv1. The initial product release builds on the R&D efforts of the Greenland ice sheet CCI project's surface elevation measurements from ERS1, ERS2, Envisat and CryoSat-2.
CDRv2. Adds the use of observations from Sentinel-3A.
CDRv3. The time series have been reprocessed, with the upgraded baselines for Envisat and CryoSat-2 data and with the inclusion of data from Sentinel-3B.
Like in the CDRv1 and v2 products, the current CDRv3 will be updated monthly as iCDRs, with a 2-month time lag between the acquisition of new measurements and their addition to the product.
Validation is performed with respect to airborne measurements collected with the Airborne Topographic Mapper, a scanning laser altimeter, during Operation IceBridge campaigns (Studinger 2014). Validation is performed annually since Arctic campaigns take place once a year during the spring season (March-May).
Specification summary:
- Sensors: ERS1 RA, ERS2 RA, Envisat RA-2 CryoSat-2 SIRAL, Sentinel-3A SRAL and Sentinel-3B SRAL
- Grid: 25km by 25km polar stereographic projection (EPSG 3413)
- Time range: 1992 to present
- Time intervals: 5-year window moving in monthly steps (3-year in the CryoSat-2 and Sentinel-3 era)
- Validation frequency: Annual
- Quality flagging: Flags are provided for steep terrain, missing data and nearest distance to the original RA-altimeter measurement.
Details of the product file format and contents are provided in the following section. Details of methodology may be found in the related document, the Algorithm Theoretical Basis Document [D1].
2.2. Target requirements
The target requirements are generic for surface elevation change measurements. Hence, the same requirements listed for the Antarctic surface elevation change in Table 1 apply to Greenland too.
2.3. Data usage information
2.3.1. Product data format and content
The product is provided as a NetCDF file containing stacked grids of the surface elevation change rate and associated uncertainty, and validity flags. Given in a north polar stereographic projection, with centre longitude 45W and latitude 70N, the grids cover the Greenland ice sheet at 25km x 25km. Solutions for the Greenland surface elevation change are given at a temporal resolution of one month, which combined with the 25 km grid gives a parameter-dimension of (65x123x340)1. The change rate and its uncertainty are given in m/year. Missing data is indicated by a floating-point NaN (not a number) value. Single-layer flag grids are provided for surface type (0 for Land/Ocean and 1 for an ice cover of more than 95%) and high slope (< 2⁰, 2⁰ <= slope <= 5⁰, > 5⁰). All single-layer flags are gridded using the same north polar stereographic projection as the data. The main variables are listed in Table 4.
Table 4: Greenland surface elevation change main data variables.
Variable name | Variable description | Type |
---|---|---|
x | Centre of grid cell on X axis, in m | 32-bit float |
y | Centre of grid cell on Y axis, in m | 32-bit float |
lon | Longitude of grid cell centre, in degrees east | 32-bit float |
lat | Latitude of grid cell centre, in degrees north | 32-bit float |
time | Central time of surface elevation change rate derivation, in hours since 1990.0 | 32-bit float |
dhdt | Surface elevation change rate, in m/year | 32-bit float |
dhdt_uncert | Uncertainty on surface elevation change rate, in m/yr | 32-bit float |
dhdt_ok | Validity flag for surface elevation change rate | Byte |
surface_type | Flag for geographical surface type in cell | Byte |
high_slope | Flag for geographical slope class (i.e. low/medium/high) in cell | Byte |
dist | Distance to the nearest observational node, in m | int |
2.3.2. Product known limitations
Please note: The gridded data and its uncertainties are well-defined at the pixel level (25x25 km). If they are to be combined to drainage basin-level the user should be sure to consider the varying terrain and ice dynamics within the basin - a simple mean value will not be representative. The performed kriging procedure has the capability of extrapolating data over undesired distances and the distance-flag should be consulted before any averaging of elevation change is performed.
2.3.3. Product nomenclature
The dataset filename as delivered to EODC is of the format
C3S_GrIS_RA_SEC_ggkm_vv_yyyy-mm-dd.nc, where:
- C3S is the overall project
- GrIS indicates the Greenland ice sheet region
- RA indicates that the data source is radar altimetry
- SEC is the ECV addressed, i.e. surface elevation change
- gg the grid size in km
- vv is the version number of the processing code
- yyyy-mm-dd is the creation date of the dataset
The filename is changed within the CDS to fit their information storage system.
The latest file created should always be used, as the files are cumulative – each one contains all previous data as well as its monthly updates.
2.3.4. Structure of NetCDF files
The header data for an example NetCDF data file is given Table 5.
Table 5: Example of the NetCDF file structure for the Greenland Ice Sheet Surface elevation change based on the November 2020 version 3 test dataset.
netcdf C3S_GrIS_RA_SEC_25km_3.0-test_2020-10-27.nc { |
3. Data access information
3.1. All products
Data will be made available through the Copernicus Climate Data Store (CDS), which is the sole data distributor. Registration (free) is required to access the CDS and its toolbox software suite. The CDS is a web-based service, with its homepage at https://cds.climate.copernicus.eu/cdsapp#!/home
Data can be downloaded from the website and used under the License to Use Copernicus Products (included on download page). Data may also be viewed online.
All requests for information or further data should be channelled through the CDS Knowledge Base at https://confluence.ecmwf.int//display/CKB/
References
Studinger, M. (2014). IceBridge ATM L4 Surface Elevation Rate of Change, Version 299 1, Antarctica subset. N. S. a. I. D. C. D. A. A. Center. Boulder, Colorado, USA. DOI: 10.5067/BCW6CI3TXOCY
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 (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.