This data set provides measurements of sea ice freeboard and sea ice thickness for the Arctic region. The data were acquired from the Ice, Cloud, and land Elevation Satellite (ICESat) Geoscience Laser Altimeter System (GLAS) instrument, the Special Sensor Microwave/Imager (SSM/I), and climatologies of snow and drift of ice. The data span six GLAS campaigns, laser 3D through 3I, from 21 October 2005 to 05 November 2007. Data parameters include sea ice freeboard and thickness measured in meters derived from GLAS Release 28 data. The data are provided in three formats: ASCII track data derived from binary track data, binary gridded polar stereographic data derived from ASCII gridded polar stereographic files, and Portable Network Graphic (PNG) image files. Also included are mask files used in preparation of the image files, a Mapx grid definition file, and grid cell center latitude and longitude files. The ASCII track data vectors of position and ice thickness have a resolution of about 170 meters in the along-track direction. The binary gridded polar stereographic data have a resolution of 25 km. The data set is approximately 940 megabytes. The data are available via FTP.
The following example shows how to cite the use of this data set in a publication. For more information, see our Use and Copyright Web page.
The following example shows how to cite the use of this data set in a publication. List the principal investigators, year of data set release, data set title and version number, dates of the data you used (for example, March to June 2004), publisher: NSIDC, and digital media.
Yi, Donghui and Jay Zwally. 2010. Arctic Sea Ice Freeboard and Thickness. Boulder, Colorado USA: National Snow and Ice Data Center. Digital media.
| Category | Description |
|---|---|
| Data format | Space delimited ASCII text Binary gridded 32-byte little-endian |
| Spatial coverage and resolution |
Spatial Coverage: Arctic Region Southernmost Latitude: 65° N Northernmost Latitude: 86° N Westernmost Longitude: 180° W Easternmost Longitude: 180° E Spatial resolution: The centers of 70 m spots illuminated by the laser on the earth's surface are separated in the along-track direction by 170 m. |
| Temporal coverage and resolution | Laser 3D 2005-10-21 to 2005-11-24 Laser 3E 2006-02-22 to 2006-03-27 Laser 3F 2006-05-24 to 2006-06-26 Laser 3G 2006-10-25 to 2006-11-27 Laser 3H 2007-03-12 to 2007-04-14 Laser 3I 2007-10-02 to 2007-11-05 A Release Schedule lists the temporal coverage of each ICESat/GLAS product. Data are sampled 40 times per second, captured during 14.8 orbits per day. |
| Tools for accessing data | ASCII data: any Web browser, plain text display or spreadsheet software. Binary gridded data: ENVI, ArcGIS, other similar software packages. |
| Grid/projection description | Binary gridded data are in polar stereographic projection. ASCII track data are provided with latitude and longitude for each data point and are not projected. Geoid models: EGM96 global gravity field model is used below 64° N. ArcGP is used north of 64° N. |
| File naming convention | ASCII files: laser3h1341001.txt Binary gridded data: laser3d_freeboard_mskd.img laser3d_freeboard_mskd.img.hdr laser3d_freeboard_mskd.png laser3d_thickness_mskd.img laser3d_thickness_mskd.img.hdr laser3d_thickness_mskd.png Mask: gsfc_25n.msk gsfc_25n.msk.hdr Mapx grid definition file: PS25km_north.gpd Grid cell center lat / long files: PS25km_north_lon.img PS25km_north_lon.img.hdr PS25km_north_lat.img PS25km_north_lat.img.hdr |
| File size | ASCII files range from approximately 2 KB to 1 MB. Binary gridded files: ENVI header files (.hdr) are 749 bytes each. Image files (.img) are 532 KB each. Portable Network Graphics files (.png) range from approximately 31 KB to 42 KB. Mapx grid definition file: 1 KB. Grid cell center lat / long image files are 532 KB each. Grid cell center lat / long files ENVI header files are 1 KB each. Mask file is 133 KB. Mask file ENVI header is 1KB. |
| Parameter(s) | Sea ice freeboard in meters Sea ice thickness in meters |
| Metadata Access | View Metadata |
| Procedures for obtaining data | Available via FTP. |
Donghui Yi
Stinger Ghaffarian Technologies, Inc.
Cryospheric Sciences Branch
NASA Goddard Space Flight Center
Greenbelt, MD 20771 USA
NSIDC User Services
National Snow and Ice Data Center
CIRES, 449 UCB
University of Colorado
Boulder, CO 80309-0449 USA
phone: +1 303.492.6199
fax: +1 303.492.2468
form: Contact NSIDC User Services
e-mail: nsidc@nsidc.org
Each space delimited ASCII vector data file contains file headers followed by columns for latitude, longitude, freeboard, and thickness.
Longitude in WGS 84 degrees in the approximate range of 0 to 360.
Freeboard in meters above local sea level at the time of the measurement.
Thickness in meters.
Freeboard = 0 means the sea ice surface is at sea level. The average of the lowest 1percent reference elevation is used as reference sea level, so about 0.5 percent of the data are below sea level, that is, negative freeboard. Freeboard values less than zero are set to zero.
Note: Thickness values of -999 indicate that a thickness value could not be calculated and should be treated as missing data.
The freeboard and thickness binary gridded files are 304 columns x 448 rows containing little-endian floating-point values in meters. The grid is the 25 km polar stereographic grid used for SSM/I processing.
Data are available on the FTP site in the following directory:
ftp://sidads.colorado.edu/pub/DATASETS/NSIDC0393_GLAS_SI_Freeboard_v01
ASCII files are named according to the following convention and as described in Table 1:
Example file name: laser3i1343001.txt
laserLPTTTTCCC.txt
| Variable | Description |
|---|---|
| laserLP | GLAS Laser Period used to collect data. LP = 3d, 3e, 3f, 3g, 3h, or 3i |
| TTTT | Laser track number: 0001 through 0432, or 1282 through 1354 |
| CCC | Cycle of Reference Orbit: 001, 002, or 003 |
| .txt | Identifies the file as ASCII data |
Binary gridded files are named according to the following convention and as described in Table 2:
Example file name: laser3d_freeboard_mskd.img
Example file name: laser3d_freeboard_mskd.img.hdr
laserLP_nnnnnnnnn.mskd.zzz
| Variable | Description |
|---|---|
| laserLP | GLAS Laser Period used to collect data. LP = 3d, 3e, 3f, 3g, 3h, or 3i |
| _nnnnnnnnn | Indicating 'freeboard' or 'thickness' |
| _mskd | Indicating masked image |
| .zzz | Identifies file type: .img (image), .img.hdr (header), or .png (Portable Network Graphics) |
The grid cell center files provide the latitude and longitude for the center of each grid cell as a little-endian floating-point number. The files are named according to the following convention, and as described in Table 3.
Example file name: PS25km_north_lon.img
PS25km_north_mmm.zzz
| Variable | Description |
|---|---|
| PS | Polar Stereographic |
| 25km | 25 kilometer grid |
| _north | northern hemisphere |
| mmm | latitude (lat) or longitude (lon) |
| .zzz | Identifies file type: .img (image), img.hdr (header) |
The mask file consists of values 0 = water, 1 = land. The file is named gsfc_25n.msk. The corresponding header file is gsfc_25n.msk.hdr.
ASCII file sizes range from 2 K to 1 MB
The .hdr files are each 749 bytes.
The .img files are each 532 KB.
The .png files range from 31 KB to 44 KB.
Arctic region
Southernmost Latitude: 65° N
Northernmost Latitude: 86° N
Westernmost Longitude: 180° W
Easternmost Longitude: 180° E
Note: The ASCII data are reported with longitude of 0 degrees to 360 degrees.
The resolution of the gridded data is 25 km. Vectors have a resolution of about 170 meters in the along-track direction, or the distance between the centers of adjacent GLAS footprints. Several kilometers typically separate the ICESat GLAS data tracks.
Two geoid models are used in this study. North of 64° N, the Arctic Gravity Project (ArcGP) latitude limit, the ArcGP geoid is used. The Earth Gravitational Model 1996 (EGM96) geoid is used below 64° N. The difference, D, between ICESat measured sea level and the geoid was used to create a 5 km D-grid for each period. The mean of the 5 km D-grids was used as an improved geoid.
Release 28 Laser Identifier 3D: 21 Oct 2005 to 24 Nov 2005
Release 28 Laser Identifier 3E: 22 Feb 2006 to 27 Mar 2006
Release 28 Laser Identifier 3F: 24 May 2006 to 26 Jun 2006
Release 28 Laser Identifier 3G: 25 Oct 2006 to 27 Nov 2006
Release 28 Laser Identifier 3H: 12 Mar 2007 to 14 Apr 2007
Release 28 Laser Identifier 3I: 02 Oct 2007 to 05 Nov 2007
A Release Schedule lists the temporal coverage of each ICESat/GLAS product.
Each ICESat campaign period, 3D through 3I, represents one Arctic-wide ice thickness assessment. Thus the temporal resolution for this data set is approximately two to three times per year.
The parameters of this data set are sea ice freeboard, and sea ice thickness.
Sea ice freeboard is the height in meters of the sea ice above the water level. Thickness is the thickness in meters of the sea ice.
In the ASCII files, the thickness value -999 indicates missing values, possibly due to failure of the thickness algorithm for certain corresponding freeboard measurements.
The binary gridded *.img and *.png files have been masked. Table 4 shows the mask color value assignments to surface features.
| .img File Value | .png File Value | Surface Feature |
|---|---|---|
| -4 | gray = 128 | land south of 65° N |
| -3 | gray = 164 | land at or north of 65° N |
| -2 | gray = 64 | water south of 65° N |
| -1 | gray = 96 | water at or north of 65° N |
| greater than or equal to 0 | full color (rainbow) applied across parameter value range | sea ice freeboard or thickness measurements in meters |
The following sample of the laser3d0001002.txt ASCII file shows header information, and the first four records of Latitude, Longitude, Freeboard and Thickness values.

The following sample shows the laser3d_freeboard_mskd.img masked binary gridded image.

The following sample shows the laser3d_freeboard_mskd.png masked binary gridded image.

Table 5 summarizes the GLAS single-shot error budget for elevation measurements (Zwally et al. 2002).
| Error Source | Error Limit |
|---|---|
| Precision orbit determination (POD) | 5 cm |
| Precision attitude determination (PAD) | 7.5 cm |
| Atmospheric delay | 2 cm |
| Atmospheric forward scattering | 2 cm |
| Other (tides, etc.) | 1 cm |
| RSS | 13.8 cm |
For further discussion on error sources, see Zwally et al. 2002, and also refer to the Error Sources section of the NSIDC GLAS/ICESat L1 and L2 Global Altimetry Data Web page.
For further discussion on quality assessment of GLAS products, refer to the Quality Assessment section of the NSIDC GLAS/ICESat L1 and L2 Global Altimetry Data Web page.
Data are available via FTP.
The total volume for the ASCII track data files is approximately 932 megabytes.
The total volume for the binary gridded image, header and masked image files is approximately eight megabytes.
The ASCII track data may be displayed with any Web browser or plain text display software.
The binary gridded data may be displayed using ENVI, ArcGIS, or other similar software packages.
Two geoid models are used in this study. The ArcGP geoid is used north of 64° N since that is the ArcGP latitude limit. The EGM96 geoid is used below 64° N.
Elevations varying more than plus-or-minus 4 meters are not used. This condition filters out some land, island, and iceberg data. These conditions are the same as in Zwally et al. 2008.
ICESat measures a surface elevation profile referenced to an ellipsoid. Due to the limited accuracy of the geoids and ocean tide models, and poor knowledge of the dynamic topography, sea-ice surface elevation referenced to a geoid cannot be regarded as sea-ice freeboard. The information needed to calculate sea-ice freeboard is the elevation difference between the top of the snow surface, local sea levels, and snow height and density above the snow/ice interface. If the elevation difference is known, even if the absolute elevations are biased, the sea-ice freeboard can be determined. Thus, the knowledge of relative elevation is crucial while absolute elevation is less important. This is the underlying concept in the derivation of freeboard.
In this study, constant densities of ρW = 1023.9 kg m-3 and ρI = 915.1 kg m-3 are used to calculate sea ice thickness from the freeboard. There is no spatial variation of now density ρS. Snow density, including the time variation, is based on Kwok 2008. The range of the snow density is 0.16 to 0.40.
Refer to the Processing Steps section below for details on the application of theory of measurements.
The data for this data set were acquired with the GLAS instrument onboard the Ice, Cloud, and land Elevation satellite (ICESat), from the SSM/I instrument on board the Defense Meteorological Satellite Program (DMSP), and from weather stations on the sea ice.
Freeboard is measured from ICESat elevation profiles (Zwally et al. 2008). Snow depth is interpolated, both spatially and temporally, from climatology snow depth in situ measurements (Warren et al. 1999). Thickness is estimated from ICESat freeboard and climatology snow depth. SSM/I daily ice concentration data (Gloersen et al. 1992) from January 2003 to October 2008 are used to determine sea ice boundaries.
The grids of sea ice freeboard and thickness were derived from the corresponding ICESat point measurements of freeboard and thickness using the following steps:
![]() |
(Equation 1) |
| Variable | Description |
|---|---|
| w | Weight |
| d | Distance |
| r | Maximum radius from which to draw samples. To interpolate values across the polar region, r was set to 420 km. |
Sea ice freeboard and thickness are calculated from ICESat GLAS ground tracks using the following steps.
![]() |
(Equation 2) |
| Variable | Description |
|---|---|
| S | Measure of the broadening of the transmitted pulse associated with surface topography and the undesirable effects of saturation and atmospheric forward scattering |
| c | Speed of light |
| σR | Echo waveform 1-sigma pulse width |
| σT | Transmitted waveform 1-sigma pulse width |
| ΔHib = 9.948 × (Patm - P) | (Equation 3) |
| Variable | Description |
|---|---|
| ΔHib | Inverse barometer correction |
| Patm | Surface atmospheric pressure |
| P | Time varying mean of the global surface atmospheric pressure over the oceans |
| h = Hie + ΔHib + ΔHsat - hg | (Equation 4) |
| Variable | Description |
|---|---|
| h | Elevation above the geoid |
| Hie | ICESat measured surface elevation |
| ΔHib | Inverse barometer correction |
| ΔHsat | Saturation correction |
| hg | Geoid height |
| hr = h - hm | (Equation 5) |
| Variable | Description |
|---|---|
| hr | Relative elevation |
| h | Elevation above geoid |
| hm | 50 km running mean of elevation h |
| hd = hm - hs | (Equation 6) |
| Variable | Description |
|---|---|
| hd | Difference between the ICESat measured sea level and the geoid used |
| hm | 50 km running mean of elevation h |
| hs | Sea level |
| F = hr - hs | (Equation 7) |
| Variable | Description |
|---|---|
| F | Freeboard height |
| hr | Relative elevation |
| hs | Sea level |
| T = | ρW | F - | ρW - ρS | TS | (Equation 8) |
|
|
|
||||
| ρW - ρI | ρW - ρI |
| Variable | Description |
|---|---|
| T | Sea ice thickness |
| F | Freeboard height |
| TS | Snow depth |
| ρW | Water density |
| ρS | Snow density |
| ρI | Sea ice density |
Where:
| Variable | Description |
|---|---|
| Fi | freeboard |
| Fx | snow accumulation factor |
| d | ratio Fi/Fx |
| Ts | snow depth |
| Ci | bilaterally interpolated ice concentration at a point from SSMI ice concentration |
| Ti | sea ice thickness |
Dorandeu, J. and P. Y. Le Traon. 1999. Effects of Global Mean Atmospheric Pressure Variations on Mean Sea Level Changes from TOPEX/POSEIDON. Journal of Atmospheric and Oceanic Technology, 16(9): 1279-1283.
Gloersen, P., W. J. Campbell, D. J. Cavalieri, J. C. Comiso, C. L. Parkinson, and H. J. Zwally. 1992. Arctic and Antarctic Sea Ice, 1978-1987: Satellite Passive-Microwave Observations and Analysis. NASA SP-511, 290 pp.
Kwok, R., and G. F. Cunningham. 2008. ICESat Over Arctic Sea Ice: Estimation of Snow Depth and Ice Thickness. Journal of Geophysical Research 113(C08010), doi: 10.1029/2008JC004753.
Picot, N., K. Case, S. Desai, and P. Vincent. 2003. AVISO and PODAAC User Handbook, IGDR and GDR Jason Products. SMM-MU-M5-OP-13184-CN (AVISO), JPL D-21352 (PO.DAAC).
Stackpole, J. D. 1994. The WTO Format for the Storage of Weather Product Information and the Exchange of Weather Product Message in Gridded Binary Form. NOAA Office Note 388.
Warren, S. G., I. G. Rigor, N. Untersteiner, V. F. Radionov, N. N. Bryazgin, Y. I. Aleksandrov, and R. Colony. 1999. Snow Depth on Arctic Sea Ice. Journal of Climate 12: 1814 - 1829.
Zwally, H. J., D. Yi, R. Kwok, and Y. Zhao. 2008. ICESat Measurements of Sea Ice Freeboard and Estimates of Sea Ice Thickness in the Weddell Sea. Journal of Geophysical Research 113(C02S15), doi: 10.1029/2007JC004284.
Zwally, H. J., B. Schutz, W. Abdalati, et al. 2002. ICESat's Laser Measurements of Polar Ice, Atmosphere, Ocean, and Land. Journal of Geodynamics 34(3-4): 405-445.
The acronyms used in this document are listed in Table 13.
| Acronym | Description |
|---|---|
| ArcGP | Arctic Gravity Project |
| ASCII | American Standard Code for Information Interchange |
| CIRES | Cooperative Institute for Research in Environmental Science |
| DMSP | Defense Meteorological Satellite Program |
| EGM96 | Earth Gravitational Model 1996 |
| ENVI | Environment for Visualizing Images |
| FTP | File Transfer Protocol |
| GLAS | Geoscience Laser Altimeter System |
| GSFC | Goddard Space Flight Center |
| ICESat | Ice, Cloud, and land Elevation Satellite |
| NCEP | National Center for Environmental Protection |
| NSIDC | National Snow and Ice Data Center |
| PNG | Portable Network Graphics |
| PAD | Precision Attitude Determination |
| POD | Precision Orbit Determination |
| PODAAC | Physical Oceanography Distributed Active Archive Center |
| SSM/I | Special Sensor Microwave/Imager |
| URL | Uniform Resource Locator |
| WGS 84 | World Geodetic System 1984 |
September 2010
http://nsidc.org/data/docs/daac/nsidc0393_arctic_seaice_freeboard/index.html