Paleomegnetic and geochemistry of sediments from the early Aptian Oceanic Anoxic Event

The early Aptian Oceanic Anoxic Event (OAE1a, 120 Ma) represents a geologically brief time interval in the mid-Cretaceous greenhouse world that is characterized by increased organic carbon accumulation in marine sediments, sudden biotic changes, and abrupt carbon-isotope excursions indicative of significant perturbations to global carbon cycling. The brevity of these drastic environmental changes (< 10**6 year) and the typically 10**6 year temporal resolution of the available chronologies, however, represent a critical gap in our knowledge of OAE1a. We have conducted a high-resolution investigation of three widely distributed sections, including the Cismon APTICORE in Italy, Santa Rosa Canyon in northeastern Mexico, and Deep Sea Drilling Project (DSDP) Site 398 off the Iberian margin in the North Atlantic Ocean, which represent a range of depositional environments where condensed and moderately expanded OAE1a intervals are recorded. The objectives of this study are to establish orbital chronologies for these sections and to construct a common, high-resolution timescale for OAE1a. Spectral analyses of the closely-spaced (corresponding to ~5 to 10 kyr) measurements of calcium carbonate content of the APTICORE, magnetic susceptibility (MS) and anhysteretic remanent magnetization (ARM) of the Santa Rosa samples, and MS, ARM and ARM/IRM, where IRM is isothermal remanent magnetization, of Site 398 samples reveal statistically significant cycles. These cycles exhibit periodicity ratios and modulation patterns similar to those of the mid-Cretaceous orbital cycles, suggesting that orbital variations may have modulated depositional processes. Orbital control allows us to estimate the duration of unique, globally identifiable stages of OAE1a. Although OAE1a had a duration of ~1.0 to 1.3 Myr, the initial perturbation represented by the negative carbon-isotope excursion was rapid, lasting for ~27-44 kyr. This estimate could serve as a basis for constraining triggering mechanisms for OAE1a.

Chemical, and Sr and Li isotopic composition of pore fluids and sediments from ODP Leg 170 and Leg 205 sites

Pore fluid and sediment chemical and isotopic data were obtained for samples from Ocean Drilling Program (ODP) Leg 205 Sites 1253, 1254, and 1255 in the Costa Rica subduction zone. The chemical and isotopic data reported here were generated in our shore-based laboratories to complement shipboard inorganic geochemical data. Li isotopic analyses were carried out by L.-H. Chan at Louisiana State University (USA). The data reported herein include fluoride, bromide, rubidium, cesium, and barium concentrations; Li and Sr isotopic compositions in pore fluids; and Rb, Cs, and Ba concentrations in representative bulk sediments. The data also include new pore fluid fluoride and bromide concentrations from corresponding ODP Leg 170 Sites 1039, 1040, and 1043. O.M. Saether's Site 1039 and 1040 fluoride concentration data are shown for comparison. Basal sediment fluoride concentrations and Li and Sr isotope ratios at both Sites 1253 and 1039 show reversals that approach modern seawater values. Br/Cl ratios are, however, conservative throughout the sediment section at Sites 1039 and 1253. The observed sharp F and Br concentration maxima, Rb and K concentration minima, the most radiogenic 87Sr/86Sr ratios, and highest 7Li values along the décollement and fracture zone (Sites 1040, 1043, 1254, and 1255) strengthen the evidence obtained during Leg 170 that a deeply sourced fluid, originating from fluid-rock reactions at ~150°C and corresponding to between 10 and 15 km depth, is transporting solutes to the ocean.

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.

Horizontal profiles of longwave and shortwave radiation components over sea ice near Barrow, Alaska during the 2011 melt

An integrated instrument package for measuring and understanding the surface radiation budget of sea ice is presented, along with results from its first deployment. The setup simultaneously measures broadband fluxes of upwelling and downwelling terrestrial and solar radiation (four components separately), spectral fluxes of incident and reflected solar radiation, and supporting data such as air temperature and humidity, surface temperature, and location (GPS), in addition to photographing the sky and observed surface during each measurement. The instruments are mounted on a small sled, allowing measurements of the radiation budget to be made at many locations in the study area to see the effect of small-scale surface processes on the large-scale radiation budget. Such observations have many applications, from calibration and validation of remote sensing products to improving our understanding of surface processes that affect atmosphere-snow-ice interactions and drive feedbacks, ultimately leading to the potential to improve climate modelling of ice-covered regions of the ocean. The photographs, spectral data, and other observations allow for improved analysis of the broadband data. An example of this is shown by using the observations made during a partly cloudy day, which show erratic variations due to passing clouds, and creating a careful estimate of what the radiation budget along the observed line would have been under uniform sky conditions, clear or overcast. Other data from the setup's first deployment, in June 2011 on fast ice near Point Barrow, Alaska, are also shown; these illustrate the rapid changes of the radiation budget during a cold period that led to refreezing and new snow well into the melt season.

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.