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: 27/04/2021
Ref: C3S_312b_Lot4.D3.IS.7-v3.0_IV_Product_User_Guide_Specification_i1.0
Official reference number service contract: 2018/C3S_312b_Lot4_EODC/SC2
History of modifications
Issue | Date | Description of modification | Chapters / Sections |
i0.1 | 13/01/2021 | The present document is a subsection/update for Ice Velocity only, based on C3S_312b_Lot4.D3.IS.7-v2.0_202001_Product_User_Guide_i1.2. Minor revisions section 1.1. Minor revisions in section 1.2 to Table 1 and Table 2. Revision to section 1.3.3 to reflect temporal extension for CDR v3.0. Provision of updated metadata examples for CDR v3.0 in Table 4 and Table 5. Revised modification table. Updated all fields, accepted all formatting changes, updated ToC. Introduced product change log. | JW |
I0.4 | 20/04/2021 | Accepted all reviewed changes. |
List of datasets covered by this document
Deliverable ID | Product title | Product type (CDR, ICDR) | Version number | Delivery date |
D3.IS.4-v3.0 | Ice velocity | CDR | 3.0 | 31/01/2021 |
Related documents
Reference ID | Document |
D1 | Algorithm Theoretical Basis Document D1.IS.2-v3.0 |
D2 | Product Quality Assessment Report D2.IS.2-v3.0 |
Acronyms
Acronym | Definition |
AIS | Antarctic Ice Sheet |
ALOS | Advanced Land Observation Satellite |
ATBD | Algorithm Theoretical Basis Document |
CCI | Climate Change Initiative |
CDR | Climate Data Record |
CDS | Climate Data Store |
DEM | Digital Elevation Model |
ECV | Essential Climate Variable |
EODC | Earth Observation Data Centre |
EPSG | European Petroleum Survey Group map projection database |
ERS | European Remote-sensing Satellite |
ESP | ENVEO Software Package |
GCOS | Global Climate Observing System |
GMB | Gravimetric Mass Balance |
GPS | Global Positioning System |
GrIS | Greenland Ice Sheet |
IV | Ice Velocity |
PALSAR | Phased Array type L-band Synthetic Aperture Radar |
RA | Radar Altimeter |
SAR | Synthetic Aperture Radar |
SRAL | Sentinel-3 ku/c Radar Altimeter |
S1 | Sentinel-1 |
TSX | TerraSAR-X |
Scope of the document
This document is the Product User Guide and Specification for Ice Velocity (IV) as 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 provide here the description of the Ice Velocity CDR v3 with a guide to its usage and specifications.
Product Change Log
Version | Product Changes |
---|---|
v202101 | Release of data set CDR v1.3. Product algorithm updated from v202001, providing maps at higher resolution (250m instead of 500m) and uncertainty/error isprovided component wise for easting and northing. CDR v1.2 product extended to 2020-09-31, including also reprocessed CDR v1.2 dataset at 250m. |
v202001 | Release of data set CDR v1.2. Product algorithm same as v201901. CDR v1.1 product extended to include map for 2018-10-01 to 2019-09-31, including also reprocessed CDR v1.1 dataset with adjusted end month: 2018-09-31 instead of 2018-10-31. |
v201901 | First release of the dataset: CDR v1.1. Ice velocity map produced for 2017-10-01 to 2018-10-31. |
1. Greenland ice sheet velocity
1.1. Product description
The product contains annually averaged ice velocity (IV) for Greenland derived from Sentinel-1 (S1) SAR data (Figure 1). The surface velocity is derived by applying advanced iterative feature tracking techniques. The primary processor for IV generation is the ENVEO software package (ESP v2.1), see the related ATBD [D1]. The system’s core module performs coherent and incoherent offset tracking utilizing long stripes of S1 data acquired in interferometric wide (IW) swath mode. The ESP v2.1 is a state-of-the-art IV retrieval algorithm designed for SAR sensors (e.g. Sentinel-1, TerraSAR-X, ERS-1, ERS-2, ALOS PALSAR, Cosmo-SkyMed), and has been tested rigorously through intercomparisons with other packages and extensive validation efforts.
Feature tracking is a technique capable of acquiring ice flow velocity data over short (days) and longer time spans (years) and in regions with fast flow, as no coherence is required. The method uses the displacement of surface features such as crevasses or rifts and edges that move approximately with the same speed as the ice and are identifiable on two co-registered satellite images to derive velocity. The ice velocity product is derived using feature tracking on repeat pass SAR images. The velocity grid for a given file represents the average ice surface velocity over the repeat pass period. For Sentinel-1 acquisitions in Greenland the repeat pass period is 6 to 12 days. The C3S Greenland Ice Sheet Velocity combines all ice velocity maps acquired over a full year in an annually averaged product provided at 250 m grid spacing. The averaging year runs from Oct 1st until September 30th, roughly mimicking a glaciological mass balance year.
Figure 1: C3S ice velocity map of the Greenland Ice Sheet based on Sentinel-1 data acquired from October 2019 to September 2020.
1.2. Target requirements
Table 1 shows the C3S IV product achievement versus the user requirements identified through an extensive user survey within the glaciological community as part of the Ice Sheets CCI project (Hvidberg, et al, 2012). The C3S product fulfills all minimum requirements for spatial and temporal resolution and falls well within the optimum accuracy range.
Table 1: C3S IV product achievement versus Ice Sheets CCI user requirements ( *)
Requirement | C3S GrIS IV Product | Minimum Req ( *) | Optimum Req ( *) |
---|---|---|---|
Spatial resolution | 250m | 100m-1km | 50m-100m |
Temporal resolution | Annual | annual | monthly |
Accuracy | 1 cm/d | 8-27 cm/d | 3-8 cm/d |
What times are observations needed | All year | All year | - |
1.3. Data usage information
1.3.1. Product Data Format and Content
The product is provided as a NetCDF file with separate layers for the velocity components: vx, vy, vz and vv (magnitude of the horizontal components), and maps showing the valid pixel count and uncertainty (std) (Table 2). The ice velocity map is annually averaged and provided at 250m grid spacing in North Polar Stereographic projection (EPSG: 3413, defined at https://epsg.io/3413). The horizontal velocity is provided in true meters per day, towards easting (vx) and northing (vy) direction of the grid, and the vertical displacement (vz), is derived from a digital elevation model (TanDEM-X 90m DEM; Rizzoli et al., 2017). For all maps a nodata value of 3.4028235e+38 is used.
Table 2: IV main data variables
Variable name | Variable description | Type |
---|---|---|
land_ice_surface_easting_velocity | Ice velocity East component [m/day] | 32-bit floating-point |
land_ice_surface_northing_velocity | Ice velocity North component [m/day] | 32-bit floating-point |
land_ice_surface_vertical_velocity | Ice velocity Vertical component [m/day] | 32-bit floating-point |
land_ice_surface_velocity_magnitude | Ice velocity magnitude [m/day] | 32-bit floating-point |
land_ice_surface_measurement_count | Valid pixel count [#] | 32-bit integer |
land_ice_surface_easting_stddev | Standard deviation [m/day] (Easting) | 32-bit floating-point |
land_ice_surface_northing_stddev | Standard deviation [m/day] (Northing) | 32-bit floating-point |
1.3.2. Product Known Limitations
The following lists some known product limitations:
1) The IV products contain separate layers for the horizontal (Easting, Northing) and the vertical components of velocity. This is, however, not the true 3D velocity, which requires both ascending and descending image pairs acquired close in time. The vertical component is derived from the difference in height of start and end position of the displacement vector taken from a DEM.
2) The IV products do not have a time stamp for a single date, but give the average velocity over the time-period covered.
3) For various reasons, the tracking software sometimes fails to find matching features leading to gaps in the velocity fields. This can be caused by a lack of surface features or when features, for example crevasses, rapidly change due to shearing leading to low correlation. Other reasons for gaps in the IV maps can be areas affected by radar shadow or anomalous pixels that are filtered out. A simple distance-weighted averaging filter is applied to get rid of outliers and to fill small gaps in the data (<5 pixels), further filtering/gap filling is left to the user if required. The annual map has only a few gaps.
4) Due to different acquisition modes, sensor type, resolution and processing strategy there can be differences between S1 IV products and IV products derived from other sensors that complicate a direct comparison between the data sets. Because of differences in resolution, the image patches used for feature tracking have different dimensions impacting the type of features that can be resolved. S1 can for instance not capture the high velocity gradients that may be found in shear zones with the same detail as for example TerraSAR-X (TSX). On the other hand, due to the regular repeat acquisition the temporal sequence of S-1 is much higher than that of TSX and the covered area of the IV maps is much larger.
5) In-situ GPS data for validation of ice velocity are only sparsely available. Velocity products are therefore also compared with ice velocity maps retrieved from other sensors (e.g. S1 vs TSX) to estimate product performance and uncertainty. As an additional quality test, velocity results on stable terrain (rock outcrops), where no movement is expected, are analysed.
1.3.3. Product Nomenclature
The Filename of the IV Products consists of:
C3S_GrIS_IV_250m_S1_<Startdate><Enddate><Version>.<Format>
Table 3: Product nomenclature
C3S | Program |
---|---|
GrIS | Region (Greenland Ice Sheet) |
IV | Parameter (Ice Velocity) |
250m | Grid Spacing |
S1 | Platform |
<Startdate> | Date of first acquisition |
<Enddate> | Date of last acquisition |
<Version> | Version of the product in the format vMAJOR_MINOR |
<Format> | Fileformat (nc: NetCDF) |
Example: C3S_GrIS_IV_250m_S1_20191001_20200930_v1_3.nc, Explanation: Annually averaged Greenland Ice Sheet surface velocity derived from Sentinel-1 for Oct 1st 2019 – Sept 30th 2020 provided at 250 m grid spacing, processing version 1_3 for C3S.
1.3.4. Structure of NetCDF files
The product metadata is included in the NetCDF file and is shown in Table 4. The NetCDF file contains 7 sub-datasets. Table 5 shows as an example the metadata for the "land_ice_surface_easting-_velocity" sub-dataset.
Table 4: Example metadata for Greenland Ice Sheet Velocity product
NC_GLOBAL#comment=Ice velocity map of Greenland derived from Sentinel-1 SAR data acquired from 2018-10-01 to 2019-09-30. The surface velocity is derived applying feature tracking techniques. The ice velocity map is provided at 250m grid spacing in North Polar Stereographic projection (EPSG: 3413). The horizontal velocity is provided in true meters per day, towards EASTING(vx) and NORTHING(vy) direction of the grid, and the vertical displacement (vz), is derived from a digital elevation model. Provided is a NetCDF file with the velocity components: vx, vy, vz and vv (magnitude of the horizontal components), along with maps showing valid pixel count and uncertainty (stdx, stdy). The product was generated by ENVEO. NC_GLOBAL#contact=www.c3s.com NC_GLOBAL#Conventions=CF-1.7 NC_GLOBAL#creation_date=2020-06-29 NC_GLOBAL#history=Initial product version 1.3 NC_GLOBAL#institution=Copernicus Climate Change Service NC_GLOBAL#keywords=EARTH SCIENCE CLIMATE INDICATORS CRYOSPHERIC INDICATORS GLACIAL MEASUREMENTS ICE SHEET VELOCITY CRYOSPHERE GLACIERS/ICE SHEETS NC_GLOBAL#license=C3S general license NC_GLOBAL#project=C3S_312b_Lot4_ice_sheets_and_shelves NC_GLOBAL#reference=Main: Nagler, T.; Rott, H.; Hetzenecker, M.; Wuite, J.; Potin, P. The Sentinel-1 Mission: New Opportunities for Ice Sheet Observations. Remote Sens. 2015, 7, 9371-9389. NC_GLOBAL#source=Copernicus Sentinel-1A and Sentinel-1B NC_GLOBAL#summary=Ice velocity derived for Greenland Ice Sheet gridded at 250m averaged from 2018-10-01 to 2019-09-30. NC_GLOBAL#title=Ice Velocity of the Greenland Ice Sheet Subdatasets: SUBDATASET_1_NAME=NETCDF:"C3S_GrIS_IV_250m_S1_20181001_20190930_v1_3.nc":land_ice_surface_easting_velocity SUBDATASET_1_DESC=[10801x5984] land_ice_surface_easting_velocity (32-bit floating-point) SUBDATASET_2_NAME=NETCDF:"C3S_GrIS_IV_250m_S1_20181001_20190930_v1_3.nc":land_ice_surface_northing_velocity SUBDATASET_2_DESC=[10801x5984] land_ice_surface_northing_velocity (32-bit floating-point) SUBDATASET_3_NAME=NETCDF:"C3S_GrIS_IV_250m_S1_20181001_20190930_v1_3.nc":land_ice_surface_vertical_velocity SUBDATASET_3_DESC=[10801x5984] land_ice_surface_vertical_velocity (32-bit floating-point) SUBDATASET_4_NAME=NETCDF:"C3S_GrIS_IV_250m_S1_20181001_20190930_v1_3.nc":land_ice_surface_velocity_magnitude SUBDATASET_4_DESC=[10801x5984] land_ice_surface_velocity_magnitude (32-bit floating-point) SUBDATASET_5_NAME=NETCDF:"C3S_GrIS_IV_250m_S1_20181001_20190930_v1_3.nc":land_ice_surface_measurement_count SUBDATASET_5_DESC=[10801x5984] land_ice_surface_measurement_count (32-bit integer) SUBDATASET_6_NAME=NETCDF:"C3S_GrIS_IV_250m_S1_20181001_20190930_v1_3.nc":land_ice_surface_easting_stddev SUBDATASET_6_DESC=[10801x5984] land_ice_surface_easting_stddev (32-bit floating-point) SUBDATASET_7_NAME=NETCDF:"C3S_GrIS_IV_250m_S1_20181001_20190930_v1_3.nc":land_ice_surface_northing_stddev SUBDATASET_7_DESC=[10801x5984] land_ice_surface_northing_stddev (32-bit floating-point)
Table 5: Example metadata for land_ice_surface_easting_velocity subdataset.
DDriver: netCDF/Network Common Data Format Files: C3S_GrIS_IV_250m_S1_20181001_20190930_v1_3.nc Size is 5984, 10801 Coordinate System is: PROJCRS["WGS 84 / NSIDC Sea Ice Polar Stereographic North", BASEGEOGCRS["WGS 84", DATUM["World Geodetic System 1984", ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], ID["EPSG",4326]], CONVERSION["unnamed", METHOD["Polar Stereographic (variant B)", ID["EPSG",9829]], PARAMETER["Latitude of standard parallel",70 ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8832]], PARAMETER["Longitude of origin",-45, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8833]], PARAMETER["False easting",0, LENGTHUNIT["metre",1], ID["EPSG",8806]], PARAMETER["False northing",0, LENGTHUNIT["metre",1], ID["EPSG",8807]]], CS[Cartesian,2], AXIS["x",east, ORDER[1], LENGTHUNIT["metre",1]], AXIS["y",north, ORDER[2], LENGTHUNIT["metre",1]], ID["EPSG",3413]] Data axis to CRS axis mapping: 1,2 Origin = (-640000.000000000000000,-655500.000000000000000) Pixel Size = (250.000000000000000,-250.000000000000000) Metadata: crs#false_easting=0 crs#false_northing=0 crs#grid_mapping_name=polar_stereographic crs#latitude_of_projection_origin=90 crs#spatial_ref=PROJCS["WGS 84 / NSIDC Sea Ice Polar Stereographic North", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.0174532925199433, AUTHORITY["EPSG","9122"]], AUTHORITY["EPSG","4326"]], PROJECTION["Polar_Stereographic"], PARAMETER["latitude_of_origin",70], PARAMETER["central_meridian",-45], PARAMETER["scale_factor",1], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]], AXIS["X",EAST], AXIS["Y",NORTH], AUTHORITY["EPSG","3413"]] crs#standard_parallel=70 crs#straight_vertical_longitude_from_pole=-45 crs#unit=metre land_ice_surface_easting_velocity#coordinates=y x land_ice_surface_easting_velocity#description=easting ice velocity land_ice_surface_easting_velocity#grid_mapping=crs land_ice_surface_easting_velocity#units=m/day land_ice_surface_easting_velocity#_FillValue=3.4028235e+38 NC_GLOBAL#comment=Ice velocity map of Greenland derived from Sentinel-1 SAR data acquired from 2018-10-01 to 2019-09-30. The surface velocity is derived applying feature tracking techniques. The ice velocity map is provided at 250m grid spacing in North Polar Stereographic projection (EPSG: 3413). The horizontal velocity is provided in true meters per day, towards EASTING(vx) and NORTHING(vy) direction of the grid, and the vertical displacement (vz), is derived from a digital elevation model. Provided is a NetCDF file with the velocity components: vx, vy, vz and vv (magnitude of the horizontal components), along with maps showing valid pixel count and uncertainty (stdx, stdy). The product was generated by ENVEO. NC_GLOBAL#contact=www.c3s.com NC_GLOBAL#Conventions=CF-1.7 NC_GLOBAL#creation_date=2020-06-29 NC_GLOBAL#history=Initial product version 1.3 NC_GLOBAL#institution=Copernicus Climate Change Service NC_GLOBAL#keywords=EARTH SCIENCE CLIMATE INDICATORS CRYOSPHERIC INDICATORS GLACIAL MEASUREMENTS ICE SHEET VELOCITY CRYOSPHERE GLACIERS/ICE SHEETS NC_GLOBAL#license=C3S general license NC_GLOBAL#project=C3S_312b_Lot4_ice_sheets_and_shelves NC_GLOBAL#reference=Main: Nagler, T.; Rott, H.; Hetzenecker, M.; Wuite, J.; Potin, P. The Sentinel-1 Mission: New Opportunities for Ice Sheet Observations. Remote Sens. 2015, 7, 9371-9389. NC_GLOBAL#source=Copernicus Sentinel-1A and Sentinel-1B NC_GLOBAL#summary=Ice velocity derived for Greenland Ice Sheet gridded at 250m averaged from 2018-10-01 to 2019-09-30. NC_GLOBAL#title=Ice Velocity of the Greenland Ice Sheet x#axis=X x#long_name=x coordinate of projection x#standard_name=projection_x_coordinate x#units=m y#axis=Y y#long_name=y coordinate of projection y#standard_name=projection_y_coordinate y#units=m Geolocation: LINE_OFFSET=0 LINE_STEP=1 PIXEL_OFFSET=0 PIXEL_STEP=1 SRS=GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNI T["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]] X_BAND=1 X_DATASET=NETCDF:"C3S_GrIS_IV_250m_S1_20181001_20190930_v1_3.nc":x Y_BAND=1 Y_DATASET=NETCDF:"C3S_GrIS_IV_250m_S1_20181001_20190930_v1_3.nc":y Corner Coordinates: Upper Left ( -640000.000, -655500.000) ( 89d18'52.26"W, 81d33'27.74"N) Lower Left ( -640000.000,-3355750.000) ( 55d47'51.49"W, 59d11'51.92"N) Upper Right ( 856000.000, -655500.000) ( 7d33'22.15"E, 80d 4'16.22"N) Lower Right ( 856000.000,-3355750.000) ( 30d41'23.64"W, 58d47'38.46"N) Center ( 108000.000,-2005625.000) ( 41d55' 3.66"W, 71d36'42.42"N) Band 1 Block=748x1351 Type=Float32, ColorInterp=Undefined NoData Value=3.4028234663852886e+38 Unit Type: m/day Metadata: coordinates=y x description=easting ice velocity grid_mapping=crs NETCDF_VARNAME=land_ice_surface_easting_velocity units=m/day _FillValue=3.4028235e+38
2. Data access information
2.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
Hvidberg, C.S., et al., User Requirements Document for the Ice_Sheets_cci project of ESA's Climate Change Initiative, version 1.5, 03 Aug 2012. Available from: http://www.esa-icesheets-cci.org/
Rizzoli, P., Martone, M., Gonzalez, C., Wecklich, C., Tridon, D.B., Bräutigam, B., Bachmann, M., Schulze, D., Fritz, T., Huber, M. and Wessel, B., (2017). Generation and performance assessment of the global TanDEM-X digital elevation model. ISPRS Journal of Photogrammetry and Remote Sensing, 132, pp.119-139
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.