Laura Lindzey

robots, science, code

A Brief Introduction to Ice-Penetrating Radar

July 27, 2015

TOT/JKB2d/X16a

I've been working with ice-penetrating radar data for years, and I still think it's so flipping cool that we can see through kilometers of ice using radar. This is possible because ice is mostly transparent to electromagnetic energy at radar frequencies. We see reflections at the air/ice interface, within the ice where its properties change, and at the ice/rock (or ice/water) interface.

Over the last few decades1, my advisor's group has collected airborne ice-penetrating data all over Antarctica:

UTIG flights

In order to collect this data, we mount antennas and recording equipment on a modified DC3. These are unpressurized WWII-era airplanes that have been updated with new engines. The left image below shows JKB taking off from the Casey skiway; the right image shows one of my colleagues monitoring the instruments during a flight. As the plane flies, the radar transmits 6250 pulses every second. After each pulse, it listens for reflected energy; this returned signal is what you're looking at in the radargram.2

platform's antennas airplane's racks

I spend a lot of time studying radar profiles like the one below. The x-axis is distance along the flight line, and the y-axis is height.3 In effect, we're looking at a slice of the ice sheet. White corresponds to a strong reflection, and we can clearly see the surface and bed all along the transect. This particular profile is running along the Totten Ice Shelf; it is highlighted in red on the right (zoomed-in) image below. The star indicates the start of the profile. Note how the bed signal gets weaker as we cross the grounding zone - radar energy reflects more strongly from water than from rock.

TOT/JKB2d/X16a full transect

context map zoom0 context map zoom1 context map zoom2

Totten Glacier is interesting because it is losing mass relatively quickly, as determined by satellite gravity4 and altimetry.5 Among other things, radar data is used in areas like this to better determine the location of the grounding line and to calculate the volume of ice in the drainage feeding an outlet glacier.6

Multiple Reflections

When looking at radargrams, another thing to watch out for are multiple reflections, which can be confused for the bed or a layer. These most commonly occur when the radar energy bounces back off the plane, to the surface, and back to the plane. It is possible to distinguish such a multiple from a true layer because it is always the same distance below the surface as the plane is above the surface, which means that the multiple's apparent slope is twice that of the surface's. An example of this type of multiple can be found in the Totten transect above. Below, we have data from over the Ross Ice Shelf. The ice shelf is thin enough, and the bottom echo strong enough, that we see multiples due to the radar energy bouncing between the ice's interfaces, as well as back off the plane.

MIS/JKB2e/Y25 full transect

MIS zoom 0 MIS zoom 1 MIS zoom 2

Beam Pattern

Radar data can be more difficult to interpret than laser data because the beam pattern is not anywhere near as focused. This means that the energy received at a given point in space is not necessarily a reflection from directly below the airplane, as would be most convenient/intuitive to assume. On our airplane, the along-track pattern is a single lobe with half power at 152-degrees (below, left). Cross-track, the first lobe is at 22-degrees, and is only 2.4dB weaker than the main lobe directly below the airplane (below, right).

side view of JKB front view of JKB side view of antenna pattern front view of antenna pattern

(Figures from7, showing the antenna's gain in dB as a function of angle.)

Reflections from the side

The previous radargram has an undulating surface that could be confused for the bottom of an ice shelf. However, looking at a map, we can see that this region is either open water or sea ice. The location of the ice shelf's edge roughly matches the reflector's geometry, so this is actually the side of the ice shelf reflecting energy from one of the sidelobes! Echoes off to the side aren't always so easy to interpret. Here's a more confusing transect flying along a valley:

ICP2/JKB1a/F09T01a full transect

A number of research groups are actively working on how to determine whether a signal is a side echo or originates from directly below the aircraft, but it's not yet a solved problem. Techniques that have been tried include recording from multiple antennas at the same time and comparing the phase difference,8 or comparing the radar data to a known elevation map of the surface.9

Reflections from ahead/behind

The wide beam pattern along-track also means that a strong reflector will be seen long before/after we actually fly over it. These show up as hyperbolas in the radargram, as illustrated in the cartoon below:

geometry forming hyperbola

Crevasses are very good reflectors since a corner will reflect energy along the same path that it came in on. They show up well in some regions of the Totten Ice Shelf transect above. Man-made structures also show up as strong reflectors. In the radargram below, we see South Pole Station as a set of hyperbolas. Their peaks appear below the surface reflection because this transect didn't directly cross the station - it was about a kilometer to the side.

There are existing processing techniques for dealing with a wide along-track beam angle that reduce the appearance of hyperbolas. The basic idea is that you've observed the same location from a number of different positions as you're flying over it, and can use that redundancy to figure out where the "real" reflector is,10 but the details are a topic for another post.

NPX/SJB2/NPX02a

The radargram above has striking layers which are also wonderfully useful for modelling the ice sheet. A lot of work is done to correlate layers observed in radargrams to ice properties directly measured from a drilled core. This is analogous to studying tree rings, but over a much longer time frame. There are ice cores from East Antarctica that provide a climate record for the last 800,000 years.11 Once you have this correlation, and if you can trace a continuous layer across a region, you then know the age-depth relationship for the ice sheet. The deformation of layers also contains information about how the ice has been moving.

A colleague claims that no presentation about Antarctica is complete without a penguin. They're far too small to show up in radar data, but here's a picture of a pair from Shirley Island, near Casey Station:

pair of penguins

Thanks to Mindy Preston for comments and questions about an earlier version of this post!


  1. I was looking for extra storage space in the office the other day, and ran across a fantastic filing cabinet: the top three drawers had data on 3 different generations of tape: LTO3, which we're only just now phasing out for LTO6, one that looked more like a cassette tape, and one set of huge spools that were labeled "800Bytes per Inch", and dated from the early 1980's. The bottom drawer contained a few rolls of photographic paper - still smelling of darkroom chemicals! - with radar transects printed on them, dating from 1979. 

  2. We transmit a 1us chirped pulse, with a 15MHz bandwidth centered around 60MHz. We digitize the returned signal at 50MHz, and stack 32 traces before recording to disk. The recorded signal is dechirped and stacked again to create the images I'm showing here, but those are details for another post. 

  3. More precisely, both axes are in units of time - each column of the image corresponds to the information received from a single pulse of the radar, and the the rows are how much time has passed since the transmission. 

  4. Chen, J. L., Wilson, C. R., Blankenship, D., and Tapley, B. D. Accelerated Antarctic ice loss from satellite gravity measurements. Nature Geoscience 2, 12 (2009), 859–862. 

  5. Pritchard, H. D., Arthern, R. J., Vaughan, D. G., and Edwards, L. A. Extensive dynamic thinning on the margins of the Greenland and Antarctic ice sheets. Nature 461, 7266 (2009), 971–975. 

  6. Greenbaum, J. S., Blankenship, D. D., Young, D. A., Richter, T. G., Aitken, A. R. A., Legresy, B., Schroeder, D. M., Roberts, J. L., Warner, R. C., van Ommen, T. D., and Siegert, M. J. Increased Ocean Access to Totten Glacier, East Antarctica. Nature Geos (2015). 

  7. Peters, M. E., Blankenship, D. D., Carter, S. P., Kempf, S. D., Young, D. A., and Holt, J. W. Along-track focusing of airborne radar sounding data from West Antarctica for improving basal reflection analysis and layer detection. Geoscience and Remote Sensing, IEEE Transactions on 45, 9 (2007), 2725–2736. 

  8. Gogineni, S., YAN, J., Paden, J., Leuschen, C., Li, J., Rodriguez-Morales, F., Braaten, D., Purdon, K., Wang, Z., Liu, W., et al. Bed topography of Jakobshavn Isbrae, Greenland, and Byrd Glacier, Antarctica. Journal of Glaciology 60, 223 (2014), 813. 

  9. Holt, J. W., M. E. Peters, S. D. Kempf, D. L. Morse, and D. D. Blankenship (2006), Echo source discrimination in single-pass airborne radar sounding data from the Dry Valleys, Antarctica: Implications for orbital sounding of Mars, J. Geophys. Res., 111, E06S24, doi:10.1029/2005JE002525. 

  10. Peters, M. E., Blankenship, D. D., Carter, S. P., Kempf, S. D., Young, D. A., and Holt, J. W. Along-track focusing of airborne radar sounding data from West Antarctica for improving basal reflection analysis and layer detection. Geoscience and Remote Sensing, IEEE Transactions on 45, 9 (2007), 2725–2736 

  11. Lambert, Fabrice, et al. "Dust-climate couplings over the past 800,000 years from the EPICA Dome C ice core." Nature 452.7187 (2008): 616-619.