Forward modelling of the completeness and preservation of palaeoclimate signals recorded by ice-marginal moraines.

Glaciers fluctuate in response to climate change and record these changes by building sedimentary landforms, including moraines. Therefore, glacial landscapes are a potentially valuable archive of terrestrial palaeoclimate change. Typically, a cooling climate causes glaciers to expand and a warming climate causes glaciers to shrink. However, the glacier response time and the influence of mountainous topography on glacier dynamics complicates this behaviour, such that moraines are not always a straightforward indicator of glacier change in response to climate change. We used a glacial landscape evolution model to simulate the response of a hypothetical mountain glacier to simple changes in climate and the resulting formation and preservation of moraines. These results show that the rate of climate change relative to the glacier response time determines the geometry, number, and position of moraines. Glaciers can build distinct moraines in the absence of climate change. The distance from the maximum ice extent may not represent the chronological order of moraine formation. Moraines can be preserved after being overrun and eroded by subsequent glaciations, but moraine sequences may also contain gaps that are unidentifiable in the field.

Subglacial meltwater routes of the Fennoscandian Ice Sheet.

Subglacial drainage systems are crucial elements of glaciers and ice sheets because they modulate ice flow velocity. However, logistical challenges of measuring subglacial processes beneath contemporary ice and natural limitations in long-term monitoring hinder our understanding about their spatio-temporal evolution. Subglacial meltwater landforms created by palaeo-ice sheets are records of past subglacial drainage systems and offer the potential to study their large-scale development throughout deglaciation. Although collectively recording subglacial drainage, individual meltwater landforms such as eskers, meltwater channels and meltwater corridors, which comprise tunnel valleys and meltwater tracks (assemblages of landforms in broad, elongated paths with irregular surface texture), have mostly been investigated as separate entities. Using high-resolution (1-2 m) digital elevation models, we map integrated networks of subglacial meltwater landforms, herein called subglacial meltwater routes, on an ice-sheet scale in Fennoscandia. Our map provides a basis for future research on the long-term evolution of subglacial drainage networks and its effect on ice dynamics of the Fennoscandian Ice Sheet.

Ice-stream demise dynamically conditioned by trough shape and bed strength.

Ice sheet mass loss is currently dominated by fast-flowing glaciers (ice streams) terminating in the ocean as ice shelves and resting on beds below sea level. The factors controlling ice-stream flow and retreat over longer time scales (>100 years), especially the role of three-dimensional bed shape and bed strength, remain major uncertainties. We focus on a former ice stream where trough shape and bed substrate are known, or can be defined, to reconstruct ice-stream retreat history and grounding-line movements over 15 millennia since the Last Glacial Maximum. We identify a major behavioral step change around 18,500 to 16,000 years ago-out of tune with external forcing factors-associated with the collapse of floating ice sectors and rapid ice-front retreat. We attribute this step change to a marked geological transition from a soft/weak bed to a hard/strong bed coincident with a change in trough geometry. Both these factors conditioned and ultimately hastened ice-stream demise.

Numerical simulations of a kilometre-thick Arctic ice shelf consistent with ice grounding observations.

Recently obtained geophysical data show sets of parallel erosional features on the Lomonosov Ridge in the central Arctic Basin, indicative of ice grounding in water depths up to 1280 m. These features have been interpreted as being formed by an ice shelf-either restricted to the Amerasian Basin (the "minimum model") or extending across the entire Arctic Basin. Here, we use a numerical ice sheet-shelf model to explore how such an ice shelf could form. We rule out the "minimum model" and suggest that grounding on the Lomonosov Ridge requires complete Arctic ice shelf cover; this places a minimum estimate on its volume, which would have exceeded that of the modern Greenland Ice Sheet. Buttressing provided by an Arctic ice shelf would have increased volumes of the peripheral terrestrial ice sheets. An Arctic ice shelf could have formed even in the absence of a hypothesised East Siberian Ice Sheet.

Exploring Explanations of Subglacial Bedform Sizes Using Statistical Models.

Sediments beneath modern ice sheets exert a key control on their flow, but are largely inaccessible except through geophysics or boreholes. In contrast, palaeo-ice sheet beds are accessible, and typically characterised by numerous bedforms. However, the interaction between bedforms and ice flow is poorly constrained and it is not clear how bedform sizes might reflect ice flow conditions. To better understand this link we present a first exploration of a variety of statistical models to explain the size distribution of some common subglacial bedforms (i.e., drumlins, ribbed moraine, MSGL). By considering a range of models, constructed to reflect key aspects of the physical processes, it is possible to infer that the size distributions are most effectively explained when the dynamics of ice-water-sediment interaction associated with bedform growth is fundamentally random. A 'stochastic instability' (SI) model, which integrates random bedform growth and shrinking through time with exponential growth, is preferred and is consistent with other observations of palaeo-bedforms and geophysical surveys of active ice sheets. Furthermore, we give a proof-of-concept demonstration that our statistical approach can bridge the gap between geomorphological observations and physical models, directly linking measurable size-frequency parameters to properties of ice sheet flow (e.g., ice velocity). Moreover, statistically developing existing models as proposed allows quantitative predictions to be made about sizes, making the models testable; a first illustration of this is given for a hypothesised repeat geophysical survey of bedforms under active ice. Thus, we further demonstrate the potential of size-frequency distributions of subglacial bedforms to assist the elucidation of subglacial processes and better constrain ice sheet models.

Discovery of relict subglacial lakes and their geometry and mechanism of drainage.

Recent proxy measurements reveal that subglacial lakes beneath modern ice sheets periodically store and release large volumes of water, providing an important but poorly understood influence on contemporary ice dynamics and mass balance. This is because direct observations of how lake drainage initiates and proceeds are lacking. Here we present physical evidence of the mechanism and geometry of lake drainage from the discovery of relict subglacial lakes formed during the last glaciation in Canada. These palaeo-subglacial lakes comprised shallow (<10 m) lenses of water perched behind ridges orientated transverse to ice flow. We show that lakes periodically drained through channels incised into bed substrate (canals). Canals sometimes trend into eskers that represent the depositional imprint of the last high-magnitude lake outburst. The subglacial lakes and channels are preserved on top of glacial lineations, indicating long-term re-organization of the subglacial drainage system and coupling to ice flow.

Ice stream motion facilitated by a shallow-deforming and accreting bed.

Ice streams drain large portions of ice sheets and play a fundamental role in governing their response to atmospheric and oceanic forcing, with implications for sea-level change. The mechanisms that generate ice stream flow remain elusive. Basal sliding and/or bed deformation have been hypothesized, but ice stream beds are largely inaccessible. Here we present a comprehensive, multi-scale study of the internal structure of mega-scale glacial lineations (MSGLs) formed at the bed of a palaeo ice stream. Analyses were undertaken at macro- and microscales, using multiple techniques including X-ray tomography, thin sections and ground penetrating radar (GPR) acquisitions. Results reveal homogeneity in stratigraphy, kinematics, granulometry and petrography. The consistency of the physical and geological properties demonstrates a continuously accreting, shallow-deforming, bed and invariant basal conditions. This implies that ice stream basal motion on soft sediment beds during MSGL formation is accommodated by plastic deformation, facilitated by continuous sediment supply and an inefficient drainage system.