Search

Abstract Only a few localised ice streams drain most of the ice from the Greenland Ice Sheet. Thus, understanding ice stream behaviour and its temporal variability is crucially important to predict future sea-level change. The interior trunk of the 700 km-long North-East Greenland Ice Stream (NEGIS) is remarkable due to the lack of any clear bedrock channel to explain its presence. Here, we present a 3-dimensional analysis of the folding and advection of its stratigraphic horizons, which shows that the localised flow and shear margins in the upper NEGIS were fully developed only ca 2000 years ago. Our results contradict the assumption that the ice stream has been stable throughout the Holocene in its current form and show that upper NEGIS-type development of ice streaming, with distinct shear margins and no bed topography relationship, can be established on time scales of hundreds of years, which is a major challenge for realistic mass-balance and sea-level rise projections.

Abstract The dynamic mass loss of ice sheets constitutes one of the biggest uncertainties in projections of ice-sheet evolution. One central, understudied aspect of ice flow is how the bulk orientation of the crystal orientation fabric translates to the mechanical anisotropy of ice. Here we show the spatial distribution of the depth-averaged horizontal anisotropy and corresponding directional flow-enhancement factors covering a large area of the Northeast Greenland Ice Stream onset. Our results are based on airborne and ground-based radar surveys, ice-core observations, and numerical ice-flow modelling. They show a strong spatial variability of the horizontal anisotropy and a rapid crystal reorganisation on the order of hundreds of years coinciding with the ice-stream geometry. Compared to isotropic ice, parts of the ice stream are found to be more than one order of magnitude harder for along-flow extension/compression while the shear margins are potentially softened by a factor of two for horizontal-shear deformation.

During the past 20 years, multi-channel radar emerged as a key tool for deciphering an ice sheet's internal architecture. To assign ages to radar reflections and connect them over large areas in the ice sheet, the layer genesis has to be understood on a microphysical scale. Synthetic radar trace modelling based on the dielectric profile of ice cores allows for the assignation of observed physical properties’ variations on the decimetre scale to radar reflectors extending from the coring site to a regional or even whole-ice-sheet scale. In this paper we rely on the available dielectric profiling data of the northern Greenland deep ice cores: NGRIP, NEEM and EGRIP. The three records are well suited for assigning an age model to the stratigraphic radar-mapped layers, and linking up the reflector properties to observations in the cores. Our modelling results show that the internal reflections are mainly due to conductivity changes. Furthermore, we deduce fabric characteristics at the EGRIP drill site from two-way-travel-time differences of along and across-flow polarized radarwave reflections of selected horizons (below 980 m). These indicate in deeper parts of the ice column an across-flow concentrated c-axis fabric.

Crary Ice Rise formed after the Ross Ice Shelf re-grounded ~1 kyr BP. We present new ice-penetrating radar data from two systems operating at center frequencies of 7 and 750 MHz that confirm the ice rise is composed of a former ice shelf buried by subsequent accumulation. Stacks of englacial diffraction hyperbolas are present almost everywhere across the central ice rise and extend up to ~350 m above the bed. In many cases, bed reflections beneath the diffraction hyperbolas are obscured for distances up to 1 km. Waveform modeling indicates that the diffraction hyperbolas are likely caused by marine ice deposits in former basal crevasses and rifts. The in-filling of rifts and basal crevasses may have strengthened the connection between the ice rise and the surrounding ice shelf, which could have influenced local and regional ice dynamics. Three internal reflection horizons mark the upper limit of disturbed ice and diffraction hyperbolas in different sections of the ice rise, indicating at least three stages of flow stabilization across the ice rise. A surface lineation visible in MODIS imagery corresponds spatially to deepening and strong deformation of these layers, consistent with the characteristics of former grounding lines observed elsewhere in Antarctica.

In this study, our goal is to track internal ice layers on the Snow Radar data collected by NASA Operation IceBridge. We examine the application of deep learning methods on radar data gathered from polar regions. Artificial intelligence techniques have displayed impressive success in many practical fields. Deep neural networks owe their success to the availability of massive labeled data. However, in many real-world problems, even when a large dataset is available, deep learning methods have shown less success, due to causes such as lack of a large labeled dataset, presence of noise in the data or missing data. In our radar data, the presence of noise is one of the main obstacles in utilizing popular deep learning methods such as transfer learning. Our experiments show that if the neural network is trained to detect contours of objects in electro-optical imagery, it can only track a low percentage of contours in radar data. Fine-tuning and further training do not provide any better results. However, we show that selecting the right model and training it on the radar imagery from the start yields far better results.

Monitoring properties of ice sheets in polar regions is one of the main challenges in glaciology. There is a large amount of heterogeneous radar data from the polar regions that have been gathered through expensive missions. However, retrieving meaningful information from this large volume of data is still a great challenge. With the advancement of machine learning techniques in recent years, many scientists are eager to take advantage of these algorithms and techniques to explore and mine Arctic and Antarctic data. These advancements, however, have happened mainly in the area of supervised learning where the models are data hungry and require large amounts of annotated data. Generating simulated data can be an effective and inexpensive approach to provide large labeled datasets for training machine learning models. In this work, we explore two approaches to simulate arctic snow radar echogram images, namely a radar scattering physics based approach combined with some statistical measures and a purely data-driven approach based on a conditional generative adversarial network. Using several image comparison metrics, we analyze the utility of both methods for the purpose of simulating echograms. Our results show that the physics simulator generates images with good structural similarities, while the purely data-driven approach achieves better textural similarities for simulated image. Finally, we also show that by augmenting our real dataset by the simulated echograms, we can improve our deep learning model for tracking internal layers of snow.

The Northeast Greenland Ice Stream (NEGIS) is an important dynamic component for the total mass balance of the Greenland ice sheet, as it reaches up to the central divide and drains 12% of the ice sheet. The geometric boundary conditions and in particular the nature of the subglacial bed of the NEGIS are essential to understand its ice flow dynamics. We present a record of more than 8000 km of radar survey lines of multi-channel, ultra-wideband radio echo sounding data covering an area of 24 000 km2, centered on the drill site for the East Greenland Ice-core Project (EGRIP), in the upper part of the NEGIS catchment. Our data yield a new detailed model of ice-thickness distribution and basal topography in the region. The enhanced resolution of our bed topography model shows features which we interpret to be caused by erosional activity, potentially over several glacial–interglacial cycles. Off-nadir reflections from the ice–bed interface in the center of the ice stream indicate a streamlined bed with elongated subglacial landforms. Our new bed topography model will help to improve the basal boundary conditions of NEGIS prescribed for ice flow models and thus foster an improved understanding of the ice-dynamic setting.

This paper provides an update and overview of the Center for Remote Sensing of Ice Sheets (CReSIS) radars and platforms, including representative results from these systems. CReSIS radar systems operate over a frequency range of 14–38 GHz. Each radar system's specific frequency band is driven by the required depth of signal penetration, measurement resolution, allocated frequency spectra, and antenna operating frequencies (often influenced by aircraft integration). We also highlight recent system advancements and future work, including (1) increasing system bandwidth; (2) miniaturizing radar hardware; and (3) increasing sensitivity. For platform development, we are developing smaller, easier to operate and less expensive unmanned aerial systems. Next-generation platforms will further expand accessibility to scientists with vertical takeoff and landing capabilities.

We developed a high-performance, multichannel, ultra-wideband radar system for measurements of the base and interior of the East Antarctic Ice Sheet. We designed the radar to be of high power (4000-W peak) yet portable and to be able to operate with 60-MHz bandwidth at a center frequency of 200 MHz, providing high sensitivity and fine vertical resolution relative to current technology. We used the radar to perform extensive measurements as a part of a multinational collaboration. We collected data onboard a tracked vehicle outfitted with an array of high-gain antennas. We sounded 2- to 3-km thick ice near Dome Fuji. Preliminary ice thickness data match those obtained via semicoincident measurements performed with a different surface-based pulse-modulated radar system operated during the same field campaign, as well as previous airborne measurements. In addition, we mapped internal reflection horizons with fine vertical resolution from 300 m below the ice surface to ~100 m above the bed. In this article, we provide a detailed overview of the radar instrument design, implementation, and field measurement setup. We present sample data to illustrate its capabilities and discuss how the data collected with it will be valuable for the assessment of promising drilling sites to recover ice cores that are 0.9-1.5 million years old.

Climate change is extensively affecting ice sheets resulting in accelerating mass loss in recent decades. Assessment of this reduction and its causes is required to project future ice mass loss. Annual snow accumulation is an important component of the surface mass balance of ice sheets. While in situ snow accumulation measurements are temporally and spatially limited due to their high cost, airborne radar sounders can achieve ice sheet wide coverage by capturing and tracking annual snow layers in the radar images or echograms. In this paper, we use deep learning to uniquely identify the position of each annual snow layer in the Snow Radar echograms taken across different regions over the Greenland ice sheet. We train with more than 15,000 images generated from radar echograms and estimate the thickness of each snow layer within a mean absolute error of 0.54 to 7.28 pixels, depending on dataset. A highly precise snow layer thickness can help improve weather models and, thus, support glaciological studies. Such a well-trained deep learning model can be used with ever-growing datasets to aid in the accurate assessment of snow accumulation on the dynamically changing ice sheets.

Ice loss measurements around the periphery of the Greenland Ice Sheet can provide key information on the response to climate change. Here we use the excellent spatial and temporal coverage provided by the European Space Agency (ESA) CryoSat satellite, together with NASA airborne Operation IceBridge and automatic weather station data, to study the influence of changing conditions on the bias between the height estimated by the satellite radar altimeter and the ice sheet surface. Surface and near-surface conditions on the ice sheet periphery change with season and geographic position in a way that affects the returned altimeter waveform and can therefore affect the estimate of the surface height derived from the waveform. Notwithstanding the possibility of a varying bias between the derived and real surface, for the lower accumulation regions in the western and northern ice sheet periphery (

Ice thickness and bed topography of fast-flowing outlet glaciers are large sources of uncertainty for the current ice sheet models used to predict future contributions to sea-level rise. Due to a lack of coverage and difficulty in sounding and imaging with ice-penetrating radars, these regions remain poorly constrained in models. Increases in off-nadir scattering due to the highly crevassed surfaces, volumetric scattering (due to debris and/or pockets of liquid water), and signal attenuation (due to warmer ice near the bottom) are all impediments in detecting bed-echoes. A set of high-frequency (HF)/very high-frequency (VHF) radars operating at 14 MHz and 30–35 MHz were developed at the University of Kansas to sound temperate ice and outlet glaciers. We have deployed these radars on a small unmanned aircraft system (UAS) and a DHC-6 Twin Otter. For both installations, the system utilized a dipole antenna oriented in the cross-track direction, providing some performance advantages over other temperate ice sounders operating at lower frequencies. In this paper, we describe the platform-sensor systems, field operations, data-processing techniques, and preliminary results. We also compare our results with data from other ice-sounding radars that operate at frequencies both above (Center for Remote Sensing of Ice Sheets (CReSIS) Multichannel Coherent Depth Sounder (MCoRDS)) and below (Jet Propulsion Laboratory (JPL) Warm Ice Sounding Explorer (WISE)) our HF/VHF system. During field campaigns, both unmanned and manned platforms flew closely spaced parallel and repeat flight lines. We examine these data sets to determine image coherency between flight lines and discuss the feasibility of forming 2D synthetic apertures by using such a mission approach.

The Center for Remote Sensing of Ice Sheets (CReSIS) at the University of Kansas has collected approximately 1000 terabytes (TB) of radar depth sounding data over the Arctic and Antarctic ice sheets since 1993 in an effort to map the thickness of the ice sheets and ultimately understand the impacts of climate change and sea level rise. In addition to data collection, the storage, management, and public distribution of the dataset are also primary roles of the CReSIS. The Open Polar Server (OPS) project developed a free and open source infrastructure to store, manage, analyze, and distribute the data collected by CReSIS in an effort to replace its current data storage and distribution approach. The OPS infrastructure includes a spatial database management system (DBMS), map and web server, JavaScript geoportal, and MATLAB application programming interface (API) for the inclusion of data created by the cryosphere community. Open source software including GeoServer, PostgreSQL, PostGIS, OpenLayers, ExtJS, GeoEXT and others are used to build a system that modernizes the CReSIS data distribution for the entire cryosphere community and creates a flexible platform for future development. Usability analysis demonstrates the OPS infrastructure provides an improved end user experience. In addition, interpolating glacier topography is provided as an application example of the system.