Several decades of ice‐penetrating radar surveys of the Greenland and Antarctic ice sheets have observed numerous widespread internal reflections. Analysis of this radiostratigraphy has produced ...valuable insights into ice sheet dynamics and motivates additional mapping of these reflections. Here we present a comprehensive deep radiostratigraphy of the Greenland Ice Sheet from airborne deep ice‐penetrating radar data collected over Greenland by The University of Kansas between 1993 and 2013. To map this radiostratigraphy efficiently, we developed new techniques for predicting reflection slope from the phase recorded by coherent radars. When integrated along track, these slope fields predict the radiostratigraphy and simplify semiautomatic reflection tracing. Core‐intersecting reflections were dated using synchronized depth‐age relationships for six deep ice cores. Additional reflections were dated by matching reflections between transects and by extending reflection‐inferred depth‐age relationships using the local effective vertical strain rate. The oldest reflections, dating to the Eemian period, are found mostly in the northern part of the ice sheet. Within the onset regions of several fast‐flowing outlet glaciers and ice streams, reflections typically do not conform to the bed topography. Disrupted radiostratigraphy is also observed in a region north of the Northeast Greenland Ice Stream that is not presently flowing rapidly. Dated reflections are used to generate a gridded age volume for most of the ice sheet and also to determine the depths of key climate transitions that were not observed directly. This radiostratigraphy provides a new constraint on the dynamics and history of the Greenland Ice Sheet.
Key Points
Phase information predicts reflection slope and simplifies reflection tracing
Reflections can be dated away from ice cores using a simple ice flow model
Radiostratigraphy is often disrupted near the onset of fast ice flow
A large portion of the recent increase in the rate of mass loss from the Greenland ice sheet is from increased outlet glacier discharge along its southeastern margin. While previous investigations of ...the region’s two largest glaciers suggest that acceleration is a dynamic response to thinning and retreat of the calving front, it is unknown whether this mechanism can explain regional acceleration and what forcing is responsible for initiating rapid thinning and retreat. We examine seasonal and interannual changes in ice-front position, surface elevation and flow speed for 32 glaciers along the southeastern coast between 2000 and 2006. While substantial seasonality in front position and speed is apparent, nearly all the observed glaciers show net retreat, thinning and acceleration, with speed-up corresponding to retreat. The ratio of retreat to the along-flow stress-coupling length is proportional to the relative increase in speed, consistent with typical ice-flow and sliding laws. This affirms that speed-up results from loss of resistive stress at the front during retreat, which leads to along-flow stress transfer. Large retreats were often preceded by the formation of a flat or reverse-sloped surface near the front, indicating that subsequent retreats were influenced by the reversed bed slope. Many retreats began with an increase in thinning rates near the front in the summer of 2003, a year of record high coastal-air and sea-surface temperatures. This anomaly was driven in part by recent warming, suggesting that episodes of speed-up and retreat may become more common in a warmer climate.
We used satellite images to examine the calving behavior of Helheim and Kangerdlugssuaq Glaciers, Greenland, from 2001 to 2006, a period in which they retreated and sped up. These data show that many ...large iceberg‐calving episodes coincided with teleseismically detected glacial earthquakes, suggesting that calving‐related processes are the source of the seismicity. For each of several events for which we have observations, the ice front calved back to a large, pre‐existing rift. These rifts form where the ice has thinned to near flotation as the ice front retreats down the back side of a bathymetric high, which agrees well with earlier theoretical predictions. In addition to the recent retreat in a period of higher temperatures, analysis of several images shows that Helheim retreated in the 20th Century during a warmer period and then re‐advanced during a subsequent cooler period. This apparent sensitivity to warming suggests that higher temperatures may promote an initial retreat off a bathymetric high that is then sustained by tidewater dynamics as the ice front retreats into deeper water. The cycle of frontal advance and retreat in less than a century indicates that tidewater glaciers in Greenland can advance rapidly. Greenland's larger reservoir of inland ice and conditions that favor the formation of ice shelves likely contribute to the rapid rates of advance.
We assess ice loss and velocity changes between 1985 and 2014 of three tidewater and five-land terminating glaciers in Godthåbsfjord (Nuup Kangerlua), Greenland. Glacier thinning accounted for 43.8 ± ...0.2 km3 of ice loss, equivalent to 0.10 mm eustatic sea-level rise. An additional 3.5 ± 0.3 km3 was lost to the calving retreats of Kangiata Nunaata Sermia (KNS) and Narsap Sermia (NS), two tidewater glaciers that exhibited asynchronous behavior over the study period. KNS has retreated 22 km from its Little Ice Age (LIA) maximum (1761 AD), of which 0.8 km since 1985. KNS has stabilized in shallow water, but seasonally advects a 2 km long floating tongue. In contrast, NS began retreating from its LIA moraine in 2004–06 (0.6 km), re-stabilized, then retreated 3.3 km during 2010–14 into an over-deepened basin. Velocities at KNS ranged 5–6 km a−1, while at NS they increased from 1.5 to 5.5 km a−1 between 2004 and 2014. We present comprehensive analyses of glacier thinning, runoff, surface mass balance, ocean conditions, submarine melting, bed topography, ice mélange and conclude that the 2010–14 NS retreat was triggered by a combination of factors but primarily by an increase in submarine melting.
Jakobshavn Isbræ, a tidewater glacier that produces some of Greenland's largest icebergs and highest speeds, reached record-high flow rates in 2012 (Joughin and others, 2014). We use terrestrial ...radar interferometric observations from August 2012 to characterize the events that led to record-high flow. We find that the highest speeds occurred in response to a small calving retreat, while several larger calving events produced negligible changes in glacier speed. This non-linear response to calving events suggests the terminus was close to flotation and therefore highly sensitive to terminus position. Our observations indicate that a glacier's response to calving is a consequence of two competing feedbacks: (1) an increase in strain rates that leads to dynamic thinning and faster flow, thereby promoting destabilization, and (2) an increase in flow rates that advects thick ice toward the terminus and promotes restabilization. The competition between these feedbacks depends on temporal and spatial variations in the glacier's proximity to flotation. This study highlights the importance of dynamic thinning and advective processes on tidewater glacier stability, and further suggests the latter may be limiting the current retreat due to the thick ice that occupies Jakobshavn Isbræ’s retrograde bed.
Outlet glaciers undergo rapid spatial and temporal changes in flow velocity during calving events. Observing such changes requires both high temporal and high spatial resolution methods, something ...now possible with terrestrial radar interferometry. While a single such radar provides line-of-sight velocity, two radars define both components of the horizontal flow field. To assess the feasibility of obtaining the two-dimensional (2-D) flow field, we deployed two terrestrial radar interferometers at Jakobshavn Isbrae, a major outlet glacier on Greenland's west coast, in the summer of 2012. Here, we develop and demonstrate a method to combine the line-of-sight velocity data from two synchronized radars to produce a 2-D velocity field from a single (3 min) interferogram. Results are compared with the more traditional feature-tracking data obtained from the same radar, averaged over a longer period. We demonstrate the potential and limitations of this new dual-radar approach for obtaining high spatial and temporal resolution 2-D velocity fields at outlet glaciers.
Earth’s thick, polar ice sheets insulate the bedrock beneath them from the colder temperatures at the ice surface. Consequently, the small amount of geothermal heat that Earth’s bedrock releases can ...have a disproportionate role in controlling ice flow. Geothermal heat flux measurements beneath ice sheets are sparse, and models do not presently account for variable bed topography. Here, we present a simple statistical method for correcting geothermal heat flux models to make them consistent with known subglacial bed topography. This method is based on previous geothermal heat flux measurements across valleys and ridges, and it helps us explore the broader significance of geothermal heat flux variability beneath the Greenland and Antarctic ice sheets. Our method can also be used to estimate the influence of topography on geothermal heat flux in non-ice-covered areas.
Age-depth relations from internal layering reveal a large region of rapid basal melting in Greenland. Melt is localized at the onset of rapid ice flow in the large ice stream that drains north off ...the summit dome and other areas in the northeast quadrant of the ice sheet. Locally, high melt rates indicate geothermal fluxes 15 to 30 times continental background. The southern limit of melt coincides with magnetic anomalies and topography that suggest a volcanic origin.
Dynamic changes on Greenland outlet glaciers are a primary driver for increases in ice-sheet mass loss and its contribution to sea-level rise. One dramatic example of such change has been observed at ...Jakobshavn Isbræ, which has thinned, retreated and doubled in speed since the early 2000s. Complementary to large changes on decadal scales, we observe the glacier response on shorter timescales, driven by tidal forcing and calving events. During a 14 day period in August 2009, we documented changes in geometry and speed near the terminus. On this timescale, ice flow responds to forcing at the front from iceberg calving and ocean tides. We observe a step-increase in velocity near the terminus during a large calving event, with transient deceleration in the days following the event. A simple calving-response model explains 94–99% of variations in detrended positions at the four sites considered. During each day, variability due to tidal forcing covers 10–90% of the variability that remains after removing effects accounted for by the calving-response model. The influence of the tidal forcing on flow decays upstream with a characteristic length scale of 2 km, comparable with about two ice thicknesses.
Ice cores provide unique archives of past atmospheres and climate, but interpretation of trapped-gas records and their climatic significance has been hampered by a poor knowledge of the prevalence of ...air convection in the firn layer on top of polar ice sheets. In particular, the phasing of greenhouse gases and climate from ice cores has been obscured by a discrepancy between empirical and model-based estimates of the age difference between trapped gases and enclosing ice, which may be due to air convection. Here we show that deep air convection (>23m) occurs at a windy, near-zero-accumulation rate site in central Antarctica known informally as the Megadunes site (80.77914°S, 124.48796°E). Deep convection is evident in depth profiles of air withdrawn from the firn layer, in the observed pattern of the nitrogen isotope ratio 15N/14N, the argon isotope ratio 40Ar/36Ar, and in the mixing ratios of the anthropogenic halocarbons methyl chloroform (CH3CCl3) and HFC-134a (CH2FCF3). Transport parameters (diffusivities) were inferred and air was dated using measured carbon dioxide (CO2) and methane (CH4) mixing ratios, by comparing with the Law Dome atmospheric record, which shows that these are the oldest firn air samples ever recovered (CO2 mean age=1863AD). The low accumulation rate and the consequent intense metamorphism of the firn (due to prolonged exposure to seasonal temperature cycling) likely contribute to deep air convection via large grain size and vertical cracks that act as conduits for vigorous air motion. The Megadunes site provides a possible modern analog for the glacial conditions in the Vostok, Dome Fuji, and Dome C ice core records and a possible explanation for lower-than-expected 15N/14N ratios in trapped air bubbles at these times. A general conclusion is that very low accumulation rate causes deep air convection via its effect on firn structural characteristics.