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.

Characterization of late Campanian and Maastrichtian planktonic foraminiferal depth habitats and vital activities based on stable isotopes

Depth habitats of 56 late Cretaceous planktonic foraminiferal species from cool and warm climate modes were determined based on stable isotope analyses of deep-sea samples from the equatorial Pacific DSDP Sites 577A and 463, and South Atlantic DSDP Site 525A. The following conclusions can be reached: Planoglobulina multicamerata (De Klasz) and Heterohelix rajagopalani (Govindan) occupied the deepest plankton habitats, followed by Abathomphalus mayaroensis (Bolli), Globotruncanella havanensis (Voorwijk), Gublerina cuvillieri Kikoine, and Laeviheterohelix glabrans (Cushman) also at subthermocline depth. Most keeled globotruncanids, and possibly Globigerinelliodes and Racemiguembelina species, lived at or within the thermocline layer. Heterohelix globulosa (Ehrenberg) and Rugoglobigerina, Pseudotextularia and Planoglobulina occupied the subsurface depth of the mixed layer, and Pseudoguembelina species inhabited the surface mixed layer. However, depth ranking of some species varied depending on warm or cool climate modes, and late Campanian or Maastrichtian age. For example, most keeled globotruncanids occupied similar shallow subsurface habitats as Rugoglobigerina during the warm late Campanian, but occupied the deeper thermocline layer during cool climatic intervals. Two distinct types of "vital effect" mechanisms reflecting photosymbiosis and respiration effects can be recognized by the exceptional delta13C signals of some species. (1) Photosymbiosis is implied by the repetitive pattern of relatively enriched delta13C values of Racemiguembelina (strongest), Planoglobulina, Rosita and Rugoglobigerina species, Pseudoguembelina excolata (weakest). (2) Enriched respiration 12C products are recognized in A. mayaroensis, Gublerina acuta De Klasz, and Heterohelix planata (Cushman). Isotopic trends between samples suggest that photosymbiotic activities varied between localities or during different climate modes, and may have ceased under certain environmental conditions. The appearance of most photosymbiotic species in the late Maastrichtian suggests oligotrophic conditions associated with increased water-mass stratification.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.