Specific conductivity and hydrogen ions measured on two ice core (NEEM and NGRIP)

A stratigraphy-based chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core has been derived by transferring the annual layer counted Greenland Ice Core Chronology 2005 (GICC05) and its model extension (GICC05modelext) from the NGRIP core to the NEEM core using 787 match points of mainly volcanic origin identified in the electrical conductivity measurement (ECM) and dielectrical profiling (DEP) records. Tephra horizons found in both the NEEM and NGRIP ice cores are used to test the matching based on ECM and DEP and provide five additional horizons used for the timescale transfer. A thinning function reflecting the accumulated strain along the core has been determined using a Dansgaard-Johnsen flow model and an isotope-dependent accumulation rate parameterization. Flow parameters are determined from Monte Carlo analysis constrained by the observed depth-age horizons. In order to construct a chronology for the gas phase, the ice age-gas age difference (Delta age) has been reconstructed using a coupled firn densification-heat diffusion model. Temperature and accumulation inputs to the Delta age model, initially derived from the water isotope proxies, have been adjusted to optimize the fit to timing constraints from d15N of nitrogen and high-resolution methane data during the abrupt onset of Greenland interstadials. The ice and gas chronologies and the corresponding thinning function represent the first chronology for the NEEM core, named GICC05modelext-NEEM-1. Based on both the flow and firn modelling results, the accumulation history for the NEEM site has been reconstructed. Together, the timescale and accumulation reconstruction provide the necessary basis for further analysis of the records from NEEM.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Trace element concentrations and Li isotope composition in DSDP and ODP sediments from island arcs

We have studied the sedimentary and basaltic inputs of lithium to subduction zones. Various sediments from DSDP and ODP drill cores in front of the Mariana, South Sandwich, Banda, East Sunda and Lesser Antilles island arcs have been analysed and show highly variable Li contents and d7Li values. The sediment piles in front of the Mariana and South Sandwich arcs largely consist of pelagic sediments (clays and oozes). The pelagic clays have high Li contents (up to 57.3 ppm) and Li isotope compositions ranging from +1.3? to +4.1?. The oozes have lower Li contents (7.3-16 ppm) with d7Li values of the diatom oozes from the South Sandwich lower (+2.8? to +3.2?) than those of the radiolarian oozes from the Mariana arc (+8.1? to +14.5?). Mariana sediment also contains a significant portion of volcanogenic material, which is characterised by a moderate Li content (14 ppm) and a relatively heavy isotope composition (+6.4?). Sediments from the Banda and Lesser Antilles contain considerable amounts of continental detritus, and have high Li contents (up to 74.3 ppm) and low d7Li values (around 0?), caused by weathering of continental bedrock. East Sunda sediments largely consist of calcareous oozes. These carbonate sediments display intermediate to high Li contents (2.4-41.9 ppm) and highly variable d7Li values (-1.6? to +12.8?). Basaltic oceanic crust samples from worldwide DSDP and ODP drill cores are characterised by enrichment of Li compared to fresh MORB (6.6-33.1 vs. 3.6-7.5 ppm, respectively), and show a large range in Li isotope compositions (+1.7? to +11.8?). The elemental and isotopic enrichment of Li in altered basalts is due to the uptake of isotopically heavy seawater Li during weathering. However, old oceanic crust samples from Sites 417/418 exhibit lighter Li isotope compositions compared to young basaltic crust samples from Sites 332B and 504B. This lighter Li isotope signature in old crust is unexpected and further research is needed to explore this issue.

(Table S3) Benthic foraminiferal Mg/Ca, Li/Ca and stable isotope ratios of ODP Site 122-761

The sensitivities of benthic foraminiferal Mg/Ca and Li/Ca to bottom water temperature and carbonate saturation state have recently been assessed. Here we present a new approach that uses paired Mg/Ca and Li/Ca records to calculate simultaneous changes in temperature and saturation state. Using previously published records, we first use this approach to document a cooling of deep ocean waters associated with the establishment of the Antarctic ice sheet at the Eocene-Oligocene climate transition. We then apply this approach to new records of the Middle Miocene Climate Transition from ODP Site 761 to estimate variations in bottom water temperature and the oxygen isotopic composition of seawater. We estimate that the oxygen isotopic composition of seawater varied by ~1 per mil between the deglacial extreme of the Miocene Climatic Optimum and the glacial maximum following the Middle Miocene Climate Transition, indicating large amplitude variations in ice volume. However, the longer-term change between 15.3 and 12.5 Ma is marked by a ~1°C cooling of deep waters, and an increase in the oxygen isotopic composition of seawater of ~0.6 per mil. We find that bottom water saturation state increased in the lead up to the Middle Miocene Climate Transition and decreased shortly after. This supports decreasing pCO2 as a driver for global cooling and ice sheet expansion, in agreement with existing boron isotope and leaf stomatal index CO2 records but in contrast to the published alkenone CO2 records.

Pore water chemistry of sediment core GeoB15101-7 and water chemistry of CTD casts GeoB15101-1 and GeoB15101-6

Submarine brine lakes feature sharp and persistent concentration gradients between seawater and brine, though these should be smoothed out by free diffusion in open ocean settings. The anoxic Urania basin of the Eastern Mediterranean contains an ultra sulfidic, hypersaline brine of Messinian origin above a thick layer of suspended sediments. With a dual modeling approach we reconstruct its contemporary stratification by geochemical solute transport fundamentals, and show that thermal convection is required to maintain mixing in the brine and mud layer. The origin of the Urania basin stratification was dated to 1650 years before present, which may be linked to a major earthquake in the region. The persistence of the chemoclines may be key to the development of diverse and specialized microbial communities. Ongoing thermal convection in the fluid mud layer may have important, yet unresolved consequences for sedimentological and geochemical processes, also in similar environments.

Ion concentrations and extracelluar pH of Calanus glacialis in Billefjorden, Svalbard

Arctic shelf zooplankton communities are dominated by the copepod Calanus glacialis. This species feeds in surface waters during spring and summer and accumulates large amounts of lipids. Autumn and winter are spent in dormancy in deeper waters. Lipids are believed to play a major role in regulating buoyancy, however, they cannot explain fine-tuning of the depth distribution. To investigate whether ion exchange processes and acid-base regulation support ontogenetic migration as suggested for Antarctic copepods, we sampled C. glacialis in monthly intervals for 1 yr in a high-Arctic fjord and determined cation concentrations and the extracellular pH (pHe) in its hemolymph. During the winter/spring transition, prior to the upward migration of the copepods, Li+ ions were exchanged with cations (Na+, Mg2+, and Ca2+) leading to Li+ concentrations of 197 mmol/L. This likely decreased the density and promoted upward migration in C. glacialis. Our data thus suggest that Li+ has a biological function in this species. Ion and pHe regulation in the hemolymph were not directly correlated, but the pHe revealed a seasonal pattern and was low (5.5) in winter and high (7.9) in summer. Low pHe during overwintering might be related to metabolic depression and thus, support diapause.