Hostname: page-component-7c8c6479df-hgkh8 Total loading time: 0 Render date: 2024-03-29T08:05:38.207Z Has data issue: false hasContentIssue false

Sub-ice topography in Patriot Hills, West Antarctica: first results of a newly developed high-resolution FM-CW radar system

Published online by Cambridge University Press:  08 September 2017

Carlos A. Cárdenas Mansilla
Affiliation:
Dirección de Programas Antárticos, Universidad de Magallanes, Casilla 113-D, Punta Arenas, Chile E-mail: carlos.cardenas@umag.cl Centro de Estudios del Cuaternario, Casilla 113-D, Punta Arenas, Chile
Martin Jenett
Affiliation:
Hamburg University of Technology, Denickestraße 22, D-21073 Hamburg, Germany
Klaus Schünemann
Affiliation:
Hamburg University of Technology, Denickestraße 22, D-21073 Hamburg, Germany
Jürgen Winkelmann
Affiliation:
Hamburg University of Technology, Denickestraße 22, D-21073 Hamburg, Germany
Rights & Permissions [Opens in a new window]

Abstract

We present a newly developed high-resolution frequency-modulated continuous-wave (FMCW) radar system for sounding the sub-ice topography. The system, working in the frequency range from 200 to 400 MHz, was developed to measure thickness and internal layers in cold ice with a resolution better than 1 m. This system has the potential to measure accumulation rates, an important input for improved knowledge of the mass balance of polar ice sheets. First measurements for the test, calibration and optimization of the new ice sounder were made in December 2005 near the Ellsworth Mountains, specifically in Patriot Hills, West Antarctica, at 80°18′ S, 81°22′ W. The complete radar system was installed on a ski-mounted support frame, and towed by hand across the ice surface. The measurement results show the capability of this system to measure ice thickness up to 1000 m and to define internal layers within the ice body.

Type
Instruments and Methods
Copyright
Copyright © International Glaciological Society 2010

Introduction

Since the 1950s, radar systems have been an efficient tool for mapping the ice thickness of the polar ice sheets. Mainly focused on thickness measurements by means of pulse systems, previous studies (Reference NixdorfNixdorf and others, 1999; Reference Damm, Eisenburger and JenettDamm and others, 2001) with these systems suffered from low depth resolution not better than a few meters in cold ice. Former frequency-modulated continuous-wave (FM-CW) systems (Reference Forster, Davis, Rand and MooreForster and others, 1991; Reference YankielunYankielun, 1992) suffered due to heavy and expensive equipment for the measurement set-up. Reference YankielunYankielun (1992) reported on a Ka-band FM-CW system for measuring ice thickness with a depth resolution of 5 cm in ice, and Reference Forster, Davis, Rand and MooreForster and others (1991) reported on an X-band FMCW radar for snow-stratification investigation showing reflections to depths of ∼6 m. Faster digital waveform generators and increasing data acquisition capabilities led to the development of strongly improved FM-CW radar systems for the measurement of ice thickness as well as internal reflection horizons (Reference Kanagaratnam, Gogieneni, Ramasami and BraatenKanagaratnam and others, 2004).

Here we present the development of a lightweight and easy-to-handle FM-CW radar system with a bandwidth of 200 MHz for sounding polar ice caps. The objective of this system development was a depth resolution in ice better than 1 m, with the prospect of further increase in resolution. This should make the system well suited to accumulation rate measurements, a key piece of information for determining the mass balance of the polar ice caps (Reference Steinhage, Eisen and ClausenSteinhage and others, 2005).

The first test measurements were carried out at Patriot Hills, near the Ellsworth Mountains, West Antarctica, at 80°18′ S, 81°22′ W. Due to strong katabatic winds north of Patriot Hills, a large and flat blue-ice area exists, which can be used as a landing strip by large cargo aircraft to maintain logistical support of some summer camps in this region. In addition to other research groups, the Instituto Antártico Chileno (INACH) operates a summer camp at Patriot Hills, which was used as base station for the Antarctic test campaign in December 2005. The increasing ice thickness, from 0 m at the foot of the mountains to about 600 m at a distance of ∼2.5 km (Cassasa and others, 2004), and the homogeneous ice body makes this region well suited for test measurements carried out on a simple sledge with the support frame for the radar system. Several test profiles were taken to optimize the measurement capabilities of the FMCW sounder. We present technical details of the radar system, including waveform generation, broadband antennas and data acquisition. Typical measurement results of our test profiles, showing ice thickness up to 600 m and internal layers of the ice body, are provided.

FM-CW Radar System

The principle of FM-CW radar operation is depicted in Figure 1, and a block diagram of the FM-CW system is shown in Figure 2. The transmitted signal is a repetitive linear modulated frequency sweep of bandwidth B and sweep time T. The received signal, reflected from the ice surface or internal ice layers, is delayed by the two-way travel time, τ, which is usually much shorter than the sweep time. After mixing the transmitted and received signals and low-pass filtering (LP) the mixer output (Fig. 2), the intermediate frequency f IF = Δf(Fig. 1) contains the distance information of the target.

Fig. 1. Signal waveform of transmitted and received signal.

Fig. 2. Block diagram of the FM-CW system.

By means of a direct digital sythesizer (DDS), the signal generator provides a linear frequency sweep of bandwidth B and sweep time T. Both parameters can be set by means of the program of the microcontroller (μP), controlling the DDS. The ramp parameters of the frequency generator are calculated from the demands of the frequency sweep by means of respective formulas given in the data sheet of the DDS (Analog Devices, http://www.analog.com/static/imported-files/data_sheets/AD9858.pdf). A part of the filtered and amplified output signal is used as local oscillator for the subsequent mixer. After mixing the transmitted and received signals and lowpass filtering the mixer output, the intermediate frequency (IF) contains the thickness information of the ice layers.

Analog to digital (A/D) conversion of the IF signal and a subsequent fast Fourier transform (FFT) yields the spectrum of the IF signal, including distinct difference frequencies, Δf (Fig. 2), which are directly proportional to different reflector distances. The distances of the reflection horizons are determined with the well-known propagation velocity of electromagnetic waves in cold ice (v = 168 × 106 m s−1) (Reference Bogorodsky, Bentley and GudmandsenBogorodsky and others, 1985). Data acquisition, storage of data and data visualization are carried out with a laptop, running a LabVIEW® program, also setting the parameters of the A/D conversion.

The components of the FM-CW radar system used are shown in Figure 3.

Fig. 3. Components of the FM-CW radar system.

System Components

Signal generator

A block diagram of the signal generator is shown in Figure 4. The main signal source is a DDS, AD9858, which produces an analog waveform by generating a time-varying signal in digital form and then performing a digital to analog conversion (Analog Devices, http://www.analog.com/static/imported-files/data_sheets/AD9858.pdf). By this technique the DDS of the signal generator directly generates an output signal with a stepwise increase of its output frequency, thus producing a quasi-linear stepwise frequency sweep. The time-step is Δt = 7.8125 ns, and the frequency step is Δf = 7.8125 Hz, which leads to an almost continuous but highly linear frequency sweep from f min = 200 MHz to f max = 400 MHz. For the DDS, the phase-locked loop (PLL Reference BestBest, 2007) provides a stable reference frequency, f ref = 1024 MHz with a fixed phase relation to a crystal oscillator.

Fig. 4. Block diagram of the frequency generator.

From this reference clock (SYSCLK) the DDS generates all ramp parameters for the frequency-modulated output signal. With a 1024 MHz reference clock, the smallest time-step (eight clock cycles) can be calculated to 7.8125 ns, and the widths of frequency steps are determined from T and B. The frequency step, Δf, and the length of each step, Δt, are controlled by different tuning words (Analog Devices, http://www.analog.com/static/importedfiles/data_sheets/AD9858.pdf) to be calculated by the μP and written into the registers of the DDS.

The DDS increases its output frequency step by step from the starting frequency, f s = 200 MHz. After sweep time, T, also controlled by the μP, the frequency of the DDS is reset to f s, thus starting the next sweep. The output power of the DDS is 3 ± 1 dB m across the frequency range used. Typical ramp parameters of the signal generator, used for our radar system, are summarized in Table 1.

Table 1. Typical ramp parameters

Together with the bandwidth B and propagation velocity v of the radar signal in cold ice, a theoretical depth resolution can be calculated by means of the definition of radar resolution, Δr = v/2B = 0.42 m, a promising value for the investigation of internal ice layers.

Antennas

For the FM-CW radar system a broadband two-reflector antenna was developed, fed by a slotted quasi-coaxial line between both reflectors (Fig. 5). This antenna is a redesign of a similar antenna for the frequency range 100–200 MHz used for subsurface soil investigations (Reference Eisenburger, Krellmann, Lentz and TriltzschEisenburger and others, 2008). The optimization of return loss and bandwidth was done by CST MICROWAVE STUDIO®. The simulated and measured return loss of this antenna is shown in Figure 6 for comparison.

Fig. 5. Dimensions of the FM-CW antenna.

Fig. 6. Simulated and measured return loss of the FM-CW antenna.

From 230 to 370 MHz the simulated return loss (RL) is better than 15 dB, resulting in a relative bandwidth of ∼45%. The simulated gain in the band 200–400 MHz is between 8 dB at the lower end and 9 dB at the upper end of the frequency band used. There is very good agreement between simulated and measured return loss of the antenna. From the measured RL curve the relative bandwidth of the antenna at an RL level of 15 dB is calculated to 55%.

HF components

Other system components of the FM-CW system (e.g. amplifiers, bandpass filters, mixer and couplers) are commercially available components. The output power of the system is 1 W in the frequency band used. The power supply of the radar system and GPS is a 12 V, 7.5 A h battery, sufficient for ∼4 hours operating time, much longer than the operating time of the laptop used.

Data acquisition system

Sampling of the IF signal is carried out by means of a 16-bit A/D converter with a maximum sampling rate of 500 kHz. From the laptop, complete control of converter parameters and data processing is provided by the LabVIEW® program. An FFT of the time signal is performed in real time, thus giving the raw data of the measured depth profile on screen. For further data processing, the sampled time signals and the raw depth profiles are stored on the hard disk. By means of a time stamp, taken from the GPS-controlled system time of the laptop, in the trace header of each set of depth profile data, the exact position of the measurements can be determined in later post-processing of the data.

Field measurements

First measurements for the test, calibration and optimization of the new ice sounder were made at Patriot Hills (Fig. 7). Test measurements were taken from the base camp Antonio Hunneus, operated by INACH, across the blue-ice field towards the mountains, thus providing a test profile with decreasing ice depth. For comparison of the test measurements, detailed information about the ice depth was available from former measurements (Reference Casassa, Brecher, Cárdenas and RiveraCasassa and others, 1998).

Fig. 7. Geographic location of Patriot Hills.

The complete radar system was installed on a ski-mounted support frame, pulled by hand across the ice surface (Fig. 8). A well-isolated aluminum box contained the transmitter and receiver of the radar system, data acquisition and GPS receiver. The inside temperature of the box was kept well above 0°C only by the rejected heat of the laptop inside the box. Both antennas were mounted at a height of ∼2.5 m at a crossbar of the support frame. The total weight of the mounted system was ∼30 kg. With this set-up, several depth profiles of a total length of ∼15 km were measured.

Fig. 8. Radar system on the ice at Patriot Hills.

Typical results are shown in Figure 9. For these measurements the bandwidth used was 200 MHz and the sweep time was chosen to be 800 ms, resulting in ∼75 traces min−1, which is sufficient for a low ground velocity.

Fig. 9. Measured depth profile. (a) Depth profile from the camp to the mountains. (b) Zoom of the upper left corner.

The maximum measured thickness is about 600 m, with still promising signal strength of the reflected signals. The depth scale in both figures is relative to the antenna position and is calculated from the propagation velocity of the radar signal, which is well known for cold ice (v = 168 × 106 m s−1). Post-processing and verification of the field data are conducted using the commercially available REFLEXW® software.

The length of the measured profile in Figure 9a is ∼3.5 km. At the beginning of the depth profile (starting point at the base camp, ice thickness 550 m) almost to the foot of the mountains (distance 2700 m, ice thickness ∼50 m) and halfway back, the bedrock is seen.

These results agree with the ice-thickness data given by Reference Casassa, Brecher, Cárdenas and RiveraCasassa and others (1998). The reported data show a very similar depth profile across the blue-ice field towards the mountains, with a maximum ice thickness of 351 m at the center of the blue-ice field.

Figure 9b shows a zoomed cut-out of Figure 9a from the beginning of the profile. Continuous layers in the upper 100 m of the ice body clearly demonstrate a depth resolution of the radar system of ∼1 m. The onset of the reflections is seen at 60 ns, resulting from the distance between antennas and ice as well as from the additional cable length between radar system and antennas. Any interference due to the frame was not recognized in the measurement results.

Conclusions

We have presented an FM-CW radar system for investigation of polar ice sheets. In contrast to former pulse radar systems, the depth resolution of this system in cold ice is ∼1 m or better. The relative bandwidth of the antennas used is about 50% at a return loss level of 15 dB. A software-controlled DDS is used to generate the quasi-linear frequency sweep of the system. Sweep bandwidth and sweep time are easily changed by means of the software. A 16-bit A/D converter with 500 kHz sampling rate is used together with a LabVIEW® program running on a laptop for data acquisition, visualization and storage of the measured data. The measurement results show the capabilities of the system. The bedrock at 600 m depth is clearly seen, and internal layers in the upper 100 m of the ice sheet can be resolved with a resolution of ∼1 m. Further development of the system will lead to an airborne ice-sounder system, which will be flown with ‘Polar 5’, the research aircraft of the Alfred Wegener Institute for Polar and Marine Research in Bremerhaven, Germany.

Acknowledgements

Financial support from the International Bureau of the Federal Ministry of Education and Research (BMBF), Germany, under contract CHL04/Z02 is gratefully acknowledged, as well as financial support from Centro de Estudios del Cuaternario (CEQUA), Punta Arenas, Chile. We thank INACH, Punta Arenas, Chile, and Dirección de Programas Antárticos, Universidad de Magallanes, for logistical support during the test campaign in Antarctica. We also thank V. Damm and D. Eisenburger of the Federal Institute for Geosciences and Natural Resources (BGR), Hannover, Germany, for help visualizing the measurement profiles.

References

Best, R.E. 2007. Phase-locked loops: design, simulation, and applications. Sixth edition. New York, McGraw-Hill.Google Scholar
Bogorodsky, V.V., Bentley, C.R. and Gudmandsen, P.E.. 1985. Radioglaciology. Dordrecht, etc., D. Reidel Publishing Co.Google Scholar
Casassa, G., Brecher, H.H., Cárdenas, C. and Rivera, A.. 1998. Mass balance of the Antarctic ice sheet at Patriot Hills. Ann. Glaciol., 27, 130134.Google Scholar
Casassa, G., Rivera, A., Acuña, C., Brecher, H.H. and Lange, H.. 2004. Elevation change and ice flow at Horseshoe Valley, Patriot Hills, West Antarctica. Ann. Glaciol., 39, 2028.Google Scholar
Damm, V., Eisenburger, D. and Jenett, M.. 2001. Ice thickness data acquired using a helicopter-borne pulse radar system. In Proceedings of the Workshop on Remote Sensing by Low-frequency Radars, 20–21 September 2001, Naples, Italy. Paris, European Association of Remote Sensing Laboratories.Google Scholar
Eisenburger, D., Krellmann, Y., Lentz, H. and Triltzsch, G.. 2008. Stepped-frequency radar system in gating mode: an experiment as a new helicopter-borne GPR system for geological applications. In Proceedings of the International Geoscience and Remote Sensing Symposium (IGARSS 2008), 6–11 July 2008 Boston, MA. Piscataway, NJ, Institute of Electrical and Electronics Engineers, I-153–I156.Google Scholar
Forster, R.R., Davis, C.H., Rand, T.W. and Moore, R.K.. 1991. Snow-stratification investigation on an Antarctic ice stream with an X-band radar system. J. Glaciol., 37(127), 323325.Google Scholar
Kanagaratnam, P., Gogieneni, S.P., Ramasami, V. and Braaten, D.. 2004. A wideband radar for high-resolution mapping of near-surface internal layers in glacial ice. IEEE Trans. Geosci. Remote Sens., 42(3), 483490.Google Scholar
Nixdorf, U. and 6 others. 1999. The newly developed airborne radio-echo sounding system of the AWI as a glaciological tool. Ann. Glaciol., 29, 231238.CrossRefGoogle Scholar
Steinhage, D., Eisen, O. and Clausen, B.H.. 2005. Regional and temporal variation of accumulation around NorthGRIP derived from ground-penetrating radar. Ann. Glaciol., 42, 326330.Google Scholar
Yankielun, N.E. 1992. An airborne millimeter-wave FM-CW radar for thickness profiling of freshwater ice. CRREL Rep., 92–20.Google Scholar
Figure 0

Fig. 1. Signal waveform of transmitted and received signal.

Figure 1

Fig. 2. Block diagram of the FM-CW system.

Figure 2

Fig. 3. Components of the FM-CW radar system.

Figure 3

Fig. 4. Block diagram of the frequency generator.

Figure 4

Table 1. Typical ramp parameters

Figure 5

Fig. 5. Dimensions of the FM-CW antenna.

Figure 6

Fig. 6. Simulated and measured return loss of the FM-CW antenna.

Figure 7

Fig. 7. Geographic location of Patriot Hills.

Figure 8

Fig. 8. Radar system on the ice at Patriot Hills.

Figure 9

Fig. 9. Measured depth profile. (a) Depth profile from the camp to the mountains. (b) Zoom of the upper left corner.