(Table S1) Boron isotope and trace element data and reconstructed carbonate system parameters, ODP Holes 122-761B and 154-926A

The middle Miocene Climatic Optimum (17-15 Ma; MCO) is a period of global warmth and relatively high CO2 and is thought to be associated with a significant retreat of the Antarctic Ice Sheet (AIS). We present here a new planktic foraminiferal d11B record from 16.6 to 11.8 Ma from two deep ocean sites currently in equilibrium with the atmosphere with respect to CO2. These new data demonstrate that the evolution of global climate during the middle Miocene (as reflected by changes in the cyrosphere) was well correlated to variations in the concentration of atmospheric CO2. What is more, within our sampling resolution (~1 sample per 300 kyr) there is no evidence of hysteresis in the response of ice volume to CO2 forcing during the middle Miocene, contrary to what is understood about the Antarctic Ice Sheet from ice sheet modelling studies. In agreement with previous data, we show that absolute levels of CO2 during the MCO were relatively modest (350-400 ppm) and levels either side of the MCO are similar or lower than the pre-industrial (200-260 ppm). These new data imply the presence of either a very dynamic AIS at relatively low CO2 during the middle Miocene or the advance and retreat of significant northern hemisphere ice. Recent drilling on the Antarctic margin and shore based studies indicate significant retreat and advance beyond the modern limits of the AIS did occur during the middle Miocene, but the complete loss of the AIS was unlikely. Consequently, it seems that ice volume and climate variations during the middle Miocene probably involved a more dynamic AIS than the modern but also some component of land-based ice in the northern hemisphere.

Strontium and boron geochemistry of ODP sites at Hydrate Ridge, eastern Pacific Ocean

ODP Leg 204, which drilled at Hydrate Ridge, provides unique insights into the fluid regime of an accretionary complex and delineates specific sub-seafloor pathways for fluid transport. Compaction and dewatering due to smectite-illite transition increase with distance from the toe of the accretionary prism and bring up fluids from deep within the accretionary complex to sampled depths (<= 600 mbsf). These fluids have a distinctly non-radiogenic strontium isotope signature indicating reaction with the oceanic basement. Boron isotopes are also consistent with a deep fluid source that has been modified by desorption of heavy boron as clay minerals change from smectite to illite. One of three major horizons serves as conduit for the transport of mainly fluid. Our results enable us to evaluate fluid migration pathways that play important roles on massive gas hydrate accumulations and seepage of methane-rich fluids on southern Hydrate Ridge.

Chemical composition and solubility of ocean phosphatized rocks in seawater

In order to elucidate possibility of limestone phosphatization during contact with seawater two sets of experiments were carried out: that of solubility of natural phosphates in seawater and that of deposition of phosphorus onto calcareous phase. Concentration of phosphorus in seawater used for solubility experiments varied from 39 to 338 µ/l, temperature was 5.5-17.5°C and pH 7.8-7.9. The lowest solubility was characteristic of ancient crystallized samples. The deposition experiments showed that in case when concentration of dissolved phosphate in seawater reaches 3 mg/l it might be sorbed on solid CaCO3 phase without forming its own mineral. The latter is able to form rapidly but only if magnesium is not present in solution. In magnesium-free water calcium phosphates precipitate when concentration of dissolved phosphorus is higher than 0.9-1.2 mg/l. These results show that phosphatization of limestones in marine environment takes place during their contact with pore water but not with marine bottom water.

Chemistry and boron isotope ratios of Izu arc front volcanic rocks

The process of fluid release from the subducting slab beneath the Izu arc volcanic front (Izu VF) was examined by measuring B concentrations and B isotope ratios in the Neogene fallout tephra (ODP Site 782A). Both were measured by secondary ion mass spectrometry, in a subset of matrix glasses and glassy plagioclase-hosted melt inclusions selected from material previously analyzed for major and trace elements (glasses) and radiogenic isotopes (Sr, Nd, Pb; bulk tephra). These tephra glasses have high B abundances (~10-60 ppm) and heavy delta11B values (+4.5? to +12.0?), extending the previously reported range for Izu VF rocks (delta11B, +7.0? to +7.3?). The glasses show striking negative correlations of delta11B with large ion lithophile element (LILE)/Nb ratios. These correlations cannot be explained by mixing two separate slab fluids, originating from the subducting sediment and the subducting basaltic crust, respectively (model A). Two alternative models (models B and C) are proposed. Model B proposes that the inverse correlations are inherited from altered oceanic crust (AOC), which shows a systematic decrease of B and LILE with increasing depth (from basaltic layer 2A to layer 3), paralleled by an increase in delta11B (from ~ +1? to +10? to +24?). In this model, the contribution of sedimentary B is insignificant (<4% of B in the Izu VF rocks). Model C explains the correlation as a mixture of a low-delta11B (~ +1?) 'composite' slab fluid (a mixture of metasediment- and metabasalt-derived fluids) with a metasomatized mantle wedge containing elevated B (~1-2 ppm) and heavy delta11B (~ +14?). The mantle wedge was likely metasomatized by 11B-rich fluids beneath the outer forearc, and subsequently down dragged to arc front depths by the descending slab. Pb-B isotope systematics indicate that, at arc front depths, ~ 53% of the B in the Izu VF is derived from the wedge. This implies that the heavy delta11B values of Izu VF rocks are largely a result of fluid fractionation, and do not reflect variations in slab source provenance (i.e. subducting sediment vs. basaltic crust). Since the B content of the peridotite at the outer forearc (7-58 ppm B, mean 24 +/- 16 ppm) is much higher than beneath the arc front (~1-2 ppm B), the hydrated mantle wedge must have released a B-rich fluid on its downward path. This 'wedge flux' can explain (1) the across-arc decrease in B and delta11B (e.g. Izu, Kuriles), without requiring a progressive decrease in fluid flux from the subducting slab, and (2) the thermal structure of volcanic arcs, as reflected in the B and delta11B variations of volcanic arc rocks.

Summary of boron isotope and Mg/Ca data

Biological productivity and carbon export in the equatorial Atlantic are thought to have been dramatically higher during the last glacial period than during the Holocene. Here we reconstruct the pH and CO2 content of surface waters from the eastern equatorial Atlantic Ocean over the past ~30 k.y. using the boron isotope composition of Globigerinoides ruber (a mixed-layer-dwelling planktic foraminifera). Our new record, combined with previously published data, indicates that during the last glacial, in contrast to today, a strong west to east gradient existed in the extent of air:sea equilibrium with respect to pCO2 (DeltapCO2), with the eastern equatorial Atlantic acting as a significant source of CO2 (+100 µatm) while the western Atlantic remained close to equilibrium (+25 µatm). This pattern suggests that a fivefold increase in the upwelling rate of deeper waters drove increased Atlantic productivity and large-scale regional cooling during the last glacial, but the higher than modern DeltapCO2 in the east indicates that export production did not keep up with enhanced upwelling of nutrients. However, the downstream decline of DeltapCO2 provides evidence that the unused nutrients from the east were eventually used for biologic carbon export, thereby effectively negating the impact of changes in upwelling on atmospheric CO2 levels. Our findings indicate that the equatorial Atlantic exerted a minimal role in contributing to lower glacial-age atmospheric CO2.