This data set contains snow thickness and internal layer measurements taken from the Center for Remote Sensing of Ice Sheets (CReSIS) ultra wide-band snow radar over land and sea ice in Greenland, the Arctic, and Antarctica. The data were collected as part of Operation IceBridge funded campaigns.
Operation IceBridge products may include test flight data that are not useful for research and scientific analysis. Test flights usually occur at the beginning of campaigns. Users should read flight reports for the flights that collected any of the data they intend to use. Check IceBridge campaign Flight Reports for dates and information about test flights.
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.
Leuschen, Carl. 2010, updated current year. IceBridge Snow Radar L1B Geolocated Radar Echo Strength Profiles, [list dates of data used]. Boulder, Colorado USA: NASA Distributed Active Archive Center at the National Snow and Ice Data Center. Digital media. http://nsidc.org/data/irsno1b.html.
| NASA P-3 and DC-8 aircraft | |
CReSIS Snow Radar |
|
Antarctica, Greenland, and Arctic and Antarctic oceans |
|
| Varies dependent on along-track, cross-track, and aircraft height characteristics |
|
31 March 2009 to the present. |
|
| Seasonal |
|
Snow Thickness |
|
Binary |
|
Metadata Access |
|
Data Access |
Carl Leuschen
CReSIS
Nichols Hall
2335 Irving Hill Road
University of Kansas
Lawrence, Kansas 66045
John Paden, Ben Panzer
CReSIS
Nichols Hall
2335 Irving Hill Road
University of Kansas
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
Data and data products from CReSIS were generated with support from NSF grant ANT-0424589 and NASA grant NNX10AT68G. CReSIS faculty, staff, and students designed, developed, operated, and processed data from the radar systems.
The Snow Radar L1B Geolocated Radar Echo Strength Profiles data set includes echograms with measurements for time, latitude, longitude, elevation, as well as flight path charts and echogram images.
The binary files contain a vector stream of record data. Each record includes seven header values and a Fast Fourier Transform (FFT) window of snow radar data. The FFT window has dimensions N-number of range bins.
The digital image files are JPEG files. The y-axis in the FFT JPEG files shows depth relative to a range around the surface. The surface is in the center of the y-axis and the y-axis is set to a fixed range, usually from 0 m to 60 or 80 m for the land ice, and 0 m to 4 m for sea ice.
The MATLAB files are binary files produced and readable by the proprietary Matlab software or other tools such as the Octave high-level language.
The KML files are flight line browse images for each segment.
The most convenient way to browse the imagery quickly is through the JPEG files. The quickest way to plot the data set is to look at the KML browse files for the entire season.
The 2011 Greenland and Antarctica MATLAB data are divided into segments. A segment is a contiguous data set where the radar settings do not change. A day is divided into segments if the radar settings were changed, hard drives were switched, or other operational constraints required that the radar recording be turned off and on. All data from a particular segment are stored in a directory with the following nomenclature YYYYMMDD_SS where YYYY is the year, MM is the month, DD is the day, and SS is the segment. Segments are always sorted in the order in which the data were collected. Currently at NSIDC, the data directories for 2011 Antarctica and Greenland are named according to this convention, but the previous year directories are not.
Each segment is broken into frames, analogous to satellite SAR scenes, to make analyzing the data easier. Frames span 33 seconds covering 4 to 5 km dependent upon aircraft speed. The frame ID is a concatenation of the segment ID and a frame number and follows the format YYYYMMDD_SS_FFF where FFF is the frame number from 000 to 999. Generally the FFF starts with 0 or 1 and increments by 1 for each new frame, but this is not always the case; only the ordering is guaranteed to match the order of data collection.
For each data frame there is a flight path file (0map) and an echogram file (1echo).
The files are organized on the FTP site, ftp://n4ftl01u.ecs.nasa.gov/SAN2/ICEBRIDGE_FTP/, as described in Figure 1.
Figure 1. Directory Structure
The binary files are named according to the following convention and as described in Table 2:
data01.0180.bin
dataSS.xxxx.bin
Where:
| Variable | Description |
|---|---|
| data | Indicates data file |
| SS | Day segment |
| xxxx | Flight ID number |
| .bin | Indicates a binary file |
The MATLAB files for 2011 Greenland and 2011 Antarctica are named according to the following convention and as described in Table 3:
Data_20110316_01_000.mat
Data_YYYYMMDD_SS_FFF.mat
Where:
| Variable | Description |
|---|---|
| Data | Indicates data file |
| YYYY | Four-digit year |
| MM | Two-digit month |
| DD | Two-digit day |
| SS | Day segment |
| FFF | Frame number |
| .mat | Indicates a MATLAB file |
The flight track MATLAB files and GPS MATLAB files are named according to the following convention and as described in Table 4.
flight_track.mat
gps.mat
Where:
| Variable | Description |
|---|---|
| flight_track or gps | Indicates flight track file or GPS file | .mat | Indicates MATLAB file |
The FFT image JPEG files and Flight Path image JPEG files are named according to the following convention and as described in Table 5.
FFT_image.00.0001.jpg
FFT_image.SS.xxxx.jpg
Flight_Path.00.0001.jpg
Flight_Path.SS.xxxx.jpg
Where:
| Variable | Description |
|---|---|
| FFT_image or Flight_Path | FFT_image = echogram browse image Flight_Path = flight path location plot |
| SS | Day segment |
| xxxx | Flight ID number | .jpg | Indicates JPEG image file |
The flight track image JPEG files are named according to the following convention and as described in Table 6.
flight_track_image10162009.jpg
flight_track_imageMMDDYYYY.jpg
Where:
| Variable | Description |
|---|---|
| flight_track_image | Flight track image map |
| MM | Two-digit month |
| DD | Two-digit day |
| YYYY | Four-digit year | .jpg | Indicates JPEG image file |
For each data frame of the 2011 data, there is a flight path file (0map) and an echogram file (1echo). The file naming conventions are shown below and as described in Table 7.
20110316_01_000_124503_0maps.jpg
20110316_01_000_124503_1echo.jpg
YYYYMMDD_SS_FFF_HHmmss_0maps.jpg
YYYYMMDD_SS_FFF_HHmmss_1echo.jpg
Where:
| Variable | Description |
|---|---|
| YYYY | Four-digit year |
| MM | Two-digit month |
| DD | Two-digit day |
| SS | Segment number |
| FFF | Frame number |
| HHmmss | GPS time stamp for the first range line in the image where HH is 00-23 hours, mm is 00-59 minutes, and ss is 00-59 seconds. |
| 0maps | Flight path file |
| 1echo | Echogram file |
| .jpg | Indicates JPEG image file |
Browse_Data_20111012_01.kml
Browse_Data_YYYYMMDD_SS.kml
Where:
| Variable | Description |
|---|---|
| Browse_Data | Indicates browse data file |
| YYYY | Four-digit year |
| MM | Two-digit month |
| DD | Two-digit day |
| SS | Day segment | .kml | Indicates KML file |
Binary files range from approximately 9 MB to 80 MB.
FFT_image JPEG files range from approximately 31 KB to 48 KB.
Flight path JPEG files are 17 KB each.
Flight track MATLAB files are approximately 76 KB each.
Flight track image JPEG files are 544 KB each.
GPS MATLAB files are approximately 156 MB each.
KML files are approximately 1 KB to 306 KB each.
The entire data set is approximately 6.9 TB.
Spatial coverage for the IceBridge snow radar campaigns include the Arctic, Greenland, Antarctica, and surrounding ocean areas. In effect, this represents the two coverages noted below.
Arctic / Greenland:
Southernmost Latitude: 60° N
Northernmost Latitude: 90° N
Westernmost Longitude: 180° W
Easternmost Longitude: 180° E
Antarctic:
Southernmost Latitude: 90° S
Northernmost Latitude: 53° S
Westernmost Longitude: 180° W
Easternmost Longitude: 180° E
Spatial Resolution varies dependent on along-track, cross-track, and aircraft height characteristics.
Referenced to WGS-84 Ellipsoid
These data were collected as part of Operation IceBridge funded campaigns from 31 March 2009 to the present.
IceBridge campaigns are conducted on an annual repeating basis. Arctic and Greenland campaigns are conducted during March, April, and May, and Antarctic campaigns are conducted during October and November.
The Snow Radar L1B Geolocated Radar Echo Strength Profiles data set contains radar backscatter measurements sensitive to snow thickness and internal layering, collected over land ice and sea ice.
The Snow Radar MATLAB files contain fields as described in Table 9.
| Parameter | Description | Units |
|---|---|---|
| Data | Radar echogram data. The data may be truncated. Presence of the Truncate_Bins variable indicates data has been truncated. Lack of the Truncate_Bins variable indicates data has not been truncated. | Represents relative received power (Watts) |
| Time | Fast time. Zero time is approximately the beginning of the transmit event. | Seconds |
| Depth | Range axis with the origin at the median of the Surface values in this frame and assuming a dielectric of 1.53 (Depth = (Time - median(Surface)) * c/2/sqrt(1.53)) | Meters |
| Truncate_Bins | Indices into the Time and Depth vectors for which the Data is available. Only available when Data are truncated. | n/a |
| Truncate_Mean | Represents a mean of the noise power for the truncated range bins before the surface return. When no range bins were truncated before the surface return the value is NaN. Only available when Data are truncated. | n/a |
| Truncate_Median | Represents a median of the noise power for the truncated range bins before the surface return. When no range bins were truncated before the surface return the value is NaN. Only available when Data are truncated. | n/a |
| Truncate_Std_Dev | Represents a standard deviation of the noise power for the truncated range bins before the surface return. When no range bins were truncated before the surface return the value is NaN. Only available when Data are truncated. | n/a |
| GPS_Time | GPS time when data were collected, seconds since January 1, 1970 00:00:00. This is the ANSI C standard. | Seconds |
| Latitude | WGS-84 geodetic latitude coordinate where data were collected, potentially modified by motion compensation. Always referenced to North. Represents the location of the origin of the trajectory data which is generally not the radar's phase center, but some other point on the aircraft, for example the GPS antenna or the INS. | Degrees |
| Longitude | WGS-84 geodetic longitude coordinate where data were collected, potentially modified by motion compensation. Always referenced to East. Represents the location of the origin of the trajectory data which is generally not the radar's phase center, but some other point on the aircraft, for example the GPS antenna or the INS. | Degrees |
| Elevation | Elevation where data were collected, potentially modified by motion compensation. Referenced to WGS-84 ellipsoid. Positive is outward from the center of the Earth. Represents the location of the origin of the trajectory data which is generally not the radar's phase center, but some other point on the aircraft, for example the GPS antenna or the INS. | Meters |
| Elevation_Correction | Represents the number of zeros that were inserted during elevation compensation for each range line to simulate near-level flight. These zeros are not included in the truncation noise statistics. Only available when data are elevation compensated. | Range bins |
| Surface | Estimated two way propagation time to the surface from the collection platform. This uses the same frame of reference as the Time variable. This information is sometimes used during truncation to determine the range bins that can be truncated. | Seconds |
| *param* | Multiple variables with a name containing the string "param." Contains radar and processing settings, and processing software version and time stamp information. Fields of structures are not static and may change from one version to the next. | n/a |
Data are available via FTP.
MATLAB files may be opened using the NSIDC MATLAB reader, or the Octave high-level language.
JPEG files may be opened using any image viewing program that recognizes the JPEG file format.
KML files are read by GIS software packages and earth browsers such as Google Earth or Google Maps.
The CReSIS snow and kuband radars use a Frequency Modulated Continuous Wave (FMCW) architecture (Carrara et al. 1995). This is done to reduce the required sampling frequency of the Analog to Digital Converter (ADC) and is possible when the range gate is limited. Currently, the range gate is limited to low altitude flights. In the FMCW radars, a long, approximately 250 μs, chirp signal is generated which sweeps linearly in frequency from the start frequency to the stop frequency. This signal is transmitted and also fed to a mixer in the receiver to be used to demodulate the received signal. Signals outside the range gate are suppressed and aliased by the system.
The dominant scattered signal is the specular or coherent reflection from the air-snow surface and shallow layers beneath the surface. A bistatic antenna configuration is used to provide isolation between the transmit and receive paths which is important because the FMCW system receives while transmitting and too little isolation means that the direct path from the transmitter to the receiver will be too strong and will saturate the receiver. The antennas are mounted so the main beam is pointed in the nadir direction to capture the specular surface and layer reflections.
The Pulse Repetition Frequency (PRF), or along-track sampling rate, does not necessarily capture the full Doppler bandwidth for point scatterers without aliasing. However, since the target energy is mostly coherent, it occupies only a small portion of the Doppler spectrum so the undersampling in along-track is not generally a problem. Since the coherent portion of the surface and layer scattering is the primary signal of interest, presumming is used to lower the data rate which effectively low-pass-filters and decimates the Doppler spectrum.
Echograms posted include altitude correction, but the binary files do not. Correction can be applied by shifting a record from bottom to top by the altitude correction value. Altitude variations within a data file are removed by subtracting the minimum altitude from all values. The result is variation in meters from the minimum. These values are then converted to whole pixel values given the radar parameters: sampling frequency = 58.32 MHz, pulse duration, FFT length, and bandwidth. Sampling frequency after the 2009 Greenland campaign is 62.5 MHz.
The code below shows how to load the data in Matlab to generate typical figures, and the illustrations show examples of the figures that were generated from the code.



For a flat surface the range resolution is expressed by Equation 1:
(Equation 1)
Where:
| Variable | Description |
|---|---|
| kt | kt = 1.6 due to the application of a Hanning time-domain window to reduce the range sidelobes of the chirped transmit waveform |
| c | Speed of light in a vacuum |
| B | Bandwidth, nominally 4500 MHz (2 to 6.5 GHz range) |
| n | Index of refraction for the medium |
The bandwidth for a particular segment can be determined by reading the param_radar structure in the echogram data file or by looking at the parameter spreadsheet values f0, f1, and fmult and doing the calculation in Equation 2:
(Equation 2)
Where:
| Variable | Description |
|---|---|
| B | Bandwidth |
| param_radar.f1 | Stop frequency of chirp out of Direct Digital Synthesis (DDS) and into Phase-Locked Loop (PLL) |
| param_radar.f0 | Start frequency of chirp out of DDS and into PLL |
| param_radar.fmult | PLL frequency multiplication factor |
The range resolutions for several indices of refraction are shown in Table 12.
| Index of Refraction | Range Resolution (cm) | Medium |
|---|---|---|
| 1 | 5.0 | Air |
| sqrt(1.53) | 4.0 | Snow |
| sqrt(3.15) | 2.8 | Solid Ice |
The index of refraction can be approximated by Equation 3:
(Equation 3)
Where:
| Variable | Description |
|---|---|
| ρsnow | Density of the snow in grams per cm3 |
In the data, a dielectric of 1.53 is used which corresponds to a snow density of 0.3 g per cm3 (Warren 1999).
In the along-track dimension, the raw data before any hardware or software coherent averages have a resolution derived in the same manner as the cross-track resolution. However, a basic form of focusing is applied called unfocussed Synthetic Aperture Radar (SAR) processing, also known as stacking or coherent averaging. If all affects are accounted for, the data may be coherently averaged to the SAR aperture length defined by Equation 4.
(Equation 4)
Where:
| Variable | Description |
|---|---|
| H | Height above ground level |
| λc | Wavelength at the center frequency |
For H = 500 m, data may be averaged to a length of 4.3 m. The resolution turns out to be approximately equal to this with the actual resolution definition given below. However, these data are only coherently averaged 16 times which includes both hardware and software averaging, and decimated by this same amount. At a platform speed of 140 m/s this is an aperture length, L, of 1.12 m. The sample spacing is likewise 1.12 m. The actual resolution is substantially less fine. The approximation is given by Equation 5.
(Equation 5)
Where:
| Variable | Description |
|---|---|
| H | Height above ground level |
| λc | Wavelength at the center frequency |
| L | Aperture length |
For H = 500 m, the along-track resolution is 16.7 m.
A 1 range-bin by 5 along-track-range-line boxcar filter is applied to the power detected data and then decimated in the along-track by 5 so the data product has an along-track sample spacing of 5.6 m.
For a smooth or quasi-specular target, for example internal layers, the primary response is from the first Fresnel zone. Therefore, the directivity of specular targets effectively creates the appearance of a cross-track resolution equal to this first Fresnel zone. The first Fresnel zone is a circle with diameter given by Equation 6.
(Equation 6)
Where:
| Variable | Description |
|---|---|
| H | Height above the air/ice interface |
| T | Depth in ice of the target |
| λc | Wavelength at the center frequency |
Table 17 gives the cross-track resolution for this case.
| Center Frequency (MHz) | Cross-track Resolution H = 500 m T = 0 m |
|---|---|
| 3500 | 8.7 |
For a rough surface with no appreciable layover, the cross-track resolution will be constrained by the pulse-limited footprint, approximated in Equation 7.
(Equation 7)
Where:
| Variable | Description |
|---|---|
| H | Height above the air/ice interface |
| T | Depth in ice of the target |
| c | Speed of light in a vacuum |
| kt | kt = 1.5 due to the application of a hanning time-domain window to reduce the range sidelobes of the chirped transmit waveform |
| B | Bandwidth in radians |
Table 19 gives the cross-track resolution with windowing.
| Bandwidth (MHz) | Cross-track Resolution H = 500 m T = 0 m |
|---|---|
| 4500 | 14.1 |
The antenna installed in the bomb bay of the P-3 aircraft and the wing roots of the DC-8 aircraft is an ETS Lindgren 3115. The E-plane of the antenna is aligned in the along-track for the P-3 and in the cross-track for the DC-8. The approximate beamwidths are 45 degrees in along-track and 45 degrees in cross-track. The footprint is a function of range as shown in Equation 8.
(Equation 8)
Where:
| Variable | Description |
|---|---|
| H | Height above ground level. For H = 500 m, the footprint is 414 m in along-track and 414 m in cross-track. |
| β | Beamwidth in radians |
The trajectory data used for this data release was from a basic GPS receiver. Lever arm and attitude compensation has not been applied to the data.
The following processing steps are performed by the data provider.
The purpose of the elevation compensation, when applied, is to remove the large platform elevation changes to make truncation more effective. The process is not designed to perform precision elevation compensation and is probably not sufficient for scientific analysis. The following steps are performed.
As of June 10, 2011, all of the 2009 Greenland binary files were replaced. The GPS time written in the headers of the previously published binary files contained an error. The error is corrected in the replacement binary files.
The CReSIS accumulation, snow, and kuband data acquisition systems have a known issue with radar data synchronization with GPS time. When the radar system is initially turned on, the radar system acquires UTC time from the GPS National Marine Electronics Association (NMEA) string. If this is done too soon after the GPS receiver has been turned on, the NMEA string sometimes returns GPS time rather than UTC time. GPS time is 15 seconds ahead of UTC time during this field season. The corrections for the whole day must include the offset (-15 second correction). GPS corrections have been applied to all of the data using a comparison between the accumulation, snow, and kuband radars which have independent GPS receivers. A comparison to geographic features and between ocean surface radar return and GPS elevation is also made to ensure GPS synchronization. GPS time corrections are given in the vector worksheet of the parameter spreadsheet.
As described on the CReSIS Sensors Development Radar Web site, the ultra-wideband radar operates over the frequency range from 2 to 8 GHz to map near-surface internal layers in polar firn with fine vertical resolution. The radar also has been used to measure thickness of snow over sea ice. Information about snow thickness is essential to estimate sea ice thickness from ice freeboard measurements performed with satellite radar and laser altimeters. This radar has been successfully flown on NASA P-3 and DC-8 aircraft.
Carrara, W. G., R. S. Goodman, and R. M. Majewski. 1995. Spotlight Synthetic Aperture Radar: Signal Processing Algorithms, Artech House, Norwood, MA, pp. 26-31.
Kanagaratnam, P., T. Markus, V. Lytle, B. Heavey, P. Jansen, G. Prescott, P. Gogineni. 2007. Ultrawideband Radar Measurements of Thickness of Snow Over Sea Ice, IEEE Transactions on Geoscience and Remote Sensing, 45(9): 2715-2724.
Panzer, B., C. Leuschen, A. Patel, T. Markus, and P. S. Gogineni. 2010. Ultra-wideband Radar Measurements of Snow Thickness Over Sea Ice, Geoscience and Remote Sensing Symposium (IGARSS), 2010 IEEE International, 3130-3133, 25-30 July 2010, doi: 10.1109/IGARSS.2010.5654342.
Panzer, B., C. Leuschen, W. Blake, R. Crowe, A. Patel, P. S. Gogineni, and T. Markus. 2010. Wideband Radar for Airborne Measurements of Snow Thickness on Sea Ice, Abstract C21D-01, 2010 Fall Meeting, AGU, San Francisco, California, 13-17 Dec, 2010.
Kwok, R., C. Leuschen, B. Panzer, A. Patel, N. T. Kurtz, T. Markus, B. Holt, and P. S. Gogineni. 2010. Radar Surveys of Snow Depth Over Arctic Sea Ice During Operation IceBridge, Abstract C21D-02, 2010 Fall Meeting, AGU, San Francisco, Caliornia, 13-17 Dec, 2010.
Rodriguez-Morales, F., P. Gogineni, C. Leuschen, C. T. Allen, C. Lewis, A. Patel, L. Shi, W. Blake, B. Panzer, K. Byers, R. Crowe, L. Smith, and C. Gifford. 2010. Development of a Multi-Frequency Airborne Radar Instrumentation Package for Ice Sheet Mapping and Imaging, Proc. 2010 IEEE Int. Microwave Symp., Anaheim, CA, May 2010, 157-160.
Warren, S., I. Rigor, and N. Untersteiner. 1999. Snow Depth on Arctic Sea Ice, Journal of Climate, 12: 1814-1829.
The acronyms used in this document are listed in Table 21.
| Acronym | Description |
|---|---|
| ADC | Analog to Digital Converter |
| CIRES | Cooperative Institute for Research in Environmental Science |
| CReSIS | Center for Remote Sensing of Ice Sheets |
| CSV | Comma Separated Values |
| DC-8 | Douglas DC-8 aircraft |
| DDS | Direct Digital Synthesis |
| FFT | Fast Fourier Transform |
| FMCW | Frequency Modulated Continuous Wave |
| FTP | File Transfer Protocol |
| GPS | Global Positioning System |
| JPEG | Joint Photographic Experts Group |
| KML | Keyhole Markup Language |
| L1B | Processing Level 1B |
| MATLAB | MATrix LABoratory numerical computing file |
| NASA | National Aeronautics and Space Administration |
| NMEA | National Marine Electronics Association |
| NSF | National Science Foundation |
| NSIDC | National Snow and Ice Data Center |
| P-3 | Lockheed P-3B Orion aircraft |
| PLL | Phase-Locked Loop |
| URL | Uniform Resource Locator |
| UTC | Universal Time Code |
| WFF | Wallops Flight Facility |
27 March 2013
.....
http://nsidc.org/data/docs/daac/icebridge/irsno1b/index.html