Search

The Antarctic Ice Sheet (AIS) has experienced accelerated loss of ice over the last decades and could become the main contributor to sea level rise in the coming centuries. However, the associated uncertainty is very large. The main sources of this uncertainty lie in the future scenarios, the climatic forcing and, most notably, the structural uncertainty due to our lack of understanding of ice–ocean interaction processes, in particular, the representation of subshelf basal melt. In this study, we use a higher-order ice sheet model to investigate the impact of these three sources of uncertainty on the contribution of the AIS to sea level in the coming centuries in the context of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) but extending the projections to 2500. We test the sensitivity of the model to basal melting parameters using several forcings and scenarios simulated in the CMIP5 and CMIP6 ensembles. Results show a strong dependency on the values of the parameter that controls the heat exchange velocity between ice and ocean as well as the forcing and scenario. Higher values of the heat exchange parameter lead to higher sea level rise, with the contribution depending on the forcing–scenario configuration and reaching in some cases more than 3 m in sea level equivalent by the end of 2500. Idealized simulations considering the individual effects of the atmospheric and the oceanic forcing have been performed, demonstrating that the oceanic forcing plays a dominant role over the western sector of the AIS, while the atmospheric forcing is more important for the eastern sector and the interior.

A thorough understanding of ice thermodynamics is essential for an accurate description of glaciers, ice sheets and ice shelves. Yet there exists a significant gap in our theoretical knowledge of the time-dependent behaviour of ice temperatures due to the inevitable compromise between mathematical tractability and the accurate description of physical phenomena. In order to bridge this shortfall, we have analytically solved the 1D time-dependent advective–diffusive heat problem including additional terms due to strain heating and depth-integrated horizontal advection. Newton's law of cooling is applied as a Robin-type top boundary condition to consider potential non-equilibrium temperature states across the ice–air interface. The solution is expressed in terms of confluent hypergeometric functions following a separation of variables approach. Non-dimensionalization reduces the parameter space to four numbers that fully determine the shape of the solution at equilibrium: surface insulation, effective geothermal heat flow, the Péclet number and the Brinkman number. The initial temperature distribution exponentially converges to the stationary solution. Transient decay timescales are only dependent on the Péclet number and the surface insulation, so higher advection rates and lower insulating values imply shorter equilibration timescales, respectively. In contrast, equilibrium temperature profiles are mostly independent of the surface insulation parameter. We have extended our study to a broader range of vertical velocities by using a general power-law dependence on depth, unlike prior studies limited to linear and quadratic velocity profiles. Lastly, we present a suite of benchmark experiments to test numerical solvers. Four experiments of gradually increasing complexity capture the main physical processes for heat propagation. Analytical solutions are then compared to their numerical counterparts upon discretization over unevenly spaced coordinate systems. We find that a symmetric scheme for the advective term and a three-point asymmetric scheme for the basal boundary condition best match our analytical solutions. A further convergence study shows that n≥15 vertical points are sufficient to accurately reproduce the temperature profile. The solutions presented herein are general and fully applicable to any problem with an equivalent set of boundary conditions and any given initial temperature distribution.

Tipping elements, including the Antarctic Ice Sheet (AIS), are Earth system components that could reach critical thresholds due to anthropogenic emissions. Increasing our understanding of past warm climates can help to elucidate the future contribution of the AIS to emissions. The mid-Pliocene Warm Period (mPWP; ∼ 3.3–3.0 million years ago) serves as an ideal benchmark experiment. During this period, CO2 levels were similar to the present day (PD; 350–450 ppmv), but global mean temperatures were 2.5–4.0 K higher. Sea level reconstructions from that time indicate a rise of 5–25 m compared to the present, highlighting the potential crossing of tipping points in Antarctica. In order to achieve a sea level contribution far beyond 10 m, not only the West Antarctic Ice Sheet (WAIS) needs to largely decrease, but a significant response in the East Antarctic Ice Sheet (EAIS) is also required. A key question in reconstructions and simulations is therefore which of the AIS basins retreated during the mPWP. In this study, we investigate how the AIS responds to climatic and bedrock conditions during the mPWP. To this end, we use the Pliocene Model Intercomparison Project, Phase 2 (PlioMIP2), general circulation model ensemble to force a higher-order ice sheet model. Our simulations reveal that the WAIS experiences collapse with a 0.5 K oceanic warming. The Wilkes Basin shows retreat at 3 K oceanic warming, although higher precipitation rates could mitigate such a retreat. Totten Glacier shows slight signs of retreats only under high-oceanic warming conditions (greater than 4 K oceanic anomaly). If only the WAIS collapses, we simulate a mean contribution of 2.7 to 7.0 ms.l.e. (metres of sea level equivalent). If, in addition, the Wilkes Basin retreats, our simulations suggest a mean contribution of 6.0 to 8.9 ms.l.e. Besides uncertainties related to the climate forcing, we also examine other sources of uncertainty related to initial ice thickness and ice dynamics. We find that the climatologies yield a higher uncertainty than the dynamical configuration if parameters are constrained with PD observations and that starting from Pliocene reconstructions leads to smaller ice sheet configurations due to the hysteresis behaviour of marine bedrocks. Ultimately, our study concludes that marine ice cliff instability is not a prerequisite for the retreat of the Wilkes Basin. Instead, a significant rise in oceanic temperatures can initiate such a retreat.

The vast majority of ice-sheet modelling studies rely on simplified representations of the glacial isostatic adjustment (GIA), which, among other limitations, do not account for lateral variations in the lithospheric thickness and upper-mantle viscosity. In studies of the last glacial cycle using 3D GIA models, this has however been shown to have major impacts on the dynamics of marine-based sectors of Antarctica, which are likely to be the greatest contributors to sea-level rise in the coming centuries. This gap in comprehensiveness is explained by the fact that 3D GIA models are computationally expensive, rarely open-source and require a complex coupling scheme. To close this gap between “best” and “tractable” GIA models, we propose FastIsostasy here, a regional GIA model capturing lateral variations in the lithospheric thickness and mantle viscosity. By means of fast Fourier transforms and a hybrid collocation scheme to solve its underlying partial differential equation, FastIsostasy can simulate 100 000 years of high-resolution bedrock displacement in only minutes of single-CPU computation, including the changes in sea-surface height due to mass redistribution. Despite its 2D grid, FastIsostasy parameterises the depth-dependent viscosity and therefore represents the depth dimension to a certain extent. FastIsostasy is benchmarked here against analytical, as well as 1D and 3D numerical solutions, and shows good agreement with them. For a simulation of the last glacial cycle, its mean and maximal error over time and space respectively yield less than 5 % and 16 % compared to a 3D GIA model over the regional solution domain. FastIsostasy is open-source, is documented with many examples and provides a straightforward interface for coupling to an ice-sheet model. The model is benchmarked here based on its implementation in Julia, while a Fortran version is also provided to allow for compatibility with most existing ice-sheet models. The Julia version provides additional features, including a vast library of adaptive time-stepping methods and GPU support.

In the last decades, great effort has been made to reconstruct the Laurentide Ice Sheet (LIS) during the Last Glacial Maximum (LGM; ca. 21 000 years before present, 21 kyr ago). Uncertainties underlying its modelling have led to notable differences in fundamental features such as its maximum elevation, extent and total volume. As a result, the uncertainty in ice dynamics and thus in ice extent, volume and ice stream stability remains large. We herein use a higher-order three-dimensional ice sheet model to simulate the LIS under LGM boundary conditions for a number of basal friction formulations of varying complexity. Their consequences for the Laurentide ice streams, configuration, extent and volume are explicitly quantified. Total volume and ice extent generally reach a constant equilibrium value that falls close to prior LIS reconstructions. Simulations exhibit high sensitivity to the dependency of the basal shear stress on the sliding velocity. In particular, a regularised Coulomb friction formulation appears to be the best choice in terms of ice volume and ice stream realism. Pronounced differences are found when the basal friction stress is thermomechanically coupled: the base remains colder, and the LIS volume is lower than in the purely mechanical friction scenario counterpart. Thermomechanical coupling is fundamental for producing rapid ice streaming, yet it leads to a similar ice distribution overall.

Little is known about the distribution of ice in the Antarctic Ice Sheet (AIS) during the Last Glacial Maximum (LGM). Whereas marine and terrestrial geological data indicate that the grounded ice advanced to a position close to the continental-shelf break, the total ice volume is unclear. Glacial boundary conditions are potentially important sources of uncertainty, in particular basal friction and climatic boundary conditions. Basal friction exerts a strong control on the large-scale dynamics of the ice sheet and thus affects its size and is not well constrained. Glacial climatic boundary conditions determine the net accumulation and ice temperature and are also poorly known. Here we explore the effect of the uncertainty in both features on the total simulated ice storage of the AIS at the LGM. For this purpose we use a hybrid ice sheet shelf model that is forced with different basal drag choices and glacial background climatic conditions obtained from the LGM ensemble climate simulations of the third phase of the Paleoclimate Modelling Intercomparison Project (PMIP3). Overall, we find that the spread in the simulated ice volume for the tested basal drag parameterizations is about the same range as for the different general circulation model (GCM) forcings (4 to 6 m sea level equivalent). For a wide range of plausible basal friction configurations, the simulated ice dynamics vary widely but all simulations produce fully extended ice sheets towards the continental-shelf break. More dynamically active ice sheets correspond to lower ice volumes, while they remain consistent with the available constraints on ice extent. Thus, this work points to the possibility of an AIS with very active ice streams during the LGM. In addition, we find that the surface boundary temperature field plays a crucial role in determining the ice extent through its effect on viscosity. For ice sheets of a similar extent and comparable dynamics, we find that the precipitation field determines the total AIS volume. However, precipitation is highly uncertain. Climatic fields simulated by climate models show more precipitation in coastal regions than a spatially uniform anomaly, which can lead to larger ice volumes. Our results strongly support using these paleoclimatic fields to simulate and study the LGM and potentially other time periods like the last interglacial. However, their accuracy must be assessed as well, as differences between climate model forcing lead to a large spread in the simulated ice volume and extension.

We describe the physics and features of the ice-sheet model Yelmo, an open-source project intended for collaborative development. Yelmo is a thermomechanical model, solving for the coupled velocity and temperature solutions of an ice sheet simultaneously. The ice dynamics are currently treated via a “hybrid” approach combining the shallow-ice and shallow-shelf/shelfy-stream approximations, which makes Yelmo an apt choice for studying a wide variety of problems. Yelmo's main innovations lie in its flexible and user-friendly infrastructure, which promotes portability and facilitates long-term development. In particular, all physics subroutines have been designed to be self-contained, so that they can be easily ported from Yelmo to other models, or easily replaced by improved or alternate methods in the future. Furthermore, hard-coded model choices are eschewed, replaced instead with convenient parameter options that allow the model to be adapted easily to different contexts. We show results for different ice-sheet benchmark tests, and we illustrate Yelmo's performance for the Antarctic ice sheet.

The Northeast Greenland Ice Stream (NEGIS) has been suffering a significant ice mass loss during the last decades. This is partly due to increasing oceanic temperatures in the subpolar North Atlantic, which enhance submarine basal melting and mass discharge. This demonstrates the high sensitivity of this region to oceanic changes. In addition, a recent study suggested that the NEGIS grounding line was 20–40 km behind its present-day location for 15 ka during Marine Isotope Stage (MIS) 3. This is in contrast with Greenland temperature records indicating cold atmospheric conditions at that time, expected to favour ice-sheet expansion. To explain this anomalous retreat a combination of atmospheric and external forcings has been invoked. Yet, as the ocean is found to be a primary driver of the ongoing retreat of the NEGIS glaciers, the effect of past oceanic changes in their paleo evolution cannot be ruled out and should be explored in detail. Here we investigate the sensitivity of the NEGIS to the oceanic forcing during the last glacial period using a three-dimensional hybrid ice-sheet–shelf model. We find that a sufficiently high oceanic forcing could account for a NEGIS ice-margin retreat of several tens of kilometres, potentially explaining the recently proposed NEGIS grounding-line retreat during Marine Isotope Stage 3.

The last glacial period (LGP; ca. 110–10 kyr BP) was marked by the existence of two types of abrupt climatic changes, Dansgaard–Oeschger (DO) and Heinrich (H) events. Although the mechanisms behind these are not fully understood, it is generally accepted that the presence of ice sheets played an important role in their occurrence. While an important effort has been made to investigate the dynamics and evolution of the Laurentide ice sheet (LIS) during this period, the Eurasian ice sheet (EIS) has not received much attention, in particular from a modeling perspective. However, meltwater discharge from this and other ice sheets surrounding the Nordic seas is often implied as a potential cause of ocean instabilities that lead to glacial abrupt climate changes. Thus, a better comprehension of the evolution of the EIS during the LGP is important to understand its role in glacial abrupt climate changes. Here we investigate the response of the EIS to millennial-scale climate variability during the LGP. We use a hybrid, three-dimensional, thermomechanical ice-sheet model that includes ice shelves and ice streams. The model is forced off-line via a novel perturbative approach that, as opposed to conventional methods, clearly differentiates between the spatial patterns of millennial-scale and orbital-scale climate variability. Thus, it provides a more realistic treatment of the forcing at millennial timescales. The effect of both atmospheric and oceanic variations are included. Our results show that the EIS responds with enhanced ice discharge in phase with interstadial warming in the North Atlantic when forced with surface ocean temperatures. Conversely, when subsurface ocean temperatures are used, enhanced ice discharge occurs both during stadials and at the beginning of the interstadials. Separating the atmospheric and oceanic effects demonstrates the major role of the ocean in controlling the dynamics of the EIS on millennial timescales. While the atmospheric forcing alone is only able to produce modest iceberg discharges, warming of the ocean leads to higher rates of iceberg discharges as a result of relatively strong basal melting at the margins of the ice sheet. Our results clearly show the capability of the EIS to react to glacial abrupt climate changes, and highlight the need for stronger constraints on the ice sheet's glacial dynamics and climate–ocean interactions.

Temperature reconstructions from Greenland ice-sheet (GrIS) ice cores indicate the occurrence of more than 20 abrupt warmings during the last glacial period (LGP) known as Dansgaard-Oeschger (D-O) events. Although their ultimate cause is still debated, evidence from both proxy data and modelling studies robustly links these to reorganisations of the Atlantic Meridional Overturning Circulation (AMOC). During the LGP, the GrIS expanded as far as the continental shelf break and was thus more directly exposed to oceanic changes than in the present. Therefore oceanic temperature fluctuations on millennial timescales could have had a non-negligible impact on the GrIS. Here we assess the effect of millennial-scale oceanic variability on the GrIS evolution from the last interglacial to the present day. To do so, we use a three-dimensional hybrid ice-sheet–shelf model forced by subsurface oceanic temperature fluctuations, assumed to increase during D-O stadials and decrease during D-O interstadials. Since in our model the atmospheric forcing follows orbital variations only, the increase in total melting at millennial timescales is a direct result of an increase in basal melting. We show that the GrIS evolution during the LGP could have been strongly influenced by oceanic changes on millennial timescales, leading to oceanically induced ice-volume contributions above 1 m sea level equivalent (SLE). Also, our results suggest that the increased flux of GrIS icebergs as inferred from North Atlantic proxy records could have been triggered, or intensified, by peaks in melting at the base of the ice shelves resulting from increasing subsurface oceanic temperatures during D-O stadials. Several regions across the GrIS could thus have been responsible for ice mass discharge during D-O events, opening the possibility of a non-negligible role of the GrIS in oceanic reorganisations throughout the LGP.

The Antarctic Ice Sheet (AIS) is the largest ice sheet on Earth and hence a major potential contributor to future global sea-level rise. A wealth of studies suggest that increasing oceanic temperatures could cause a collapse of its marine-based western sector, the West Antarctic Ice Sheet, through the mechanism of marine ice-sheet instability, leading to a sea-level increase of 3–5 m. Thus, it is crucial to constrain the sensitivity of the AIS to rapid climate changes. The last glacial period is an ideal benchmark period for this purpose as it was punctuated by abrupt Dansgaard–Oeschger events at millennial timescales. Because their center of action was in the North Atlantic, where their climate impacts were largest, modeling studies have mainly focused on the millennial-scale evolution of Northern Hemisphere (NH) paleo ice sheets. Sea-level reconstructions attribute the origin of millennial-scale sea-level variations mainly to NH paleo ice sheets, with a minor but not negligible role of the AIS. Here we investigate the AIS response to millennial-scale climate variability for the first time. To this end we use a three-dimensional, thermomechanical hybrid, ice sheet–shelf model. Different oceanic sensitivities are tested and the sea-level equivalent (SLE) contributions computed. We find that whereas atmospheric variability has no appreciable effect on the AIS, changes in submarine melting rates can have a strong impact on it. We show that in contrast to the widespread assumption that the AIS is a slow reactive and static ice sheet that responds at orbital timescales only, it can lead to ice discharges of around 6 m SLE, involving substantial grounding line migrations at millennial timescales.

Offline forcing methods for ice-sheet models often make use of an index approach in which temperature anomalies relative to the present are calculated by combining a simulated glacial–interglacial climatic anomaly field, interpolated through an index derived from the Greenland ice-core temperature reconstruction, with present-day climatologies. An important drawback of this approach is that it clearly misrepresents climate variability at millennial timescales. The reason for this is that the spatial glacial–interglacial anomaly field used is associated with orbital climatic variations, while it is scaled following the characteristic time evolution of the index, which includes orbital and millennial-scale climate variability. The spatial patterns of orbital and millennial variability are clearly not the same, as indicated by a wealth of models and data. As a result, this method can be expected to lead to a misrepresentation of climate variability and thus of the past evolution of Northern Hemisphere (NH) ice sheets. Here we illustrate the problems derived from this approach and propose a new offline climate forcing method that attempts to better represent the characteristic pattern of millennial-scale climate variability by including an additional spatial anomaly field associated with this timescale. To this end, three different synthetic transient forcing climatologies are developed for the past 120 kyr following a perturbative approach and are applied to an ice-sheet model. The impact of the climatologies on the paleo-evolution of the NH ice sheets is evaluated. The first method follows the usual index approach in which temperature anomalies relative to the present are calculated by combining a simulated glacial–interglacial climatic anomaly field, interpolated through an index derived from ice-core data, with present-day climatologies. In the second approach the representation of millennial-scale climate variability is improved by incorporating a simulated stadial–interstadial anomaly field. The third is a refinement of the second one in which the amplitudes of both orbital and millennial-scale variations are tuned to provide perfect agreement with a recently published absolute temperature reconstruction over Greenland. The comparison of the three climate forcing methods highlights the tendency of the usual index approach to overestimate the temperature variability over North America and Eurasia at millennial timescales. This leads to a relatively high NH ice-volume variability on these timescales. Through enhanced ablation, this results in too low an ice volume throughout the last glacial period (LGP), below or at the lower end of the uncertainty range of estimations. Improving the representation of millennial-scale variability alone yields an important increase in ice volume in all NH ice sheets but especially in the Fennoscandian Ice Sheet (FIS). Optimizing the amplitude of the temperature anomalies to match the Greenland reconstruction results in a further increase in the simulated ice-sheet volume throughout the LGP. Our new method provides a more realistic representation of orbital and millennial-scale climate variability and improves the transient forcing of ice sheets during the LGP. Interestingly, our new approach underestimates ice-volume variations on millennial timescales as indicated by sea-level records. This suggests that either the origin of the latter is not the NH or that processes not represented in our study, notably variations in oceanic conditions, need to be invoked to explain millennial-scale ice-volume fluctuations. We finally provide here both our derived climate evolution of the LGP using the three methods as well as the resulting ice-sheet configurations. These could be of interest for future studies dealing with the atmospheric or/and oceanic consequences of transient ice-sheet evolution throughout the LGP and as a source of climate input to other ice-sheet models.

Observations suggest that during the last decades the Greenland Ice Sheet (GrIS) has experienced a gradually accelerating mass loss, in part due to the observed speed-up of several of Greenland's marine-terminating glaciers. Recent studies directly attribute this to warming North Atlantic temperatures, which have triggered melting of the outlet glaciers of the GrIS, grounding-line retreat and enhanced ice discharge into the ocean, contributing to an acceleration of sea-level rise. Reconstructions suggest that the influence of the ocean has been of primary importance in the past as well. This was the case not only in interglacial periods, when warmer climates led to a rapid retreat of the GrIS to land above sea level, but also in glacial periods, when the GrIS expanded as far as the continental shelf break and was thus more directly exposed to oceanic changes. However, the GrIS response to palaeo-oceanic variations has yet to be investigated in detail from a mechanistic modelling perspective. In this work, the evolution of the GrIS over the past two glacial cycles is studied using a three-dimensional hybrid ice-sheet–shelf model. We assess the effect of the variation of oceanic temperatures on the GrIS evolution on glacial–interglacial timescales through changes in submarine melting. The results show a very high sensitivity of the GrIS to changing oceanic conditions. Oceanic forcing is found to be a primary driver of GrIS expansion in glacial times and of retreat in interglacial periods. If switched off, palaeo-atmospheric variations alone are not able to yield a reliable glacial configuration of the GrIS. This work therefore suggests that considering the ocean as an active forcing should become standard practice in palaeo-ice-sheet modelling.

Despite large-scale infrastructure development, deforestation, mining and petroleum exploration in the Amazon Basin, relatively little attention has been paid to the management scale required for the protection of wetlands, fisheries and other aspects of aquatic ecosystems. This is due, in part, to the enormous size, multinational composition and interconnected nature of the Amazon River system, as well as to the absence of an adequate spatial model for integrating data across the entire Amazon Basin. In this data article we present a spatially uniform multi-scale GIS framework that was developed especially for the analysis, management and monitoring of various aspects of aquatic systems in the Amazon Basin. The Amazon GIS-Based River Basin Framework is accessible as an ESRI geodatabase at doi:10.5063/F1BG2KX8.

Studies addressing climate variability during the last millennium generally focus on variables with a direct influence on climate variability, like the fast thermal response to varying radiative forcing, or the large-scale changes in atmospheric dynamics (e.g. North Atlantic Oscillation). The ocean responds to these variations by slowly integrating in depth the upper heat flux changes, thus producing a delayed influence on ocean heat content (OHC) that can later impact low frequency SST (sea surface temperature) variability through reemergence processes. In this study, both the externally and internally driven variations of the OHC during the last millennium are investigated using a set of fully coupled simulations with the ECHO-G (coupled climate model ECHAMA4 and ocean model HOPE-G) atmosphere–ocean general circulation model (AOGCM). When compared to observations for the last 55 yr, the model tends to overestimate the global trends and underestimate the decadal OHC variability. Extending the analysis back to the last one thousand years, the main impact of the radiative forcing is an OHC increase at high latitudes, explained to some extent by a reduction in cloud cover and the subsequent increase of short-wave radiation at the surface. This OHC response is dominated by the effect of volcanism in the preindustrial era, and by the fast increase of GHGs during the last 150 yr. Likewise, salient impacts from internal climate variability are observed at regional scales. For instance, upper temperature in the equatorial Pacific is controlled by ENSO (El Niño Southern Oscillation) variability from interannual to multidecadal timescales. Also, both the Pacific Decadal Oscillation (PDO) and the Atlantic Multidecadal Oscillation (AMO) modulate intermittently the interdecadal OHC variability in the North Pacific and Mid Atlantic, respectively. The NAO, through its influence on North Atlantic surface heat fluxes and convection, also plays an important role on the OHC at multiple timescales, leading first to a cooling in the Labrador and Irminger seas, and later on to a North Atlantic warming, associated with a delayed impact on the AMO.

The study of Greenland ice cores revealed two decades ago the abrupt character of glacial millennial-scale climate variability. Several triggering mechanisms have been proposed and confronted against growing proxy-data evidence. Although the implication of North Atlantic deep water (NADW) formation reorganisations in glacial abrupt climate change seems robust nowadays, the final cause of these reorganisations remains unclear. Here, the role of CO2 and Southern Ocean winds is investigated using a coupled model of intermediate complexity in an experimental setup designed such that the climate system resides close to a threshold found in previous studies. An initial abrupt surface air temperature (SAT) increase over the North Atlantic by 4 K in less than a decade, followed by a more gradual warming greater than 10 K on centennial timescales, is simulated in response to increasing atmospheric CO2 levels and/or enhancing southern westerlies. The simulated peak warming shows a similar pattern and amplitude over Greenland as registered in ice core records of Dansgaard-Oeschger (D/O) events. This is accompanied by a strong Atlantic meridional overturning circulation (AMOC) intensification. The AMOC strengthening is found to be caused by a northward shift of NADW formation sites into the Nordic Seas as a result of a northward retreat of the sea-ice front in response to higher temperatures. This leads to enhanced heat loss to the atmosphere as well as reduced freshwater fluxes via reduced sea-ice import into the region. In this way, a new mechanism that is consistent with proxy data is identified by which abrupt climate change can be promoted.

Heinrich events, identified as enhanced ice-rafted detritus (IRD) in North Atlantic deep sea sediments (Heinrich, 1988; Hemming, 2004) have classically been attributed to Laurentide ice-sheet (LIS) instabilities (MacAyeal, 1993; Calov et al., 2002; Hulbe et al., 2004) and assumed to lead to important disruptions of the Atlantic meridional overturning circulation (AMOC) and North Atlantic deep water (NADW) formation. However, recent paleoclimate data have revealed that most of these events probably occurred after the AMOC had already slowed down or/and NADW largely collapsed, within about a thousand years (Hall et al., 2006; Hemming, 2004; Jonkers et al., 2010; Roche et al., 2004), implying that the initial AMOC reduction could not have been caused by the Heinrich events themselves. Here we propose an alternative driving mechanism, specifically for Heinrich event 1 (H1; 18 to 15 ka BP), by which North Atlantic ocean circulation changes are found to have strong impacts on LIS dynamics. By combining simulations with a coupled climate model and a three-dimensional ice sheet model, our study illustrates how reduced NADW and AMOC weakening lead to a subsurface warming in the Nordic and Labrador Seas resulting in rapid melting of the Hudson Strait and Labrador ice shelves. Lack of buttressing by the ice shelves implies a substantial ice-stream acceleration, enhanced ice-discharge and sea level rise, with peak values 500–1500 yr after the initial AMOC reduction. Our scenario modifies the previous paradigm of H1 by solving the paradox of its occurrence during a cold surface period, and highlights the importance of taking into account the effects of oceanic circulation on ice-sheets dynamics in order to elucidate the triggering mechanism of Heinrich events.