Search

During the concluding phase of the NASA Operation IceBridge (OIB), we successfully completed two airborne measurement campaigns (in 2018 and 2021, respectively) using a compact S and C band radar installed on a Single Otter aircraft and collected data over Alaskan mountains, ice fields, and glaciers. This paper reports seasonal snow depths derived from radar data. We found large variations in seasonal radar-inferred depths with multi-modal distributions assuming a constant relative permittivity for snow equal to 1.89. About 34 % of the snow depths observed in 2018 were between 3.2 and 4.2 m, and close to 30 % of the snow depths observed in 2021 were between 2.5 and 3.5 m. We observed snow strata in ice facies, combined percolation and wet-snow facies, and dry-snow facies from radar data and identified the transition areas from wet-snow facies to ice facies for multiple glaciers based on the snow strata and radar backscattering characteristics. Our analysis focuses on the measured strata of multiple years at the caldera of Mount Wrangell (K'elt'aeni) to estimate the local snow accumulation rate. We developed a method for using our radar readings of multi-year strata to constrain the uncertain parameters of interpretation models with the assumption that most of the snow layers detected by the radar at the caldera are annual accumulation layers. At a 2004 ice core and 2005 temperature sensor tower site, the locally estimated average snow accumulation rate is ∼2.89 m w.e. a−1 between the years 2003 and 2021. Our estimate of the snow accumulation rate between 2005 and 2006 is 2.82 m w.e. a−1, which matches closely to the 2.75 m w.e. a−1 inferred from independent ground-truth measurements made the same year. The snow accumulation rate between the years 2003 and 2021 also showed a linear increasing trend of 0.011 m w.e. a−2. This trend is corroborated by comparisons with the surface mass balance (SMB) derived for the same period from the regional atmospheric climate model MAR (Modèle Atmosphérique Régional). According to MAR data, which show an increase of 0.86 ∘C in this area for the period of 2003–2021, the linear upward trend is associated with the increase in snowfall and rainfall events, which may be attributed to elevated global temperatures. The findings of this study confirmed the viability of our methodology, as well as its underlying assumptions and interpretation models.

The North East Greenland Ice Stream clearly stands out in the surface velocity field of the ice flow of Greenland, with its sharp and narrow shear margins visible in the flow field almost up to the central divide. While the current extent and strength of the streaming can be determined from remotely sensed velocities of the ice surface, it is less known how the ice stream is affecting the deeper layers of ice in its catchment area, and how it may have evolved over time. The deformation of the ice due to streaming can be made visible by mapping the distortion of the isochronous stratigraphy of the ice. This has been done by an airborne radar survey centering on the location of the EGRIP drilling camp, carried out with the ultra wide band radar system (AWI UWB). The dense grid of profiles arranged mainly perpendicular to the ice flow reveals the imprint that the strong shearing leaves within the layering of the ice. Although the layers are tightly folded and distorted within the shear zones, it is possible to continuously trace reflections within the upper half of the ice column throughout the entire survey area. It can be shown that the intensity of the folding is linked to the strain rate field derived from the surface velocities, and that the deformation history of the ice is preserved in the folded layers, even after it is no longer affected by high strain rates. The advection patterns of the mapped stratigraphic features reveal how the streaming of the ice and the resulting local changes of surface topography may have affected the inflow into the stream and the position of the shear margins over time.

Crystal anisotropy of ice causes slight birefringence for electromagnetic waves. At the same time, the mechanical anisotropy amounts to several orders of magnitude, thus making fabric properties highly-relevant for internal deformation. To date, bulk anisotropy of glaciers and ice sheets can be determined by geophysical methods, such as polarimetric radar, or direct sampling from ice cores. A shortcoming has been so far that changes of bulk anisotropy could mainly be inferred at single point observations, but less so as continuous profiles. Here, we exploit the effect of birefringence caused by bulk anisotropy in co-polarized airborne radar data to determine the horizontal anisotropy across the North-East Greenland Ice Stream. We base our analysis on the fact that birefringence causes a second-order effect on radar amplitudes, which leads to a beat frequency in the low and medium frequency range (O(100 kHz)), which is proportional to the horizontal anisotropy. Complementing our radar analysis with direct fabric and dielectric property observations we can constrain the range of all three fabric eigenvalues as a function of space across and along the ice stream. Finally, we assess the effect of the inferred fabric distribution on the overall ice rheology in the context of ice stream dynamics and compare it with numerical model results. Our overall approach has the advantage that it can be applied to co-polarized radar systems, as commonly used in profiling surveys, and does not require dedicated cross-polarized radar set-up. This provides the opportunity to revisit older data, especially from Greenland and Antarctica, to map fabric anisotropy in ice-dynamically interesting regions.

We present an ice-penetrating radar data set which consists of 26 internal reflecting horizons (IRHs) that cover the entire Dome C area of the East Antarctic plateau, the most extensive to date in the region. This data set uses radar surveys collected over the space of 10 years, starting with an airborne international collaboration in 2008 to explore the region, up to the detailed ground based surveys in support of the Beyond EPICA – Oldest Ice (BE-OI) European Consortium. Through direct correlation with the EPICA-DC ice core, we date 19 IRHs that span the past four glacial cycles, from the beginning of the Holocene to over 350 ka. We indirectly date and provide stratigraphic information for seven older IRHs using an 1-D ice flow inverse model, going back to an estimated 700 ka. Depth and age uncertainties are quantified for all IRHs and provided as part of the data set. The IRH data set presented in this study is available at the U.S. Antarctic Program Data Center (USAP-DC) (https://doi.org/10.15784/601411, Cavitte et al., 2020) and represents a contribution to the SCAR AntArchitecture program.

Ice shelf buttressing plays a critical role in the long-term stability of ice sheets. The underlying bathymetry and cavity thickness therefore is a key to accurate models of future ice sheet evolution. However, direct observation of sub-ice shelf bathymetry is time consuming, logistically risky, and in some areas simply not possible, meaning there is a blind-spot in our understanding of this key system. Here we use airborne gravity anomaly data to provide new estimates of sub-ice shelf bathymetry outboard of the rapidly changing West Antarctic Thwaites Glacier, and beneath the adjacent Dotson and Crosson Ice Shelves. These regions are of especial interest as the low-lying inland reverse slope of the Thwaites glacier system makes it vulnerable to collapse through marine ice sheet instability, with rapid grounding-line retreat observed since 1993 suggesting this process may be underway. Our results confirm a major marine channel > 800 m deep extends to the front of Thwaites Glacier, while the adjacent ice shelves are underlain by more complex bathymetry. Comparison of our new bathymetry with ice shelf draft reveals that ice shelves formed since 1993 comprise a distinct population where the draft conforms closely to the underlying bathymetry, unlike the older ice shelves which show a more uniform depth of the ice base. This indicates that despite rapid basal melting in some areas, these “new” ice shelves are not yet in equilibrium with the underlying ocean system. We propose qualitative models of how this transient ice-shelf configuration may have developed, but further investigation is required to constrain the longevity and full impact of these newly recognised systems.

Ice shelves play a critical role in the long-term stability of ice sheets through their buttressing effect. The underlying bathymetry and cavity thickness are key inputs for modelling future ice sheet evolution. However, direct observation of sub-ice-shelf bathymetry is time-consuming, logistically risky, and in some areas simply not possible. Here we use new compilations of airborne and marine gravity, radar depth sounding, and swath bathymetry to provide new estimates of sub-ice-shelf bathymetry outboard of the rapidly changing West Antarctic Thwaites Glacier and beneath the adjacent Dotson and Crosson ice shelves. This region is of special interest, as the low-lying inland reverse slope of the Thwaites Glacier system makes it vulnerable to marine ice sheet instability, with rapid grounding line retreat observed since 1993 suggesting this process may be underway. Our results confirm a major marine channel >800 m deep extends tens of kilometres to the front of Thwaites Glacier, while the adjacent ice shelves are underlain by more complex bathymetry. Comparison of our new bathymetry with ice shelf draft reveals that ice shelves formed since 1993 comprise a distinct population where the draft conforms closely to the underlying bathymetry, unlike the older ice shelves, which show a more uniform depth of the ice base. This indicates that despite rapid basal melting in some areas, these recently floated parts of the ice shelf are not yet in dynamic equilibrium with their retreated grounding line positions and the underlying ocean system, a factor which must be included in future models of this region's evolution.

Since 2009, the ultra-wideband snow radar on Operation IceBridge (OIB; a NASA airborne mission to survey the polar ice covers) has acquired data in annual campaigns conducted during the Arctic and Antarctic springs. Progressive improvements in radar hardware and data processing methodologies have led to improved data quality for subsequent retrieval of snow depth. Existing retrieval algorithms differ in the way the air–snow (a–s) and snow–ice (s–i) interfaces are detected and localized in the radar returns and in how the system limitations are addressed (e.g., noise, resolution). In 2014, the Snow Thickness On Sea Ice Working Group (STOSIWG) was formed and tasked with investigating how radar data quality affects snow depth retrievals and how retrievals from the various algorithms differ. The goal is to understand the limitations of the estimates and to produce a well-documented, long-term record that can be used for understanding broader changes in the Arctic climate system. Here, we assess five retrieval algorithms by comparisons with field measurements from two ground-based campaigns, including the BRomine, Ozone, and Mercury EXperiment (BROMEX) at Barrow, Alaska; a field program by Environment and Climate Change Canada at Eureka, Nunavut; and available climatology and snowfall from ERA-Interim reanalysis. The aim is to examine available algorithms and to use the assessment results to inform the development of future approaches. We present results from these assessments and highlight key considerations for the production of a long-term, calibrated geophysical record of springtime snow thickness over Arctic sea ice.

Subglacial roughness can be determined at a variety of length scales from radio-echo sounding (RES) data either via statistical analysis of topography or inferred from basal radar scattering. Past studies have demonstrated that subglacial terrain exhibits self-affine (power law) roughness scaling behaviour, but existing radar scattering models do not take this into account. Here, using RES data from northern Greenland, we introduce a self-affine statistical framework that enables a consistent integration of topographic-scale roughness with the electromagnetic theory of radar scattering. We demonstrate that the degree of radar scattering, quantified using the waveform abruptness (pulse peakiness), is topographically controlled by the Hurst (roughness power law) exponent. Notably, specular bed reflections are associated with a lower Hurst exponent, with diffuse scattering associated with a higher Hurst exponent. Abrupt waveforms (specular reflections) have previously been used as a RES diagnostic for basal water, and to test this assumption we compare our radar scattering map with a recent prediction for the basal thermal state. We demonstrate that the majority of thawed regions (above pressure melting point) exhibit a diffuse scattering signature, which is in contradiction to the prior approach. Self-affine statistics provide a generalised model for subglacial terrain and can improve our understanding of the relationship between basal properties and ice-sheet dynamics.

We present a compilation of radio-echo sounding (RES) measurements of five radar systems (AWI, BAS, CReSIS, INGV and UTIG) around the EPICA Dome C (EDC) drill site, East Antarctica. The aim of our study is to investigate the differences of the various systems in their resolution of internal reflection horizons (IRHs) and bedrock topography, penetration depth, and quality of imaging the basal layer. We address the questions of the compatibility of existing radar data for common interpretation, and the suitability of the individual systems for Oldest Ice reconnaissance surveys. We find that the most distinct IRHs and IRH patterns can be identified and transferred between most data sets. Considerable differences between the RES systems exist in range resolution and depiction of the basal layer. Considering both aspects, which we judge as crucial factors in the search for old ice, the CReSIS and the UTIG systems are the most valuable ones. In addition to the RES data set comparison we calculate a synthetic radar trace from EDC density and conductivity profiles. We identify ten common IRHs in the measured RES data and the synthetic trace. The reflection-causing conductivity sections are determined by sensitivity studies with the synthetic trace. In this way, we accomplish an accurate two-way travel time to depth conversion for the reflectors, without having to use a precise velocity-depth function that would accumulate depth uncertainties with increasing depth. The identified IRHs are assigned with the AICC2012 time scale age. Due to the isochronous character of these conductivity-caused IRHs, they are a means to extend the Dome C age structure by tracing the IRHs along the RES profiles.

Radar-inference of the bulk properties of glacier beds, most notably identifying basal melting, is, in general, derived from the basal reflection coefficient. On the scale of an ice-sheet, unambiguous determination of basal reflection is primarily limited by uncertainty in the englacial attenuation of the radio wave, which is an exponential function of temperature. Most existing radar algorithms assume stationarity in the attenuation rate, which is not feasible at an ice-sheet wide scale. Here we introduce a new framework for deriving englaical attenuation and basal reflection, and, to demonstrate its efficacy, we apply it to the Greenland Ice-Sheet. A central feature is the use of a prior Arrhenius temperature model to estimate the spatial variation in englaical attenuation as a first guess input for the radar algorithm. We demonstrate regions of solution convergence for two input temperature fields, and for independently analysed field campaigns. The coverage achieved is a trade-off with uncertainty and we propose that the algorithm can be "tuned" for discrimination of basal melt (attenuation loss uncertainty ∼ 5 dB). This is supported by our physically realistic (∼ 20 dB) range for the basal reflection coefficient. Finally, we show that the attenuation solution can be used to predict the temperature bias of thermomechanical ice-sheet models.

The catastrophic collapses of Larsen A and B ice shelves on the eastern Antarctic Peninsula have caused their tributary glaciers to accelerate, contributing to sea-level rise and freshening the Antarctic Bottom Water formed nearby. The surface of Larsen C Ice Shelf (LCIS), the largest ice shelf on the peninsula, is lowering. This could be caused by unbalanced ocean melting (ice loss) or enhanced firn melting and compaction (englacial air loss). Using a novel method to analyse eight radar surveys, this study derives separate estimates of ice and air thickness changes during a 15 year period. The uncertainties are considerable, but the primary estimate is that the surveyed lowering (0.066 ± 0.017 m yr−1) is caused by both ice loss (0.28 ± 0.18 m yr−1) and firn air loss (0.037 ± 0.026 m yr−1). Though the ice loss is much larger, ice and air loss contribute approximately equally to the lowering. The ice loss could be explained by high basal melting and/or ice divergence, and the air loss by low surface accumulation or high surface melting and/or compaction. The primary estimate therefore requires that at least two forcings caused the surveyed lowering. Mechanisms are discussed by which LCIS stability could be compromised in future, suggesting destabilisation timescales of a few centuries. The most rapid pathways to collapse are offered by a flow perturbation arising from the ungrounding of LCIS from Bawden Ice Rise, or ice-front retreat past a "compressive arch" in strain rates.

Increased surface melt over the Greenland Ice Sheet (GrIS) is now estimated to account for half or more of the ice sheet's total mass loss. Here, we show that some meltwater is stored, over winter, in buried supraglacial lakes. We use airborne radar from Operation IceBridge between 2009 and 2012 to detect buried supraglacial lakes, and we find that they were distributed extensively around the GrIS margin through that period. Buried supraglacial lakes can persist through multiple winters and are, on average, ~ 1.9 + 0.2 m below the surface. Most buried supraglacial lakes exist with no surface expression of their occurrence in visible imagery. The few buried supraglacial lakes that do exhibit surface expression have a unique visible signature associated with a darker blue color where subsurface water is located. The volume of retained water in the buried supraglacial lakes is likely insignificant compared to the total mass loss from the GrIS, but the water may have important implications locally for the development of the englacial hydrologic system and ice temperatures. Buried supraglacial lakes represent a small but year-round source of meltwater in the GrIS hydrologic system.

We present a new bed elevation dataset for Greenland derived from a combination of multiple airborne ice thickness surveys undertaken between the 1970s and 2011. Around 344 000 line kilometres of airborne data were used, with the majority of this having been collected since the year 2000, when the last comprehensive compilation was undertaken. The airborne data were combined with satellite-derived elevations for non glaciated terrain to produce a consistent bed digital elevation model (DEM) over the entire island including across the glaciated/ice free boundary. The DEM was extended to the continental margin with the aid of bathymetric data, primarily from a compilation for the Arctic. Ice shelf thickness was determined where a floating tongue exists, in particular in the north. The across-track spacing between flight lines warranted interpolation at 1 km postings near the ice sheet margin and 2.5 km in the interior. Grids of ice surface elevation, error estimates for the DEM, ice thickness and data sampling density were also produced alongside a mask of land/ocean/grounded ice/floating ice. Errors in bed elevation range from a minimum of ±6 m to about ±200 m, as a function of distance from an observation and local topographic variability. A comparison with the compilation published in 2001 highlights the improvement in resolution afforded by the new data sets, particularly along the ice sheet margin, where ice velocity is highest and changes most marked. We use the new bed and surface DEMs to calculate the hydraulic potential for subglacial flow and present the large scale pattern of water routing. We estimate that the volume of ice included in our land/ice mask would raise eustatic sea level by 7.36 m, excluding any solid earth effects that would take place during ice sheet decay.