A new global ice sheet reconstruction for the past 80 000 years

Global sea-level rise in the early Holocene revealed from North Sea peats

Rates of relative sea-level rise during the final stage of the last deglaciation, the early Holocene, are key to understanding future ice melt and sea-level change under a warming climate1. Data about these rates are scarce2, and this limits insight into the relative contributions of the North American and Antarctic ice sheets to global sea-level rise during the early Holocene. Here we present an early Holocene sea-level curve based on 88 sea-level data points (13.7-6.2 thousand years ago (ka)) from the North Sea (Doggerland3,4). After removing the pattern of regional glacial isostatic adjustment caused by the melting of the Eurasian Ice Sheet, the residual sea-level signal highlights two phases of accelerated sea-level rise. Meltwater sourced from the North American and Antarctic ice sheets drove these two phases, peaking around 10.3 ka and 8.3 ka with rates between 8 mm yr-1 and 9 mm yr-1. Our results also show that global mean sea-level rise between 11 ka and 3 ka amounted to 37.7 m (2σ range, 29.3-42.2 m), reconciling the mismatch that existed between estimates of global mean sea-level rise based on ice-sheet reconstructions and previously limited early Holocene sea-level data. With its broad spatiotemporal coverage, the North Sea dataset provides critical constraints on the patterns and rates of the late-stage deglaciation of the North American and Antarctic ice sheets, improving our understanding of the Earth-system response to climate change.

The ghost of ice ages past: Impact of Last Glacial Maximum landscapes on modern biodiversity

Modeled modern and Last Glacial Maximum (LGM) climate ranges for 47 genetically confirmed small Holarctic land snails documented profound landscape dynamism over the last 21,000 years. Following deglaciation, range areas tended to increase by 50% while isolating barrier widths were cut in half. At the same time, the nature of isolating barriers underwent profound change, with the North American continental ice sheet becoming as important in the LGM as the Atlantic Ocean is today in separating Nearctic and Palearctic faunas. Because appropriate modern climate occurs for these species throughout the Holarctic, with no clear barriers being present-especially for such efficient passive dispersers-the current >90% turnover observed between Eurasian and North American species pools appears at least in part related to the LGM landscape. Understanding current and predicting potential future biodiversity patterns thus requires consideration of the landscape template across at least 15,000 years time scales.

Mid-Pleistocene climate transition triggered by Antarctic Ice Sheet growth

Despite extensive investigation, the nature and causes of the Mid-Pleistocene Transition remain enigmatic. In this work, we assess its linkage to asynchronous development of bipolar ice sheets by synthesizing Pleistocene mid- to high-latitude proxy records linked to hemispheric ice sheet evolution. Our results indicate substantial growth of the Antarctic Ice Sheets (AISs) at 2.0 to 1.25 million years ago, preceding the rapid expansion of Northern Hemisphere Ice Sheets after ~1.25 million years ago. Proxy-model comparisons suggest that AIS and associated Southern Ocean sea ice expansion can induce northern high-latitude cooling and enhanced moisture transport to the Northern Hemisphere, thus triggering the Mid-Pleistocene Transition. The dynamic processes involved are crucial for assessing modern global warming that is already inducing asynchronous bipolar melting of ice sheets.

Paleoclimate data assimilation with CLIMBER-X: An ensemble Kalman filter for the last deglaciation

Using the climate model CLIMBER-X, we present an efficient method for assimilating the temporal evolution of surface temperatures for the last deglaciation covering the period 22000 to 6500 years before the present. The data assimilation methodology combines the data and the underlying dynamical principles governing the climate system to provide a state estimate of the system, which is better than that which could be obtained using just the data or the model alone. In applying an ensemble Kalman filter approach, we make use of the advances in the parallel data assimilation framework (PDAF), which provides parallel data assimilation functionality with a relatively small increase in computation time. We find that the data assimilation solution depends strongly on the background evolution of the decaying ice sheets rather than the assimilated temperatures. Two different ice sheet reconstructions result in a different deglacial meltwater history, affecting the large-scale ocean circulation and, consequently, the surface temperature. We find that the influence of data assimilation is more pronounced on regional scales than on the global mean. In particular, data assimilation has a stronger effect during millennial warming and cooling phases, such as the Bølling-Allerød and Younger Dryas, especially at high latitudes with heterogeneous temperature patterns. Our approach is a step toward a comprehensive paleo-reanalysis on multi-millennial time scales, including incorporating available paleoclimate data and accounting for their uncertainties in representing regional climates.

Genomic Analyses Capture the Human-Induced Demographic Collapse and Recovery in a Wide-Ranging Cervid

The glacial cycles of the Quaternary heavily impacted species through successions of population contractions and expansions. Similarly, populations have been intensely shaped by human pressures such as unregulated hunting and land use changes. White-tailed and mule deer survived in different refugia through the Last Glacial Maximum, and their populations were severely reduced after the European colonization. Here, we analyzed 73 resequenced deer genomes from across their North American range to understand the consequences of climatic and anthropogenic pressures on deer demographic and adaptive history. We found strong signals of climate-induced vicariance and demographic decline; notably, multiple sequentially Markovian coalescent recovers a severe decline in mainland white-tailed deer effective population size (Ne) at the end of the Last Glacial Maximum. We found robust evidence for colonial overharvest in the form of a recent and dramatic drop in Ne in all analyzed populations. Historical census size and restocking data show a clear parallel to historical Ne estimates, and temporal Ne/Nc ratio shows patterns of conservation concern for mule deer. Signatures of selection highlight genes related to temperature, including a cold receptor previously highlighted in woolly mammoth. We also detected immune genes that we surmise reflect the changing land use patterns in North America. Our study provides a detailed picture of anthropogenic and climatic-induced decline in deer diversity and clues to understanding the conservation concerns of mule deer and the successful demographic recovery of white-tailed deer.

Dansgaard-Oeschger cycles of the penultimate and last glacial period recorded in stalagmites from Türkiye

The last glacial period is characterized by abrupt climate oscillations, also known as Dansgaard-Oeschger (D-O) cycles. However, D-O cycles remain poorly documented in climate proxy records covering the penultimate glacial period. Here we present highly resolved and precisely dated speleothem time series from Sofular Cave in northern Türkiye to provide clear evidence for D-O cycles during Marine Isotope Stage (MIS) 6 as well as MIS 2-4. D-O cycles are most clearly expressed in the Sofular carbon isotope time series, which correlate inversely with regional sea surface temperature (SST) records from the Black Sea. The pacing of D-O cycles is almost twice as long during MIS 6 compared to MIS 2-4, and could be related to a weaker Atlantic Meridional Overturning Circulation (AMOC) and a different mean climate during MIS 6 compared to MIS 2-4, leading most likely to a higher threshold for the occurrence of D-O cycles.

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.

The ROCEEH Out of Africa Database (ROAD): A large-scale research database serves as an indispensable tool for human evolutionary studies

Large scale databases are critical for helping scientists decipher long-term patterns in human evolution. This paper describes the conception and development of such a research database and illustrates how big data can be harnessed to formulate new ideas about the past. The Role of Culture in Early Expansions of Humans (ROCEEH) is a transdisciplinary research center whose aim is to study the origins of culture and the multifaceted aspects of human expansions across Africa and Eurasia over the last three million years. To support its research, the ROCEEH team developed an online tool named the ROCEEH Out of Africa Database (ROAD) and implemented its web-based applications. ROAD integrates geographical data as well as archaeological, paleoanthropological, paleontological and paleobotanical content within a robust chronological framework. In fact, a unique feature of ROAD is its ability to dynamically link scientific data both spatially and temporally, thereby allowing its reuse in ways that were not originally conceived. The data stem from published sources spanning the last 150 years, including those generated by the research team. Descriptions of these data rely on the development of a standardized vocabulary and profit from online explanations of each table and attribute. By synthesizing legacy data, ROAD facilitates the reuse of heritage data in novel ways. Database queries yield structured information in a variety of interoperable formats. By visualizing data on maps, users can explore this vast dataset and develop their own theories. By downloading data, users can conduct further quantitative analyses, for example with Geographic Information Systems, modeling programs and artificial intelligence. In this paper, we demonstrate the innovative nature of ROAD and show how it helps scientists studying human evolution to access datasets from different fields, thereby connecting the social and natural sciences. Because it permits the reuse of "old" data in new ways, ROAD is now an indispensable tool for researchers of human evolution and paleogeography.

Geographic Genetic Structure of Alectoris chukar in Türkiye: Post-LGM-Induced Hybridization and Human-Mediated Contaminations

Türkiye is considered an important evolutionary area for Chukar partridge (Alectoris chukar), since it is both a potential ancestral area and a diversification center for the species. Using 2 mitochondrial (Cty-b and D-loop) and 13 polymorphic microsatellite markers, we investigated the geographic genetic structure of A. chukar populations to determine how past climatic fluctuations and human activities have shaped the gene pool of this species in Türkiye. Our results indicate, firstly, that only A. chukar of the genus Alectoris is present in Türkiye (Anatolia and Thrace), with no natural or artificial gene flow from congenerics. Secondly, the geographic genetic structure of the species in Türkiye has been shaped by topographic heterogeneity, Pleistocene climatic fluctuations, and artificial transport by humans. Third, there appears to be three genetic clusters: Thracian, Eastern, and Western. Fourth, the post-LGM demographic expansion of the Eastern and Western populations has formed a hybrid zone in Central Anatolia (~8 kyBP). Fifth, the rate of China clade-B contamination in Türkiye is about 8% in mtDNA and about 12% in nuDNA, with the Southeastern Anatolian population having the highest contamination. Sixth, the Thracian population was the most genetically distinct, with the lowest genetic diversity and highest level of inbreeding and no China clad-B contamination. These results can contribute to the conservation regarding A. chukar populations, especially the Thracian population.

Ice and ocean constraints on early human migrations into North America along the Pacific coast

Founding populations of the first Americans likely occupied parts of Beringia during the Last Glacial Maximum (LGM). The timing, pathways, and modes of their southward transit remain unknown, but blockage of the interior route by North American ice sheets between ~26 and 14 cal kyr BP (ka) favors a coastal route during this period. Using models and paleoceanographic data from the North Pacific, we identify climatically favorable intervals when humans could have plausibly traversed the Cordilleran coastal corridor during the terminal Pleistocene. Model simulations suggest that northward coastal currents strengthened during the LGM and at times of enhanced freshwater input, making southward transit by boat more difficult. Repeated Cordilleran glacial-calving events would have further challenged coastal transit on land and at sea. Following these events, ice-free coastal areas opened and seasonal sea ice was present along the Alaskan margin until at least 15 ka. Given evidence for humans south of the ice sheets by 16 ka and possibly earlier, we posit that early people may have taken advantage of winter sea ice that connected islands and coastal refugia. Marine ice-edge habitats offer a rich food supply and traversing coastal sea ice could have mitigated the difficulty of traveling southward in watercraft or on land over glaciers. We identify 24.5 to 22 ka and 16.4 to 14.8 ka as environmentally favorable time periods for coastal migration, when climate conditions provided both winter sea ice and ice-free summer conditions that facilitated year-round marine resource diversity and multiple modes of mobility along the North Pacific coast.

Terrestrial amplification of past, present, and future climate change

Terrestrial amplification (TA) of land warming relative to oceans is apparent in recent climatic observations. TA results from land-sea coupling of moisture and heat and is therefore important for predicting future warming and water availability. However, the theoretical basis for TA has never been tested outside the short instrumental period, and the spatial pattern and amplitude of TA remain uncertain. Here, we investigate TA during the Last Glacial Maximum (LGM; ~20 thousand years) in the low latitudes, where the theory is most applicable. We find remarkable consistency between paleotemperature proxies, theory, and climate model simulations of both LGM and future climates. Paleoclimate data thus provide crucial new support for TA, refining the range of future low-latitude, low-elevation TA to [Formula: see text] (95% confidence interval), i.e., land warming ~40% more than oceans. The observed data model theory agreement helps reconcile LGM marine and terrestrial paleotemperature proxies, with implications for equilibrium climate sensitivity.

Beringia and the peopling of the Western Hemisphere

Did Beringian environments represent an ecological barrier to humans until less than 15 000 years ago or was access to the Americas controlled by the spatial-temporal distribution of North American ice sheets? Beringian environments varied with respect to climate and biota, especially in the two major areas of exposed continental shelf. The East Siberian Arctic Shelf ('Great Arctic Plain' (GAP)) supported a dry steppe-tundra biome inhabited by a diverse large-mammal community, while the southern Bering-Chukchi Platform ('Bering Land Bridge' (BLB)) supported mesic tundra and probably a lower large-mammal biomass. A human population with west Eurasian roots occupied the GAP before the Last Glacial Maximum (LGM) and may have accessed mid-latitude North America via an interior ice-free corridor. Re-opening of the corridor less than 14 000 years ago indicates that the primary ancestors of living First Peoples, who already had spread widely in the Americas at this time, probably dispersed from the NW Pacific coast. A genetic 'arctic signal' in non-arctic First Peoples suggests that their parent population inhabited the GAP during the LGM, before their split from the former. We infer a shift from GAP terrestrial to a subarctic maritime economy on the southern BLB coast before dispersal in the Americas from the NW Pacific coast.

The Bering Strait was flooded 10,000 years before the Last Glacial Maximum

The cyclic growth and decay of continental ice sheets can be reconstructed from the history of global sea level. Sea level is relatively well constrained for the Last Glacial Maximum (LGM, 26,500 to 19,000 y ago, 26.5 to 19 ka) and the ensuing deglaciation. However, sea-level estimates for the period of ice-sheet growth before the LGM vary by > 60 m, an uncertainty comparable to the sea-level equivalent of the contemporary Antarctic Ice Sheet. Here, we constrain sea level prior to the LGM by reconstructing the flooding history of the shallow Bering Strait since 46 ka. Using a geochemical proxy of Pacific nutrient input to the Arctic Ocean, we find that the Bering Strait was flooded from the beginning of our records at 46 ka until [Formula: see text] ka. To match this flooding history, our sea-level model requires an ice history in which over 50% of the LGM's global peak ice volume grew after 46 ka. This finding implies that global ice volume and climate were not linearly coupled during the last ice age, with implications for the controls on each. Moreover, our results shorten the time window between the opening of the Bering Land Bridge and the arrival of humans in the Americas.

Pleistocene sea level changes and crocodile population histories on the Isthmus of Panama: A comment on Avila-Cervantes et al. (2020)

Avila-Cervantes et al. proposed that glacial-interglacial sea level changes played an important role in the evolutionary and demographic histories of the crocodile Crocodylus acutus on the Isthmus of Panama. However, the study used erroneous sea level proxy data that produced flawed paleogeographic reconstructions. We present new paleogeographic reconstructions and review the timing of, and proposed mechanisms behind, the demographic events estimated by Avila-Cervantes et al.. With the data currently available, we find little evidence to support the hypothesis that sea level changes drove population demographic events in crocodiles on the Isthmus. Alternative hypotheses, including changing climate and habitat suitability, are equally valid and should be considered along with well-supported sea level models.

Sea ice fluctuations in the Baffin Bay and the Labrador Sea during glacial abrupt climate changes

Sea ice decline in the North Atlantic and Nordic Seas has been proposed to contribute to the repeated abrupt atmospheric warmings recorded in Greenland ice cores during the last glacial period, known as Dansgaard-Oeschger (D-O) events. However, the understanding of how sea ice changes were coupled with abrupt climate changes during D-O events has remained incomplete due to a lack of suitable high-resolution sea ice proxy records from northwestern North Atlantic regions. Here, we present a subdecadal-scale bromine enrichment (Br enr ) record from the NEEM ice core (Northwest Greenland) and sediment core biomarker records to reconstruct the variability of seasonal sea ice in the Baffin Bay and Labrador Sea over a suite of D-O events between 34 and 42 ka. Our results reveal repeated shifts between stable, multiyear sea ice (MYSI) conditions during cold stadials and unstable, seasonal sea ice conditions during warmer interstadials. The shift from stadial to interstadial sea ice conditions occurred rapidly and synchronously with the atmospheric warming over Greenland, while the amplitude of high-frequency sea ice fluctuations increased through interstadials. Our findings suggest that the rapid replacement of widespread MYSI with seasonal sea ice amplified the abrupt climate warming over the course of D-O events and highlight the role of feedbacks associated with late-interstadial seasonal sea ice expansion in driving the North Atlantic ocean–climate system back to stadial conditions.

Terrestrial evidence for ocean forcing of Heinrich events and subglacial hydrologic connectivity of the Laurentide Ice Sheet

During the last glacial period, the Laurentide Ice Sheet (LIS) underwent episodes of rapid iceberg discharge, recorded in ocean sediments as "Heinrich events" (HEs). Two competing models attempt to describe the stimulus for HEs via either internal ice sheet oscillations or external ocean-climate system forcing. We present a terrestrial record of HEs from the northeastern LIS that strongly supports ocean-climate forcing. Subglacial carbonate precipitates from Baffin Island record episodes of subglacial melting coincident with the three most recent HEs, resulting from acceleration of nearby marine-terminating ice streams. Synchronized ice stream acceleration over Baffin Island and Hudson Strait is inconsistent with internal ice sheet oscillations alone and indicates a shared ocean-climate stimulus to coordinate these different glaciological systems. Isotopic compositions of these precipitates record widespread subglacial groundwater connectivity beneath the LIS. Extensive basal melting and flushing of these aquifers during the last HE may have been a harbinger for terminal deglaciation.

Glacial isostatic adjustment: physical models and observational constraints

By far the most prescient insights into the interior structure of the planet have been provided on the basis of elastic wave seismology. Analysis of the travel times of shear or compression wave phases excited by individual earthquakes, or through analysis of the elastic gravitational free oscillations that individual earthquakes of sufficiently large magnitude may excite, has been the central focus of Earth physics research for more than a century. Unfortunately, data provide no information that is directly relevant to understanding the solid state 'flow' of the polycrystalline outer 'mantle' shell of the planet that is involved in the thermally driven convective circulation that is responsible for powering the 'drift' of the continents and which controls the rate of planetary cooling on long timescales. For this reason, there has been an increasing focus on the understanding of physical phenomenology that is unambiguously associated with mantle flow processes that are distinct from those directly associated with the convective circulation itself. This paper reviews the past many decades of work that has been invested in understanding the most important of such processes, namely that which has come to be referred to as 'glacial isostatic adjustment' (GIA). This process concerns the response of the planet to the loading and unloading of the high latitude continents by the massive accumulations of glacial ice that have occurred with almost metronomic regularity over the most recent million years of Earth history. Forced by the impact of gravitationaln-body effects on the geometry of Earth's orbit around the Sun through the impact upon the terrestrial regime of received solar insolation, these surface mass loads on the continents have left indelible records of their occurrence in the 'Earth system' consisting of the oceans, continents, and the great polar ice sheets on Greenland and Antarctica themselves. Although this ice-age phenomenology has been clearly recognized since early in the last century, it was for over 50 years considered to be no more than an interesting curiosity, the understanding of which remained on the periphery of the theoretical physics of the Earth. This was the case in part because no globally applicable theory was available that could be applied to rigorously interpret the observations. Equally important to understanding the scientific lethargy that held back the understanding of this phenomenon involving mantle flow processes was the lack of appreciation of the wide range of observations that were in fact related to GIA physics. This paper is devoted to a review of the global theories of the GIA process that have since been developed as a means of interpreting the extensive variety of observations that are now recognized as being involved in the response of the planet to the loading and unloading of its surface by glacial ice. The paper will also provide examples of the further analyses of Earth physics and climate related processes that applications of the modern theoretical structures have enabled.

Norway spruce postglacial recolonization of Fennoscandia

Contrasting theories exist regarding how Norway spruce (Picea abies) recolonized Fennoscandia after the last glaciation and both early Holocene establishments from western microrefugia and late Holocene colonization from the east have been postulated. Here, we show that Norway spruce was present in southern Fennoscandia as early as 14.7 ± 0.1 cal. kyr BP and that the millennia-old clonal spruce trees present today in central Sweden likely arrived with an early Holocene migration from the east. Our findings are based on ancient sedimentary DNA from multiple European sites (N = 15) combined with nuclear and mitochondrial DNA analysis of ancient clonal (N = 135) and contemporary spruce forest trees (N = 129) from central Sweden. Our other findings imply that Norway spruce was present shortly after deglaciation at the margins of the Scandinavian Ice Sheet, and support previously disputed finds of pollen in southern Sweden claiming spruce establishment during the Lateglacial.

Thermochronologic constraints on the origin of the Great Unconformity

The Great Unconformity involves a common gap of hundreds of millions to billions of years in the geologic record. The cause of this missing time has long eluded explanation, but recently two opposing hypotheses claim either a glacial or a plate tectonic origin in the Neoproterozoic. We provide thermochronologic evidence of rock cooling and multiple kilometers of exhumation in the Cryogenian Period in support of a glacial origin for erosion contributing to the composite basement nonconformity found across the North American interior. The broad synchronicity of this cooling signal at the continental scale can only be readily explained by glacial denudation.

The origin of the phenomenon known as the Great Unconformity has been a fundamental yet unresolved problem in the geosciences for over a century. Recent hypotheses advocate either global continental exhumation averaging 3 to 5 km during Cryogenian (717 to 635 Ma) snowball Earth glaciations or, alternatively, diachronous episodic exhumation throughout the Neoproterozoic (1,000 to 540 Ma) due to plate tectonic reorganization from supercontinent assembly and breakup. To test these hypotheses, the temporal patterns of Neoproterozoic thermal histories were evaluated for four North American locations using previously published medium- to low-temperature thermochronology and geologic information. We present inverse time–temperature simulations within a Bayesian modeling framework that record a consistent signal of relatively rapid, high-magnitude cooling of ∼120 to 200 °C interpreted as erosional exhumation of upper crustal basement during the Cryogenian. These models imply widespread, synchronous cooling consistent with at least ∼3 to 5 km of unroofing during snowball Earth glaciations, but also demonstrate that plate tectonic drivers, with the potential to cause both exhumation and burial, may have significantly influenced the thermal history in regions that were undergoing deformation concomitant with glaciation. In the cratonic interior, however, glaciation remains the only plausible mechanism that satisfies the required timing, magnitude, and broad spatial pattern of continental erosion revealed by our thermochronological inversions. To obtain a full picture of the extent and synchroneity of such erosional exhumation, studies on stable cratonic crust below the Great Unconformity must be repeated on all continents.

Impact of paleoclimate on present and future evolution of the Greenland Ice Sheet

Using transient climate forcing based on simulations from the Alfred Wegener Institute Earth System Model (AWI-ESM), we simulate the evolution of the Greenland Ice Sheet (GrIS) from the last interglacial (125 ka, kiloyear before present) to 2100 AD with the Parallel Ice Sheet Model (PISM). The impact of paleoclimate, especially Holocene climate, on the present and future evolution of the GrIS is explored. Our simulations of the past show close agreement with reconstructions with respect to the recent timing of the peaks in ice volume and the climate of Greenland. The maximum and minimum ice volume at around 18-17 ka and 6-5 ka lag the respective extremes in climate by several thousand years, implying that the ice volume response of the GrIS strongly lags climatic changes. Given that Greenland's climate was getting colder from the Holocene Thermal Maximum (i.e., 8 ka) to the Pre-Industrial era, our simulation implies that the GrIS experienced growth from the mid-Holocene to the industrial era. Due to this background trend, the GrIS still gains mass until the second half of the 20th century, even though anthropogenic warming begins around 1850 AD. This is also in agreement with observational evidence showing mass loss of the GrIS does not begin earlier than the late 20th century. Our results highlight that the present evolution of the GrIS is not only controlled by the recent climate changes, but is also affected by paleoclimate, especially the relatively warm Holocene climate. We propose that the GrIS was not in equilibrium throughout the entire Holocene and that the slow response to Holocene climate needs to be represented in ice sheet simulations in order to predict ice mass loss, and therefore sea level rise, accurately.

Landscape evolution under the southern Laurentide Ice Sheet

Subglacial landscapes, revealed in regions of recent ice-sheet retreat, provide a window into ice-sheet dynamics and interactions with evolving subglacial topography. Here, we document landscape evolution beneath the southern Laurentide Ice Sheet of North America since the end of the Pliocene, 2.6 million years (Ma) ago, by reconstructing the isostatically adjusted preglacial surface and modern bedrock topography at 250 m horizontal resolution. We use flow routing to reconstruct drainage networks and river longitudinal profiles, revealing the pattern and extent of their glacially forced reorganization. The overall mean Quaternary (2.6 Ma ago to present) erosion rate is 27 m/Ma, rising within ice-streaming corridors to 35 m/Ma (and locally reaching 400 m/Ma) and falling to 22 m/Ma in non–ice-streaming regions. Our results suggest that subglacial erosion was sufficient to lower the southern Laurentide Ice Sheet into warmer environments, thereby enhancing ablation and reducing ice-sheet extent over time.