Topographic Survey

During December 1993, topographic and geodetic measurements were carried out. The centre of the topographic dome obtained by ERS-1 was monumented and located in absolute geographic coordinates utilizing the DORIS system (Vincent, 1994). Around this point a global positioning system (GPS) survey was carried out in order to obtain a detailed topographic map of the summital area (Cefalo and others, 1994). The measurements were carried out in kinematic mode, using two Trimble 4000SEs, the first installed on a snowcat (rover), the second one (master) located close to the DORIS beacon, in order to have a good reference position for differential processing of data. To avoid passive trajectory, the survey was carried out on a “star” grid, centred on the DORIS point, composed of six rectangular isosceles triangles, with cathetus length of

20 km, rotated from each other by 60. The total length of the GPS profiles was about 400 km. In 1995, another topographic map of a larger area (contour line interval 1 m), obtained by a newly available series of ERS-1 data which was merged with 1993 GPS data, was worked out by F. Rémy and R. Cefalo (Fig. 1); this map shows the asymmetry of the dome structure and suggests a regional ire divide oriented northeast -southwest.

Radar Survey

The radar equipment was designed and built by M. R. Gorman of SPRI, Cambridge, U.K., in a development funded by the Italian Antarctic Programme (PNRA); some minor modifications and adjustments (balun, antennas, addition of GPS system, radio-frequency pulse alimentation) were carried out by an Italian team (F. Corbelli and A. Passerini). The radar utilises a digital pulse generator to generate a 60 MHz pulse of switchable duration (0.3 or 1.0 μs), at about 1 mW peak power. This is amplified by a solid-state power amplifier to more than 2 kW peak power. The receiver consists of a solid-slate low-noise preamplifier, followed by a switchable-bandwidth successive-detection logarithmic amplifier. The dynamic range of this amplifier is 80dB, and the output is digitised to 8 bit resolution, acquiring 512 samples spaced by 100 ns, for a total time range of. 51.2 μs. The digitised data pass via an RS232 serial link to the control computer (IBM-compatible 80386SX, 40 MB hard disk).

The calibration test of the overall system shows a system performance of 150-160dB (not including antenna gain), depending on pulse length and bandwidth settings. The same radar equipment was deployed both by aircraft and by snowcat. For airborne survey we employed a Twin Otter aircraft; a GPS Trimble 4000SSE system (L1 and L2 frequency), with geodetic antenna mounted on the fuselage, was installed and linked to the radar. The aircraft radar antennas used were folded dipoles under each wing, with one used for transmit, the other for receive. Antenna gain is estimated at 8.2 dB. For ground survey we used a snowcat, modified by a PNRA team (R. Buccolini and L. Blasi) and suited for GPS antenna and radar installation. The radar antenna was installed on a non-metallic sledge (designed and built by P. Godon, Institut Français pour la Recherche et la Technologie Polaire, Plouzane) to be towed by the snow-eat. A master GPS station, synchronised with the rover, was installed at Dome C camp to enable differential correction of data.

Radar measurements were carried oui on three grids (Fig. 2). The “big grid”, used for airborne survey, is an 80 km ? 120 km rectangular grid, with a spacing of 10 km between flight-lines, centred on the summit point, with the longer side parallel to the apparent dome axis. The “small grid”, co-centred with the previous one, is a 50 km ? 50 km square grid (line spacing 10 km), shifted by 5 km from the lines of the “big grid” in order to obtain in the central part of the area a final coverage with a spacing of 5 km between the flight-lines. The geometry of the “star grid” for ground-based measurements was similar to that used during the 1993 topographic survey. The total airborne survey length was about 2800 km, the ground-based one about 300 km. The grids were designed to have a large number of crossover check points between all the radar lines; we totalled about 240 crossover points.

Fig. 1.

Surface topography by ERS-1 and GPS.

Radar Data and Processing

The radar survey consisted of 34 profiles by air and 21 by ground. The airborne survey was achieved in four flights, with a cruising speed of 100-120 knots (185-220 km h^-1), at 1000 ft (305 m) average elevation over the surface (controlled by radar altimeter). Data were acquired using a pulse length of 1 μs, with an acquisition rate of 1 record/3 s, that is, 1 record/150—185 m. About 14 000 records were acquired by air. The ground-based survey was carried out with a travelling speed of 4-10 knots (7.5-18.5 km h^-1), with an acquisition rate of 1 record/6-15 m; in total, about 30 000 records were acquired. All the radar data were located by a synchronised GPS system with a precision in x, y coordinates of better than ± 20 m (only pseudo-range differential corrections were done).

All the radar profiles were processed with a software package specially designed for the new radio-echo sounder. A constant electromagnetic-wave propagation velocity of 168 m μs^-1was assumed (Robin, 1975; Patterson, 1981; Bogodorosky and others, 1985). Thickness value was determined to the nearest value of digit sample on the record, with a digitising error of ±100 ns for both surface and bed reflections. This corresponds to an rms uncertainty of about ±12 m in thickness. No corrections for higher velocity in the firn were calculated, which implies that all thickness values will be underestimated by very approximately 30 m. It is not possible to improve this estimate without information on firn thickness and the density profile. Good-quality radar records were obtained over all the grids, with strong and unambiguous bottom returns; measured ice thickness was 2500-4000 in. Internal layering is apparent everywhere, mostly from 800 to 2800 m depth below surface. Normally the layering follows the bed topography with reduced amplitude; a featureless lower zone of about 700-800 m occurs over the deeper bed areas. The featureless zone is reduced to less than 300 m where the bed topography is higher. Four subglacial lakes were detected. One of them, at 74.91°S, 124.65°E, corresponds to Lake No. 25 of the inventory produced by Siegert and others (1996); the lake closest to the dome summit (74.95°S, 123.76° E) is about 20 km northeast. Examples of airborne radar lines are shown in Figure 3a-c. A ground-based radar line, located near the centre, is shown in Figure 4. Here the flat bed morphology appears to vary within a range of less than 100 m, at the mean ice depth of 3250 m. Horizontal and continuous layering is observed down to 2200 m depth; below this, until about 2500 m, the layering appears more irregular, and below this depth the internal layering disappears.

Comparative ice-thickness values at each crossover point of all the airborne and ground-based radar profiles (more than 220 points) were calculated; they presented differences of less than 10 m, i.e. less than the estimated accuracy in thickness determination. A first grid with averaged values of crossing points was calculated. Over the whole area, 3800 ice-thickness values were calculated, about 1800 (spacing 1500 m) on the flight-lines and about 2000 (spacing 150 m) on ground-based ones. Different gridding methods were used; each gridding result was compared with the crossover-point dataset. The best results were obtained by the minimum curvature method that gave minimum differences with the crossover-point depth.

In order to calculate the bed topography, a new detailed surface map of the “big grid” area (minimum curvature gridding) was worked out (Cefalo and others, 1996). This map was obtained by merging the 1994 ERS-1 data with all the GPS and DORIS values. Slight differences with the Cefalo 1994 map were observed and the new summit of the dome was fixed at 75°06′06″ S, 123°23′42″ E. By subtracting thickness values from the surface topography a bed map was obtained (Fig. 5). A tridimensional view of the area with surface and bed topography is presented in Figure 6. The maps are World Geodetic System 1984 ellipsoid (WGS84) referenced, and elevation values refer to height above ellipsoid (HAE).

Conclusions

With regard to the main aim of the work, i.e. the drilling-point selection, we can conclude:

Bed topography

The plateau located over the central part of the area (at a mean elevation corresponding to zero ?A?) is bounded eastward by a well-defined system composed of a deep valley and a parallel sharp ridge oriented north-south. South and southwest it is surrounded by a complex chain of mountains reaching a maximum altitude of about +500 to +700 m. Northward the plateau is bordered by a system of smooth hills (reaching a maximum altitude of +100 m), then gently deepens to -300

m without prominent obstacles.

Fig. 2.

Location of GPS tracks of radar grids. (a) GPS path “big grid”; (b) “small grid”; (c) “star grid”.

Fig. 3.

Examples of airborne radar records. For flight location see Figure 2a.

Fig. 4.

Example of ground-based radar records. Fur profile location see Figure 2c.

Fig. 5.

Bed topography. Height above WGS84.

Fig. 6.

Tridimensional view of surface and bed topography. Height above WGS 84 (bold line an bed topography indicates ellipsoid elevation).

The north south valley, with a width of about 15 km, deepens northward from -200 to -600 m; its eastern side is steeper than the western one. The altitude of the main four summits of the north-south ridge ranges from +400 to +700 m. The relative elevations of valley plateau and valley ridge are respectively about 400 and 1000 m; the morphologic appearance suggests a tectonic-fault angle depression.