Behavior Of Ice Sheets Research Articles

Ice sheet model dependency of the simulated Greenland Ice Sheet in the mid-Pliocene

Abstract. The understanding of the nature and behavior of ice sheets in past warm periods is important for constraining the potential impacts of future climate change. The Pliocene warm period (between 3.264 and 3.025 Ma) saw global temperatures similar to those projected for future climates; nevertheless, Pliocene ice locations and extents are still poorly constrained. We present results from the efforts to simulate mid-Pliocene Greenland Ice Sheets by means of the international Pliocene Ice Sheet Modeling Intercomparison Project (PLISMIP). We compare the performance of existing numerical ice sheet models in simulating modern control and mid-Pliocene ice sheets with a suite of sensitivity experiments guided by available proxy records. We quantify equilibrated ice sheet volume on Greenland, identifying a potential range in sea level contributions from warm Pliocene scenarios. A series of statistical measures are performed to quantify the confidence of simulations with focus on inter-model and inter-scenario differences. We find that Pliocene Greenland Ice Sheets are less sensitive to differences in ice sheet model configurations and internal physical quantities than to changes in imposed climate forcing. We conclude that Pliocene ice was most likely to be limited to the highest elevations in eastern and southern Greenland as simulated with the highest confidence and by synthesizing available regional proxies; however, the extent of those ice caps needs to be further constrained by using a range of general circulation model (GCM) climate forcings.

Seismic wave propagation in anisotropic ice – Part 1: Elasticity tensor and derived quantities from ice-core properties

Abstract. A preferred orientation of the anisotropic ice crystals influences the viscosity of the ice bulk and the dynamic behaviour of glaciers and ice sheets. Knowledge about the distribution of crystal anisotropy is mainly provided by crystal orientation fabric (COF) data from ice cores. However, the developed anisotropic fabric influences not only the flow behaviour of ice but also the propagation of seismic waves. Two effects are important: (i) sudden changes in COF lead to englacial reflections, and (ii) the anisotropic fabric induces an angle dependency on the seismic velocities and, thus, recorded travel times. A framework is presented here to connect COF data from ice cores with the elasticity tensor to determine seismic velocities and reflection coefficients for cone and girdle fabrics. We connect the microscopic anisotropy of the crystals with the macroscopic anisotropy of the ice mass, observable with seismic methods. Elasticity tensors for different fabrics are calculated and used to investigate the influence of the anisotropic ice fabric on seismic velocities and reflection coefficients, englacially as well as for the ice–bed contact. Hence, it is possible to remotely determine the bulk ice anisotropy.

Recharge of a subglacial lake by surface meltwater in northeast Greenland.

In a warming climate, surface meltwater production on large ice sheets is expected to increase. If this water is delivered to the ice sheet base it may have important consequences for ice dynamics. For example, basal water distributed in a diffuse network can decrease basal friction and accelerate ice flow, whereas channelized basal water can move quickly to the ice margin, where it can alter fjord circulation and submarine melt rates. Less certain is whether surface meltwater can be trapped and stored in subglacial lakes beneath large ice sheets. Here we show that a subglacial lake in Greenland drained quickly, as seen in the collapse of the ice surface, and then refilled from surface meltwater input. We use digital elevation models from stereo satellite imagery and airborne measurements to resolve elevation changes during the evolution of the surface and basal hydrologic systems at the Flade Isblink ice cap in northeast Greenland. During the autumn of 2011, a collapse basin about 70 metres deep and about 0.4 cubic kilometres in volume formed near the southern summit of the ice cap as a subglacial lake drained into a nearby fjord. Over the next two years, rapid uplift of the floor of the basin (which is approximately 8.4 square kilometres in area) occurred as surface meltwater flowed into crevasses around the basin margin and refilled the subglacial lake. Our observations show that surface meltwater can be trapped and stored at the bed of an ice sheet. Sensible and latent heat released by this trapped meltwater could soften nearby colder basal ice and alter downstream ice dynamics. Heat transport associated with meltwater trapped in subglacial lakes should be considered when predicting how ice sheet behaviour will change in a warming climate.

Marine ice-sheet profiles and stability under Coulomb basal conditions

AbstractThe behavior of marine-terminating ice sheets, such as the West Antarctic ice sheet, is of interest due to the possibility of rapid grounding-line retreat and consequent catastrophic loss of ice. Critical to modeling this behavior is a choice of basal rheology, where the most popular approach is to relate the ice-sheet velocity to a power-law function of basal stress. Recent experiments, however, suggest that near-grounding line tills exhibit Coulomb friction behavior. Here we address how Coulomb conditions modify ice-sheet profiles and stability criteria. The basal rheology necessarily transitions to Coulomb friction near the grounding line, due to low effective stresses, leading to changes in ice-sheet properties within a narrow boundary layer. Ice-sheet profiles ‘taper off’ towards a flatter upper surface, compared with the power-law case, and basal stresses vanish at the grounding line, consistent with observations. In the Coulomb case, the grounding-line ice flux also depends more strongly on flotation ice thickness, which implies that ice sheets are more sensitive to climate perturbations. Furthermore, with Coulomb friction, the ice sheet grounds stably in shallower water than with a power-law rheology. This implies that smaller perturbations are required to push the grounding line into regions of negative bed slope, where it would become unstable. These results have important implications for ice-sheet stability in a warming climate.

A way forward to discover Antarctica's past

Fifty-four participants attended the Antarctic Geologic Drilling Workshop (AGDW) to discuss science objectives and develop key projects. The goal of the NSF-sponsored AGDW was specifically to discuss the interests of US-based scientists in Antarctic and Southern Ocean projects, foster interactions within the Antarctic geoscience community, and discuss top-priority scientific questions and technological requirements to advance outstanding scientific goals. Antarctica, with its thick shroud of glacial ice and fringed belt of floating oceanic ice, offers challenges to unraveling the history of one of Earth’s last frontiers. Rocks deposited in and around Antarctica, reachable through a wide-range of innovative geological drilling approaches, formed during times that witnessed climate and ice-sheet changes. Recovery of rocks through drilling expands the understanding of the interplay of Earth’s dynamic processes that control and respond to the Antarctic cryosphere. Much of our knowledge of past climate changes, and inferred ice-sheet history, has been obtained from drill cores taken in low-latitude settings. Such far-field proxies offer an outline of ice-sheet behavior, but cannot show which part of the ice sheet changed or what the ocean currents, temperatures, or other controlling parameters were in ice-proximal settings. Drilling in Antarctica can yield samples of rock that were influenced directly by glacial processes and which provide access to Antarctica’s ice-covered geology. Limited outcrops, short gravity cores, and drill cores with partial recovery have been studied from many locations, but such records cannot give the continuous temporal record needed to determine the timing and rates of ice-sheet change or boundary conditions controlling that behavior. Spatially distributed records, including transects from onshore to distal records from each major ice drainage basin, are needed to resolve the individual histories of each area. Seismic data linking drill sites can increase the impact of individual sites by extending the details over a broader area and tightening the time constraints at each site. New over-ice seismic data acquisition through systems like Vibroseis will identify new subglacial drilling targets. New methods to sample bedrock and measure conditions beneath the ice will help refine basal-bed boundary conditions that are vital to reconstructions of icesheet behavior. Numerical modeling can test data-driven hypotheses and evaluate forcing mechanisms under different atmospheric boundary conditions. Within the workshop discussions, two general themes rose to highest priority. One is on late Quaternary interglacials, when Earth and ocean conditions were similar to today and ice retreated landward of its current position. The other priority is the study of mid-Cenozoic ice-sheet history during times when boundary conditions of high atmospheric CO2 approached those estimated for the next century, but when Earth and ocean conditions were different from today. Combining records of ice-sheet behavior during these two time intervals and these two sets of conditions that are converging in our future will yield information needed to advance the development of computer simulations for studies of near-term future behavior of Antarctica’s ice sheets. The full report is available at http://agdw.uh.edu/.

Hydrostatic grounding line parameterization in ice sheet models

Abstract. Modeling of grounding line migration is essential to accurately simulate the behavior of marine ice sheets and investigate their stability. Here, we assess the sensitivity of numerical models to the parameterization of the grounding line position. We run the MISMIP3D benchmark experiments using the Ice Sheet System Model (ISSM) and a two-dimensional shelfy-stream approximation (SSA) model with different mesh resolutions and different sub-element parameterizations of grounding line position. Results show that different grounding line parameterizations lead to different steady state grounding line positions as well as different retreat/advance rates. Our simulations explain why some vertically depth-averaged model simulations deviate significantly from the vast majority of simulations based on SSA in the MISMIP3D benchmark. The results reveal that differences between simulations performed with and without sub-element parameterization are as large as those performed with different approximations of the stress balance equations in this configuration. They also demonstrate that the reversibility test is passed at relatively coarse resolution while much finer resolutions are needed to accurately capture the steady-state grounding line position. We conclude that fixed grid SSA models that do not employ such a parameterization should be avoided, as they do not provide accurate estimates of grounding line dynamics, even at high spatial resolution. For models that include sub-element grounding line parameterization, in the MISMIP3D configuration, a mesh resolution finer than 2 km should be employed.

A temperate former West Antarctic ice sheet suggested by an extensive zone of subglacial meltwater channels

Several recent studies predict that the West Antarctic Ice Sheet will become increasingly unstable under warmer conditions. Insights on such change can be assisted through investigations of the subglacial landscape, which contains imprints of former ice-sheet behavior. Here, we present radio-echo sounding data and satellite imagery revealing a series of ancient large sub-parallel subglacial bed channels preserved in the region between the Möller and Foundation Ice Streams, West Antarctica. We suggest that these newly recognized channels were formed by significant meltwater routed along the ice-sheet bed. The volume of water required is likely substantial and can most easily be explained by water generated at the ice surface. The Greenland Ice Sheet today exemplifies how significant seasonal surface melt can be transferred to the bed via englacial routing. For West Antarctica, the Pliocene (2.6–5.3 Ma) represents the most recent sustained period when temperatures could have been high enough to generate surface melt comparable to that of present-day Greenland. We propose, therefore, that a temperate ice sheet covered this location during Pliocene warm periods.

Probabilistic calibration of a Greenland Ice Sheet model using spatially resolved synthetic observations: toward projections of ice mass loss with uncertainties

Abstract. Computer models of ice sheet behavior are important tools for projecting future sea level rise. The simulated modern ice sheets generated by these models differ markedly as input parameters are varied. To ensure accurate ice sheet mass loss projections, these parameters must be constrained using observational data. Which model parameter combinations make sense, given observations? Our method assigns probabilities to parameter combinations based on how well the model reproduces the Greenland Ice Sheet profile. We improve on the previous state of the art by accounting for spatial information and by carefully sampling the full range of realistic parameter combinations, using statistically rigorous methods. Specifically, we estimate the joint posterior probability density function of model parameters using Gaussian process-based emulation and calibration. This method is an important step toward calibrated probabilistic projections of ice sheet contributions to sea level rise, in that it uses data–model fusion to learn about parameter values. This information can, in turn, be used to make projections while taking into account various sources of uncertainty, including parametric uncertainty, data–model discrepancy, and spatial correlation in the error structure. We demonstrate the utility of our method using a perfect model experiment, which shows that many different parameter combinations can generate similar modern ice sheet profiles. This result suggests that the large divergence of projections from different ice sheet models is partly due to parametric uncertainty. Moreover, our method enables insight into ice sheet processes represented by parameter interactions in the model.

Enthalpy balance methods versus temperature models in ice sheets

In this paper we propose and numerically solve an original enthalpy formulation for the problem governing the thermal behaviour of polythermal ice sheets. Although the modelling follows some ideas introduced in Aschwanden and Blatter (2009), nonlinear basal boundary conditions in both cold and temperate regions are also considered, thus including the sliding effects in the frame of a fully coupled shallow ice approximation (SIA) model. One of the main novelties of this work comes from the introduction of the Heaviside multivalued operator to take into account the discontinuity of the thermal diffusion function at the cold-temperate transition surface (CTS) free boundary. Moreover, we propose a duality method for maximal monotone operators to solve simultaneously the nonlinear diffusive term and the free boundary. Some numerical simulation examples with real data from Antarctica are presented and illustrate the small differences between the computed results from the enthalpy formulation here proposed and the alternative formulation in terms of the temperature (Calvo et al., 2001).

Evidence for the asynchronous retreat of large outlet glaciers in southeast Greenland at the end of the last glaciation

Recent rapid changes in the marine-terminating sectors of the Greenland Ice Sheet (GrIS) have prompted concerns about the future stability of the ice sheet. Long-term records of ice sheet behaviour provide valuable context to assess the magnitude of current change and may help resolve the mechanisms driving deglaciation. We report 23 new 10Be exposure ages which constrain the deglacial history of two large fjord systems in southeast (SE) Greenland. We compare our chronologies with existing data from the centre of the sector to examine the timing and style of deglaciation at a regional-scale. Glacial erratic 10Be exposure ages demonstrate that Kangerdlugssuaq Fjord deglaciated at ∼11.8 ka at the end of the Younger Dryas (YD – 12.8–11.6 ka). Retreat at Kangerdlugssuaq Fjord coincided with known incursion of the warm Irminger Current (IC) onto the continental shelf; this is inferred to have initiated retreat. Comparison with recently published results from Sermilik Fjord and new 10Be ages from Bernstorffs Fjord indicates deglaciation occurred ∼1 ka later in the south of the SE region. Sermilik Fjord (∼10.9 ka) and Bernstorffs Fjord (∼10.4 ka) deglaciated later; retreat likely occurred in response to dramatic climatic amelioration at the termination of the YD stadial. We suggest that the disparate timing of deglaciation across the SE region may be primarily explained by the varying influence of the warm IC; glaciers in southern SE Greenland were isolated from warm Atlantic waters during the YD by complex shelf bathymetry. In all fjord settings ice retreat was rapid and persistent, consistent with the absence of geomorphological evidence for stillstand or readvance events. Ice retreat was accompanied by rapid thinning and likely continued to well within present-day ice sheet margins. Glacial erratic 10Be age determinations and geomorphological observations show no evidence for Holocene readvance events prior to the Little Ice Age (LIA).

Neoproterozoic ice sheets and olistoliths: multiple glacial cycles in the Kingston Peak Formation, California

The Kingston Peak Formation is a diamictite-bearing succession that crops out in the Death Valley region, California, USA. An exceptionally thick (>1.5 km) outcrop belt in its type area (the Kingston Range) provides clear insights into the dynamics of mid-Cryogenian (‘Sturtian’) ice sheets in Laurentia. Seven detailed logs allow the lateral and vertical distribution of facies associations to be assessed. We recognize (1) diamictite facies association (ice-proximal glacigenic debris flows), (2) lonestone-bearing facies association (ice-marginal hemipelagic deposits and low-density gravity flows with iceberg rafting), (3) pebble to boulder conglomerate facies association (ice-proximal cogenetic glacigenic debris flows and high-density turbidites), (4) megaclast facies association (olistostrome and hemipelagic sediments subject to ice-rafting), and (5) interbedded heterolithics facies association (low-density turbidites and hemipelagic deposits). The stratigraphic motif allows three glacial cycles to be inferred across the range. Ice-minimum conditions interrupting the Kingston Peak succession are associated with the development of an olistostrome complex, succeeded by a thick accumulation of boulder conglomerates deposited during ice readvance. The data testify to a strong glacial influence on sedimentation within this ancient subaqueous succession, and to highly dynamic ice sheet behaviour with clear glacial cycles during the Sturtian glaciation.

The subglacial geology of Wilkes Land, East Antarctica

Wilkes Land is a key region for studying the configuration of Gondwana and for appreciating the role of tectonic boundary conditions on East Antarctic Ice Sheet (EAIS) behavior. Despite this importance, it remains one of the largest regions on Earth where we lack a basic knowledge of geology. New magnetic, gravity, and subglacial topography data allow the region's first comprehensive geological interpretation. We map lithospheric domains and their bounding faults, including the suture between Indo-Antarctica and Australo-Antarctica. Furthermore, we image subglacial sedimentary basins, including the Aurora and Knox Subglacial Basins and the previously unknown Sabrina Subglacial Basin. Commonality of structure in magnetic, gravity, and topography data suggest that pre-EAIS tectonic features are a primary control on subglacial topography. The preservation of this relationship after glaciation suggests that these tectonic features provide topographic and basal boundary conditions that have strongly influenced the structure and evolution of the EAIS.

Testing techniques to quantify drumlin height and volume: synthetic DEMs as a diagnostic tool

ABSTRACTGlacial bedform height (H) and volume (V) likely preserve important information about the behaviour of former ice sheets. However, large systematic errors exist in the measurement of H and V. Three semi‐automated methods to isolate drumlins from other components of the landscape (e.g. trees, hills) as portrayed by NEXTMap have recently been devised; however, it is unclear which is most accurate. This paper undertakes the first quantitative comparison of such readily implementable methods, illustrating the use of statistically representative ‘synthetic landscapes’ as a diagnostic tool. From this analysis, guidelines for quantifying the 3D attributes of drumlins are proposed. Specifically, to avoid obtaining incorrect estimates caused by substantial systematic biases, interpreters should currently take three steps: declutter the digital elevation model for estimating H but not V; remove height data within the drumlin; then interpolate across the resultant hole to estimate a basal surface using Delaunay triangulation. Results are demonstrated through analysis of drumlins in an area in western central Scotland. The guidance arguably represents the best current advice for subglacial bedforms in general, highlighting the need for more studies into the quality of mapped data using synthetic landscapes. Copyright © 2014 John Wiley & Sons, Ltd.

Sequencing the Sturtian icehouse: dynamic ice behaviour in South Australia

The Cryogenian record of South Australia includes the type region of the Sturtian glaciation, the oldest of three pan-global icehouse intervals during the Neoproterozoic. Data are presented from previously little described sections at Holowilena Creek, Oladdie Creek and Hillpara Creek in the central and southern Flinders Ranges, where five facies associations are recognized. These are (1) diamictite and conglomerate, (2) interbedded heterolithic deposits, (3) hummocky cross-stratified sandstone, (4) lonestone-bearing siltstone, and (5) ferruginous siltstone and sandstone. The succession reveals significant lateral and vertical facies variation, which is linked to a complex inherited palaeotopography and distance from the sediment source. Repeated stratigraphic occurrences of striated clasts and abundant ice-rafted debris strongly support recurrent glacial influence on sedimentation. The intercalation of gravitationally reworked diamictites, dropstone-bearing siltstone and dropstone-free siltstone testifies to dynamic sedimentation within a periodically glacially influenced subaqueous environment. Sequence stratigraphic analysis identifies four glacial advance systems tracts, separated by three glacial retreat systems tracts, wherein hummocky cross-stratified sandstones attest to open water conditions. These findings support dynamic ice sheet behaviour in South Australia, and provide clear evidence for repeated intra-Sturtian ice sheet recession.

Greenland ice sheet contribution to future global sea level rise based on CMIP5 models

Sea level rise (SLR) is one of the major socioeconomic risks associated with global warming. Mass losses from the Greenland ice sheet (GrIS) will be partially responsible for future SLR, although there are large uncertainties in modeled climate and ice sheet behavior. We used the ice sheet model SICOPOLIS (SImulation COde for POLythermal Ice Sheets) driven by climate projections from 20 models in the fifth phase of the Coupled Model Intercomparison Project (CMIP5) to estimate the GrIS contribution to global SLR. Based on the outputs of the 20 models, it is estimated that the GrIS will contribute 0–16 (0–27) cm to global SLR by 2100 under the Representative Concentration Pathways (RCP) 4.5 (RCP 8.5) scenarios. The projected SLR increases further to 7–22 (7–33) cm with 2×basal sliding included. In response to the results of the multimodel ensemble mean, the ice sheet model projects a global SLR of 3 cm and 7 cm (10 cm and 13 cm with 2×basal sliding) under the RCP 4.5 and RCP 8.5 scenarios, respectively. In addition, our results suggest that the uncertainty in future sea level projection caused by the large spread in climate projections could be reduced with model-evaluation and the selective use of model outputs.

Landscape imprints of changing glacial regimes during ice-sheet build-up and decay: a conceptual model from Svalbard

The behaviour of ice sheets and their geologic imprints in fjord regions are often multifaceted. Fjords, which were temporarily occupied by fast flowing ice-streams during major glaciations, and inter-fjord areas, which were covered by less active ice, show different signatures of past glaciations. The land and marine records of glaciations over the western Svalbard fjord region have been extensively studied during the last few decades. We have re-examined ice-flow records from stratigraphic and geomorphic settings, and propose a succession of ice-flow events that occurred repeatedly over glacial cycles: the maximum, the transitional, and the local flow style. The differently topographically constrained segments of the ice-sheet switched behaviour as glacial dynamics developed through each glacial cycle. These segments, as well as the different flow styles, are reflected differently in the offshore stratigraphic record. We propose that the glacial geomorphological signatures in the inter ice-stream areas mostly developed under warm-based conditions during a late phase of the glaciations, and that the overall glacial imprints in the landscape are strongly biased towards the youngest events.

Rapid marine deglaciation: asynchronous retreat dynamics between the Irish Sea Ice Stream and terrestrial outlet glaciers

Abstract. Understanding the retreat behaviour of past marine-based ice sheets provides vital context for accurate assessments of the present stability and long-term response of contemporary polar ice sheets to climate and oceanic warming. Here new multibeam swath bathymetry data and sedimentological analysis are combined with high resolution ice-sheet modelling to reveal complex landform assemblages and process dynamics associated with deglaciation of the Celtic ice sheet within the Irish Sea Basin. Our reconstruction indicates a non-linear relationship between the rapidly receding Irish Sea Ice Stream and the retreat of outlet glaciers draining the adjacent, terrestrially based ice cap centred over Wales. Retreat of Welsh ice was episodic; superimposed over low-order oscillations of its margin are asynchronous outlet readvances driven by catchment-wide mass balance variations that are amplified through migration of the ice cap's main ice divide. Formation of large, linear ridges which extend at least 12.5 km offshore (locally known as sarns) and which dominate the regional bathymetry are attributed to repeated frontal and medial morainic deposition associated with the readvancing phases of these outlet glaciers. Our study provides new insight into ice-sheet extent, dynamics and non-linear retreat across a major palaeo-ice stream confluence zone, and has ramifications for the interpretation of recent fluctuations observed by satellites over short timescales across marine sectors of the Greenland and Antarctic ice sheets.

A 2 million year glacial chronology of the Hatherton Glacier, Antarctica and implications for the size of the East Antarctic Ice Sheet at the Last Glacial Maximum

A series of distinct glacial deposits flanking the margins of the upper Hatherton Glacier, an outlet glacier in the central Transantarctic Mountains, are used to constrain the behaviour of the Antarctic ice-sheets. Cosmogenic exposure ages of 18 erratics from four glacial drifts covering the ice free Dubris and Bibra valleys, range in age from 5 to 1997 ka. Our results document four glacial advance and retreat events superimposed on an overall long-term ice thickness reduction of about 500 m since the mid-Pleistocene. The lack of field evidence and absence of LGM exposure ages in the glacial deposits of the Hatherton Glacier supports our conclusion that at the LGM the East Antarctic Ice Sheet was of similar size, or may have been slightly smaller, than present. Minimum exposure ages from the oldest two glacial events, represented by the Isca and Danum drifts, are ∼1–2 Ma and ∼0.5 Ma respectively. The Britannia-II Drift, previously assumed to mark the maximum extent of the Last Glacial Maximum advance, has a mean 10Be age of 126 ± 3.2 ka (n = 5). Ages from the younger Britannia-I Drift suggest that since the mid-Holocene (6.5 ± 1.2 ka, n = 5), approximately 200 m of additional ice has been lost.

Evidence for multiple Quaternary ice advances and fan development from the Amundsen Gulf cross-shelf trough and slope, Canadian Beaufort Sea margin

The seismic stratigraphy and sedimentary architecture of the Amundsen Gulf Trough and adjacent slope, Canadian Beaufort Sea margin, are investigated using a grid of 2-D seismic reflection data. The inner-shelf of the Amundsen Gulf Trough is interpreted to be composed predominantly of exposed or near-surface bedrock, overlain by a spatially-discontinuous veneer of glacimarine to open-marine sediment. There is a seaward transition from exposed bedrock on the inner-shelf to a thick (up to 500 m) outer-shelf prograding wedge of acoustically semi-transparent sediment. Eight seismic sequences, divided into four megasequences, are described from the outer-shelf stratigraphy. Eight till sheets are identified from Megasequences A to C, providing evidence for at least eight Quaternary ice-stream advances through the Amundsen Gulf Trough to the shelf break. A trough-mouth fan with a minimum volume of about 10,000 km3 is present on the adjacent slope. The Amundsen Gulf ice stream probably represented the most northwesterly marine-terminating ice stream of the Laurentide Ice Sheet through much of the Quaternary, providing a major route for ice and sediment transfer to the Arctic Ocean. The youngest till sheet within the Amundsen Gulf Trough, Megasequence D, was probably deposited by a subsidiary ice stream, the Anderson ice stream, subsequent to retreat of the last, Late Wisconsinan Amundsen Gulf ice stream. This provides evidence of dynamic ice-sheet behaviour and the reorganisation of the northwest Laurentide Ice Sheet margin during the last deglaciation. A number of buried glacigenic landforms, including palaeo-shelf break gullies and a grounding-zone wedge with a volume of 90 km3, are described from the Amundsen Gulf Trough stratigraphy. Lateral grounding-zone wedges are identified at the northern and southern lateral margins of the Amundsen Gulf and M'Clure Strait troughs, respectively, and are interpreted to have been formed roughly contemporaneously by ice streams in Amundsen Gulf and M'Clure Strait.

Hot rocks in a cold place: high sub-glacial heat flow in East Antarctica

Numerical models are the primary predictive tools for understanding the dynamic behaviour of the Antarctic ice sheet. However, a key boundary parameter, sub-glacial heat flow, remains poorly constrained. We show that variations in abundance and distribution of heat-producing elements within the Antarctic continental crust result in greater and more variable regional sub-glacial heat flows than currently assumed in ice modelling studies. Such elevated heat flows would have a fundamental effect on ice sheet behaviour and highlight that geological controls on heat flow must be considered to obtain more accurate and refined predictions of ice mass balance and sea-level change. Supplementary materials: Heat flow and heat production definitions, heat production data, and details of the 2D model of the gross geometry of the continental lithosphere for a section through Prydz Bay are available at www.geolsoc.org.uk/SUP18690 .