Description: We present an ensemble of last glacial inception (LGI) simulations for the Northern Hemisphere that captures a significant fraction of inferred ice volume changes within proxy uncertainties. This ensemble was performed with LCice 1.0, a coupled ice sheet and climate model, varying parameters of both climate and ice sheet components, as well as the coupling between them. Certain characteristics of the spatiotemporal pattern of ice growth and subsequent retreat in both North America (NA) and Eurasia (EA) are sensitive to parameter changes while others are not. We find that the initial inception of ice over NA and EA is best characterized by the nucleation of ice at high-latitude and high-elevation sites. Subsequent spreading and merger along with large-scale conversion of snowfields dominate in different sectors. The latter plays an important role in the merging of eastern and western ice regions in NA. The inception peak ice volume in the ensemble occurs approximately at 111 ka and therefore lags the summer 60∘ N insolation minimum by more than 3 kyr. Ice volumes consistently peak earlier over EA than NA. The inception peak in North America is characterized by a merged Laurentide and Cordilleran ice sheet, with the Davis Strait covered in ice in ∼80 % of simulations. Ice also bridges Greenland and Iceland in all runs by 114 ka and therefore blocks the Denmark Strait. This latter feature would thereby divert the East Greenland Current and Denmark Strait overflow with a potentially significant impact on ocean circulation. The Eurasian ice sheet at its inception peak varies across ensemble runs between a continuous ice sheet and multiple smaller ice caps. In both continents, the colder high latitudes (i.e. Ellesmere and Svalbard) tend to grow ice through the entire simulation (to 102 ka), while lower latitudes lose ice after ∼110 ka. We find temperature decreases over the initial phases of the inception lead to the expansion of NA ice sheet area and that subsequent precipitation increases contribute to its thickening. EA ice sheet area also expands with decreasing temperatures, but sea ice limits any increases in precipitation, leading to an earlier retreat away from the EA maximum ice sheet volume. We also examine the extent to which the capture of both LGI ice growth and retreat constrains the coupled ice–climate model sensitivity to changing atmospheric pCO2. The 55-member sub-ensemble that meets our criteria for “acceptable” ice growth and retreat has an equilibrium climate sensitivity lower bound that is 0.3 ∘C higher than that of the full ensemble. This suggests some potential value of fully coupled ice–climate modelling of the last glacial inception to constrain future climate change.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Love, R., Andres, H. J., Condron, A., & Tarasov, L. Freshwater routing in eddy-permitting simulations of the last deglacial: the impact of realistic freshwater discharge. Climate of the Past, 17(6), (2021): 2327–2341. https://doi.org/10.5194/cp-17-2327-2021.

Description: Freshwater, in the form of glacial runoff, is hypothesized to play a critical role in centennial- to millennial-scale climate variability, such as the Younger Dryas and Dansgaard–Oeschger events, but this relationship is not straightforward. Large-scale glacial runoff events, such as Meltwater Pulse 1a (MWP1a), are not always temporally proximal to subsequent large-scale cooling. Moreover, the typical design of hosing experiments that support this relationship tends to artificially amplify the climate response. This study explores the impact that limitations in the representation of runoff in conventional “hosing” simulations has on our understanding of this relationship by examining where coastally released freshwater is transported when it reaches the ocean. We particularly focus on the impact of (1) the injection of freshwater directly over sites of deep-water formation (DWF) rather than at runoff locations (i.e. hosing), (2) excessive freshwater injection volumes (often by a factor of 5), and (3) the use of present-day (rather than palaeo) ocean gateways. We track the routing of glaciologically constrained freshwater volumes from four different inferred injection locations in a suite of eddy-permitting glacial ocean simulations using the Massachusetts Institute of Technology General Circulation Model (MITgcm) under both open and closed Bering Strait conditions. Restricting freshwater forcing values to realistic ranges results in less spreading of freshwater across the North Atlantic and indicates that the freshwater anomalies over DWF sites depend strongly on the geographical location of meltwater input. In particular, freshwater released into the Gulf of Mexico generates a very weak freshwater signal over DWF regions as a result of entrainment by the turbulent Gulf Stream. In contrast, freshwater released into the Arctic with an open Bering Strait or from the Eurasian ice sheet is found to generate the largest salinity anomalies over DWF regions in the North Atlantic and GIN (Greenland–Iceland–Norwegian) seas region respectively. Experiments show that when the Bering Strait is open, the Mackenzie River source exhibits more than twice as much freshening of the North Atlantic deep-water formation regions as when the Bering Strait is closed. Our results illustrate that applying freshwater hosing directly into the North Atlantic with even “realistic” freshwater amounts still overestimates the amount of terrestrial runoff reaching DWF regions. Given the simulated salinity anomaly distributions and the lack of reconstructed impact on deep-water formation during the Bølling–Allerød, our results support that the majority of the North American contribution to MWP1a was not routed through the Mackenzie River.

Description: This research has been supported by the Bundesministerium für Bildung und Forschung (Research for Sustainability initiative, FONA, through the PalMod project), the Natural Sciences and Engineering Research Council of Canada Discovery Grant, and by the NSF (grant no. OCE-1903427).