The build-up, configuration, and dynamical sensitivity of the Eurasian ice-sheet complex to Late Weichselian climatic and oceanic forcing

The Eurasian ice-sheet complex (EISC) was the third largest ice mass during the Last Glacial Maximum (LGM), after the Antarctic and North American ice sheets. Despite its global significance, a comprehensive account of its evolution from independent nucleation centres to its maximum extent is conspicuously lacking. Here, a first-order, thermomechanical model, robustly constrained by empirical evidence, is used to investigate the dynamics of the EISC throughout its build-up to its maximum configuration. The ice flow model is coupled to a reference climate and applied at 10 km spatial resolution across a domain that includes the three main spreading centres of the Celtic, Fennoscandian and Barents Sea ice sheets. The model is forced with the NGRIP palaeo-isotope curve from 37 ka BP onwards and model skill is assessed against collated flowsets, marginal moraines, exposure ages and relative sea-level history. The evolution of the EISC to its LGM configuration was complex and asynchronous; the western, maritime margins of the Fennoscandian and Celtic ice sheets responded rapidly and advanced across their continental shelves by 29 ka BP, yet the maximum aerial extent (5.48 × 106 km2) and volume (7.18 × 106 km3) of the ice complex was attained some 6 ka later at c. 22.7 ka BP. This maximum stand was short-lived as the North Sea and Atlantic margins were already in retreat whilst eastern margins were still advancing up until c. 20 ka BP. High rates of basal erosion are modelled beneath ice streams and outlet glaciers draining the Celtic and Fennoscandian ice sheets with extensive preservation elsewhere due to frozen subglacial conditions, including much of the Barents and Kara seas. Here, and elsewhere across the Norwegian shelf and North Sea, high pressure subglacial conditions would have promoted localised gas hydrate formation

Sensitivity of the Antarctic ice sheets to the warming of marine isotope substage 11c

Funding: This research has been supported by the Swedish Research Council (grant no. 2016-04422), the German Research Foundation (grant no. 1158-365737614), the US National Science Foundation (grant no. PLR-1542930), and the Norwegian Polar Institute/NARE (grant no. 2015/38/7/NK/ihs). Jorge Bernales has been supported by the MAGIC-DML project through DFG SPP 1158 (RO 4262/1-6). The article processing charges for this open-access publication were covered by Stockholm University.Studying the response of the Antarctic ice sheets during periods when climate conditions were similar to the present can provide important insights into current observed changes and help identify natural drivers of ice sheet retreat. In this context, the marine isotope substage 11c (MIS11c) interglacial offers a suitable scenario, given that during its later portion orbital parameters were close to our current interglacial. Ice core data indicate that warmer-than-present temperatures lasted for longer than during other interglacials. However, the response of the Antarctic ice sheets and their contribution to sea level rise remain unclear. We explore the dynamics of the Antarctic ice sheets during this period using a numerical ice sheet model forced by MIS11c climate conditions derived from climate model outputs scaled by three glaciological and one sedimentary proxy records of ice volume. Our results indicate that the East and West Antarctic ice sheets contributed 4.0–8.2 m to the MIS11c sea level rise. In the case of a West Antarctic Ice Sheet collapse, which is the most probable scenario according to far-field sea level reconstructions, the range is reduced to 6.7–8.2 m independently of the choices of external sea level forcing and millennial-scale climate variability. Within this latter range, the main source of uncertainty arises from the sensitivity of the East Antarctic Ice Sheet to a choice of initial ice sheet configuration. We found that the warmer regional climate signal captured by Antarctic ice cores during peak MIS11c is crucial to reproduce the contribution expected from Antarctica during the recorded global sea level highstand. This climate signal translates to a modest threshold of 0.4 °C oceanic warming at intermediate depths, which leads to a collapse of the West Antarctic Ice Sheet if sustained for at least 4000 years.Publisher PDFPeer reviewe

The protocataclasite dilemma: in situ 36Cl and REE-Y lessons from an impure limestone fault scarp at Sparta, Greece

Reconstructions of palaeoseismicity are useful for understanding and mitigating seismic hazard risks. We apply cosmogenic 36Cl exposure-age dating and measurements of rare-earth elements and yttrium (REE-Y) concentrations to the palaeoseismic history of the Sparta Fault, Greece. Bayesian-inference Markov chain Monte Carlo (MCMC) modelling of 36Cl concentrations along a 7.2 m long vertical profile on the Sparta Fault scarp at Anogia indicate an increase in the average slip rate of the scarp from 0.8–0.9 mm yr−1 6.5–7.7 kyr ago to 1.1–1.2 mm yr−1 up to the devastating 464 BCE earthquake. The average exhumation of the entire scarp up to the present day is 0.7–0.8 mm yr−1. Modelling does not indicate additional exhumation of the Sparta Fault after 464 BCE. The Sparta Fault scarp is composed of fault breccia, containing quartz and clay-lined pores, in addition to host-rock-derived clasts of calcite and microcrystalline calcite cement. The impurities control the distribution of REE-Y in the fault scarp surface and contribute spatial variation to 36Cl concentrations, which precludes the identification of individual earthquakes that have exhumed the Sparta Fault scarp from either of these data sets. REE-Y may illustrate processes that localize slip to a discrete fault plane in the Earth's near-surface, but their potential use in palaeoseismicity would benefit from further evaluation.</p

Nunataks as barriers to ice flow : implications for palaeo ice sheet reconstructions

Funding: This research has been supported by the Vetenskapsrådet (grant no. 2016-04422), the Deutsche Forschungsgemeinschaft (grant no. 1158-365737614), the National Science Foundation (grant no. OPP-1542930), and the Norsk Polarinstitutt (grant no. 2015/38/7/NK/ihs).Numerical models predict that discharge from the polar ice sheets will become the largest contributor to sea-level rise over the coming centuries. However, the predicted amount of ice discharge and associated thinning depends on how well ice sheet models reproduce glaciological processes, such as ice flow in regions of large topographic relief, where ice flows around bedrock summits (i.e. nunataks) and through outlet glaciers. The ability of ice sheet models to capture long-term ice loss is best tested by comparing model simulations against geological data. A benchmark for such models is ice surface elevation change, which has been constrained empirically at nunataks and along margins of outlet glaciers using cosmogenic exposure dating. However, the usefulness of this approach in quantifying ice sheet thinning relies on how well such records represent changes in regional ice surface elevation. Here we examine how ice surface elevations respond to the presence of strong topographic relief that acts as an obstacle by modelling ice flow around and between idealised nunataks during periods of imposed ice sheet thinning. We find that, for realistic Antarctic conditions, a single nunatak can exert an impact on ice thickness over 20 km away from its summit, with its most prominent effect being a local increase (decrease) of the ice surface elevation of hundreds of metres upstream (downstream) of the obstacle. A direct consequence of this differential surface response for cosmogenic exposure dating is a delay in the time of bedrock exposure upstream relative to downstream of a nunatak summit. A nunatak elongated transversely to ice flow is able to increase ice retention and therefore impose steeper ice surface gradients, while efficient ice drainage through outlet glaciers produces gentler gradients. Such differences, however, are not typically captured by continent-wide ice sheet models due to their coarse grid resolutions. Their inability to capture site-specific surface elevation changes appears to be a key reason for the observed mismatches between the timing of ice-free conditions from cosmogenic exposure dating and model simulations. We conclude that a model grid refinement over complex topography and information about sample position relative to ice flow near the nunatak are necessary to improve data–model comparisons of ice surface elevation and therefore the ability of models to simulate ice discharge in regions of large topographic relief.Publisher PDFPeer reviewe

Deglaciation of Fennoscandia

To provide a new reconstruction of the deglaciation of the Fennoscandian Ice Sheet, in the form of calendar-year time-slices, which are particularly useful for ice sheet modelling, we have compiled and synthesized published geomorphological data for eskers, ice-marginal formations, lineations, marginal meltwater channels, striae, ice-dammed lakes, and geochronological data from radiocarbon, varve, optically-stimulated luminescence, and cosmogenic nuclide dating. This 25 is summarized as a deglaciation map of the Fennoscandian Ice Sheet with isochrons marking every 1000 years between 22 and 13 cal kyr BP and every hundred years between 11.6 and final ice decay after 9.7 cal kyr BP. Deglaciation patterns vary across the Fennoscandian Ice Sheet domain, reflecting differences in climatic and geomorphic settings as well as ice sheet basal thermal conditions and terrestrial versus marine margins. For example, the ice sheet margin in the high-precipitation coastal setting of the western sector responded sensitively to climatic variations leaving a detailed record of prominent moraines and ice-marginal deposits in many fjords and coastal valleys. Retreat rates across the southern sector differed between slow retreat of the terrestrial margin in western and southern Sweden and rapid retreat of the calving ice margin in the Baltic Basin. Our reconstruction is consistent with much of the published research. However, the synthesis of a large amount of existing and new data support refined reconstructions in some areas. For example, we locate the LGM extent of the ice sheet in northwestern Russia further east than previously suggested and conclude that it occurred at a later time than the rest of the ice sheet, at around 17-15 cal kyr BP, and propose a slightly different chronology of moraine formation over southern Sweden based on improved correlations of moraine segments using new LiDAR data and tying the timing of moraine formation to Greenland ice core cold stages. Retreat rates vary by as much as an order of magnitude in different sectors of the ice sheet, with the lowest rates on the high-elevation and maritime Norwegian margin. Retreat rates compared to the climatic information provided by the Greenland ice core record show a general correspondence between retreat rate and climatic forcing, although a close match between retreat rate and climate is unlikely because of other controls, such as topography and marine versus terrestrial margins. Overall, the time slice reconstructions of Fennoscandian Ice Sheet deglaciation from 22 to 9.7 cal kyr BP provide an important dataset for understanding the contexts that underpin spatial and temporal patterns in retreat of the Fennoscandian Ice Sheet, and are an important resource for testing and refining ice sheet models

Changes in vertical ice extent along the East Antarctic Ice Sheet margin in western Dronning Maud Land – initial field and modelling results of the MAGIC-DML collaboration

Constraining numerical ice sheet models by comparison with observational data is crucial to address the interactions between cryosphere and climate at a wide range of scales. Such models are tested and refined by comparing model predictions of past ice geometries with field-based reconstructions from geological, geomorphological, and ice core data. However, for the East Antarctic Ice sheet, there is a critical gap in the empirical data necessary to reconstruct changes in ice sheet geometry in the Dronning Maud Land (DML) region. In addition, there is poor control on the regional climate history of the ice sheet margin, because ice-core locations, where detailed reconstructions of climate history exist, are located on high inland domes. This leaves numerical models ofregional glaciation history in this near-coastal area largely unconstrained. MAGIC-DML is an ongoing Swedish-US-Norwegian-German-UK collaboration with a focus on improvingice sheet models of the western DML margin by combining advances in modeling with filling critical data gaps regarding the timing and pattern of ice-surface changes. A combination of geomorphological mapping using remote sensing data, field observations, cosmogenic nuclide surface exposure dating, and numerical ice sheetmodeling are being used in an iterative manner to produce a comprehensive reconstruction of the glacial historyof western DML. Here, we present an overview of the project, field evidence for formerly higher ice surfaces and in-situ cosmogenic nuclide measurements from the 2016/17 expedition. Preliminary field evidence indicate that interior sectors of DML have experienced a general decrease in ice sheet thickness since the late Miocene, with potential episodes of increasing thickness in the late Pleistocene (700-300 ka, 250-75 ka). To aid in interpreting these field data, new high-resolution ice sheet model reconstructions, constraining ice sheet configurations during key episodes, are presented

Mid-Pleistocene ice sheet fluctuations from cosmogenic nuclide field constraints in western Dronning Maud Land, Antarctica

The East Antarctic Ice Sheet (EAIS) is generally assumed to have been relatively insensitive to Quaternary climat echange. However, recent studies have shown potential instabilities in coastal, marine sectors of the EAIS. In addition, long-term climate reconstructions and modelling experiments indicate the potential for significant changes in ice volume and ice sheet configuration since the Pliocene. Hence, more empirical evidence for ice surface and ice volume changes is required to discriminate between contrasting inferences. MAGIC-DML is an ongoing Swedish-US-Norwegian-German-UK collaboration focused on improving ice sheetm odels by filling critical data gaps that exist in our knowledge of the timing and pattern of ice surface changes along the western Dronning Maud Land (DML) margin and combining this with advances in numerical techniques. As part of the project, field studies in the 2016/17 and 2017/18 austral summers targeted selected sites spanning accessible altitudes in the Heimefrontfjella, Vestfjella, Ahlmannryggen, Borgmassivet, and Kirwanveggen nunatakranges for in situcosmogenic nuclide sampling. Comparing concentrations of nuclides with widely differing half-lives in bedrock and erratics from a range of altitudes above modern ice surfaces can provide information on ice sheet fluctuations and complex burial and exposure histories, and thus, past configurations of non-erosive ice. Quartz-bearing rock types were sampled and analyzed for 10Be (t1/21.4 My),14C (t1/25.7 ky),26Al (t1/2705ky), and 21Ne (stable), and mafic lithologies for36Cl (t1/2301 ky). Results thus far for 3210Be and 26Al isotope pairs complemented with seven21Ne measurements have yielded some consistent patterns of paleoglaciation for the western DML margin. Eight out of fourteen bedrock samples from high-elevation (1700-2238 m a.s.l.) ridges and summits return some of the oldest exposure ages in Antarctica and have consistent 10Be,26Al, and 21Ne minimum apparent exposure ages of 1.8-4.1 Ma. Initial results therefore indicate that parts of the ice sheet marginal to the Antarctic plateau, along the Heimefrontfjella range, generally have experienced a decrease in ice thickness since the late Miocene. Another six bedrock samples (1556-1732 ma.s.l.) fall in the 300-700 ka range, and they all show significant burial. At face value, perhaps this indicates aregional ice sheet surface above 1700 m a.s.l. for much of the Plio-early Pleistocene. All other samples analyzedto date are erratics from lower elevation and more coastal sites (10 from nunataks at 553-1400 m a.s.l., and 6 froma surface moraine at 1385 m a.s.l.), exhibiting ages between 59 and 275 ka, save for two (4 and 6 ka). Whereas almost all of the nunatak erratics (including the young ones) show significant burial durations, five of the six surface moraine samples do not. These 2016/17 field samples are not yet leading to conclusive age constraints but already start to paint a picture of the western DML margin being relatively stable although there was possibly one or more episodes of relatively limited ice thickening during the last 700 ka

DECIPHERING THE FINEST IMPRINT OF GLACIAL EROSION: OBJECTIVE ANALYSIS OF STRIAE PATTERNS ON BEDROCK

The aim of this study is to compare the efficiency of different mathematical and statistical geometrical methods applied to characterise the orientation distribution of striae on bedrock for deciphering the finest imprint of glacial erosion. The involved methods include automatic image analysis techniques of Fast Fourier Transform (FFT), and the experimental investigations by means of Saltikov's directed secants analysis (rose of intersection densities), applied to digital and analogue images of the striae pattern, respectively. In addition, the experimental data were compared with the modelling results made on the basis of Underwood's concept of linear systems in a plane. The experimental and modelling approaches in the framework of stereology yield consistent results. These results reveal that stereological methods allow a reliable and efficient delineation of different families of glacial striae from a complex record imprinted in bedrock