Last Glacial Maximum Ice Research Articles

Pronounced changes of subterranean biodiversity patterns along a Late Pleistocene glaciation gradient

Understanding spatial patterns of biodiversity within the context of long‐term climatic shifts is of high importance, particularly in the face of contemporary climate change. In comparison to aboveground taxa, subterranean organisms respond to changing climates with generally much lower dispersal and recolonization potential, yet possible persistence in refugial groundwater habitats under ice‐shields. However, knowledge on general and geographically large‐scale effects of glaciation on contemporary groundwater biodiversity patterns is still very limited. Here, we tested how Late Pleistocene glaciation influenced the diversity and distribution of 36 groundwater amphipod species in Alpine and peri‐Alpine regions, characterized by extensive glaciation cycles, and how its legacy explains contemporary diversity patterns. We based our analysis on an unprecedented density of ~ 1000 systematic sampling sites across Switzerland. Using presence–absence data, we assessed biodiversity and species' ranges, and calculated for each site within‐catchment distance to the Last Glacial Maximum (LGM) glacier extent. We then applied a sliding window approach along the obtained distance gradient from LGM ice‐covered to ice‐free sites to compute biodiversity indices reflecting local richness, regional richness, and differentiation, respectively. We found a strong signal of the LGM ice extent on the present‐day distribution of groundwater amphipods. Our findings revealed pronounced species turnover and spatial envelopes of individual species' occurrences in formerly ice‐covered, ice‐free, or transitional zones, respectively. While local richness remained constant and low along the LGM distance gradient, groundwater communities in LGM ice‐covered areas were more similar to each other and had lower gamma diversities and decreased occurrence probabilities per sliding window compared to communities in Pleistocene ice‐free areas. These results highlight the significant impact of Pleistocene glaciation on shaping biodiversity patterns of subterranean communities and imprinting contemporary distribution of groundwater organisms.

Dynamical response of the southwestern Laurentide Ice Sheet to rapid Bølling–Allerød warming

Abstract. The shift in climate that occurred between the Last Glacial Maximum (LGM) and the Early Holocene (ca. 18–12 kyr BP) displayed rates of temperature increase similar to present-day warming trends. The most rapid recorded changes in temperature occurred during the abrupt climate oscillations known as the Bølling–Allerød interstadial (14.7–12.9 kyr BP) and the Younger Dryas stadial (12.9–11.7 kyr BP). Reconstructing ice sheet dynamics during these climate oscillations provides the opportunity to assess long-term ice sheet evolution in reaction to a rapidly changing climate. Here, we use glacial geomorphological inversion methods (flowsets) to reconstruct the ice flow dynamics and the marginal retreat pattern of the southwestern sector of the Laurentide Ice Sheet (SWLIS). We combine our reconstruction with a recently compiled regional deglaciation chronology to depict ice flow dynamics that encompass the time period from pre-LGM to the Early Holocene. Our reconstruction portrays three macroscale reorganizations in the orientation and dynamics of ice streaming followed by regional deglaciation associated with rapid warming during the Bølling–Allerød interstadial. Initial westward flow is documented, likely associated with an early set of ice streams that formed during the advance to the LGM. During the LGM ice streaming displays a dominant north to south orientation. Ice sheet thinning at ∼15 ka is associated with a macroscale reorganization in ice stream flow, with a complex of ice streams recording south-eastward flow. A second macroscale reorganization in ice flow is then observed at ∼14 ka, in which southwestern ice flow is restricted to the Hay, Peace, Athabasca, and Churchill river lowlands. Rates of ice sheet retreat then slowed considerably during the Younger Dryas stadial; at this time, the ice margin was situated north of the Canadian Shield boundary and ice flow continued to be sourced from the northeast. Resulting from these changes in ice sheet dynamics, we recognize a three-part pattern of deglacial landform zonation within the SWLIS characterized by active ice margin recession, stagnation and downwasting punctuated by local surging (terrestrial ice sheet collapse): the outer deglacial zone contains large recessional moraines aligned with the direction of active ice margin retreat; the intermediate deglacial zone contains large regions of hummocky and stagnation terrain, in some areas crosscut by the signature of local surges, reflecting punctuated stagnation and downwasting; and the inner deglacial zone contains inset recessional moraines demarcating progressive regional ice margin retreat. We attribute these macroscale changes in ice flow geometry and associated deglacial behaviour to external climatic controls during the Bølling–Allerød and Younger Dryas but also recognize the role of internal (glaciological, lithological and topographic) controls in SWLIS dynamics.

Modeling the timing of Patagonian Ice Sheet retreat in the Chilean Lake District from 22–10 ka

Abstract. Studying the retreat of the Patagonian Ice Sheet (PIS) during the last deglaciation represents an important opportunity to understand how ice sheets outside the polar regions have responded to deglacial changes in temperature and large-scale atmospheric circulation. At the northernmost extension of the PIS during the Last Glacial Maximum (LGM), the Chilean Lake District (CLD) was influenced by the southern westerly winds (SWW), which strongly modulated the hydrologic and heat budgets of the region. Despite progress in constraining the nature and timing of deglacial ice retreat across this area, considerable uncertainty in the glacial history still exists due to a lack of geologic constraints on past ice margin change. Where the glacial chronology is lacking, ice sheet models can provide important insight into our understanding of the characteristics and drivers of deglacial ice retreat. Here we apply the Ice Sheet and Sea-level System Model (ISSM) to simulate the LGM and last deglacial ice history of the PIS across the CLD at high spatial resolution (450 m). We present a transient simulation of ice margin change across the last deglaciation using climate inputs from the National Center for Atmospheric Research Community Climate System Model (CCSM3) Trace-21ka experiment. At the LGM, the simulated ice extent across the CLD agrees well with the most comprehensive reconstruction of PIS ice history (PATICE). Coincident with deglacial warming, ice retreat ensues after 19 ka, with large-scale ice retreat occurring across the CLD between 18 and 16.5 ka. By 17 ka, the northern portion of the CLD becomes ice free, and by 15 ka, ice only persists at high elevations as mountain glaciers and small ice caps. Our simulated ice history agrees well with PATICE for early deglacial ice retreat but diverges at and after 15 ka, where the geologic reconstruction suggests the persistence of an ice cap across the southern CLD until 10 ka. However, given the high uncertainty in the geologic reconstruction of the PIS across the CLD during the later deglaciation, this work emphasizes a need for improved geologic constraints on past ice margin change. While deglacial warming drove the ice retreat across this region, sensitivity tests reveal that modest variations in wintertime precipitation (∼10 %) can modulate the pacing of ice retreat by up to 2 ka, which has implications when comparing simulated outputs of ice margin change to geologic reconstructions. While we find that TraCE-21ka simulates large-scale changes in the SWW across the CLD that are consistent with regional paleoclimate reconstructions, the magnitude of the simulated precipitation changes is smaller than what is found in proxy records. From our sensitivity analysis, we can deduce that larger anomalies in precipitation, as found in paleoclimate proxies, may have had a large impact on modulating the magnitude and timing of deglacial ice retreat. This fact highlights an additional need for better constraints on the deglacial change in strength, position, and extent of the SWW as it relates to understanding the drivers of deglacial PIS behavior.

Millennial-scale variability of Greenland dust provenance during the last glacial maximum as determined by single particle analysis

Greenland ice core records exhibited 100-fold higher dust concentrations during the Last Glacial Maximum (LGM) than during the Holocene, and dust input temporal variability corresponded to different climate states in the LGM. While East Asian deserts, the Sahara, and European loess have been suggested as the potential source areas (PSAs) for Greenland LGM dust, millennial-scale variability in their relative contributions within the LGM remains poorly constrained. Here, we present the morphological, mineralogical, and geochemical characteristics of insoluble microparticles to constrain the provenance of dust in Greenland NEEM ice core samples covering cold Greenland Stadials (GS)-2.1a to GS-3 (~ 14.7 to 27.1 kyr ago) in the LGM. The analysis was conducted on individual particles in microdroplet samples by scanning electron microscopy with energy dispersive X-ray spectroscopy and Raman microspectroscopy. We found that the kaolinite-to-chlorite (K/C) ratios and chemical index of alteration (CIA) values were substantially higher (K/C: 1.4 ± 0.7, CIA: 74.7 ± 2.9) during GS-2.1a to 2.1c than during GS-3 (K/C: 0.5 ± 0.1, CIA: 65.8 ± 2.8). Our records revealed a significant increase in Saharan dust contributions from GS-2.1a to GS-2.1c and that the Gobi Desert and/or European loess were potential source(s) during GS-3. This conclusion is further supported by distinctly different carbon contents in particles corresponding to GS-2.1 and GS-3. These results are consistent with previous estimates of proportional dust source contributions obtained using a mixing model based on Pb and Sr isotopic compositions in NEEM LGM ice and indicate millennial-scale changes in Greenland dust provenance that are probably linked to large-scale atmospheric circulation variabilities during the LGM.

Relative importance of the mechanisms triggering the Eurasian ice sheet deglaciation in the GRISLI2.0 ice sheet model

Abstract. The last deglaciation (21 to 8 ka) of the Eurasian ice sheet (EIS) is thought to have been responsible for a sea level rise of about 20 m. While many studies have examined the timing and rate of the EIS retreat during this period, many questions remain about the key processes that triggered the EIS deglaciation 21 kyr ago. Due to its large marine-based parts in the Barents–Kara (BKIS) and British Isles sectors, the BKIS is often considered to be a potential analogue of the current West Antarctic ice sheet (WAIS). Identifying the mechanisms that drove the EIS evolution might provide a better understanding of the processes at play in the West Antarctic destabilization. To investigate the relative impact of key drivers on the EIS destabilization, we used the three-dimensional ice sheet model GRISLI (GRenoble Ice Shelf and Land Ice) (version 2.0) forced by climatic fields from five Paleoclimate Modelling Intercomparison Project phases 3 and 4 (PMIP3, PMIP4) Last Glacial Maximum (LGM) simulations. In this study, we performed sensitivity experiments to test the response of the simulated Eurasian ice sheets to surface climate, oceanic temperatures (and thus basal melting under floating ice tongues), and sea level perturbations. Our results highlight that the EIS retreat simulated with the GRISLI model is primarily triggered by atmospheric warming. Increased atmospheric temperatures further amplify the sensitivity of the ice sheets to sub-shelf melting. These results contradict those of previous modelling studies mentioning the central role of basal melting on the deglaciation of the marine-based Barents–Kara ice sheet. However, we argue that the differences with previous works are mainly related to differences in the methodology followed to generate the initial LGM ice sheet. Due to the strong sensitivity of EIS to the atmospheric forcing highlighted with the GRISLI model and the limited extent of the confined ice shelves during the LGM, we conclude by questioning the analogy between EIS and the current WAIS. However, because of the expected rise in atmospheric temperatures, the risk of hydrofracturing is increasing and could ultimately put the WAIS in a configuration similar to the past Eurasian ice sheet.

Unraveling the complexities of the Last Glacial Maximum climate: the role of individual boundary conditions and forcings

Abstract. In order to quantify the relative importance of individual boundary conditions and forcings, including greenhouse gases, ice sheets, and Earth's orbital parameters, on determining Last Glacial Maximum (LGM) climate, we have performed a series of LGM experiments using a state-of-the-art climate model AWI-ESM, in which different combinations of boundary conditions and forcings have been applied following the protocol of Paleoclimate Modelling Intercomparison Project phase 4 (PMIP4). In good agreement with observational proxy records, a general colder and drier climate is simulated in our full-forced LGM experiment as compared to the present-day simulation. Our simulated results from non-full-forced sensitivity simulations reveal that both the greenhouse gases and ice sheets play a major role in defining the anomalous LGM surface temperature compared to today. Decreased greenhouse gases in LGM as compared to present day leads to a non-uniform global cooling with polar amplification effect. The presence of LGM ice sheets favors a warming over the Arctic and northern Atlantic oceans in boreal winter, as well as a cooling over regions with the presence of ice sheets. The former is induced by a strengthening in the Atlantic meridional overturning circulation (AMOC) transporting more heat to high latitudes, whilst the latter is due to the increased surface albedo and elevation of ice sheets. We find that the Northern Hemisphere monsoon precipitation is influenced by the opposing effects of LGM greenhouse gases and ice sheets. Specifically, the presence of ice sheets leads to significant drying in the Northern Hemisphere monsoon regions, while a reduction in greenhouse gases results in increased monsoon rainfall. Based on our model results, continental ice sheets exert a major control on atmospheric dynamics and the variability of El Niño–Southern Oscillation (ENSO). Moreover, our analysis also implies a nonlinearity in climate response to LGM boundary conditions and forcings.

Limited extension of the MIS 2 proglacial lake in the Severnaya Dvina valley, south‐eastern margin of the last Scandinavian Ice Sheet

ABSTRACTThe Severnaya Dvina River valley crosses the former south‐eastern margin of the last Scandinavian Ice Sheet. Despite a long research history, there remains considerable controversy about the maximum ice‐sheet extent and the expansion of proglacial lakes within the Severnaya Dvina fluvial system. The goal of this study was to address these issues using new material from the valleys of the Severnaya Dvina and the lower Vychegda, thereby contributing to an understanding of the history of the south‐eastern sector of the last Scandinavian Ice Sheet and its periglacial areas. We studied a number of geological sections using radiocarbon and optically stimulated luminescence (OSL) dating, pollen and carpological analyses, and found that the Last Glacial Maximum (LGM) proglacial lake occupied the Severnaya Dvina valley between 20–19 ka and about 15.5 ka. The lake was localized within the proglacial isostatic depression and occupied only the Severnaya Dvina valley extending no further than the confluence of the Vychegda River, no more than 110 km from the edge of the ice sheet. The lake formation did not cause any drainage reorganization within and outside the Vychegda fluvial system. Given the small size of the lake and presence of the oldest OSL ages along the entire length of the former lake, we suggest that the onset of the proglacial lake marks the maximum extent of the ice sheet in the area. Consequently, the onset of the local LGM may be dated to 20–19 ka, somewhat earlier than assumed by most previous researchers. The LGM ice sheet boundary was located in the Severnaya Dvina valley downstream from the Vychegda confluence and did not extend into the Vychegda valley, despite previous suggestions. During deglaciation, the lake disappearance together with crustal rebound caused an incision episode in the Vychegda – Severnaya Dvina system and the formation of the Lateglacial alluvial terrace with relative height rising downstream due to the uneven rates of the postglacial uplift.

Morphometric analysis of ice scour lakes in Iceland: A proxy for ice sheet dynamics

AbstractGlacial erosion rates on the Iceland landscape throughout the Pleistocene are not well quantified. Ice scour lakes provide an opportunity to investigate glacial erosive activity and relative intensity because they commonly form in areas beneath ice sheets due to intense quarrying. This study evaluates ice scour lake morphology and density as a potential proxy for paleo‐ice flow direction and basal thermal regime for parts of the Iceland Ice Sheet. Using GIS analysis, properties of ice scour lakes are examined in regions that experienced different rates of ice flow during and following the Last Glacial Maximum (LGM), as interpreted from streamlined subglacial landforms. The primary input datasets were topographic, hydrologic and bedrock data. Lake distribution and morphology were quantified for all natural lakes of Iceland (n = 30 169), with close examination of three 400 km2 sub‐regions in Vestfirðir (n = 904), Húnaflói (n = 69) and Bakkaheiði (n = 232). Lake distribution parameters include density, packing and centroid elevation. PolyMorph‐2D was used to quantify lake morphology including orientation, area and elongation ratio. Ice scour lake density, packing and elevation were all greatest on Vestfirðir. The high density and packing of lakes across Vestfirðir suggests unique ice dynamics relative to the other two sub‐regions, and multidirectional flow of lake axes supports the presence of an independent ice cap covering Vestfirðir. Differential intensity of glacial erosion interpreted from regions of high and low lake density on Bakkaheiði supports a proposed ice divide in Northeast Iceland. The orientation of lakes align with the flow directions of proposed LGM ice streams in Northern and Northeast Iceland. Ice scour lakes were most elongate on average on Bakkaheiði, which supports fast ice flow. Improving understanding of former ice sheet basal thermal regime and paleo‐ice flow velocity serves to inform modern controls on ice sheet stability and ice sheet basal processes.

Subglacial hydrology from high-resolution ice-flow simulations of the Rhine Glacier during the Last Glacial Maximum: a proxy for glacial erosion

Abstract. At the Last Glacial Maximum (LGM), the Rhine Glacier complex (Rhine and Linth glaciers) formed large piedmont lobes extending north into the Swiss and German Alpine forelands. Numerous overdeepened valleys there were formed by repeated glaciations. A characteristic of these overdeepened valleys is their location close to the LGM ice margin, away from the Alps. Numerical models of ice flow of the Rhine Glacier indicate a poor fit between the sliding distance, a proxy for glacial erosion, and the location of these overdeepenings. Calculations of the hydraulic potential based on the computed time-dependent ice surface elevations of the Rhine Glacier lobe obtained from a high-resolution thermo-mechanically coupled Stokes flow model are used to estimate the location of subglacial water drainage routes. Results indicate that the subglacial water discharge is high and focused along glacial valleys and overdeepenings when water pressure is equal to the ice overburden pressure. These conditions are necessary for subglacial water to remove basal sediments, expose fresh bedrock, and favor further erosion by quarrying and abrasion. Knowledge of the location of paleo-subglacial water drainage routes may be useful to understand patterns of subglacial erosion beneath paleo-ice masses that do not otherwise relate to the sliding of ice. Comparison of the erosion pattern from subglacial meltwater with those from quarrying and abrasion shows the importance of subglacial water flow in the formation of distal overdeepenings in the Swiss lowlands.

A steady-state model reconstruction of the patagonian ice sheet during the last glacial maximum

During the Last Glacial Maximum (LGM), the Patagonian Ice Sheet (PIS) was the largest Quaternary ice mass in the Southern Hemisphere outside of Antarctica. Although the margins of the LGM ice sheet are now well established through end-moraine mapping and dating, apart from a few modelling and empirical studies, there remains a lack of constraint on its thickness and three-dimensional configuration. Here, we provide a high-resolution steady-state model reconstruction of the PIS at its maximum - LGM - extent applied using Nye's perfect-plastic ice rheology. The yield-strength parameter for the perfect-plastic flow model was calibrated against independent empirical reconstructions of the Lago Pueyrredón Glacier, where the former vertical extent of this major outlet glacier is well constrained by cosmogenically-dated trimlines and lateral and end-moraine limits. Using this derived yield-strength parameter, the perfect-plastic model is then applied to multiple flowlines demarking each outlet across the entirety of the PIS in a GIS framework. Our results reveal that the area of the PIS was ∼504,500 km2 (±8.5%) with a corresponding modelled ice volume of ∼554,500 km3 (±10%), equivalent to ∼1.38 m (±10%) of eustatic sea-level lowering at the LGM. Maximum surface elevation was at least 3500m asl although the majority of the ice sheet surface was below 2500 m asl. We find that our ice sheet reconstruction is in good general agreement with previous estimates of net PIS volume derived from transient modelling studies. We attribute the slightly lower aspect-ratio of our ice sheet (and its concomitant 5% reduction in volume and sea-level equivalent) to the lower yield strength applied, based on more temperate and dynamic ice sheet conditions.

Paleoclimate Changes in the Pacific Northwest Over the Past 36,000 Years From Clumped Isotope Measurements and Model Analysis

AbstractSince the last glacial period, North America has experienced dramatic changes in regional climate, including the collapse of ice sheets and changes in precipitation. We use clumped isotope (∆47) thermometry and carbonate δ18O measurements of glacial and deglacial pedogenic carbonates from the Palouse Loess to provide constraints on hydroclimate changes in the Pacific Northwest. We also employ analysis of climate model simulations to help us further provide constraints on the hydroclimate changes in the Pacific Northwest. The coldest clumped isotope soil temperatures T(47) (13.5 ± 1.9°C to 17.1 ± 1.7°C) occurred ∼34,000–23,000 years ago. Using a soil‐to‐air temperature transfer function, we estimate Last Glacial Maximum (LGM) mean annual air temperatures of ∼−5.5°C and warmest average monthly temperatures (i.e., mean summer air temperatures) of ∼4.4°C. These data indicate a regional warming of 16.4 ± 2.6°C from the LGM to the modern temperatures of 10.9°C, which was about 2.5–3 times the global average. Proxy data provide locality constraints on the boundary of the cooler anticyclone induced by LGM ice sheets, and the warmer cyclone in the Eastern Pacific Ocean. Climate model analysis suggests regional amplification of temperature anomalies is due to the proximal location of the study area to the Laurentide Ice Sheet margin and the impact of the glacial anticyclone on the region, as well as local albedo. Isotope‐enabled model experiments indicate variations in water δ18O largely reflect atmospheric circulation changes and enhanced rainout upstream that brings more depleted vapor to the region during the LGM.

Last Glacial Maximum glacier fluctuations on the northern Alpine foreland: Geomorphological and chronological reconstructions from the Rhine and Reuss glacier systems

Fresh glacial landforms of the Alpine forelands evidence the presence and extent of large piedmont glaciers during the Last Glacial Maximum (LGM) and yield valuable insights into LGM glacier dynamics. This study assesses widespread ice marginal landforms preserved within the limit of the former LGM Rhine glacier and the eastern lobes of the LGM Reuss glacier system by means of geomorphological mapping. Timing of formation of the studied ice margins in Rhine and Reuss systems are chronologically framed by new 10Be and 36Cl surface exposure ages, as well as new radiocarbon dates. This includes redating of radiocarbon samples first determined in the 1980s. Results of this and an earlier study focussing on outwash deposits downstream of the outer LGM ice margin, suggest that the Rhine glacier advanced to and reached its LGM maximum between ca. 26–22 ka, thereby forming a broad >100 km wide foreland piedmont lobe and constructing prominent and largely continuous chains of frontal moraines (outer Schaffhausen moraines). The eastern lobes of the Reuss glacier likely advanced to their LGM maximum position by 25/24 ± 2 ka as indicated by published luminescence dates. Stabilization of the corresponding Untertannwald ice margin and the slightly internal yet more prominent Mellingen moraines (Reuss glacier system) occurred no later than 22 ± 1 ka and 21 ± 1 ka, respectively. Glacier oscillations following the LGM maximum position are evidenced in both Rhine and Reuss systems but show varying degrees of preservation. Largely contemporaneous, late LGM readvances of Rhine (Stein am Rhein stadial) and Reuss (Bremgarten stadial) glaciers occurred after 20.6 ± 1.7 ka and 20.8 ± 1.3 ka, respectively. Absence of upstream moraines suggests rapid ice decay without marked stabilization, thereafter. Despite clear differences in the size and nature of their foreland lobes, with Rhein glacier as a broad piedmont lobe and narrow valley glacier like lobes characterising the eastern Reuss system, a remarkable similarity in timing is shown between the two.

Reconstructing the advance and retreat dynamics of the central sector of the last Cordilleran Ice Sheet

The advance of the Cordilleran Ice Sheet (CIS) towards its Last Glacial Maximum (LGM) configuration and its subsequent retreat remain poorly understood. We use the glacial landform record to determine ice dynamics for the central sector of the CIS in northern British Columbia, Canada, beneath the LGM ice divide. We classify seventy ice-flow indicator flowsets based on morphology, elevation, orientation and cross-cutting relationships into one of three stages, whereby stage 1 is oldest and stage 3 youngest. Combined with ice-contact geomorphology, our reconstruction highlights complex changes in ice flow over time as a result of ice divide migrations through the LGM and deglacial phases. The orientation and distribution of landforms indicates active post-LGM ice retreat westward through the Cassiar and Omineca mountains. We map the regional distribution of independent mountain glaciers, ice caps, and ice fields that regrew during a cooling event in the Late Glacial and show that some of these readvance glaciers were subsequently overrun by advancing outlet glaciers of the CIS. We use the cross-cutting relationship between readvance glaciers and CIS outlet glaciers and available chronological data to reconstruct the eastern CIS margin during the Late Glacial for the first time.

Three-dimensional gravity modelling of a Quaternary overdeepening fill in the Bern area of Switzerland discloses two stages of glacial carving

The geometry of glacial overdeepenings on the Swiss Plateau close to Bern was inferred through a combination of gravity data with a 3D gravity modelling software. The target overdeepenings have depths between 155 and > 270 m and widths between 860 and 2400 m. The models show incisions characterized by U-shaped cross-sectional geometries and steep to over-steepened lateral flanks. Existing stratigraphic data reveals that the overdeepenings were formed and then filled during at least two glacial stages, which occurred during the Last Glacial Maximum (LGM) within the Marine Isotope Stage (MIS) 2, and possibly MIS 6 or before. The U-shaped cross-sectional geometries point towards glacial erosion as the main driver for the shaping of the overdeepenings. The combination of the geometries with stratigraphic data suggests that the MIS 6 (or older) glaciers deeply carved the bedrock, whereas the LGM ice sheet only widened the existing valleys but did not further deepen them. We relate this pattern to the different ice thicknesses, where a thicker MIS 6 ice was likely more powerful for wearing down the bedrock than a thinner LGM glacier. Gravity data in combination with forward modelling thus offers robust information on the development of a landscape formed through glaciers.

LGM ice extent and deglaciation history in the Gurktal and Lavantal Alps (eastern European Alps): first constraints from 10Be surface exposure dating of glacially polished quartz veins

AbstractCompared with the western European Alps, the ice extent during the Last Glacial Maximum (LGM) and the subsequent deglaciation history of the eastern Alps east of the Tauern Window remain less well constrained. Also, considerable discrepancies exist between the mapped LGM ice margin and the ice extent predicted by ice‐sheet models. Here we present the first 10Be surface exposures ages from two regions east of the Tauern Window (the Gurktal and Lavantal Alps), which provide constraints on the LGM ice extent and the deglaciation history. Our results show that the deglaciation of the Gurktal Alps occurred between 16 and 14 ka, which agrees with the predictions from ice‐sheet models. In contrast, the 10Be ages from the Lavantal Alps located farther east are either LGM in age or predate the LGM, indicating that these regions were ice free or only partially covered by LGM ice. This finding suggests that ice‐sheet models may have overestimated the LGM ice extent in the easternmost Alps. In conclusion, our study highlights the need for more age data from the eastern Alps to refine the location of the LGM ice margin and the deglaciation history, which is also crucial for climate‐evolution and postglacial‐rebound models.

Unraveling Glacial Hydroclimate in the Indo‐Pacific Warm Pool: Perspectives From Water Isotopes

AbstractThe Indo‐Pacific Warm Pool (IPWP) is home to the warmest sea surface temperatures in the world oceans, favoring strong tropospheric convection and heavy rainfall. The mechanisms controlling long‐term change in the region's hydroclimate are still uncertain. Here, we present a 450,000‐year record of precipitation δD from southern Sumatra that records a consistent pattern of glacial isotopic enrichment and interglacial depletion. We synthesize existing paleo‐indicators of precipitation δD and δ18O in the IPWP and compare results with water isotope‐enabled climate simulations of the Last Glacial Maximum (LGM). The simulations show glacial isotopic enrichment over the eastern Indian Ocean extending into the southern IPWP and isotopic depletion over Southeast Asia, the west Pacific, and Australia. The pattern of simulated LGM isotopic change agrees generally well with our proxy synthesis. We conclude that reorganization of regional circulation under glacial conditions controls precipitation isotope variability in the IPWP: Low‐level tropospheric convergence dominates the signal in the north/east, whereas divergence controls the response in the south/west. Additional sensitivity simulations suggest that the LGM ice sheets and the associated lowering in sea level, rather than decreased greenhouse gases, are responsible for the distinctive spatial pattern in glacial changes of precipitation isotopes and hydroclimate across the IPWP.

Climatic patterns over the European Alps during the LGM derived from inversion of the paleo-ice extent

Quaternary climate has been dominated by alternating glacial and interglacial periods. While the timing and extent of past ice caps are well documented, local variations in temperature and precipitation as a response to cyclic glaciations are not resolved. Resolving these issues is necessary for understanding regional and global climate circulation. In particular, the impact of the cold high-pressure zone above the Fennoscandian ice cap on the position of the jet stream in Europe, and a possible change in the direction and the source of moisture flow are still debated. Here we reconstruct climate conditions that led to the observed ice extent in the European Alps during the last glacial maximum (LGM). Using a new inverse method to reconstruct the spatially variable position of the equilibrium line altitude (ELA), we investigate whether the Alpine LGM ice cap dominantly received precipitation in the south due to a strong southward shift of the westerlies in midlatitudes. We report inversion results that enable us to estimate the role of climate, inversion method parameters, ice dynamics, flexure and topography in modulating the inferred equilibrium line altitude. Our main finding is the increase of the position of the ELA from west to east and from north to south of the mountain range, which suggest dominant moisture delivery over the Alps during the LGM relatively similar to today.

Coupled climate–carbon cycle simulation of the Last Glacial Maximum atmospheric CO<sub>2</sub> decrease using a large ensemble of modern plausible parameter sets

Abstract. During the Last Glacial Maximum (LGM), atmospheric CO2 was around 90 ppmv lower than during the pre-industrial period. The reasons for this decrease are most often elucidated through factorial experiments testing the impact of individual mechanisms. Due to uncertainty in our understanding of the real system, however, the different models used to conduct the experiments inevitably take on different parameter values and different structures. In this paper, the objective is therefore to take an uncertainty-based approach to investigating the LGM CO2 drop by simulating it with a large ensemble of parameter sets, designed to allow for a wide range of large-scale feedback response strengths. Our aim is not to definitely explain the causes of the CO2 drop but rather explore the range of possible responses. We find that the LGM CO2 decrease tends to predominantly be associated with decreasing sea surface temperatures (SSTs), increasing sea ice area, a weakening of the Atlantic Meridional Overturning Circulation (AMOC), a strengthening of the Antarctic Bottom Water (AABW) cell in the Atlantic Ocean, a decreasing ocean biological productivity, an increasing CaCO3 weathering flux and an increasing deep-sea CaCO3 burial flux. The majority of our simulations also predict an increase in terrestrial carbon, coupled with a decrease in ocean and increase in lithospheric carbon. We attribute the increase in terrestrial carbon to a slower soil respiration rate, as well as the preservation rather than destruction of carbon by the LGM ice sheets. An initial comparison of these dominant changes with observations and paleoproxies other than carbon isotope and oxygen data (not evaluated directly in this study) suggests broad agreement. However, we advise more detailed comparisons in the future, and also note that, conceptually at least, our results can only be reconciled with carbon isotope and oxygen data if additional processes not included in our model are brought into play.

Glacial geology and cosmogenic-nuclide exposure ages from the Tucker Glacier - Whitehall Glacier confluence, northern Victoria Land, Antarctica

We describe glacial-geological observations and cosmogenic-nuclide exposure ages from the vicinity of the present grounding line of Tucker Glacier, a large alpine glacier flowing from the mountains of northern Victoria Land, East Antarctica, into the outer Ross Sea. These data are relevant for constraining the extent of ice sheet expansion and retreat in the Ross Sea, and associated eustatic sea level impact, between the Last Glacial Maximum (LGM) and the present. In addition, a terrestrial geological record of ice thickness change from this region could provide evidence for or against the hypothesis that rapid eustatic sea-level rise during meltwater pulse 1A (“MWP-1A”) at 14.6 ka was in part the result of rapid, large-scale thinning or breakup of a marine-based portion of the LGM ice sheet in the outer Ross Sea. Glacial-geological observations, exposure ages on glacial deposits, and a novel application of in-situ-produced cosmogenic 14C in quartz-bearing bedrock to identify the limits of LGM ice cover in the absence of direct geomorphic evidence, show that Tucker Glacier near its present grounding line was 300 to 350 m thicker than present during the LGM and thinned steadily between 17 to 5 ka. The largest possible rapid thickness change in the time period 14 to 15 ka that could be accommodated by the exposure-age data is ∼50 m, which is a small fraction of that predicted for the western Ross Sea by model simulations of the Antarctic contribution to MWP-1A. There do exist possible scenarios in which hypothesized marine ice sheet collapse in the outer Ross Sea during MWP-1A might not be recorded by ice thickness changes at Tucker Glacier. However, our record of ice thickness changes spanning this time period is the closest such record to the outer Ross Sea that is likely to exist, and it agrees with all exposure-age deglaciation chronologies from other regions of the Ross embayment in providing no evidence for such an event.

Two AMOC States in Response to Decreasing Greenhouse Gas Concentrations in the Coupled Climate Model MPI-ESM

This study analyzes the response of the Atlantic meridional overturning circulation (AMOC) to different CO2 concentrations and two ice sheet configurations in simulations with the coupled climate model MPI-ESM. With preindustrial (PI) ice sheets, there are two different AMOC states within the studied CO2 range: one state with a strong and deep upper overturning cell at high CO2 concentrations and one state with a weak and shallow upper cell at low CO2 concentrations. Changes in AMOC variability with decreasing CO2 indicate two stability thresholds. The strong state is stable above the first threshold near 217 ppm, and the weak state is stable below the second threshold near 190 ppm. Between the two thresholds, both states are marginally unstable, and the AMOC oscillates between them on millennial time scales. The weak AMOC state is stable when Antarctic Bottom Water becomes dense and salty enough to replace North Atlantic Deep Water (NADW) in the deep North Atlantic and when the density gain over the North Atlantic becomes too weak to sustain continuous NADW formation. With Last Glacial Maximum (LGM) ice sheets, the density gain over the North Atlantic and the northward salt transport are enhanced with respect to the PI ice sheet case. This enables active NADW formation and a strong AMOC for the entire range of studied CO2 concentrations. The AMOC variability indicates that the simulated AMOC is far away from a stability threshold with LGM ice sheets. The nonlinear relationship among AMOC, CO2, and prescribed ice sheets provides an explanation for the large intermodel spread of AMOC states found in previous coupled LGM simulations.