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.

Sea-surface reconstruction of the east-equatorial South Atlantic

In order to reconstruct Late Quatemary variations of surface oceanography in the eastequatorial South Atlantic, time series of sea-surface temperatures (SST) and paleoproductivity were established from cores recovered in the Guinea and Angola Basins, and at the Walvis Ridge. These records, based on sedimentary alkenone and organic carbon concentrations, reveal that during the last 350,000 years surface circulation and productivity changes in the east-equatorial South Atlantic were highiy sensitive to climate forcing at 23- and 100-kyr periodicities. Covarying SST and paleoproductivity changes at the equator and at the Walvis Ridge appear to be driven by variations in zonal trade-wind intensity, which forces intensification or reduction of coastal and equatorial upwelling, as well as enhanced Benguela cold water advection from the South. Phase relationships of precessional variations in the paleoproductivity and SST records from the distinct sites were evaluated with respect to boreal summer insolation over Africa, movements of southem ocean thermal fronts, and changes in global ice volume. The 23-kyr phasing implies a sensitivity of eastem South Atlantic surface water advection and upwelling to West African monsoon intensity and to changes in the position ofthe subtropical high pressure cell over the South Atlantic, both phenomena which modulate zonal strength of southeasterly trades. SST and productivity changes north of 20°S lack significant variance at the 41-kyr periodicity; and at the Walvis Ridge and the equator lead changes in ice volume. This may indicate that obliquity-driven clirnate change, characteristic for northem high latitudes, e.g fluctuations in continental ice masses, did not substantially influence subtropical and tropical surface circulation in the South Atlantic. At the 23-kyr cycle SST and productivity changes in the eastern Angola Basin lag those in the equatorial Atlantic and at the Walvis Ridge by about 3500 years. This lag is explained by variations in cross-equatorial surface water transport and west-east countercurrent retum flow modifying precessional variations of SST and productivity in the eastem Angola Basin relative to those in the mid South Atlantic area under the central field of zonal trade winds. Sea level-related shifts of upwelling cells in phase with global clirnate change may be also recorded in SST and productivity variability along the continental margin off Southwest Africa. They may account for the delay of the paleoceanogreaphic signal from continental margin sites with respect to that from the pelagic sites at the equator and the Walvis Ridge.

Heterocyst glycolipids in surface sediments and water samples of the Amazon shelf

Marine endosymbiotic heterocystous cyanobacteria make unique heterocyst glycolipids (HGs) containing pentose (C5) moieties. Functionally similar HGs with hexose (C6) moieties found in free-living cyanobacteria occur in the sedimentary record, but C5 HGs have not been documented in the natural environment. Here we developed a high performance liquid chromatography multiple reaction monitoring (MRM) mass spectrometry (HPLC-MS2) method specific for trace analysis of long chain C5HGs and applied it to cultures of Rhizosolenia clevei Ostenfeld and its symbiont Richelia intracellularis which were found to contain C5 HGs and no C6 HGs. The method was then applied to suspended particulate matter (SPM) and surface sediment from the Amazon plume region known to harbor marine diatoms carrying heterocystous cyanobacteria as endosymbionts. C5 HGs were detected in both marine SPM and surface sediments, but not in SPM or surface sediment from freshwater settings in the Amazon basin. Rather, the latter contained C6 HGs, established biomarkers for free-living heterocystous cyanobacteria. Our results indicate that the C5 HGs may be potential biomarkers for marine endosymbiotic heterocystous cyanobacteria.

Sea surface temperatur reconstruction for ODP Site 108-662

The record of planktonic foraminifer abundances at Site 662 during the late Pliocene (~1.7-2.1 Ma) was examined to determine variations in estimated sea-surface temperature (SST). We compared the results to SST estimates from a late Pleistocene record (~1.5-200 ka) from nearby piston core RC24-7. Within the primary orbital band (~20-100 k.y.), the cold-season responses of both equatorial Atlantic records are dominated by the precessional period, and the computed range of variability is quite similar. This is in contrast to the evolution of the dominant climatic response from 41 to 100 k.y. at high northern latitudes between the late Pliocene and the late Pleistocene. The orbital-band SST response in this region of greatest divergence in the equatorial Atlantic has not changed appreciably between the late Pliocene and the late Pleistocene.

Ocean temperature response to idealized Gleissberg and de Vries solar cycles in a comprehensive climate model

The ~90-year Gleissberg and ~200-year de Vries cycles have been identified as two distinctive quasi-periodic components of Holocene solar activity. Evidence exists for the impact of such multi-decadal to centennial-scale variability in total solar irradiance (TSI) on climate, but concerning the ocean, this evidence is mainly restricted to the surface response. Here we use a comprehensive global climate model to study the impact of idealized solar forcing, representing the Gleissberg and de Vries cycles, on global ocean potential temperature at different depth levels, after a recent proxy record indicates a signal of TSI anomalies in the northeastern Atlantic at mid-depth. Potential impacts of TSI anomalies on deeper oceanic levels are climatically relevant due to their possible effect on ocean circulation by altering water mass characteristics. Simulated solar anomalies are shown to penetrate the ocean down to at least deep-water levels. Despite the fact that the two forcing periods differ only by a factor of ~2, the spatial pattern of response is significantly distinctive between the experiments, suggesting different mechanisms for solar signal propagation. These are related to advection by North Atlantic Deep Water flow (200-year forcing), and barotropic adjustment in the South Atlantic in response to a latitudinal shift of the westerly wind belt (90-year forcing).

Soil, snow, weather, and sub-surface storage data from a mountain catchment in the rain-snow transition zone

A comprehensive hydroclimatic data set is presented for the 2011 water year to improve understanding of hydrologic processes in the rain-snow transition zone. This type of dataset is extremely rare in scientific literature because of the quality and quantity of soil depth, soil texture, soil moisture, and soil temperature data. Standard meteorological and snow cover data for the entire 2011 water year are included, which include several rain-on-snow events. Surface soil textures and soil depths from 57 points are presented as well as soil texture profiles from 14 points. Meteorological data include continuous hourly shielded, unshielded, and wind corrected precipitation, wind speed, air temperature, relative humidity, dew point temperature, and incoming solar and thermal radiation data. Sub-surface data included are hourly soil moisture data from multiple depths from 7 soil profiles within the catchment, and soil temperatures from multiple depths from 2 soil profiles. Hydrologic response data include hourly stream discharge from the catchment outlet weir, continuous snow depths from one location, intermittent snow depths from 5 locations, and snow depth and density data from ten weekly snow surveys. Though it represents only a single water year, the presentation of both above and below ground hydrologic condition makes it one of the most detailed and complete hydro-climatic datasets from the climatically sensitive rain-snow transition zone for a wide range of modeling and descriptive studies.

Response of deep-sea benthic foraminifera to Late Quaternary climate changes, southeast Indian Ocean, offshore Western Australia

The Late Quaternary benthic foraminifera of four deep-sea cores off Western Australia (ODP 122-760A, ODP 122-762B, BMR96GC21 and RC9-150) have been examined for evidence of increased surface productivity to explain the anomalously low sea-surface paleotemperatures inferred by planktic foraminifera for the last and penultimate glaciations. The delta13C trends of Cibicidoides wuellerstorfi, and differences between the delta13C trends of planktics (Globigerinoides sacculifer) and benthics (C. wuellerstorfi) in the four cores indicate that during stage 6 bottom waters were significantly depleted in delta13C, and strong delta13C gradients were established in the water column, while during stage 2 and the Last Glacial Maximum, delta13C trends did not differ greatly from that of the Holocene. Two main assemblages of benthic foraminifera were identified by principal component analyses: one dominated by Uvigerina peregrina, another dominated by U. proboscidea. Abundance of these Uvigerinids, and of taxa preferring an infaunal microhabitat, and of Epistominella exigua and Bulimina aculeata indicate that episodes of high influx of particulate organic matter were established in most sites during glacial episodes, and particularly so during stage 6, while evidence for upwelling during the Last Glacial Maximum is less strong. The Penultimate Glaciation upwellings were established within the areas of low sea-surface paleotemperature indicated by planktic foraminifera. During the Last Interglacial Climax, upwelling appears to have been established in an isolated region offshore from a strengthened Leeuwin Current off North West Cape. Last Glacial Maximum delta13C values of C. wuellerstorfi at waterdepths of less than 2000 m show smaller than global mean glacial-interglacial changes suggesting the development of a deep hydrological front. A similar vertical stratification/bathyal front was also established during the Penultimate Glaciation.

Pore-water chemistry of abyssal Arabian Sea surface sediments

Benthic fluxes and pore-water compositions of silicic acid, nitrate and phosphate were investigated for surface sediments of the abyssal Arabian Sea during four cruises (1995-1998). Five sites located in the northern (NAST), western (WAST), central (CAST), eastern (EAST), and southern (SAST) Arabian Sea were revisited during intermonsoonal periods after the NE- and SW-Monsoon. At these sites, benthic fluxes of remineralized nutrients from the sediment to the bottom water of 36-106, 102-350 and 4-16 mmol/m**2/yr were measured for nitrate, silicic acid and phosphate, respectively. The benthic fluxes and pore-water compositions showed a distinct regional pattern. Highest fluxes were observed in the western and northern region of the Arabian Sea, whereas decreasing fluxes were derived towards the southeast. At WAST, the general temporal pattern of primary production, related to the NE- and SW-Monsoon, is reflected by benthic fluxes. In contrast, at sites NAST, SAST, CAST, and EAST a temporal pattern of fluxes in response to the monsoon is not obvious. Our results reveal a clear coupling between the general regional pattern of production in surface waters and the response of the benthic environment, as indicated by the flux of remineralized nutrients, though a spatially differing degree of decoupling during transport and remineralization of particulate organic matter and biogenic opal was observed. This has to be taken into account regarding budget calculations and paleoceanographic topics.

At surface behaviour data of juvenile southern elephant seals from Marion Island between 2001 and 2006 with links to datasets

Marine mammals forage in dynamic environments characterized by variables that are continuously changing in relation to large-scale oceanographic processes. In the present study, behavioural states of satellite-tagged juvenile southern elephant seals (n = 16) from Marion Island were assessed for each reliable location, using variation in turning angle and speed in a state-space modelling framework. A mixed modelling approach was used to analyse the behavioural response of juvenile southern elephant seals to sea-surface temperature and proximity to frontal and bathymetric features. The findings emphasised the importance of frontal features as potentially rewarding areas for foraging juvenile southern elephant seals and provided further evidence of the importance of the area west of Marion Island for higher trophic-level predators. The importance of bathymetric features during the transit phase of juvenile southern elephant seal migrations indicates the use of these features as possible navigational cues.

Response of marine palynomorphs to Neogene climate cooling in the Iceland Sea, ODP Hole 151-907A

The present study on ODP Leg 151 Hole 907A combines a detailed analysis of marine palynomorphs (dinoflagellate cysts, prasinophytes, and acritarchs) and a low-resolution alkenone-based sea-surface temperature (SST) record for the interval between 14.5 and 2.5 Ma, and allows to investigate the relationship between palynomorph assemblages and the paleoenvironmental evolution of the Iceland Sea. A high marine productivity is indicated in the Middle Miocene, and palynomorphs and SSTs both mirror the subsequent long-term Neogene climate deterioration. The diverse Middle Miocene palynomorph assemblages clearly diminish towards the impoverished assemblages of the Late Pliocene; parallel with a somewhat gradual decrease of SSTs being as high as 20 °C at ~13.5 Ma to around 8 °C at ~3 Ma. Superimposed, palynomorph assemblages not only reflect Middle to Late Miocene climate variability partly coinciding with the short-lived global Miocene isotope events (Mi-events), but also the initiation of a proto-thermohaline circulation across the Middle Miocene Climate Transition, which led to increased meridionality in the Nordic Seas. Last occurrences of species cluster during three events in the Late Miocene to Early Pliocene and are ascribed to the progressive strengthening and freshening of the proto-East Greenland Current towards modern conditions. A significant high latitude cooling between 6.5 and 6 Ma is depicted by the supraregional "Decahedrella event" coeval with lowest Miocene productivity and a SST decline. In the Early Pliocene, a transient warming is accompanied by surface water stratification and increased productivity that likely reflects a high latitude response to the global biogenic bloom. The succeeding crash in palynomorph accumulation, and a subsequent interval virtually barren of marine palynomorphs may be attributed to enhanced bottom water oxygenation and substantial sea ice cover, and indicates that conditions seriously affecting marine productivity in the Iceland Sea were already established well before the marked expansion of the Greenland Ice Sheet at 3.3 Ma.

Consistent increase in dimethyl sulfide (DMS) in response to high CO2 in five shipboard bioassays from contrasting NW European waters

The ubiquitous marine trace gas dimethyl sulfide (DMS) comprises the greatest natural source of sulfur to the atmosphere and is a key player in atmospheric chemistry and climate. We explore the short-term response of DMS production and cycling and that of its algal precursor dimethyl sulfoniopropionate (DMSP) to elevated carbon dioxide (CO2) and ocean acidification (OA) in five 96 h shipboard bioassay experiments. Experiments were performed in June and July 2011, using water collected from contrasting sites in NW European waters (Outer Hebrides, Irish Sea, Bay of Biscay, North Sea). Concentrations of DMS and DMSP, alongside rates of DMSP synthesis and DMS production and consumption, were determined during all experiments for ambient CO2 and three high-CO2 treatments (550, 750, 1000 µatm). In general, the response to OA throughout this region showed little variation, despite encompassing a range of biological and biogeochemical conditions. We observed consistent and marked increases in DMS concentrations relative to ambient controls (110% (28-223%) at 550 µatm, 153% (56-295%) at 750 µatm and 225% (79-413%) at 1000 µatm), and decreases in DMSP concentrations (28% (18-40%) at 550 µatm, 44% (18-64%) at 750 µatm and 52% (24-72%) at 1000 µatm). Significant decreases in DMSP synthesis rate constants (µDMSP /d) and DMSP production rates (nmol/d) were observed in two experiments (7-90% decrease), whilst the response under high CO2 from the remaining experiments was generally indistinguishable from ambient controls. Rates of bacterial DMS gross consumption and production gave weak and inconsistent responses to high CO2. The variables and rates we report increase our understanding of the processes behind the response to OA. This could provide the opportunity to improve upon mesocosm-derived empirical modelling relationships and to move towards a mechanistic approach for predicting future DMS concentrations.

Diatom studies in surface and downcore sediments off southern Sumatra and Java

A wealth of sedimentary records aimed at reconstructing late Quaternary changes in productivity and temperature have been devoted to understanding linkages between the Indo-Pacific Warm Pool (IPWP) and other distant oceanic areas. Most of these reconstructions are based, however, on biogeochemical and sedimentological proxies, with comparatively less attention devoted to microfossils. A high-resolution (<1 ka) study of diatom concentrations and the community at site GeoB10038-4, recovered off southern Sumatra (ca. 6°S, 103°E), closely tracks the variations of diatom concentrations in the westernmost IPWP during the last glacial-interglacial cycle. The diatom record provides evidence that diatom paleoproductivity was highest during interglacials, primarily due to the input of lithogenics and nutrients following the rise in sea level after full glacials. In addition, the co-variation of total diatom concentration and Northern Hemisphere forcing for Marine Isotope Stage 5 suggests a direct response of diatom productivity and upwelling intensity to boreal summer insolation. Temporal shifts of the diverse diatom community at site GeoB10038-4 correspond well with the present-day seasonal monsoon pattern and the strengthening and weakening phases of upwelling along the southern coast of Sumatra. Resting spores of Chaetoceros, typical of nutrient-rich waters, were dominant during periods of highest diatom paleoproductivity and responded to the strengthening of the SE monsoon, while diatoms of oligotrophic to mesotrophic waters characterized intermonsoon periods. The close correspondence between the dominance of upwelling diatoms and the boreal summer insolation resembles the present-day dynamics of diatom production. The observed interglacial highs and glacial lows of diatom productivity at site GeoB10038-4 is a unique pattern in the late Quaternary tropics.

Stable isotope ratios of foraminifera and sea surface temperature estimation of sediment core W2804A-14

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C (epsilon p) in a central equatorial Pacific sediment core that spans the last ~255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic compostion of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon p, derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 (c e) based on a new, site-specific calibration of the relationship between epsilon p and c e. The calibration was based on reassessment of existing epsilon p versus c e data, which support a physiologically based model in which epsilon p is inversely related to c e. Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index UK 37. Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon p and 1/c e. These are discussed in detail and it is concluded that the observed record of epsilon p most probably reflects significant variations in Delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from ~110 µatm during glacial intervals (ocean > atmosphere) to ~60 µatm during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p < 0.01) inverse correlations of Delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.