Using ground-penetrating radar, topography and classification of vegetation to model the sediment and active layer thickness in a periglacial lake catchment, Western Greenland, link to shapefiles

The geometries of a catchment constitute the basis for distributed physically based numerical modeling of different geoscientific disciplines. In this paper results from ground-penetrating radar (GPR) measurements, in terms of a 3D model of total sediment thickness and active layer thickness in a periglacial catchment in western Greenland, is presented. Using the topography, thickness and distribution of sediments is calculated. Vegetation classification and GPR measurements are used to scale active layer thickness from local measurements to catchment scale models. Annual maximum active layer thickness varies from 0.3 m in wetlands to 2.0 m in barren areas and areas of exposed bedrock. Maximum sediment thickness is estimated to be 12.3 m in the major valleys of the catchment. A method to correlate surface vegetation with active layer thickness is also presented. By using relatively simple methods, such as probing and vegetation classification, it is possible to upscale local point measurements to catchment scale models, in areas where the upper subsurface is relatively homogenous. The resulting spatial model of active layer thickness can be used in combination with the sediment model as a geometrical input to further studies of subsurface mass-transport and hydrological flow paths in the periglacial catchment through numerical modelling.

High-resolution ice thickness and bed topography of a land-terminating section of the Greenland Ice Sheet, with links to GeoTIFFs and text files

We present ice thickness and bed topography maps with a high spatial resolution (250-500 m) of a land-terminating section of the Greenland Ice Sheet derived from ground-based and airborne radar surveys. The data have a total area of ~12 000 km^2 and cover the whole ablation area of the outlet glaciers of Isunnguata Sermia, Russell, Leverett, Ørkendalen and Isorlersuup up to the long-term mass balance equilibrium line altitude at ~1600 m above sea level. The bed topography shows highly variable subglacial trough systems, and the trough of Isunnguata Sermia Glacier is overdeepened and reaches an elevation of ~500 m below sea level. The ice surface is smooth and only reflects the bedrock topography in a subtle way, resulting in a highly variable ice thickness. The southern part of our study area consists of higher bed elevations compared to the northern part. The compiled data sets of ground-based and airborne radar surveys cover one of the most studied regions of the Greenland Ice Sheet and can be valuable for detailed studies of ice sheet dynamics and hydrology.

31 ground-penetrating radar lines of Area 3 - Beach ridge plain (high wave energy) from Maxwell Bay, King George Island, Antarctic Peninsula

During a field campaign in the Austral spring 2012 the sedimentary architecture of a polar gravel-beach system at the southern coast of Potter Peninsula (Area 3) was revealed using ground-penetrating radar (GPR, Geophysical Survey Systems, Inc. SIR-3000). 31 profiles were collected using a mono-static 200 MHz antenna operated in common offset mode. Trace increment was set to 0.05 m. A differential global-positioning system (dGPS, Leica GS09) was used to obtain topographical information along the GPR lines. GPR data are provided in RADAN-Format, dGPS coordinates are provided in ascii format; projection is UTM (WGS 84, zone 21S).

Radar profiles from South Shetland Islands/Western Antarctic Peninsula

The sedimentary architecture of polar gravel-beach ridges is presented and it is shown that ridge internal geometries reflect past wave-climate conditions. Ground-penetrating radar (GPR) data obtained along the coasts of Potter Peninsula (King George Island) show that beach ridges unconformably overlie the prograding strand plain. Development of individual ridges is seen to result from multiple storms in periods of increased storm-wave impact on the coast. Strand-plain progradation, by contrast, is the result of swash sedimentation at the beach-face under persistent calm conditions. The sedimentary architecture of beach ridges in sheltered parts of the coast is characterized by seaward-dipping prograding beds, being the result of swash deposition under stormy conditions, or aggrading beds formed by wave overtopping. By contrast, ridges exposed to high-energy waves are composed of seaward- as well as landward-dipping strata, bundled by numerous erosional unconformities. These erosional unconformities are the result of sediment starvation or partial reworking of ridge material during exceptional strong storms. The number of individual ridges which are preserved from a given time interval varies along the coast depending on the morphodynamic setting: sheltered coasts are characterized by numerous small ridges, whereas fewer but larger ridges develop on exposed beaches. The frequency of ridge building ranges from decades in the low-energy settings up to 1600 years under high-energy conditions. Beach ridges in the study area cluster at 9.5, 7.5, 5.5, and below 3.5 m above the present-day storm beach. Based on radiocarbon data, this is interpreted to reflect distinct periods of increased storminess and/or shortened annual sea-ice coverage in the area of the South Shetland Islands for the times around 4.3, c. 3.1, 1.9 ka cal BP, and after 0.65 ka cal BP. Ages further indicate that even ridges at higher elevations can be subject to later reactivation and reworking. A careful investigation of the stratigraphic architecture is therefore essential prior to sampling for dating purposes.

Present-day King George Island ice cap geometry measured by ground penetrating radar in 2008/09

We present a consistent data set for the ice thickness, the bedrock topography and the ice surface topography of the King George Island ice cap (Arctowski Icefield and the adjacent central part). The data set is composed of groundbased and airborne Ground Penetrating Radar (GPR) and differential GPS (DGPS) measurements, obtained during several field campaigns. The data set incorporates groundbased measurements in the safely accessible inner parts and airborne measurements in the heavily crevassed coastal areas of the ice cap. In particular, the inclusion of airborne GPR measurements with the 30MHz BGR-P30-System developed at the Institute of Geophysics (University of Münster) completes the picture of the ice geometry substantially. The compiled digital elevation model of the bedrock shows a rough, highly variable topography with pronounced valleys, ridges, and troughs. The mean ice thickness is approx. 238m, with a maximum value of approx. 400m in the surveyed area. Noticeable are bounded areas in the bedrock topography below sea level where marine based ice exists.

21 ground-penetrating radar lines of Area 1 - Beach ridge plain (low wave energy) from Potter Cove, King George Island, Antarctic Peninsula

During two field campaigns (Austral springs 2011 and 2012) the sedimentary architecture of a polar gravel-beach system at the western coast of Potter Peninsula (Area 1) was revealed using ground-penetrating radar (GPR, Geophysical Survey Systems, Inc. SIR-3000). 21 profiles were collected using a mono-static 200 MHz antenna operated in common offset mode. Trace increment was set to 0.05 m. A differential global-positioning system (dGPS, Leica GS09) was used to obtain topographical information along the GPR lines. GPR data are provided in RADAN-Format, dGPS coordinates are provided in ascii format; projection is UTM (WGS 84, zone 21S).

49 ground-penetrating radar lines of Area 2 - Gravel cuspate spit with gravel beach ridges from Potter Cove, King George Island, Antarctic Peninsula

During two field campaigns (Austral springs 2011 and 2012) the sedimentary architecture of a polar gravel-beach system at the southwestern coast of Potter Peninsula (Area 2) was revealed using ground-penetrating radar (GPR, Geophysical Survey Systems, Inc. SIR-3000). 49 profiles were collected using a mono-static 200 MHz antenna operated in common offset mode. Trace increment was set to 0.05 m. A differential global-positioning system (dGPS, Leica GS09) was used to obtain topographical information along the GPR lines. GPR data are provided in RADAN-Format, dGPS coordinates are provided in ascii format; projection is UTM (WGS 84, zone 21S).

(Supplement 3) Radar sounding observations including bedrock and surface elevation and ice thickness in the ASB region, Antarctica

The first Cenozoic ice sheets initiated in Antarctica from the Gamburtsev Subglacial Mountains and other highlands as a result of rapid global cooling ~34 million years ago. In the subsequent 20 million years, at a time of declining atmospheric carbon dioxide concentrations and an evolving Antarctic circumpolar current, sedimentary sequence interpretation and numerical modelling suggest that cyclical periods of ice-sheet expansion to the continental margin, followed by retreat to the subglacial highlands, occurred up to thirty times. These fluctuations were paced by orbital changes and were a major influence on global sea levels. Ice-sheet models show that the nature of such oscillations is critically dependent on the pattern and extent of Antarctic topographic lowlands. Here we show that the basal topography of the Aurora Subglacial Basin of East Antarctica, at present overlain by 2-4.5 km of ice, is characterized by a series of well-defined topographic channels within a mountain block landscape. The identification of this fjord landscape, based on new data from ice-penetrating radar, provides an improved under¬standing of the topography of the Aurora Subglacial Basin and its surroundings, and reveals a complex surface sculpted by a succession of ice-sheet configurations substantially different from today's. At different stages during its fluctuations, the edge of the East Antarctic Ice Sheet lay pinned along the margins of the Aurora Subglacial Basin, the upland boundaries of which are currently above sea level and the deepest parts of which are more than 1 km below sea level. Although the timing of the channel incision remains uncertain, our results suggest that the fjord landscape was carved by at least two ice- flow regimes of different scales and directions, each of which would have over-deepened existing topographic depressions, reversing valley floor slopes.

Greenland and Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo-2), links to NetCDF files

The ocean plays an important role in modulating the mass balance of the polar ice sheets by interacting with the ice shelves in Antarctica and with the marine-terminating outlet glaciers in Greenland. Given that the flux of warm water onto the continental shelf and into the sub-ice cavities is steered by complex bathymetry, a detailed topography data set is an essential ingredient for models that address ice-ocean interaction. We followed the spirit of the global RTopo-1 data set and compiled consistent maps of global ocean bathymetry, upper and lower ice surface topographies and global surface height on a spherical grid with now 30-arc seconds resolution. We used the General Bathymetric Chart of the Oceans (GEBCO, 2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the Interna- tional Bathymetric Chart of the Southern Ocean (IBCSO) version 1. While RTopo-1 primarily aimed at a good and consistent representation of the Antarctic ice sheet, ice shelves and sub-ice cavities, RTopo-2 now also contains ice topographies of the Greenland ice sheet and outlet glaciers. In particular, we aimed at a good representation of the fjord and shelf bathymetry sur- rounding the Greenland continent. We corrected data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ and Sermilik Fjord assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79°N, we incorporated a high-resolution digital bathymetry model considering original multibeam survey data for the region. Radar data for surface topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centers of Technical University of Denmark (DTU), Operation Icebridge (NASA/NSF) and Alfred Wegener Institute (AWI). For the Antarctic ice sheet/ice shelves, RTopo-2 largely relies on the Bedmap-2 product but applies corrections for the geometry of Getz, Abbot and Fimbul ice shelf cavities.