NSIDC DAAC

ICESat/GLAS Data at NSIDC

Home  |   Data Summaries  |   Order Data  |   Data Releases  |   Images  |   Tools  |   News  |   FAQs  |   Research

Data Releases

Page last updated: 05/09/07

This page summarizes changes and known limitations associated with each release of ICESat/GLAS data. Changes in input data and algorithms result in a new release, set by the GLAS Software Development Team. Products generated with a new release may be processed forward in time or reprocessed from earlier data. If reprocessed, NSIDC either deletes or hides data granules from previous releases, once they have been replaced by the new release. The first data release was Release-12; there were no products with release numbers lower than 12. We recommend you work with the latest release of data (the highest release number), which include improvements to product algorithms. See the following for more information.

Laser Operational Periods
(on the Data Summaries page)
Important information about laser operating periods and data releases, including metadata table The Attributes for ICESat Laser Operations Periods
ICESat Science Investigator-led Processing System (I-SIPS) Release Information Table summarizing release information for all ICESat/GLAS data distributed by NSIDC
NASA Wallops Flight Facility's ICESat/GLAS Web site Links to release notes on Wallops Flight Facility's Web site

Click on a release number below to view a summary of its release notes.

Release-28
Release-26
Release-24
Release-22
Release-21
Release-19
Release-18
Release-17 Atmospheric Products (GLA02, 07-09)
Release-14 Land Products (GLA14)
Release-13 and -14 Period Data (Oct./Nov. 2003)
Release-12 Cal/Val Period Data


Release-28

This release incorporates significant changes to the Level-1B and Level-2A product parameters. There are now improved corrections for saturation effects and improved flags to aid in data selection. Note that the saturation elevation correction needs to be added to the elevations by the user, but the correction is not generally recommended for land surfaces (GLA14). Many parameters have been improved due to updated algorithms, and some are now being computed for the first time. The data dictionary and flag descriptions have been updated. A significant effort has been made to improve descriptions of the altimetry products. As a result of applying the updated oscillator frequency to the initial time tag estimates, the record index may have changed from the previous release by one or two counts. All Release-28 products have consistent time tags. The sample times have not changed. The surface slope and roughness are output as invalid values for this release.

A new release number designation has been adopted (YXX Release Numbers) to denote accuracy of orbit and laser pointing knowledge, which contribute to elevation accuracy. Release-28 products generally have YXX release number = 428 in the new designation, with the exception of the Laser 1 operation period, which is designated as YXX=228. Use caution when comparing 228 data to 428 data.

Notable changes are summarized on this page. See the following for more information:

File naming has changed. Please see the following for more information:

Altimetry Products (GLA05, 06, 12, 13, 14, and 15)

Please read the Notes About Saturation Elevation Correction below. The changes to altimetry products include:

Notes About Saturation Elevation Correction

The Release-28 products contain an improved correction to the surface elevation for saturated returns. The corrections are not applied to the elevations on the products. Detector saturation, caused by higher than predicted energy returns, occurred frequently and is most pervasive in the early operation periods of each laser (Laser 1, Laser 2A, Laser 2B, Laser 3A, and Laser 3B). The saturation elevation correction mitigates this problem. Ice sheet crossovers for the Laser 3A and Laser 3B operation periods show a reduction in the ice sheet crossover standard deviation from 20 cm down to below 10 cm. The user is encouraged to apply the correction by adding it to the surface elevation when the correction value is valid. Unsaturated returns will have a correction value of zero.

The saturation elevation correction has two known problems in the Release-28 data.

The effect of both of these issues is that the correction will not be “defined” for all possible saturated returns. It is still recommended that the correction be used when valid. These issues will be corrected in an upcoming release.

Atmospheric Products (GLA07, 08, 09, 10 and 11)

The changes to atmospheric products include:

Known Problems with Release-28


Release-26

The data dictionaries and flag descriptions were significantly updated for this release. There were no L1A changes.

Notable changes are summarized here. See a list of product format and definition changes associated with Release-26. Please refer to the NASA Wallops Flight Facility's GSAS v5.2 Release Notes for detailed release notes for each product in Release-26.

Altimetry Products

The changes to altimetry products are:

Atmospheric Products

Release-26 incorporates two significant additions for GLAS 532 nm data products. First, there are improvements to parameter tables used to calculate the aerosol optical depth and cloud multiple scattering factors; these improvements considerably increase accuracy. Second, Release-26 is the first to include data products from the Laser 2B (February 2004) operating period.

Other issues for Release-26 are.

Known Problems with Release-26


Release-24

This is a major release. Most changes affect the Level-1B and Level-2 products. Only minor changes were made in Level-1A processing (correction made to the EU conversion of GLA04 IST i_VTBoreH and i_VTBoreV). The following products were updated to Release-24: GLA01 - GLA15.

Notable changes are summarized here. See a list of product format and definition changes associated with Release-24. Please refer to the NASA Wallops Flight Facility's GSAS v5.0 Release Notes for detailed release notes for each product in Release-24.

Altimetry Products

The key change made to the altimetry products is the lengthening of all of the Level-1B and Level-2 products allowing space for new variables and additional growth. These new variables include:

Changes have been made to both the transmitted and received energy calculations, as suggested by the instrument team. For details of the energy computation, see the waveform ATBD.

Two minor changes have been made to the standard parameterization waveform fitting. The standard fit now uses the absolute peak location change (instead of a percentage change) as the convergence criteria in the Gaussian fitting. An error in the computation of the standard fit sigma has also been corrected. Neither of these changes has significantly affected the elevations on the data.

In the alternate waveform fitting, an option to normalize the waveform based on the peak amplitude was implemented. An error in reporting the alternate fit standard deviation was corrected and the precision was changed on the product to accommodate the smaller numbers. The alternate fit parameters were changed to perform normalization by peak amplitude before fitting, the minimum distance between selected peaks was set to be greater than 15 ns, and the weight of the sigma widths was tuned to provide the best alternate fits on selected land waveforms.

Atmospheric Products

Through I-SIPS Release-23, cloud detection using the GLAS 1064 nm atmospheric channel was limited to 4 and 1 second resolution. This was due mainly to the low inherent signal-to-noise of the 1064 data. Over approximately the last 6 months, we have been striving to produce 1064 cloud top height on a shot-to-shot basis (40 Hz). The new algorithms have been tested and are included in this release of GSAS. The limited amount of testing done to date has shown that we can definitely detect clouds at this resolution, though the exact limits of what can be detected (in terms of optical depth) have not yet been quantified. An assessment will be performed after the laser 2A data has been reprocessed using this release of the GSAS software, by comparing the 1064 nm cloud detection with that of the 532 nm channel. A part of the cloud detection algorithm is based on the integrated 1064 nm signal, which is also stored as a separate quantity on the Release-24 GLA09 product. Testing has shown that when the magnitude of the integrated signal exceeds a certain threshold, it is highly likely (as corroborated by the 532 channel) that there is a cloud present even if the threshold algorithm did not detect a cloud. In this case, the cloud height is set to 10 km and a flag is set to indicate this fact. The algorithm also interrogates the ground return signal width in an effort to detect very low clouds. Over flat terrain such as ice sheets, the 1064 ground return signal has a characteristic width. When this signal is unusually wide, and clouds have not been detected higher up, the cloud top height is set to 100 m and a flag is set to indicate the likely presence of low clouds. While these approaches improve the overall cloud detection, they also have their limits. Extensive testing is needed to determine those limits.

Known Problems with Release-24


Release-22

The following products were updated to Release-22: GLA01, GLA02, GLA03, GLA04, GLA05, GLA07, GLA08, and GLA09. This release contains significant improvements to the waveform alternate fitting process, to energy/reflectance calculations, and to atmospheric products. Notable changes are as follows:

Altimetry Products

The energy calculation improved for GLA01. For GLA05, energy/reflectance calculations improved and a saturation index was added. The alternate fit uses the normalized waveform based on area before functional fit. (Standard fit uses raw voltage waveform.) The alternate fit now also uses up to five largest peaks plus the peak nearest the ground. The algorithm keeps all six peaks, but amplitudes are allowed to go to zero. The alternate fit least-squares exit criteria changed; maximum iterations are now significantly less frequent. If the maximum iteration flag is set, data should not be used (or used with caution). See a list of product format and definition changes associated with release-22.

Known Problems with Release-22 Altimetry Products

Atmospheric Products

The 532 nm cloud layer detection improved. The 40 Hz cloud search is now executed for every shot, independent of the results of cloud searches at lower resolutions. Also, the starting height for the cloud search is now 10 km, which means that clouds are now reported when they are found up to 10 km altitude. (In prior releases, 4 km was the maximum altitude). These changes mean that a cloud could be reported for one or more 40 Hz shots of a given second, while no cloud was reported at 1 second or 5 Hz resolution.

The major problems with these data pertain to the 532 channel only. The 1064 backscatter data (GLA07) appear well-calibrated and the 1064 cloud heights are robust both day and night.

The 532 channel had only 5 mJ of energy and a fractured far-field laser beam pattern that reduced the signal strength considerably. As a result, the calibration for daytime 532 data is often poor. The calibration of nighttime data is generally acceptable. The 532 channel signal levels are very low, which severely affects the daytime measurements.

Because of the poor 532 channel signal, elevated aerosol and PBL height data (GLA08) will be unusable during day and only marginal at night. Thicker aerosol and PBL layers can still be detected, but only at night. Cloud height data (GLA09) are good during night, but only thick clouds can be detected during the day.

Optical depth and extinction data (GLA10 and GLA11) are not included in this release but attempts will be made to release nighttime only data in the future.


Release-21

This release contains the following improvements to the elevation products only (GLA06, GLA12-15). No product format changes were implemented.

Known Problems with Release-21

Known problems are the same as Release-19, plus flagging does not correctly indicate when GLA05 pass-through PAD is used in place of data from ANC09.


Release-19

Notable improvements include the following:

Altimetry products

Atmospheric products

Known Problems with Release-19


Please refer to the NASA Wallops Flight Facility's GSAS v4.1 Release Notes for detailed release notes for each product in Release-19.


Release-18

Improvements include the following:

Known Problems in Release-18 Altimetry Products

ICESat/GLAS Release-18 and earlier data products error in received energy computation:

The 40 Hz received gain values i_gainSet1064 reported on GLA01, and i_gval_rcv reported on GLA05, are offset from the shot to which they actually pertain. The gain used onboard for shot N+1 is contained in the telemetry for shot N. Ground processing should therefore use the gain reported with for shot N in the energy computations for shot N+1.

For Release-18 and earlier ICESat/GLAS data products, this shift was not made; therefore, all parameters that used the gain on GLA01, 5, 6, 12, 13, 14, and 15 have some error. The energy parameters are in error by an amount proportional to the shot-to-shot change in the received gain setting. No error is induced in the elevation measurements from this wrong application of the gain setting.

Parameters that are known to be affected are:

i_RecNrgAll_EU
i_RecNrgLast_EU
i_reflctUncorr
i_reflctuncmxpk

Gain is also used to determine the waveform saturations flags:

i_WFqual (gwi_satFlagLo)
i_WFqual (gwi_satFlagHi)
i_WFqual (gwi_satFlagHiF).

The gain on GLA01 will be shifted to the correct shot in a future release of the data. This will fix all the problems in the ground processing and correct the parameters on the data products.

Known Problems in Release-18 Atmospheric Products

GLA08 and GLA09: Cloud aerosol discrimination usually results in layers misidentified as aerosol when they should in fact be cloud. The problem is most severe in the high-latitude regions and over the poles, and in regions of multiple cloud layers where the signal attenuation from the upper layer makes the lower layer "look" like aerosol.

GLA07: Background computation problem in areas of very bright light cases signal to be poorly calibrated in some instances where the background is high.


Please refer to the NASA Wallops Flight Facility's GSAS v4.0 Release Notes for detailed release notes for each product in Release-18.


Release-17 Atmospheric Products (GLA02, 07-09)

Release-17 data include protection for floating point exceptions, corrections to science software documentation, and improved Quality Assurance Products (QAPs) and browse products. Release-17 products include data from GLAS laser #2. Atmospheric products now have substantially improved PBL height retrieval (GLA08), improved cloud layer detection from 1064 channel, and a new time-dependent calibration for 1064 nm.

The GLAS data product Release-17 is the initial production run of data products for the GLAS 532 nm atmospheric channel. Preliminary data quality and validation activities have been preformed, but additional testing and improvements will be incorporated in further releases.

GLAS is a two channel atmospheric lidar and was operated for two periods in 2003. The instrument was run with only the less sensitive 1064nm channel from February 21 to March 29. GLAS was first fully turned on at the end of September. There is data from both 1064 and 532 nm channels from Sept. 25 to Nov. 18. From Sept. 25 to 28 the filter and alignment for the 532 nm channel were not yet set, and there were also subsequent short periods where the 532 nm data quality was variable when the alignments were reset. For most of the fall period of operation, we obtained very good data.

The GLAS instrument was launched with a number of hardware level problems that effect atmospheric data quality. These include some signal non-linearity and variation of offset under certain conditions. The large majority of errors have been corrected or mitigated in the data processing, but some known issues require further attention.

GLA02 Lidar Signal

The 1064 nm detector circuits introduce “signal droop” resulting in too low values following strong signals. The effect is corrected in the GLA02 software for low and medium strength signals, but due to saturation and other effects, is not fully corrected for strong return signals.

GLA07 Calibrated Observed (Attenuated) Backscatter Cross Section

The 532 backscatter cross section calibration is very good at night and with daytime low to medium backgrounds. At higher background levels (> 10 photons/bin), the detector response varies somewhat with range, causing problems in the background subtraction. which in turn causes the 532 profile to not follow Rayleigh with high accuracy. The problem is estimated to be in about 10% of the daytime data. The 1064 calibration must be based on comparison to the 532 nm data where signals overlap and ground data intercomparisons and is not as accurate as the 532 nm calibration. The 1064 nm calibration in release 17 are now thought to give backscatter values about 5-10% low (calibration constant about 5-10% high), based on further analysis and results from 532 calibration transfers.

GLA08 Planetary Boundary Layer Heights

Planetary boundary layer heights are in general very good. Some problems are:

1. Over polar regions some low cirrus type clouds may at times be identified as PBL.

2. In regions where the gradient of scattering at PBL top is very small, a PBL height may be assigned even though it is likely in error. This problem can be mitigated by checking the value of the PBL quality flag and eliminating those PBL heights with a low quality flag.

3. A few of the elevated aerosol layers reported in GLA08 are actually clouds. This is especially a problem in polar regions and multi layer clouds where attenuation has made the lower layers appear as lower cross sections.

GLA09 Cloud/Aerosol Heights Including Multi Layers

1. Under certain attenuation conditions in geometrically thick but optically rarefied layers, the analyzed layer bottom is too high.

2. Some false bottoms are generated in the daytime signal which sometimes leads to multiple layer results where only one layer exists.

3. PBL heights sometimes supercede elevated aerosol layer results and cloud results.

4. Layer cloud/aerosol determination depends upon magnitude of backscatter coefficient and vertical gradient of magnitude. Two known situations where clouds are designated to be aerosol layers are: a) when a vertically extensive but tenuous cloud layer exists at a low altitude, as those found over Antarctica; b) when attenuation from an overlying layers reduces the signal of a lower layer to the aerosol category. Another known situation where aerosol is designated as cloud is some cases of smoke injected into higher regions in the troposphere. (These false results may be remedied by using a color ratio parameter in the discrimination algorithm to be completed.)

5. The results of the layer boundary detection function for nighttime (constant threshold) have not been quality checked.

6. Diminishing of ground signal caused by steep terrestrial slopes can cause false negative ground signal detection results and false positive lower level layer results.

7. In cases where no layers are detected at 1 Hz resolution, the ground location detected and stored in the GLA09 product fails to get transferred to the GLA10 and GLA11 products.

GLA10 Extinction and Attenuation Corrected Backscatter Cross Sections
GLA11 Thin Cloud and Aerosol Optical Depth

Due to a coding error in the latitudes and longitude which was found after release 17 processing, these products will not be available until Release-18.


Please refer to the NASA Wallops Flight Facility's GSAS v3.9 Release Notes for detailed release notes for each product in Release-17.


Release-13 and -14 Period Data (October/November 2003)

Release-13 and Release-14 use data from GLAS laser #2, and include improved waveform fitting and precision-pointing data for altimetry products (GLA05, GLA06, GLA12-15).

Release-13 (GLA01) and Release-14 (GLA05, GLA06, GLA12-15) elevation data products covering October/November 2003 have not been fully calibrated. The fundamental elevation product contains: time tag, geodetic latitude, longitude and height (elevation) above the TOPEX reference ellipsoid. These products result from a combination of round-trip laser pulse travel time, precision orbit determination (POD), and precision attitude determination (PAD). Contributions from various sources lead to the following absolute accuracy assessment of the fundamental elevation products.

a. The time tag, when converted to laser pulse "ground bounce time", has been validated with an accuracy of about one microsec.

b. The geodetic latitude and longitude provide the horizontal location of the laser spot on the Earth's surface. The accuracy of this location is assessed to be 15 meters (1-sigma), resulting from about 5 arc-sec accuracy in the laser-pointing knowledge.

c. The geodetic height may be biased by approximately 10 cm, but this bias is also dependent on the surface slope. For example, on 1 degree slopes with the 5 arc sec horizontal position accuracy, there will be an apparent 15 cm bias in range and hence in elevation. The precision of the height will also be dependent on surface characteristics, with smooth flat surfaces producing 2-3 cm precision, while rough surfaces will be decimeters.

In some occasional cases, off-nadir pointing was used at some locations, which introduces an apparent surface slope of up to 5 degrees, depending on the angle used.

Comments on Release-14 Land Products (GLA14)

ICESat products for land areas (all continents including Antarctica, and ocean islands, including Greenland) are distributed in the GLA14 data product. The GLA14 results are produced using processing results designated as alternate in the GLA05 waveform processing data product, in contrast to the standard processing results used to create the GLA12 ice sheet products. ICESat products have been produced in a number of releases, with improving data acquisition, calibration, and processing techniques. Data covering an 8-day period in March, 2003 were distributed as Release 12. Data covering a 33-day period in October and November, 2003 are distributed as Release 14.

The primary difference between the Release 12 and Release 14 versions of the GLA14 land products is a substantially reduced frequency of saturated returns, achieved by modifying the onboard gain algorithm to better accommodate rapid changes in received signal amplitude caused by changes in atmospheric transmissivity and surface reflectance. The modified gain algorithm was implemented for the fall, 2003 operation period.

An additional difference between the Release 12 and 14 data sets is the threshold level used in ground processing to detect the presence of a return signal. A threshold 5 times the standard deviation of the background noise level was used in Release 12. This threshold detected some weak returns from clouds that were thus included in the GLA14 product. To reduce the number of cloud false-alarm returns, a higher threshold, 7.5 times the background noise standard deviation, was used in Release 14. However, the higher threshold also causes some weak returns from the land surface beneath clouds to be excluded from the GLA14 product.

In addition to detecting the presence of returns, the threshold is used to define waveform extent based on the start and end of the signal corresponding to the first and last crossing of the threshold. The higher 7.5x threshold used in Release 14 underestimates waveform extent as compared to the 5x threshold used in Release 12. The Release 12 designations of signal start and end are a more accurate measure of the highest and lowest detected surfaces within the laser footprint. Furthermore, the elevation reported in the GLA14 product corresponds to the centroid of the received signal between signal start and end. The elevation is therefore sensitive to the threshold used.

Because of the above considerations, future releases of the GLA14 products will revert to the lower 5x threshold in order to detect weak surface returns and more accurately define waveform extent.

Future releases of GLA14 products will include several additional improvements. Fitting of Gaussian distributions to waveform peaks will be refined to more fully characterize multi-peaked distributions. Waveform fitting results in Release 14 are known to have some deficiencies and should be used only as a provisional product. Additionally, Releases 12 and 14 of GLA14 include the location of the signal centroid and end, but not signal start, as defined by the threshold above noise. The signal start currently included in GLA14 is based on a different threshold defined as a percentage of the maximum signal amplitude. This signal start will be replaced by the noise threshold signal start in future releases, so as to be compatible with the centroid and signal end designations. The position of the noise threshold signal start is available in the Release 12 and 14 data sets in the GLA05 waveform processing product.

Finally, for some returns in the data acquired during 2003 the signal extent exceeds the maximum width of the waveform and the upper end of signal is truncated. This occurs for footprint locations with large relief caused by steep topographic slopes and/or tall forest canopies. In mountainous regions, approximately 10% of the waveforms are truncated in this way. During 2003 operations the sampled waveform covers an extent of 81.6 m (544 waveform bins x 0.15 m per bin). Where the signal exceeds this maximum extent, signal start is designated as the first (highest) waveform bin. For operations in 2004 an onboard compression algorithm has been implemented whereby the waveform extent is increased to 150 m, thus minimizing occurrences of signal truncation.


Please refer to the NASA Wallops Flight Facility's GSAS v3.7 Release Notes for detailed release notes for each product in Releases 13 and 14.


Release-12 Cal/Val Period Data

Release-12 is the original, at-launch version of all products, using data from GLAS laser #1.

The data calibration and validation activities by the ICESat/GLAS Science Team were interrupted by the premature end of laser-1 operation after about 36 days on March 28, 2003, and were resumed with operation of laser-2 beginning on September 25. The results from ICESat demonstrate the unique capabilities of the laser measurements with unprecedented accuracy. However, the cal/val period data, in particular, should be regarded as not fully calibrated and validated (i.e. with some uncertain errors). Some limitations may be described by flags on the data products. One important limitation (to be corrected in future releases) is the laser-attitude accuracy (at the 10 to 15 arcsec level), which affects the accuracy of the surface elevations (cms to meters depending on surface slope) and geo-locations of the measurements (tens of meters).

Another limitation is the distortion of the received-pulse shape due to "saturation" of the detector amplifier. Type-1 saturation, which non-linearly distorts the top and trailing edge of the waveform, affects much of the data over ice and other strong reflectors (tens of cm range effect) and may be correctable with future algorithms. Type-2 saturation, for which the gain goes to a maximum and stays there awhile, causes an effect in the meters range. Type 2 has been corrected for laser 2 data by changing the instrument gain algorithm. Lidar data on clouds and aerosols from the 1064 nm channel is included, but the backscatter cross-section magnitudes have a large offset that will be corrected. The more-sensitive 532 nm lidar channel was not yet activated, so there are no optical-depth or extinction cross-section products. Other limitations may be discovered as analysis and calibration/validation continues.

Please refer to the NASA Wallops Flight Facility's GSAS v3.6 Release Notes for detailed release notes for each product in Release-12.