(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.

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.

(Table 1) Boron concentration and d11B ratios of ODP Hole 130-803D

The pH of the surface ocean is a sensitive function of its alkalinity and total inorganic carbon concentration, properties which also control the partial pressure of atmospheric carbon dioxide (Broecker and Peng, 1982). Thus, an accurate proxy for past ocean pH could yield information about variations in atmospheric CO2. Recently, it has been suggested that the boron isotopic composition of foraminiferal tests depends on the pH of sea water as well as its isotopic composition (Vengosh et al., 1991, doi:10.1016/0016-7037(91)90139-V; Hemming and Hanson, 1992, doi:10.1016/0016-7037(92)90151-8). Here we present boron isotope and elemental data for sedimentary pore fluids and isotope data for bulk foraminiferal samples from a deep-sea sediment core. The composition of the pore waters implies that sea water boron concentrations and isotopic composition have been constant during the past 21 Myr, allowing us to reconstruct past ocean pH directly from the foraminiferal isotope data. We find that 21 Myr ago, surface ocean pH was only 7.4 ±0.2, but it then increased to 8.2 ±0.2 (roughly the present value) about 7.5 Myr ago. This is consistent with suggestions (Popp et al., 1989; Cerling, 1991; Arthur et al., 1991) that atmospheric CO2 concentrations may have been much higher 21 Myr ago than today.

(Table 3) Boron and chlorine in atmospheric moisture and rain water collected on the shore of the Golubaya Bay in 1970

Boron and chlorine were determined in rain water and in atmospheric moisture condensed in a "Saratov" refrigerator. Ocean is the main source of boron on the earth surface. Boron evaporates from the ocean and enriches atmospheric precipitation: B/Cl ratio of ocean water (0.00024) increases by factor of 10-15. Assuming that the average Cl content in global river runoff is 7.8 mg/l and boron content 0.013 mgl, B/Cl ratio in this runoff is 0.0017. The average B/Cl ratio in rain water of the Golubaya (Blue) Bay (Gelendzhik, Black Sea region) is 0.0026 and in condensates of atmospheric moisture during onshore and offshore winds in the same region it averages from 0.0029 to 0.0033. The maximum boron content in the condensates of this region during onshore winds was 0.032 mg/l and the minimum during offshore winds, 0.004 mg/l. /Cl ratio in sea water over the Atlantic Ocean and in the Gelendzhik area of the Black Sea varied within narrow range, mostly from 0.0025 to 0.0035. Similar B/Cl ratio (0.0024) was found for atmospheric precipitation on the slope of the Terskei Ala-Tau near the Issyk-Kul Lake in 1969. Thus, although chemistries of boron and chlorine (in chlorides) are very different, the B/Cl ratio in the atmosphere is fairly constant. This can be taken as a confirmation of an assumption that salt composition of sea water passes into the atmosphere in molecularly dispersed state. Supposing that the ocean-atmosphere system is in equilibrium as regards to the boron budget, it can be assumed that the same amount of boron passes from the ocean into bottom sediments and from lithosphere rocks and soils into the hydrosphere.

(Table 1, 2) Foraminifera biometric data, and boron isotopic composition of Globigerinoides sacculifer of sediments from the Ontong-Java Plateau

Sediment samples from the Ontong-Java Plateau in the Pacific and the 90° east ridge in the Indian Ocean were used to investigate whether shell size and early diagenesis affect d11B of the symbiont-bearing planktonic foraminifer Globigerinoides sacculifer. In pristine shells from both study locations we found a systematic increase of d11B and Mg/Ca with shell size. Shells in the sieve size class 515-865 µm revealed d11B values +2.1 to +2.3 per mil higher than shells in the 250-380 µm class. This pattern is most likely due to differences in symbiont photosynthetic activity and its integrated effect on the pH of the foraminiferal microenvironment. We therefore suggest smaller individuals must live at approximately 50-100 m water depth where ambient light levels are lower. Using the empirical calibration curve for d11B in G. sacculifer, only shells larger than 425 µm reflect surface seawater pH. Partial dissolution of shells derived from deeper sediment cores was determined by shell weight analyses and investigation of the shell surface microstructure by scanning electron microscopy. The d11B in partially dissolved shells is up to 2 per mil lower relative to pristine shells of the same size class. In agreement with a relatively higher weight loss in smaller shells, samples from the Ontong-Java Plateau show a more pronounced dissolution effect than larger shells. On the basis of the primary size effect and potential postdepositional dissolution effects, we recommend the use of shells that are visually pristine and, in the case of G. sacculifer, larger than 500 ?m for paleoreconstructions.

(Table 1) Boron concentration and d11B values of interstitial water from IODP Holes 308-U1319A, U1320A, U1322B and U1324B

The Integrated Ocean Drilling Program Expedition 308 (IODP308) drilled normal-pressured sediments from the Brazos-Trinity Basin IV and over-pressured sediments from the Ursa Basin on the northern slope of the Gulf of Mexico. The interstitial water samples from the normal-pressured basin show B concentrations and B isotopic compositions ranging from 255 to 631 µM (0.6 to 1.5 times of seawater value) and from +29.1 to +42.7 per mil (relative to NIST SRM 951), respectively. A wider range is observed both for B concentrations (292 to 865 µM, 0.7 to 2.1 times of seawater value) and d11B values (+25.5 to +43.2 per mil) of the interstitial water in the over-pressured basin. The down-core distribution of B concentrations and d11B values in the interstitial waters are sensitive tracers for assessing various processes occurring in the sediment column, including boron adsorption/desorption reactions involving clay minerals and organic matter in sediments as well as fluid migration and mixing in certain horizons and in the sediment column. In the normal-pressured basin adsorption/desorption reactions in shallow sediments play the major role in controlling the B content and B isotopic composition of the interstitial water. In contrast, multiple processes affect the B content and d11B of the interstitial water in the over-pressured Ursa Basin. There, the stratigraphic level of the maxima of B and d11B correspond to seismic reflectors. The intruded fluids along the seismic reflector boundary from high to low-topography mix with local interstitial water. Fluid flow is inferred in the Blue Unit (a coarse sandstone layer, connecting the high- to low-pressured region) from the freshening of interstitial water in Ursa Basin Site U1322, and upward flow by the overpressure expels fluid from the overburden above the Blue Unit.

(Table S1) Globigerinoides sacculifer Boron isotope ratios of ODP Hole 165-999A

everal hypotheses have been put forward to explain the onset of intensive glaciations on Greenland, Scandinavia, and North America during the Pliocene epoch between 3.6 and 2.7 million years ago (Ma). A decrease in atmospheric CO2 may have played a role during the onset of glaciations, but other tectonic and oceanic events occurring at the same time may have played a part as well. Here we present detailed atmospheric CO2 estimates from boron isotopes in planktic foraminifer shells spanning 4.6-2.0 Ma. Maximal Pliocene atmospheric CO2 estimates gradually declined from values around 410 µatm to early Pleistocene values of 300 ?atm at 2.0 Ma. After the onset of large-scale ice sheets in the Northern Hemisphere, maximal pCO2 estimates were still at 2.5 Ma +90 µatm higher than values characteristic of the early Pleistocene interglacials. By contrast, Pliocene minimal atmospheric CO2 gradually decreased from 310 to 245 µatm at 3.2 Ma, coinciding with the start of transient glaciations on Greenland. Values characteristic of early Pleistocene glacial atmospheric CO2 of 200 ?atm were abruptly reached after 2.7 Ma during the late Pliocene transition. This trend is consistent with the suggestion that ocean stratification and iron fertilization increased after 2.7 Ma in the North Pacific and Southern Ocean and may have led to increased glacial CO2 storage in the oceanic abyss after 2.7 Ma onward.