Search

The last Fennoscandian Ice Sheet provides a valuable scenario for testing and evaluating numerical ice sheet models with a large amassed database of landform, stratigraphic and dated evidence of ice sheet activity. In contrast to the core shield area (Norway, Sweden and Finland) of the ice sheet, fewer investigations beyond the shield (Denmark, Germany, Poland, Lithuania, Latvia, Estonia and Russia) attempt to gather local to regional information into ice sheet wide syntheses of ice margin and lobe dynamics. For example, many detailed investigations across these countries remain disconnected with adjacent areas applying varying methods and naming schemes making it difficult to reconcile at the ice sheet scale.Here we present a systematic and spatially coherent reconstruction of ice margin dynamics for the whole southern and eastern margin, from Denmark to arctic Russia. The landform to reconstruction method allows for a consistent approach to be applied to the 1.2+ million km2 mapping area despite a wide range of glaciological landform and data variability (DEM vary in resolution from 0.4 m-25 m) found in the 1.2+ million km2 study area. We propose this reconstruction as a first-order framework of ice marginal dynamics that can be used to develop second-order and more detailed knowledge of fluctuations when more closely connected to stratigraphic and geochronometric investigations. Rather than a simple concentric retreat pattern often envisaged the landform record and its frequent overprinting forces a solution of complexity with lobe interactions and readvances.

The Laurentide Ice Sheet (LIS) was the largest of the ephemeral Pleistocene ice sheets in the Northern Hemisphere, with a Last Glacial Maximum (LGM) ice volume similar to the modern Antarctic ice sheet. A recent inventory of paleo-ice streams across the LIS shows many similarities with present-day ice streaming in Antarctica, where ice streams account for approx. 90% of mass loss. However, in the Mackenzie Lowlands of the Northwest Territories, Canada, the paleo-ice stream record is enigmatic. Previous work has identified a number of large paleo-ice streams, including the Mackenzie Trough, Anderson, Bear Lake and Fort Simpson ice streams, however, their extent, configuration, temporal relationship to each other and spatial evolution over time remains poorly constrained. Consequently, their impact on the rate and style of deglaciation of the northwestern sector of the LIS is poorly understood.Here we utilise the newly available high resolution Arctic DEM (0.5 m resolution) to re-map glacial landforms across the Mackenzie Lowlands in greater detail (area >800,000 km2). We then use this landform record to reconstruct the ice dynamics in this region following the well-established approaches of flowset mapping and the glacial inversion method. The high resolution data allow us to present a detailed reconstruction of LGM ice flow over the Mackenzie Lowlands and resolve the configuration and evolution of ice streams over time. The landform record suggests that the ice streams operated time-transgressively during deglaciation, switching on and off at different times. While ice contact landforms, such as moraines, lateral and submarginal meltwater channels and ice-contact deltas, show the overall retreat of the LIS towards the Keewatin Dome in the east, in several regions the ice retreat record is complex, suggesting interlobate ice configurations with multiple ice retreat directions.

Glaciers and ice sheets fluctuate in response to climatic change and often record these changes by building ice-marginal (terminal and lateral) moraines. Therefore, glacial landscapes are a potentially valuable archive of terrestrial palaeoclimate change. Typically, a cooling climate causes glaciers to expand and warming causes glaciers to shrink. However, the influence of high-relief mountainous topography on glacier dynamics complicates this behaviour, such that ice-marginal moraines are not always a straightforward palaeoclimate indicator. We used a higher-order ice-flow model to simulate change in glacier erosion, extent, and thickness in the response to climatic change and the resulting formation and preservation of moraines in a synthetic mountain landscape. Our results show that the rate of palaeoclimatic change relative to the glacier’s response time determines the geometry, number and position of ice-marginal moraines. However, glaciers can build distinct moraines in the absence of climate change, and the distance from the glacial maximum may not represent the chronological order of moraine formation. While moraines can be preserved despite erosion during subsequent glaciations, moraine sequences frequently contain gaps that could be misinterpreted as representing more stable palaeoclimates. These results provide theoretical understanding for the interpretation of glacial landforms both in the field and from satellite data (e.g. digital terrain models) to understand Quaternary climate change.

Over the three last decades, great efforts have been undertaken by the glaciological community to characterize the behaviour of ice streams and better constrain the dynamics of ice sheets. Studies of modern ice stream beds reveal crucial information on ice-meltwater-till-bedrock interactions, but are restricted to punctual observations limiting the understanding of ice stream dynamics as a whole. Consequently, theoretical ice stream landsystems derived from geomorphological and sedimentological observations were developed to provide wider constraints on those interactions on palaeo-ice stream beds. Within these landsystems, the spatial distribution and formation processes of subglacial periodic bedforms transverse to the ice flow direction – ribbed bedforms – remain unclear. The purpose of this study is (i) to explore the conditions under which these ribbed bedforms develop and (ii) to constrain their spatial organisation along ice stream beds. We performed physical experiments with silicon putty (to simulate the ice), water (to simulate the meltwater) and sand (to simulate a soft sedimentary bed) to model the dynamics of ice streams and produce analog subglacial landsystems. We compare the results of these experiments with the distribution of ribbed bedforms on selected examples of palaeo-ice stream beds of the Laurentide Ice Sheet. Based on this comparison, we can draw several conclusions regarding the significance of ribbed bedforms in ice stream contexts:Ribbed bedforms tend to form where the ice flow undergoes high velocity gradients and the ice-bed interface is unlubricated. Where the ribs initiate, we hypothesize that high driving stresses generate high basal shear stresses, accommodated through bed deformation of the active uppermost part of the bed. Ribbed bedforms can develop subglacially from a flat sediment surface beneath shear margins (i.e., lateral ribbed bedforms) and stagnant lobes (i.e., submarginal ribbed bedforms) of ice streams, while they do not develop beneath surging lobes. The orientation of ribbed bedforms reflects the local stress state along the ice-bed interface, with transverse bedforms formed by compression beneath ice lobes and oblique bedforms formed by transgression below lateral shear margins. The development of ribbed bedforms where the ice-bed interface is unlubricated reveals distinctive types of discontinuous basal drainage systems below shear and lobe margins: linked-cavities and efficient meltwater channels respectively. Ribbed bedforms could thus constitute convenient geomorphic markers for the reconstruction of palaeo-ice stream margins, palaeo-ice flow dynamics and palaeo-meltwater drainage characteristics.

Meltwater exerts an important influence on ice sheet dynamics and has attracted an increasing amount of attention over the last 20 years. However, the active subglacial environment remains difficult to study mainly due to its inaccessibility. Understanding of the dimensions, pattern, and extent of subglacial meltwater conduits at the ice sheet scale is limited to relatively sparse observations. We address this gap by using the geomorphological record of Quaternary ice sheets as a proxy to quantify the dimensions and pattern of subglacial conduits at the ice sheet scale. We present the results of a high-resolution (2 m), large sample (n>50,000) study of three-dimensional esker morphometry at sample locations in SW Finland and Nunavut, Canada. Detailed mapping of esker crestlines and outlines permits the quantification of a number of parameters, including: length, width, height, cross-sectional area, volume, sinuosity, cross-sectional symmetry, and uphill/downhill trends. Whilst the dimensions of eskers reflect depositional processes as well as simply the size of the parent conduit, they nevertheless offer a powerful tool for understanding the size and shape of meltwater conduits and the configuration of subglacial drainage systems across large areas (entire ice sheets), and over long periods of time (from years to thousands of years) in both high spatial and temporal resolution. The results may be used to: (1) inform numerical models of subglacial meltwater drainage, (2) inform process models of esker formation, and (3) provide a dataset of esker morphometry against which other features may be compared (e.g. sinuous ridges on Mars).

The Western Gneiss Region (WGR), Norway is an archetypal continental ultrahigh-pressure (U)HP terrane with an extensive metamorphic history, recording the subduction and subsequent exhumation of continental crust to depths exceeding 120 km. The vast bulk of past work within the WGR has focused on mafic eclogites. In this study, data from rare garnet-kyanite metapelites in (UHP) domains of the WGR is presented. U–Pb geochronology and trace element compositions in zircon, monazite, apatite, rutile and garnet were acquired, and P–T conditions were calculated by mineral equilibria forward modelling and Zr-in-rutile thermometry. The Ulsteinvik metapelite defines a prograde path that traverses through ~600–710 °C and ~11–14 kbar. Minimum peak conditions are ~750 °C and ~2.9 GPa in an inferred garnet-kyanite-coesite-omphacite-muscovite-rutile-quartz-H2O assemblage. Plagioclase-biotite-quartz intergrowths developed after omphacite-phengite-rutile breakdown on the early retrograde path, followed by cordierite-spinel-plagioclase symplectites after garnet-kyanite-biotite, defining a retrograde P–T point at ~740 °C and ~7 kbar. Late Ordovician-Early Silurian (~470–440 Ma) zircon and rutile age data in Ulsteinvik pre-dates the major Scandian UHP subduction episode in the WGR, interpreted as recording early Caledonian subduction within the Blåhø nappe. Monazite and apatite U-Pb geochronology and trace element data suggest exhumation occurred at ~400 Ma. The Fjørtoft metapelite is a constituent of the Blåhø nappe. Minimum peak P–T conditions are ~1.8 GPa and ~750 °C, with poor peak mineral fidelity attributed to extensive retrograde deformation. Negative Eu anomalies in ~423 Ma monazite suggest retrograde conditions were reached [RJT1] by ~423 Ma. Ulsteinvik and Fjørtoft may have experienced pre-Scandian subduction together within the Blåhø nappe, but record dissimilar histories after this. Two potential scenarios are presented: (1) Ulsteinvik resided within the mantle for 20 million-years longer than Fjørtoft during Scandian subduction, or (2), the samples were exhumed at different times during pre-Scandian subduction of the Blåhø nappe. The preservation of prograde zoning within Ulsteinvik garnets precludes a long-term residence within the mantle and suggests the latter option. In this scenario, the subducting Blåhø nappe experienced a degree of slab tear and partial underplating of the upper plate during the early stages of continental underthrusting. Discrete pieces may have later reattached to the lower plate at different times, partially exhumed, and then subducted to mantle-depths during the Scandian.

The Polish sector of the last Scandinavian Ice Sheet is a key area for studying ice sheet drainage and decay from its local Last Glacial Maximum (LGM) extent, as it is located at the terrestrial terminus of the large and dynamic Baltic Ice Stream Complex. Geomorphological mapping, based on a 0.4 m LIDAR digital elevation model, revealed about 940 streamlined bedforms, many of which are shown for the first time and consisting of mega-scale glacial lineations and drumlins. The lineation flow-sets together with associated landforms were used to identify seventeen ice streams, occupying 80% of the study area. We demonstrated that subtle topographic variations played an important role in influencing ice sheet dynamics. Variations in ice dynamics were a response to external climatic forcing that controlled deglaciation at the ice sheet scale as well as internal reorganisation due to the influence of topography, subglacial hydrology and glacier thermal regime. During the local LGM, the southern sector of the Scandinavian Ice Sheet in Poland was dominated by four simultaneously operating ice streams, likely active for several millennia, followed by fast active recession interrupted by three main periods of ice stream stagnation. Increased ice flowdynamics during the period of the Young Baltic advances is suggested to be caused by variations in subglacial hydrology and the polythermal structure of the ice sheet.

We identify and map chains of esker beads (series of aligned mounds) up to 15 m high and on average ∼ 65 m wide in central Nunavut, Canada, from the high-resolution (2 m) ArcticDEM. Based on the close 1 : 1 association with regularly spaced, sharp-crested ridges interpreted as De Geer moraines, we interpret the esker beads to be quasi-annual ice-marginal deposits formed time-transgressively at the mouth of subglacial conduits during deglaciation. Esker beads therefore preserve a high-resolution record of ice-margin retreat and subglacial hydrology. The well-organised beaded esker network implies that subglacial channelised drainage was relatively fixed in space and through time. Downstream esker bead spacing constrains the typical pace of deglaciation in central Nunavut between 8.1 and 6.8 cal kyr BP to 165–370 m yr−1, although with short periods of more rapid retreat (> 400 m yr−1). Under our time-transgressive interpretation, the lateral spacing of the observed eskers provides a true measure of subglacial conduit spacing for testing mathematical models of subglacial hydrology. Esker beads also record the volume of sediment deposited from conduits in each melt season, thus providing a minimum bound on annual sediment fluxes, which is in the range of 103–104 m3 yr−1 in each 6–10 km wide subglacial conduit catchment. We suggest that the prevalence of esker beads across this predominantly marine-terminating sector of the Laurentide Ice Sheet is a result of sediment fluxes that were unable to backfill conduits at a rate faster than ice-margin retreat. Conversely, we hypothesise that esker ridges form when sediment backfilling of the subglacial conduit outpaced retreat, resulting in headward esker growth close to but behind the margin. The implication, in accordance with recent modelling results, is that eskers in general record a composite signature of ice-marginal drainage rather than a temporal snapshot of ice-sheet-wide subglacial drainage.

Geochronological data in zircon from Archaean tonalite–trondhjemite–tonalite (TTG) gneisses is commonly difficult to interpret. A notable example are TTG gneisses from the Lewisian Gneiss Complex (LGC), northwest Scotland, which have metamorphic zircon ages that define a more-or-less continuous spread through the Neoarchaean, with no clear relationship to zircon textures. These data are generally interpreted to record discrete high-grade events at c. 2.7 Ga and c. 2.5 Ga, with intermediate ages reflecting variable Pb-loss. Although ancient diffusion of Pb is commonly invoked to explain such protracted age spreads, trace element data in zircon may permit identification of otherwise cryptic magmatic and metamorphic episodes. Although zircons from the TTG gneiss analyzed here show a characteristic spread of Neoarchaean ages, they exhibit subtle but key step changes in trace element compositions that are difficult to ascribe to diffusive resetting, but which are consistent with emplacement of regionally-extensive bodies of mafic magma. These data suggest suprasolidus metamorphic temperatures persisted for 200 Myr or more during the Neoarchaean. Such long-lived high-grade metamorphism is supported by data from zircon grains from a nearby monzogranite sheet. These preserve distinctive trace element compositions suggesting derivation from a mafic source, and define a well-constrained U–Pb zircon age of c. 2.6 Ga that is intermediate between the two previously proposed discrete metamorphic episodes. The persistence for hundreds of millions of years of melt-bearing lower crust was probably the norm during the Archaean.

The configuration of subglacial drainage systems has a major impact on the dynamics of ice sheets. However, the logistical challenges of measuring subglacial processes beneath contemporary ice sheets hinder our understanding about the spatio-temporal evolution of subglacial drainage systems. Furthermore, today’s observations on contemporary ice sheets are inherently limited to a short period within the process of deglaciation. Landforms generated by the flow of meltwater at the ice-bed interface offer the potential to study both large-scale (103-106 km2) and long-term (103-105 a) developments of subglacial drainage networks beneath past ice sheets. Despite collectively recording subglacial drainage, individual meltwater landform types such as eskers, meltwater channels and tunnel valleys, and hummock corridors have mostly been considered as separate entities. Using high-resolution (1-2 m) DEMs, we summarise the suite of interconnected subglacial meltwater landforms into a common drainage signature herein called a subglacial drainage route. Our integrated map of subglacial meltwater landforms presents the large-scale distribution of major subglacial drainage routes across Scandinavia and provides a basis for future research about the long-term evolution of subglacial drainage networks and its effect on ice dynamics of the Scandinavian Ice Sheet.

Despite being widely studied, subglacial meltwater landforms are typically mapped and investigated individually, thus the drainage system as a whole remains poorly understood. Here, we identify and map all visible traces of subglacial meltwater flow across the Keewatin sector of the former Laurentide Ice Sheet from the ArcticDEM, generating significant new insights into the connectedness of the drainage system.Due to similarities in spacing, morphometry and spatial location, we suggest that the 100s-1000s m wide features often flanking and connecting sections of eskers (i.e. tunnel valleys, meltwater tracks and esker splays) are varying expressions of the same phenomena and collectively term these features ‘meltwater corridors’. Based on observations from contemporary ice masses, we propose a new formation model based on the pressure fluctuations surrounding a central conduit, in which the esker records the imprint of the central conduit and the wider meltwater corridors the interactions with the surrounding distributed drainage system, or variable pressure axis (VPA).We suggest that the widespread aerial coverage of meltwater corridors across the Keewatin sector provides constraints on the extent of basal uncoupling induced by basal water pressure fluctuation and variations in spatial distribution and evolution of the subglacial drainage system, which have important implications for ice sheet dynamics.

Uncertainties in future sea level projections are dominated by our limited understanding of the dynamical processes that control instabilities of marine ice sheets. The last deglaciation of the British–Irish Ice Sheet offers a valuable example to examine these processes. The Minch Ice Stream, which drained a large proportion of ice from the northwest sector of the British–Irish Ice Sheet during the last deglaciation, is constrained with abundant empirical data which can be used to inform, validate, and analyse numerical ice sheet simulations. We use BISICLES, a higher-order ice sheet model, to examine the dynamical processes that controlled the retreat of the Minch Ice Stream. We perform simplified experiments of the retreat of this ice stream under an idealised climate forcing to isolate the effect of marine ice sheet processes, simulating retreat from the continental shelf under constant “warm” surface mass balance and sub-ice-shelf melt. The model simulates a slowdown of retreat as the ice stream becomes laterally confined at the mouth of the Minch strait between mainland Scotland and the Isle of Lewis, resulting in a marine setting similar to many large tidewater glaciers in Greenland and Antarctica. At this stage of the simulation, the presence of an ice shelf becomes a more important control on grounded ice volume, providing buttressing to upstream ice. Subsequently, the presence of a reverse slope inside the Minch strait produces an acceleration in retreat, leading to a “collapsed” state, even when the climate returns to the initial “cold” conditions. Our simulations demonstrate the importance of the marine ice sheet instability and ice shelf buttressing during the deglaciation of parts of the British–Irish Ice Sheet. We conclude that geological data could be applied to further constrain these processes in ice sheet models used for projecting the future of contemporary ice sheets.

Uncertainties in future sea level projections are dominated by our limited understanding of the dynamical processes that control instabilities of marine ice sheets. A valuable case to examine these processes is the last deglaciation of the British-Irish Ice Sheet. The Minch Ice Stream, which drained a large proportion of ice from the northwest sector of the British-Irish Ice Sheet during the last deglaciation, is well constrained, with abundant empirical data which could be used to inform, validate and analyse numerical ice sheet simulations. We use BISICLES, a higher-order ice sheet model, to examine the dynamical processes that controlled the retreat of the Minch Ice Stream. We simulate retreat from the shelf edge under constant "warm" surface mass balance and subshelf melt, to isolate the role of internal ice dynamics from external forcings. The model simulates a slowdown of retreat as the ice stream becomes laterally confined at a "pinning-point" between mainland Scotland and the Isle of Lewis. At this stage, the presence of ice shelves became a major control on deglaciation, providing buttressing to upstream ice. Subsequently, the presence of a reverse slope inside the Minch Strait produces an acceleration in retreat, leading to a "collapsed" state, even when the climate returns to the initial "cold" conditions. Our simulations demonstrate the importance of the Marine Ice Sheet Instability and ice shelf buttressing during the deglaciation of parts of the British-Irish Ice Sheet. Thus, geological data could be used to constrain these processes in ice sheet models used for projecting the future of our contemporary ice sheets.