Scientists
Research Grants
Scientific Publications
Informatics Research
NSIDC Scientific Expeditions
Scientists blog from Antarctica and provide a glimpse of what it's like to do research in the field. 
What's hot in the news around climate and sea ice and what are scientists talking about now? 
A Digital Elevation Model (DEM), ice thickness grid, and bedrock elevation grid of Greenland acquired as part of the PARCA program are available in ASCII text format at a 5 km grid spacing in a polar stereographic projection. DEM data are a combination of ERS-1 and Geosat satellite radar altimetry data, Airborne Topographic Mapper (ATM) data, and photogrammetric digital height data. Ice thickness data are based on approximately 700,000 data points collected in the 1990s from a University of Kansas airborne ice penetrating radar (IPR). Nearly 30,000 data points were collected in the 1970s from a Technical University of Denmark (TUD) airborne echo sounder. Bamber subtracted the ice thickness grid from the DEM to produce a grid of bedrock elevation values. Applications include studies of gravitational driving stress and ice volume (mass balance) of the Greenland Ice Sheet. Each of the three grids is approximately 1.5 MB. Data are available via FTP. Data access is unrestricted, but we recommend that users register with us. Registered users automatically receive e-mail notification of product updates and changes to processing.
Get Data
DOWNLOADING DATA VIA HTTPS
To learn more about Earthdata Login and to register for an account, please see How to Register with Earthdata Login. Once you have logged in, data can be downloaded via a Web browser, command line, or client. For help with downloading data, please see Options Available for Bulk Downloading Data from HTTPS with Earthdata Login.
Greenland 5 km DEM, Ice Thickness, and Bedrock Elevation Grids, Version 1
Geographic Coverage
| Parameter(s): |
|
|---|---|
| Spatial Coverage: |
|
| Spatial Resolution: | Not Specified |
| Temporal Coverage: |
|
| Temporal Resolution: | Not specified |
| Data Format(s): |
|
| Platform(s) | AIRCRAFT, GEOSAT |
| Sensor(s): | ALTIMETERS, ATM, RA, RADAR ALTIMETERS, RADAR ECHO SOUNDERS |
| Version: | V1 |
| Data Contributor(s): | Jonathan Bamber |
| Metadata XML: | View Metadata Record |
Data Citation
As a condition of using these data, you must cite the use of this data set using the following citation. For more information, see our Use and Copyright Web page.
Bamber, J. 2001. Greenland 5 km DEM, Ice Thickness, and Bedrock Elevation Grids, Version 1. [Indicate subset used]. Boulder, Colorado USA. NASA National Snow and Ice Data Center Distributed Active Archive Center. doi: https://doi.org/10.5067/01A10Z9BM7KP. [Date Accessed].Literature Citation
As a condition of using these data, we request that you acknowledge the author(s) of this data set by referencing the following peer-reviewed publication.
Bamber, J., R. L. Layberry, and S. P. Gogenini. 2001. A new ice thickness and bed data set for the Greenland ice sheet 1: Measurement, data reduction, and errors., Journal of Geophysical Research. 106. 33,773-33,780. https://doi.org/10.1029/2001JD900054
Layberry, R. L. and J. Bamber. 2001. A New Ice Thickness and Bed Data Set for the Greenland Ice Sheet: 2. Relationship between Dynamics and Basal Topography, Journal of Geophysical Research. 106. 33,781 - 33,788. https://doi.org/10.1029/2001JD900053
A Digital Elevation Model (DEM), ice thickness grid, and bedrock elevation grid of Greenland acquired as part of the PARCA program are available in ASCII text format at a 5 km grid spacing in a polar stereographic projection. DEM data are a combination of ERS-1 and Geosat satellite radar altimetry data, Airborne Topographic Mapper (ATM) data, and photogrammetric digital height data. Ice thickness data are based on approximately 700,000 data points collected in the 1990s from a University of Kansas airborne ice penetrating radar (IPR). Nearly 30,000 data points were collected in the 1970s from a Technical University of Denmark (TUD) airborne echo sounder. Bamber subtracted the ice thickness grid from the DEM to produce a grid of bedrock elevation values. Applications include studies of gravitational driving stress and ice volume (mass balance) of the Greenland Ice Sheet. Each of the three grids is approximately 1.5 MB.
Detailed Data Description
This data set consists of a Digital Elevation Model (DEM) and ice thickness grid, at 5 km resolution. It also includes a bedrock elevation grid, computed as the difference between the DEM and the ice thickness grid. DEM data are based on a combination of ERS-1 satellite altimetry data, Airborne Topographic Mapper (ATM) data, and digital height data from Kort & Matrikelstyrelsen (KMS) and Geologiske Undersogelser (GEUS). Ice thickness data were collected with airborne ice penetrating radar (IPR) in the 1990s by the Remote Sensing Laboratory at the University of Kansas and with an airborne radar echo sounder in the 1970s by the Technical University of Denmark.
Usage
Ice thickness data were subtracted from a Digital Elevation Model (DEM) of Greenland to produce a bedrock elevation data set. Both data sets are useful in studying a variety of ice sheet characteristics. Bedrock elevation data contribute to a more detailed analysis of topographic patterns such as channels and basins that affect ice velocity rates. Ice thickness profile data, when combined with sea water density and ice column density data, are used to compute flotation thickness for an ice column. Ice thickness data are also used in estimating gravitational driving stress. In addition, recent elevation and ice thickness data provide scientists with an updated estimate of Greenland ice sheet volume to compare with previous studies.
The DEM, ice thickness grid, and bedrock grid are each provided as fixed-width, tab-delimited ASCII text grids (310 columns x 561 rows) of floating-point data. ArcInfo users need to flip the grid vertically after importing the data, so the grid displays properly in ArcInfo.
The DEM, ice thickness grid, and bedrock grid are each a separate granule (file).
In June 2001, NSIDC obtained corrected copies of all three grids from Bamber. See Known Problems with the Data for more information.
surface_5km_corrected: DEM
thick_5km_corrected: ice thickness grid
bed_5km_corrected: bedrock elevation grid
Each of the three files is 1.5 MB.
Figure 4 shows the extent of flight line coverage for the 1970s and 1990s airborne radar surveys of the Greenland Ice Sheet, which define the spatial coverage of this data set. The 1990's data are plotted in black, and the 1970's data are plotted in grey. (Bamber 2000a).
Spatial Resolution
The DEM, ice thickness grid, and bedrock elevation grid are interpolated to 5 km grid resolution. However, the true horizontal resolution varies according to slope and surface characteristics. The DEM originally had a horizontal spatial resolution of 1 km, with accuracy ranging from 20 m to 200 m over bare rock areas, depending on the source data. Horizontal resolution of the ice thickness grid is between 5 km and 50 km, depending on flight line coverage.
Projection
All grids are in a polar stereographic projection with reference to the WGS84 ellipsoid. The standard parallel is 71°N, and the projection parallel (latitude of the center of the projection) is 90°N. The North Pole is the origin of the projection. The central meridian is 39°W.
Grid Description
Each grid consists of 310 columns and 561 rows, with data points gridded every 5 km and the origin at the North Pole. Corner points are as follows:
| Longitude | Latitude | x (km) | y (km) | |
|---|---|---|---|---|
| Upper Left | 92.1301024° W | 80.8152653° N | -800 | -600 |
| Lower Left | 52.2405199° W | 58.6292691° N | -800 | -3400 |
| Upper Right | 10.3987054° E | 81.5294837° N | 700 | -600 |
| Lower Right | 27.3663660° W | 58.8136131° N | 700 | -3400 |
| Longitude | Latitude | x (km) | y (km) | |
|---|---|---|---|---|
| Upper Left | 92.3305365° W | 80.8106358° N | -802.5 | -597.5 |
| Lower Left | 52.2710201° W | 58.6035491° N | -802.5 | -3402.5 |
| Upper Right | 10.6177120° E | 81.5269394° N | 702.5 | -597.5 |
| Lower Right | 27.3342981°W | 58.7883271° N | 702.5 | -3402.5 |
The DEM ice thickness grid and bedrock grid are each a separate granule (file).
In June 2001, NSIDC obtained corrected copies of all three grids from Bamber. See Known Problems with the Data for more information.
surface_5km_corrected DEM
thick_5km_corrected ice thickness grid
bed_5km_corrected bedrock elevation grid
Data used to produce the DEM were collected at different times, beginning in March 1985 through September 1986 with the Geosat geodetic phase, and ending in April 1994 through April 1995 with the ERS-1 336-day repeat cycle. Ice thickness data from the University of Kansas IPR were collected from 1993 to 1999. Echo sounding data from the Technical University of Denmark were collected in the late 1970s.
Parameter Description
Figure 5 is a visual representation of all the parameters. Bedrock elevation is computed as the difference of DEM surface elevation and ice thickness values. For example, if the surface elevation (a) is 400 m above sea level, and ice thickness (b) is 100 m, the bedrock elevation (c) is 300 m above sea level (c = a - b).
Unit of Measurement
meters [m]
Parameter Source
- DEM data are based on a combination of the following data sources:
- ERS-1 (geodetic phase) satellite radar altimetry data
- Geosat satellite radar altimetry data
- Airborne Topographic Mapper (ATM) data
- Digital height data from photogrammetric analysis of aerial photos
- Recently unclassified DEMs provided by U.S. National Imaging and Mapping Agency (NIMA)
The primary data set used to compute ice thickness for Greenland was a collection of ice penetrating radar (IPR) data from the University of Kansas, Remote Sensing Laboratory. Beginning in 1993, a 150 MHz IPR was used for yearly airborne surveys across Greenland, which yielded nearly 700,000 ice thickness data points. A secondary source of ice thickness data came from the Technical University of Denmark (TUD) in the late 1970s. The TUD airborne radar collected nearly 30,000 data points throughout Greenland.
Parameter Range
DEM: -0.1 m to 3278.300 m
Ice thickness: 0 m to 3366.500 m
Bedrock elevation: -963.100 m to 3239.000 m
DEM values of -0.1 m represent ocean, so the true minimum elevation value is 0 m. Bedrock elevation values are listed with respect to the WGS84 ellipsoid.
Sample Data Record
Grid dimensions for the DEM, ice thickness grid, and bedrock elevation grid are 310 columns by 561 rows in fixed-width, tab-delimited ASCII text format. However, the format of the data file is actually 10 columns by 17391 rows. In other words, one row of 310 data values is represented by 31 lines in the ASCII data file. Following is sample output from the file 'thick_5km' (ice thickness grid):
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 9.105 48.894
130.515 228.119 346.149 472.796 604.376 736.562 855.897 923.472 974.034 1006.809
1053.530 1083.509 1129.979 1184.289 1242.173 1312.345 1392.936 1471.406 1554.642 1612.586
1666.496 1715.926 1762.610 1808.412 1837.086 1854.994 1859.129 1876.878 1877.157 1870.814
1851.615 1836.129 1824.894 1822.442 1828.953 1829.900 1803.543 1712.215 1658.382 1606.686
1545.734 1516.651 1551.700 1684.686 1808.865 1855.575 2056.482 2014.391 1979.463 1962.578
2033.146 1988.963 1733.050 1654.831 1769.080 1810.623 1881.694 1926.600 1964.864 1865.733
1637.862 1437.682 1343.301 1332.935 1322.639 1303.178 1124.136 1331.643 1206.667 892.235
767.588 589.441 504.706 573.235 642.080 669.964 610.454 435.203 253.982 91.241
9.533 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
|
Software and Tools
Validation by Source
For DEM validation, reference points derived from aerial photogrammetry along the rock outcrops serve as a measure of comparison for surface elevation values. They have a root-mean-square (rms) accuracy of about 10 m. Elevation values computed from aerial photogrammetry have an rms of 25-30 m. For coastal areas where elevation points were digitized from maps, the error is much larger (200 m to 300 m rms). The most accurate and extensive validation data are the ATM measurements. Bamber compared satellite radar altimetry data and ATM data by estimating the DEM elevation at each ATM point using bilinear interpolation. Areas around the Thule and Petermann glaciers showed a significant mismatch between the aerial photogrammetry and satellite radar altimetry data. Although these two areas encompass ten percent of the high slope ATM data coverage, they are poorly mapped compared to the rest of the DEM. Consequently, the ATM data for these specific areas were excluded from the comparison. Results indicate that the overall bias between the DEM and ATM data is small, but the standard deviation for all points is only 6.96 m. This is a significant improvement on a previous high-resolution Greenland DEM where the standard deviation was 10 to 11 m. The random error in the DEM is a factor of two less. This is attributed primarily to the slope correction scheme used in this data set, which uses a relocation method rather than the simpler, but less accurate, direct approach (Bamber, Ekholm, and Krabill 2000, Bamber 1994a).
Bamber used cross-over analysis as the primary method for assessing the accuracy of the ice thickness measurements collected by airborne radar echo-sounders. The average cross-over difference is 130 m for the 1970s data, which is primarily due to navigation data errors. The average cross-over difference is 15 m for the 1990s data. Bamber removed some 1990s tracks due to problems with matching radar data to navigation data. Standard deviations in the data are a combination of the random error in the ice thickness measurements and local variability in bed topography. Rougher areas have a higher standard deviation even though the accuracy of each data point is the same. The average standard deviation of ice thickness is 50 m, with the highest values concentrated around the eastern ice sheet margin where the topography is highly variable. Along-track standard deviations for unaveraged ice thickness data shows a strong correlation with actual ice thickness, with thin ice showing a much greater standard deviation.
Bamber constructed a histogram of the standard deviations for different thickness ranges. He removed data points with anomalously high standard deviations from the data set using a thickness-dependent threshold where the distribution was one percent of its maximum value. He removed two percent of the gridded data with this method (Bamber, Layberry, and Gogenini 2000).
Confidence Level or Accuracy Judgment
Horizontal resolution varies according to slope and surface characteristics. The DEM originally had a horizontal spatial resolution of 1 km, with accuracy ranging from 20 m to 200 m over bare rock areas, depending on the source data. DEM accuracy varies from -1.04 m (± 1.98 m) to -0.06 m (± 14.33 m) for slopes of 0-1°. Mean accuracy over the whole ice sheet was -0.33 m (± 6.97 m) (Bamber, Ekholm, and Krabill 2000). Horizontal resolution of the ice thickness grid is between 5 km and 50 km, and vertical accuracy is between 10 m and 100 m, depending on flight line coverage. Bogorodskiy reports that a typical radar-sounding survey has an inherent uncertainty of about 15 m for ice depth measurements (1985).
Measurement Error for Parameters
Bamber calculated errors in the estimate of derived ice thickness measurements by comparing the radar-based measurements with those from the Greenland Ice Sheet Project (GISP) and Greenland Ice Core Project (GRIP) ice cores. Results showed that radar-based measurements were within ± 10 m of ice core measurements. He made no corrections for firn effect.
Since elevation and ice thickness data from the 1970s were much less accurate than the 1990s data, errors in the merged grid were most significant where 1970s and 1990s tracks were in close proximity. Since errors in the 1970s data did not cause anomalously high ice thickness gradients at a distance of three grid points (15 km) away from the 1990s data, Bamber removed gridded data points closer than three grid cells to the 1990s data (Bamber, Layberry, and Gogenini 2000).
Limitations of the Data
Accuracy of surface elevation and ice thickness data is questionable in the ice sheet margin and coastal areas, due to the highly variable topography. Also, the improved resolution of the DEM does not significantly affect ablation estimates (Bamber, Ekholm, and Krabill 2000).
Known Problems with the Data
In June 2001, NSIDC discovered errors in the data values for 'surface_5km' (DEM). Because the ice thickness grid is subtracted from the DEM to produce the bedrock elevation grid, the incorrect DEM data resulted in inconsistent values for 'thick_5km' (ice thickness grid) and 'bed_5km' (bedrock elevation grid). NSIDC obtained corrected copies of all three grid files from Bamber to ensure consistency and accuracy among all grids.
Future Modifications and Plans
Bamber plans to update the ice thickness and bedrock grids with newly acquired radar echo sounding data within the next two to three years. The new data should improve coverage and accuracy, particularly near the margins. Also, he will update the DEM with ICESat data about 6-12 months after the launch of ICESat. This will only make a marginal improvement in accuracy over most of the ice sheet but will, in particular, improve the accuracy in the marginal areas on steeper slopes, as well as some of the non-ice covered areas that currently only have digitized cartographic coverage
Data Acquisition and Processing
Radar Altimetry
Satellite radar altimeters measure the time it takes an electromagnetic signal to travel from the altimeter antenna to the ice sheet surface and back to the altimeter's receiver. This "range measurement" allows investigators to determine the satellite's height above the ice sheet. NSIDC's Radar Altimeter document describes the instrument and how it works.
Corrections are applied to the range measurements to account for the fact that returns are different over ice than over the ocean. The correction for interpreting the data over ice is called "retracking," and is described in papers by Davis and Zwally (1993), and Zwally, et al. (1983).
Over sloping terrain, the radar altimeter measurement needs to be corrected because the return comes from a point not directly below the satellite, but to the uphill side. The elevation indicated by the return time in this case is higher than that directly below the satellite. The data are slope-corrected to reduce the errors, using the slope correction algorithm from Brenner et al. (1983).
Echo Sounding
Radar echo sounding is a method of active remote sensing that measures ice thickness. An antenna mounted on a platform moving over the ice surface emits a short electromagnetic pulse which penetrates the surface and is reflected by inhomogeneities in the ice and underlying bedrock, sending an echo back to the antenna. The pulse duration depends on the system resolution requirements and the distance between the measured layers. Ice thickness is determined by analysis of the pulse delay time in the ice. The electromagnetic pulse loses energy as it travels through the ice, causing the returning echo to be smaller than the original pulse. The change in signal strength during the pulse's travel to the bedrock and back to the receiver is described by the following equation (Bogorodskiy 1985):
Nsum = NG + NB+ NF+ ND + NA + NP
Where:
Nsum = change in signal strength
NG = geometrical spreading losses
NB = losses due to reflections from interfaces
NF = focusing factor, i.e., changes in signal strength due to focusing effects (also called refraction gain)
ND = losses due to scattering
NA = losses due to signal absorption in the ice
NP = apparent losses due to rotation in the polarization of the received signal relative to the orientation of the receiving antenna
Characteristics of the electromagnetic signal are dependent on the physical properties of the ice. Glacial ice has an inhomogeneous structure caused by variations in the ice, air bubbles, and the presence of water. The variation of electrical parameters between different layers and the boundary conditions of those layers are important factors in effective radar sounding (Bogorodskiy 1985).
See ERS-1 (PDF file) and Geosat platform documents for characteristics of the two satellite radar altimeters which were used to collect surface elevation data.
Satellites transmitted radar signals to the ground in the form of radar echograms. Scientists estimated ice thickness by tracking the peaks of surface and ice-bedrock signals and computing the number of range cells between the peaks. They multiplied the result by 4.494, which was obtained by dividing free-space range-cell dimension by the refractive index of ice, to estimate ice thickness. They processed the data further by coherently integrating two to four samples and incoherently integrating 10 to 20 samples to reduce fading. Because of increased absorption loss for warm ice over areas in the south and central parts of the ice sheet, scientists conducted additional processing using the f-k migration algorithm to extract ice thickness. This algorithm was essentially a correlation processor that reversely migrated the received signals back to their sources within the ice. Processing consisted of a variable that described the migration of the received signal back to the air-ice interface, and a second variable that described the propagation through the ice layer (Gogineni 2000).
Processing Steps
Surface Elevation (DEM)
Bamber attempted to improve the accuracy of surface elevation estimates of variable topographic surfaces with minimal biases using the following steps (Bamber 1994b):
-
Implementation of a range-estimate refinement (waveform tracking) using the offset center of gravity method for calculating the waveform amplitude, with a power threshold of 25 percent.
-
Application of a slope correction using a variation of the relocation method (Bamber 1994b).
-
Implementation of various data filtering techniques, including tests applied to the return-echo waveform shape, backscatter coefficient, and retracking correction value for each altimeter height estimate. He removed approximately 27 percent of the data during this procedure.
-
Removal of anomalous orbits by comparing one track with another where they cross (cross-over analysis).
-
Calculation of standard deviations for each filtered group. A significant bias exists in sampling of the surface, from clustering of satellite radar altimetry data points on topographic peaks after using a relocation method for slope correction. To determine the relationship between the local surface slope and the bias, Bamber ran ATM along-track profiles through a high-pass filter with a 30 km threshold. Bias and standard deviation between the satellite radar altimetry data and ATM heights decreased from 3.10 m to -0.68 m, and from 7.15 to 6.68 m, respectively (Bamber, Ekholm, and Krabill 2000).
-
Application of a least squares collocation technique combined with a second order Markov model (that defines the covariances between observed and modeled elevation) to generate the DEM from the various data sources. He selected the data used to estimate each modeled height by defining a north-south and east-west axis through each model point. He then used the three closest data points from each quadrant. He interpolated the model to a Cartesian coordinate grid with a spacing of 0.01° latitude and 0.025° longitude, which corresponds to approximately 1 km (Bamber, Ekholm, and Krabill 2000). To compute a bedrock elevation data grid from the 5 km ice thickness grid, he subsequently regridded the DEM to a 5 km resolution.
Ice Thickness and Bedrock Elevation
Bamber averaged ice thickness data from University of Kansas and TUD radar sounders in a polar stereographic projection using a linear distance weighting method to produce a 310 x 561 grid of georeferenced thickness values with 5 km spacing between data points. He calculated the thickness at each grid point as the weighted average of the data within the search radius of the cell center and the mean latitude and longitude values. He gave data from the 1990s a weighting factor value of ten higher relative to the 1970s data based on cross-over analysis. This initial averaging of the data reduced the number of points used in the final interpolation, reduced the disparity between along- and across-track spacing, and allowed filtering on the basis of the standard deviations of the mean ice thicknesses. Initial filtering involved the removal of gridded data points with only one contributing depth estimate. He implemented quality control procedures by using a histogram to filter anomalous data values (see Quality Assessment section). He appended the coordinates of non-glacier points near the edge of the ice sheet to the mean data values. These were derived from a high resolution (2 km) land ice mask (Weng 1995) to provide a boundary condition at the ice margins with zero ice thickness.
Bamber interpolated the resulting grid of mean ice thickness onto a 310 x 561 grid using a kriging procedure which was suitable for the high range of spatial density of data points. The kriging generates an estimate of the variance of fit of the data to the model it generates (Figure 1). Minimum values of uncertainty lie along the 1990s tracks and have an error of about 10 m. Maximum values lie in areas furthest from a track that have high thickness variability, and have an error of about 100 m. Figure 2 is an image derived from ice thickness values for Greenland.
Finally, Bamber subtracted the ice thickness grid from a DEM (Bamber, Ekholm, and Krabill 2000) to produce a grid of bedrock elevation values (Figure 3).
Calculated Variables
- Surface elevation
- Ice thickness
- Bedrock elevation
Errors and Limitations
Errors in satellite radar altimetry data come from several sources including geographically correlated orbit errors, errors in slope correction procedures, and non-uniform spatial sampling (Bamber, Ekholm, and Krabill 2000). Uncertainties in airborne radar sounding relate to the magnitude of the refractive indices of solid ice, uncertainty in the profile of refractive index in the firn, signal contamination by thermal noise in the equipment, insufficient clarity of the oscilloscope signal, and imperfect signal definition on the recording medium. Large errors in ice thickness measurements can result if pulses reflected from ice inhomogeneities are mistaken for bottom echoes (Bogorodskiy 1985).
The ERS-1 and Geosat satellite radar altimeters measure surface elevation. See ERS-1 (PDF file) and Geosat platform documents for characteristics of the two satellite radar altimeters. Airborne radar echo sounders measure ice thickness. See the Data Collection System section in this document for details about the University of Kansas IPR and TUD radar echo sounder.
Source or Platform Collection Environment
Collection platforms were primarily satellite and aircraft.
Source or Platform Mission Objectives
The primary objective of the missions that collected surface elevation and ice thickness data was to use airborne or satellite-based altimetry in support of scientific research of the mass balance on the Greenland ice sheet.
Source or Platform Program Management
Bamber obtained DEM data from the following sources and programs:
- The European Space Agency (ESA) manages the European Remote Sensing Satellite (ERS-1) (PDF file).
- The U.S. Navy launched and managed the Geosat satellite.
- Kort & Matrikelstyrelsen (KMS) and Geologiske Undersogelser (GEUS) provided digital height data from photogrammetric analysis of aerial photos.
- The U.S. National Imaging and Mapping Agency (NIMA)provided additional DEM data of Greenland.
Ice thickness data were obtained from two different sources:
- The University of Kansas, Remote Sensing Laboratory managed the 1990s IPR system.
- The Technical University of Denmark managed the 1970s airborne echo sounding data.
Coverage Information
Following the ERS-1 altimetry mission, DEM data were available for nearly 80% of Greenland. Ice thickness data coverage varied according to aircraft flight lines. The 1970s and 1990s flight lines are shown in Figure 4. The along-track sampling from the 1970s was much coarser (1 km) than the 1990s data (150 m average sampling). The across-track spacing is highly variable but averages 20 km. Tracks were dense near the airfields used for field surveys (especially Thule and Søndre Strømfjord), while areas near the southern tip of the ice sheet have gaps up to 120 km between tracks. Ice thickness measurements could not be obtained over the deep, narrow channel beneath the Jakobshavn glacier and near the margins of the ice sheet. Researchers collected much of the radar data over the Jakobshavn area during 1997, where post-processing could not be performed to improve the signal-to-noise ratio (Gogineni 2000).
Data Collection System
See ERS-1 (PDF file) and Geosat platform documents for characteristics of the two satellite radar altimeters, which collected surface elevation data.
The University of Kansas IPR system for collecting ice thickness data transmitted a 150 MHz chirped, 1.6 µs pulse with a peak power of 200 W. The radars used Surface Acoustic Wave (SAW) dispersive delay lines to generate the chirp signal and to compress the received signal. Antennae consisted of two four-element dipole arrays mounted under each wing of the aircraft - one for transmission and the other for reception. The radar transmitter and receiver were mounted in a rack inside the aircraft and connected to their respective antennae with cables and a feed network. The effective transmit power at the antennae was about 100 W. A low noise receiver amplified and detected the received signal, at a compressed pulse width of about 60 ns. The digital signal processor digitized the in-phase and quadrature output signals from the detector and integrated them.
Researchers collected ice thickness data from 1998 and 1999 with an improved radar system which consisted of a Digital Signal Processor (DSP), 12-bit analog/digital (A/D) converters, and high-speed Field-Programmable Gate Arrays (FPGAs). One of the data acquisition cards digitized the data and integrated the data signal, while another card generated the Pulse Repetition Frequency (PRF), the clock signal for the A/D converters, and the delayed PRF signal to start digitization. This system could operate with a maximum PRF of 18.4 kHz, in which case the number of range bins had to be decreased to 800. This corresponded to a maximum ice thickness of 3 km for radar operating at an altitude of 500 m over the ice surface. The normal operating mode was a PRF of 9.2 kHz with the integrator set to sum 256 values coherently. With the aircraft flying at an altitude of about 450 m and the airspeed at 130 m/s, the distance traveled was 3.6 m per 256 samples. The radar also employed a kinematic Global Positioning System (GPS) (Gogineni 2000).
List of Sensors or Instruments
- ERS-1 and Geosat radar altimeters
- University of Kansas 150 MHz ice penetrating radar (IPR)
- Technical University of Denmark (TUD) 60 MHz radar echo sounder
References and Related Publications
Document Information
DOCUMENT CREATION DATE
16 January 2001
DOCUMENT REVISION DATE
08 April 2005
