Modelling feedbacks between the Northern Hemisphere ice sheets and climate during the last glacial cycle

During the last glacial cycle (LGC), ice sheets covered large parts of Eurasia and North America, which resulted in ∼120 m of sea level change. Ice sheet-climate interactions have considerable influence on temperature and precipitation patterns and therefore need to be included when simulating this time period. Ideally, ice sheet-climate interactions are simulated by a high-resolution Earth system model. While these models are capable of simulating climates at a certain point in time, such as the pre-industrial (PI) or the Last Glacial Maximum (LGM; 21 000 years ago), a full transient glacial cycle is currently computationally unfeasible as it requires a too-large amount of computation time. Nevertheless, ice sheet models require forcing that captures the gradual change in climate over time to calculate the accumulation and melt of ice and its effect on ice sheet extent and volume changes. Here we simulate the LGC using an ice sheet model forced by LGM and PI climates. The gradual change in climate is modelled by transiently interpolating between pre-calculated results from a climate model for the LGM and the PI. To assess the influence of ice sheet-climate interactions, we use two different interpolation methods: the climate matrix method, which includes a temperature-albedo and precipitation-topography feedback, and the glacial index method, which does not. To investigate the sensitivity of the results to the prescribed climate forcing, we use the output of several models that are part of the Paleoclimate Modelling Intercomparison Project Phase III (PMIP3). In these simulations, ice volume is prescribed, and the climate is reconstructed with a general circulation model (GCM). Here we test those models by using their climate to drive an ice sheet model over the LGC. We find that the ice volume differences caused by the climate forcing exceed the differences caused by the interpolation method. Some GCMs produced unrealistic LGM volumes, and only four resulted in reasonable ice sheets, with LGM Northern Hemisphere sea level contribution ranging between 74-113 m with respect to the present day. The glacial index and climate matrix methods result in similar ice volumes at the LGM but yield a different ice evolution with different ice domes during the inception phase of the glacial cycle and different sea level rates during the deglaciation phase. The temperature-albedo feedback is the main cause of differences between the glacial index and climate matrix methods

Benchmarking the vertically integrated ice-sheet model IMAU-ICE (version 2.0)

Ice-dynamical processes constitute a large uncertainty in future projections of sea-level rise caused by anthropogenic climate change. Improving our understanding of these processes requires ice-sheet models that perform well at simulating both past and future ice-sheet evolution. Here, we present version 2.0 of the ice-sheet model IMAU-ICE, which uses the depth-integrated viscosity approximation (DIVA) to solve the stress balance. We evaluate its performance in a range of benchmark experiments, including simple analytical solutions and both schematic and realistic model intercomparison exercises. IMAU-ICE has adopted recent developments in the numerical treatment of englacial stress and sub-shelf melt near the grounding line, which result in good performance in experiments concerning grounding-line migration (MISMIP, MISMIP+) and buttressing (ABUMIP). This makes it a model that is robust, versatile, and user-friendly, which will provide a firm basis for (palaeo-)glaciological research in the coming years

Simulated land ice: latitudinal changes (0-2.6 Ma) and importance of CO2-glaciation divergence during times of decreasing obliquity (0-800 ka)

Following Milankovitch's theory the incoming insolation or summer energy at 65°N is typically analysed to predict the waxing or waning of land ice. We here use a model-based deconvolution of the LR04 benthic-d18O stack into land ice distribution (de Boer et al., 2014, Köhler et al., 2015) to verify if the latitudinal focal point of land ice dynamics has changed over the last 2 Myr or whether this choice of 65°N in orbital data is indeed well justified. We find that the 5°-latitudinal band which contributes most to land ice albedo radiative forcing (ΔR_[LI]) is 70-75°N between 2.0-1.5 Myr, which is then until 1.0 Myr gradually substituted by 65-70°N. During the last 1 Myr both 60-65°N and 65-70°N dominate ΔR_[LI] and contribute approximately the same amount, while the relative importance of 70-75°N is shrinking. Our analyses illustrates that the choice of 65°N seems for the last 1 Myr to be well justified, while for earlier parts of the last 2 Myr the dominant land ice changes seems to happen up to 10° further to the north. Focusing on the last 800 kyr (the time for which precise data on atmospheric CO2 concentration exists) we furthermore find that the multi-millennial land ice growth and proxy-based reconstruction of global cooling (= the glaciation) appear synchronously to each other and to decreasing obliquity, but diverge from CO2. This suggests that the global cooling associated with Earth's way into an ice age as deduced in the reconstructions has to be mainly caused by the land ice albedo feedback, and is not dominated by the CO2 greenhouse forcing. One way of perceiving this CO2-glaciation divergence in reconstructions is that the reduced incoming insolation at high latitudes causes land ice growth and cooling, while there is a coexisting process that keeps CO2 at a relatively constant level. Solid Earth modeling experiments have indicated that falling sea level might lead to enhanced magma and CO2 production at mid-ocean ridges. Hasenclever et al. (2017) suggested that the combination of marine volcanism at mid-ocean ridges and at hot spot island volcanoes might react to decreasing sea level and be a potential cause for this CO2-glaciation divergence. This CO2-glaciation divergence needs to be considered, when using paleo data to quantify paleoclimate sensitivity: periods with diverging CO2 and global temperature change should be filtered out when approximating the relationship between global temperature rise and CO2 concentrations (Köhler et al., 2018). References: de Boer et al. (2014). https://doi.org/10.1038/ncomms3999. Köhler et al. (2015). https://doi.org/10.5194/cp-11-1801-2015. Hasenclever et al. (2017). https://doi.org/10.1038/ncomms15867. Köhler et al. (2018). https://doi.org/10.1029/2018GL077717

A State-Dependent Quantification of Climate Sensitivity Based On Paleodata of the Last 2.1 Million Years

The evidence from both data and models indicates that specific equilibrium climate sensitivity S[X]—the global annual mean surface temperature change (ΔTg) as a response to a change in radiative forcing X (ΔR[X])—is state dependent. Such a state dependency implies that the best fit in the scatterplot of ΔTg versus ΔR[X] is not a linear regression but can be some nonlinear or even nonsmooth function. While for the conventional linear case the slope (gradient) of the regression is correctly interpreted as the specific equilibrium climate sensitivity S[X], the interpretation is not straightforward in the nonlinear case. We here explain how such a state-dependent scatterplot needs to be interpreted and provide a theoretical understanding—or generalization—how to quantify S[X] in the nonlinear case. Finally, from data covering the last 2.1 Myr we show that—due to state dependency—the specific equilibrium climate sensitivity which considers radiative forcing of CO2 and land ice sheet (LI) albedo, math formula, is larger during interglacial states than during glacial conditions by more than a factor 2

The effect of obliquity-driven changes on paleoclimate sensitivity during the late Pleistocene

We reanalyze existing paleodata of global mean surface temperature ΔTg and radiative forcing ΔR of CO2 and land ice albedo for the last 800,000 years to show that a state‐dependency in paleoclimate sensitivity S, as previously suggested, is only found if ΔTg is based on reconstructions, and not when ΔTg is based on model simulations. Furthermore, during times of decreasing obliquity (periods of land‐ice sheet growth and sea level fall) the multi‐millennial component of reconstructed ΔTg diverges from CO2, while in simulations both variables vary more synchronously, suggesting that the differences during these times are due to relatively low rates of simulated land ice growth and associated cooling. To produce a reconstruction‐based extrapolation of S for the future we exclude intervals with strong ΔTg‐CO2 divergence and find that S is less state‐dependent, or even constant (state‐independent), yielding a mean equilibrium warming of 2–4 K for a doubling of CO2

Anti-Phased Miocene Ice Volume and CO2 Changes by Transient Antarctic Ice Sheet Variability

Geological evidence indicates large continental‐scale Antarctic ice volume variations during the early and mid‐Miocene. On million‐year timescales, these variations can largely be explained by equilibrium Antarctic ice sheet (AIS) simulations. In contrast, on shorter orbital timescales, the AIS needs not be in equilibrium with the forcing and ice volume variations may be substantially different. Here, we introduce a conceptual model, based on ice dynamical model results, to investigate the difference between transient variability and equilibrium differences of the Miocene AIS. In our model, an ice sheet will grow (shrink) by a specific rate when it is smaller (larger) than its equilibrium size. We show that phases of concurrent ice volume increase and rising CO2 levels are possible, even though the equilibrium ice volume decreases monotonically with CO2. When the AIS volume is out of equilibrium with the forcing climate, the ice sheet can still be adapting to a relatively large equilibrium size, although CO2 is rising after a phase of decrease. A delayed response of Antarctic ice volume to (covarying) solar insolation and CO2 concentrations can cause discrepancies between Miocene solar insolation and benthic δ18O variability. Increasing forcing frequency leads to a larger disequilibrium and consequently larger CO2‐ice volume phase differences. Furthermore, an amplified forcing amplitude causes larger amplitude ice volume variability, because the growth and decay rates depend on the forcing. It also leads to a reduced average ice volume, resulting from the growth rates generally being smaller than the decay rates

A Highly Sensitive ESIPT-Based Ratiometric Fluorescence Sensor for Selective Detection of Al³⁺

An excited-state intramolecular proton transfer (ESIPT)-based highly sensitive ratiometric fluorescence sensor, <b>1</b>H was developed for selective detection of aluminum (Al³⁺) in acetonitrile as well as in 90% aqueous system. Single-crystal X-ray diffraction analysis reveals almost planar and conjugated structure of <b>1</b>H. Photophysical properties of the sensor as well as its selectivity toward Al³⁺ are explored using UV–visible, steady-state, and time-resolved fluorescence spectroscopic studies. The bright cyan (λ_em = 445 nm) fluorescence of <b>1</b>H in acetonitrile turns into deep blue (λ_em = 412 nm) with ∼2.3-fold enhancement in emission intensity, in the presence of parts per billion level Al³⁺ (detection limit = 0.5 nM). Interestingly, the probe <b>1</b>H exhibits increased selectivity toward Al³⁺ in H₂O/acetonitrile (9:1 v/v) solvent system with a change in fluorescence color from pale green to deep blue associated with ca. sixfold enhancement in emission intensity. Density functional theoretical (DFT) calculations provide the ground- and excited-state energy optimized structures and properties of the proposed aluminum complex [Al­(<b>1</b>) (OH)]₂²⁺, which is in harmony with the solution-state experimental findings and also supports the occurrence of ESIPT process in <b>1</b>H. The ESIPT mechanism was also ascertained by comparing the basic photophysical properties of <b>1</b>H with a similar O-methylated analogue, <b>1</b>′Me