Decoupled carbonate chemistry controls on the incorporation of boron into Orbulina universa

In order to fully constrain paleo-carbonate systems, proxies for two out of seven parameters, plus temperature and salinity, are required. The boron isotopic composition (δ11B) of planktonic foraminifera shells is a powerful tool for reconstructing changes in past surface ocean pH. As B(OH)−4 is substituted into the biogenic calcite lattice in place of CO2− 3 , and both borate and carbonate ions are more abundant at higher pH, it was suggested early on that B /Ca ratios in biogenic calcite may serve as a proxy for [CO2− 3 ]. Although several recent studies have shown that a direct connection of B /Ca to carbonate system parameters may be masked by other environmental factors in the field, there is ample evidence for a mechanistic relationship between B /Ca and carbonate system parameters. Here, we focus on investigating the primary relationship to develop a mechanistic understanding of boron uptake. Differentiating between the effects of pH and [CO2− 3 ] is problematic, as they co-vary closely in natural systems, so the major control on boron incorporation remains unclear. To deconvolve the effects of pH and [CO2− 3 ] and to investigate their impact on the B /Ca ratio and δ11B, we conducted culture experiments with the planktonic foraminifer Orbulina universa in manipulated culture media: constant pH (8.05), but changing [CO2− 3 ] (238, 286 and 534 μmol kg −1 CO2− 3 ) and at constant [CO2− 3 ] (276± 19.5 μmol kg −1) and varying pH (7.7, 7.9 and 8.05). Measurements of the isotopic composition of boron and the B /Ca ratio were performed simultaneously using a femtosecond laser ablation system coupled to a MC-ICP-MS (multiple-collector inductively coupled plasma mass spectrometer). Our results show that, as expected, δ11B is controlled by pH but it is also modulated by [CO2− 3 ]. On the other hand, the B /Ca ratio is driven by [HCO−3 ], independently of pH. This suggests that B /Ca ratios in foraminiferal calcite can possibly be used as a second, independent, proxy for complete paleo-carbonate system reconstructions. This is discussed in light of recent literature demonstrating that the primary relationship between B /Ca and [HCO−3 ] can be obscured by other environmental parameters.


Introduction
Before the Anthropocene, the atmospheric CO 2 concentration was governed by the surface ocean [CO 2 ], simply because the carbon content of the ocean is 65 times larger than that of the atmosphere (Siegenthaler and Sarmiento, 1993).Hence, understanding the global carbon cycle and the evolu-Published by Copernicus Publications on behalf of the European Geosciences Union.tion of atmospheric pCO 2 in Earth history requires knowledge of the dynamics of the oceanic carbonate chemistry.Since the industrial revolution, the unprecedented magnitude and rate of carbon emissions has caused both warming and acidification of the oceans (Bijma et al., 2013;Ciais et al., 2013;Gattuso and Hansson, 2011;Gattuso et al., 2015;Rhein et al., 2013).As a consequence, the interest in the reconstruction of seawater carbonate chemistry to identify ocean acidification in Earth history experienced an impetus (Hönisch et al., 2012;Martínez-Botí et al., 2015a).
Reconstruction of the full oceanic carbonate chemistry requires proxies of at least two independent parameters of the carbonate system, in addition to temperature and salinity.However, to date, all reconstructions are based on the analysis of δ 11 B of biogenic carbonates alone with assumptions regarding a secondary parameter.In these reconstructions, total alkalinity (A T ) or [CO 2− 3 ] was estimated from modern ocean conditions or from reconstructions of the carbonate compensation depth (CCD).Total alkalinity is a conservative parameter, meaning that A T is linearly correlated with salinity (Dickson, 1981(Dickson, , 1992;;Wolf-Gladrow et al., 1999, 2007).Therefore, if it is assumed that the modern salinity-A T relationship was constant over time, A T can be estimated from reconstructions of salinity using sea-level records (Foster, 2008;Hönisch et al., 2009).However, salinity and alkalinity may be decoupled in space and time through weathering and changes in riverine alkalinity input.In addition, reliable proxies for regional salinity reconstructions have yet to be developed.Another approach is based on the assumption that seawater [Ca 2+ ] has remained proportional to A T over time so that A T can be adjusted in a way that the water column is exactly saturated with respect to calcite at the lysocline (∼ 500 m above the CCD; Pearson and Palmer, 2000).Surface A T can be estimated by assuming that increases in A T with depth were the same as in the modern ocean.The CCD, however, is not uniform through space and time (Van Andel, 1975), calling into question these approaches for estimating past A T .Pearson and Palmer (2000) note that "the CCD record for the Palaeogene Pacific Ocean is relatively poorly constrained".
Although δ 11 B has proven to be a reliable proxy for pH and one can argue that ocean pH is the main driver of past atmospheric CO 2 , it is important to remember that changes in past glacial interglacial atmospheric pCO 2 can be achieved via two end-member scenarios (e.g.Lea et al., 1999;Sanyal and Bijma, 1999).In the first scenario, changes in carbonate chemistry are brought about by changes in total dissolved inorganic carbon (C T ) only.This is equivalent to varying the response of the biological pump as a reaction to variations in the nutrient content of the surface ocean.In the second scenario, changes in carbonate chemistry are solely controlled by addition (due to dissolution in sediments) or removal (due to production) of calcium carbonate.The change in surface ocean carbonate chemistry is very different in these two scenarios because the ratio of carbonate ion increase in relation to pCO 2 decrease, depending on surface ocean alkalinity (Lea et al., 1999).A smaller change is associated with the drawdown of C T under conditions of unchanging alkalinity (e.g.strengthening the biological pump without calcite compensation).The change in surface water [CO 2 ] is twice as much when the same atmospheric pCO 2 is reached solely via a change in alkalinity as in the coral reef hypothesis (Lea et al., 1999).These dependencies are nicely demonstrated in Fig. 1.1.3 of Zeebe and Wolf-Gladrow (2001) and Fig. 1 of Foster and Rae (2016).
The real ocean operates somewhere between these end member scenarios and, basically, depends on the relative delivery rates of calcium carbonate and particulate organic carbon (the CaCO 3 : POC "rain ratio") and the sensitivity of calcium carbonate preservation in deep ocean sediments.This demonstrates that a second, pH-independent, parameter could reduce the uncertainty in CO 2 estimates.On the other hand, the propagated uncertainty in the second parameter, reconstructed from an independent proxy (taking into account measurement and calibration uncertainty) might not be much lower than the margin of error that is garnered using assumptions around, for example, total alkalinity.
The boron isotope pH proxy in foraminifera has recently been reviewed by Foster and Rae (2016) and we refer to that for further reading.Here, we will briefly explain the boron systematics.Boron exists in seawater primarily in the form of two species, boric acid (B(OH) 3 ) and borate ion (B(OH) − 4 ; Fig. 1a).As for all weak acids, the relative abundance between these two species is controlled by pH (Dickson, 1990;DOE, 1994).At low pH (< 7), nearly all boron is present in the form of boric acid, whereas at high pH (> 10), boron primarily exists as borate.Because of the isotopic fractionation between the two aqueous species (Fig. 1b; ), the boron isotopic composition of each species is also pH-dependent (Hemming and Hanson, 1992;Palmer et al., 1987;Sanyal et al., 1996Sanyal et al., , 2000)).B(OH) 3 is enriched in the stable isotope 11 B compared to B(OH) − 4 , with a constant isotopic fractionation of 27.2 ‰ between the two boron species (Klochko et al., 2009(Klochko et al., , 2006)).Consequently, as the relative concentration of the dissolved species changes with pH, so does their isotopic composition.Because it is assumed that only the charged species, borate, is 4 and B(OH) 3 with thermodynamic fractionation factor (α 3−4 ) = 1.0272 (Klochko et al., 2006).4 ] / C T ) in Globigerinoides ruber from culture experiments.However, this relationship was completely lost in the plankton tow samples and the sediments they analysed.While they explicitly tested for a carbonate chemistry control on B / Ca, they found a strong relationship to [PO − 4 ] and neither a correlation with carbonate system parameters nor a covariation of phosphate with carbonate system parameters.They concluded that apparently B / Ca in G. ruber is controlled by [PO −  4 ].We will discuss why we believe that the primary (mechanistic) relationship explaining B / Ca is probably still controlled by carbonate chemistry parameters in the ambient environment of the foraminifer, but that it may be masked in the field and decoupled from the bulk seawater carbonate chemistry.
Here, we are specifically focussing on the primary controls of boron uptake and conducted experiments with the planktonic foraminifer Orbulina universa and decoupled pH and [CO 2− 3 ] in the same way as Kaczmarek et al. (2015b).We show, in principle, that combined measurements of δ 11 B calcite and B / Ca on single shells of planktonic foraminifera might be used to fully constrain the carbonate chemistry in downcore records.However, based on recent publications (Allen et al., 2012;Babila et al., 2014;Henehan et al., 2015;Salmon et al., 2016), it becomes increasingly clear that in the field and downcore, B / Ca may not be a very robust carbonate system proxy (at least in some species) as the primary relationship can be masked by other environmental factors.

Collection and culturing
Living specimens of O. universa were collected daily using a 57 cm diameter WP2 plankton net (200 µm mesh size), between July and September 2012 at Point B, Villefranchesur-Mer, France (43.41 • N, 7.19 • E), and maintained until gametogenesis in laboratory cultures at the Laboratoire d'Oceanographie de Villefranche.Established procedures for maintaining planktonic foraminifera in laboratory culture were used (Bemis et al., 1998;Bijma et al., 1998;Spero and Lea, 1993).Briefly, specimens were identified, the diameters measured with a light microscope, and they were then transferred to 0.2 µm-filtered seawater, whose carbonate chemistry was accurately determined and subsequently modified.Specimens were maintained individually in airtight 100 mL acid-washed SCHOTT DURAN ® bottles that were sealed without an air space and placed upside down into thermostated water baths maintained at a temperature of 23 • C (±0.2 • C).Light was provided by four 39 W fluorescent tubes (JBL Solar Ultra Marin Day), with reflectors (at a distance of ca. 15 cm from the water surface), with a 12:12 h L:D photoperiod.The average irradiance, measured with a LI-193 sensor (LiCOR) in the culture jars was about 290 µmol photons m −2 s −1 .
The foraminifers were fed a one-day-old brine shrimp Artemia nauplius every second day until gametogenesis.The brine shrimp were hatched in modified seawater from the same batch as used for culturing the foraminifera.Just prior to feeding, hatched nauplii were transferred once again to fresh medium from the same batch.After feeding, culture jars were topped up with medium from the same batch to prevent the formation of a headspace.Empty shells were collected within 24 h after successful gametogenesis, rinsed in deionized water and archived in covered micro-paleo-slides for later analysis.Prior to analysis, specimens were harvested, bleached in NaOCl (active chlorine: 4.6 %) for 6 h, rinsed four times using deionized water, and dried for 12 h at 50 • C. Approximately 35 tests were grown for each experimental treatment.Culture water samples were collected at the start and end of the experiments to verify the boron concentration, its isotopic composition and the carbonate system parameters.

Modified seawater chemistry
The objective of these experiments was to decouple seawater pH and [CO 2− 3 ] and create treatments with a constant pH and varying carbonate ion concentration and treatments with a constant carbonate ion concentration but varying pH.To decouple the effects of pH T and [CO 2− 3 ], seawater carbonate chemistry was modified by manipulating pH T , using NaOH and HCl, and dissolved inorganic carbon (C T ) by adding gravimetrically carbonate and bicarbonate or bub-bling with CO 2 .Calculations were made using csys_vari.m(Zeebe et al., 2001) with carbonic acid dissociation constants of Mehrbach et al. (1973).Temperature (23 • C) and salinity (38.0) were kept constant (Table 1).
To enable single-shell analysis by LA-MC-ICP-MS, the boron concentration was increased to 10 times the concentration of natural seawater by adding boric acid to the culture water (see Sanyal et al., 2001Sanyal et al., , 2000)).The pH T and C T were then modified via titration with boron-free NaOH (1N) and HCl (1N) to bring the experimental pH to desired levels of 7.70 ± 0.03, 7.90 ± 0.02, and 8.05 ± 0.05.Culture water samples collected at the start and at the end of each experiment showed that pH remained nearly constant throughout each experiment.The boron isotopic composition of each culture treatment is provided in Table 1.The pH of the culture solutions was measured using a Metrohm, 826 mobile pH meter with a glass electrode (Metrohm, electrode plus) calibrated to the total scale using TRIS and 2-aminopyridine buffer solutions (Dickson et al., 2007) adjusted to a salinity of 38.0.Total alkalinity (A T ) samples (150 mL) were filtered on GF/F and measured potentiometrically using a Metrohm Tritando 80 titrator and a Metrohm, electrode plus glass electrode (Dickson et al., 2007).60 mL samples were also taken at the start and end of incubations and poisoned with 10 µL of saturated HgCl 2 pending determination of dissolved inorganic carbon (C T ).Samples were measured using an AIR-ICA (Marianda, Kiel) fitted with a Licor 6262 infrared gas analyser.All parameters of the carbonate system were calculated from A T and pH T (Hoppe et al., 2012) using the R package seacarb (Lavigne and Gattuso, 2013).

Culture water analysis
Boron isotopic composition of the culture media were analysed by means of a Thermo ® Element XR, a single collector, sector field, high-resolution inductively coupled plasma mass spectrometer, fitted with a high-sensitivity interface pump (Jet pump) as described in Misra et al. (2014).Boron isotopic composition is reported as per mil (‰) deviation from NIST SRM 951a ( 11 B / 10 B = 4.04362 ± 0.00137) (Catanzaro et al., 1970) where: Boron isotope analyses were made following a Sample-Standard Bracketing (SSB) technique.NIST 951a was used as the standard and samples were concentration-matched, typically at 5 %, with the standard and were analysed in quintuplicate.The accuracy and precision of the analytical method was assessed by comparing 11 B measurements of seawater (from the Atlantic Ocean) and secondary boron standards (AE 120,121,122) with published (accepted) results.Our estimates of δ 11 B SW of 39.8 ± 0.4 ‰ (2 SE, n = 30) are independent of sample size and are in agreement  (Vogl and Rosner, 2012).Information about sample preparation for analysis can be found in the Supplement provided in Kaczmarek et al. (2015a).

Analysis of O. universa
For simultaneous determination of the B isotopic composition and its concentration, a Fiber Optics Spectrometer (Maya2000 Pro, Ocean Optics) was connected to the torch of a Thermo Finnigan Neptune multiple-collector inductively coupled plasma mass spectrometer (MC-ICP-MS) at the Leibniz University of Hannover.Laser ablation on reference material and samples was performed by an in-house-built UV-femtosecond laser ablation system based on a regenerative one-box femtosecond laser (Solstice Newport/Spectra Physics).A detailed description of the method used for the simultaneous determination of B concentration and δ 11 B of O. universa can be found in Kaczmarek et al. (2015a).A summary of the procedure is given below.

Simultaneous determination of B concentration and δ 11 B
The B intensity of a reference material corresponds to its known B concentration.Based on this relationship, the unknown B concentration of a sample can be calculated.However, our measurements of the reference material (NIST SRM 610) and samples were not performed at the same laser repetition rate: hence, their B ratios are not proportional.Because Ca concentrations in the reference material and in the sample are known (NIST SRM 610: 8.45 %, CaCO 3 : 40 %) a correction for different laser repetition rates was realized by the analysis of calcium using the optical spectrometer.More information on this procedure is provided by Longerich et al. (1996).

Calcium analysis
The Maya2000 Pro is a high-sensitivity fiber optical spectrometer.It has a measuring range between 250 and 460 nm with a resolution of 0.11 nm covering the first order emission lines of Mg II, Ca II, Sr II, Ba II and Li II.It is equipped with a back-thinned 2-D FFT-CCD detector, and a grating with a groove density of 1200 lines mm −1 .The optical fiber used is 2 m long (attenuation of the photon flux is length-dependent), connecting the spectrometer with the coupling lens at the end of the plasma torch of the MC-ICP-MS (Thermo Finnigan Neptune).Ca II ion lines were measured at a wavelength of 393.48 and 396.86 nm.At these wavelengths the Ca spectrum shows no detectable interferences for the matrices used.The acquisition parameters were set to acquire 220 cycles per analysis with an integration time of 1 s for each cycle.For the first 40 cycles, only background (BG) signal was detected prior to measuring the sample.The BG signal detected at the start of the analysis was later used for correcting sample measurements by subtracting BG intensity from the intensity of the reference and the sample material.

Boron isotope analysis -194 nm femtosecond laser ablation
The in-house-built laser ablation system is based on a 100 femtosecond Ti-sapphire regenerative amplifier system operating at a fundamental wavelength of 777 nm in the infrared spectrum.Subsequent harmonic generations produce the wavelengths 389 nm in the second, 259 nm in the third and 194 nm in the fourth harmonic.The pulse energies measured with a pyroelectric sensor (Molectron, USA) are 3.2 mJ pulse −1 at 777 nm, 0.7 mJ pulse −1 at 259 nm, and 0.085 mJ pulse −1 at 194 nm.After the fourth harmonic generation stage, the 194 nm beam is steered by eight dichronic mirrors into an 8 × objective (NewWave Research, USA) and focussed onto the outside of the sample.Spot size was set to 50 µm for the reference material and the samples.Within this spot, an energy density of ∼ 2 J cm −2 is maintained.Reference material measurements were performed in raster mode (100 µm × 100 µm) at 10 Hz and samples were ablated at 8-50 Hz depending on B concentration.
It should be noted that the fs laser ablation process is fundamentally different from ns laser ablation.When the pulse length is shorter than 1 ps (Hergenröder et al., 2006), the laser energy can be deposited into the material before it can thermally equilibrate.Femtosecond ablation also provides smaller aerosol particle sizes.Due to the short pulse length, fs laser ablation shows no detectable matrix dependency (e.g.Chmeleff et al., 2008;Horn et al., 2006;Kaczmarek et al., 2015a;Lazarov and Horn, 2015;Oeser et al., 2014;Schuessler and von Blanckenburg, 2014), i.e. it does not require a matrix-matched standard and therefore permits the use of NIST SRM 610 (a glass) as a reference for carbonates.As boron concentrations differ between sample and standard, and different matrices require more or less energy for ablation, the repetition rate was chosen such that the signal of sample and standard at the ion counters was comparable.This is important for normalization of the sample to the known δ 11 B of the standard and also accounts for the imprecision of the determined detector dead time.
Most previous publications on boron isotopes have used "wet chemistry" for which NIST SRM 951 is a perfect standard.We have also used this standard for the analysis of the culture waters.The foraminiferal shells, however, were referenced against NIST SRM 610.As shown by several studies (Fietzke et al., 2010;Kasemann et al., 2001;Le Roux et al., 2004), both standards are, within analytical uncertainty, isotopically equal.Hence, for comparison between δ 11 B of O. universa and δ 11 B of B(OH) − 4 , the isotopic difference between the two standards can be neglected and it does not make a difference if values are reported vs. one or the other standard.

Boron isotope analysis -acquisition parameters
All measurements are carried out in low mass resolution ( m/m = 350 where m is the mass of the ion of interest and m is the mass difference between its 5 and 95 % peak height).Compact discrete dynode multipliers (CDDs, Thermo) are attached to Faraday cups at the low site on L4 and the high site on H4.The low-resolution mode is sufficient to resolve potential interferences from doubly charged ions due to the intrinsic high resolution in the low mass region.Possible interferences are the clusters of 40 Ar 4+ or 20 Ne 2+ , which are well resolved to the background level.The instrument was tuned prior to each analytical session for optimal peak shape.Instrumental operating conditions are reported in Table 2.All measurements were performed at plateau voltage of the CDDs, which was checked prior to every analytical session.Before the beginning of sample analysis, measurements of NIST SRM 610 were continued until instrumental drift (due to warm-up) was less than 200 ppm over a bracketing sequence duration of twelve minutes.Boron signal intensities of NIST SRM 610 and samples were matched within 10 % in signal intensity by adapting the laser repetition rate.The acquisition parameters in static mode for analysis of NIST SRM 610 and samples were set to acquire 200 cycles of 1 s integrations each.During the first 40 cycles the background signal was acquired, whereas the remaining cycles represent the sum of the background and the reference material, or the background and the sample signals.A complete measurement consisting of 200 cycles of a single reference material or sample took 4 min before the next sample was introduced.For analysis we adopted the standard sample bracketing procedure and the B isotopic composition is reported using the delta notation: where NIST 610 − 1 and NIST 610 + 1 refer to the analysis of the reference material before and after the sample.The uncertainty of the samples was calculated according to: where 11/10 B ratios represent mean values of the reference material and the sample calculated from one measurement (based on 160 cycles) and SE represents the standard error of the 11/10 B ratios.Due to the natural inhomogeneity of the samples, the analytical uncertainty is represented best by repeated measurements of the homogenous reference material given by: where the measurements of the ( 11/10 B) −1 and ( 11/10 B) +1 ratios of NIST 610 were performed before and after the measurement of ( 11/10 B) 0 , respectively.For the determination of the analytical uncertainty and external reproducibility, all measurements of NIST 610 performed between each sample measurement were taken into account.On average, the analytical uncertainty and external reproducibility is 0.66 ‰.
2.9 Conversion of δ 11 B O. universa to natural seawater Due to the additional B addition to our culture media the δ 11 B seawater shifted from 37.63 (Mediterranean) to, on average, 4.66 ‰ (Table 1).Therefore, the δ 11 B O. universa shifted accordingly.In order to compare our O. universa data to published values (Fig. 3a), the measured δ 11 B from each experiment was normalized to natural seawater using the following (Zeebe and Wolf-Gladrow, 2001): where ε is (α sw−msw − 1) × 1000, δ 11 B c represents the converted δ 11 B for the measured value (δ 11 B m ), α sw−msw is the fractionation factor expressing the difference between the natural seawater and manipulated seawater: (5)

Statistics
Lamtool (a modified Excel spreadsheet, initially programmed by Jan Kosler, University of Bergen, Norway) was used for analysis and background correction of the δ 11 B data.All other statistics were carried out using R (R Core Team, 2008).Error bars represent ±2σ errors, correlations were calculated by linear regression.The procedures for data evaluation, background correction and uncertainty calculations for boron concentration and isotopes are extensively described in Kaczmarek et al. (2015a).
In contrast to "wet chemical" analysis, laser ablation (LA) records the inhomogeneous boron distribution ("boron banding", see Branson et al., 2015) within a specimen and individual shell analysis captures inter-specimen differences.Sadekov et al. (2016) demonstrated that the variability in both B / Ca and δ 11 B recurs in each chamber and, therefore, represents real data of high quality.This is supported by the fact that the values of the averaged laser data are very close to wet chemical analyses where multiple specimens are dissolved and the intra-and inter-variability is "averaged" before the analysis.The intra-specimen δ 11 B variability in Cibicidoides wuellerstorfi is up to ca. 10 ‰ (Sadekov et al., 2016), while the inter-specimen δ 11 B variability of Amphistegina lessonii from the same treatment is ca.6 ‰ (Kaczmarek et al., 2015b).Histograms of single-foram δ 11 B measurements from each of our pH treatments (Fig. S1 in the Supplement) show that the LA data is normally distributed (p values from Shapiro-Wilk tests are all higher than 0.05).This is confirmed by the box plots (Fig. S1) where the average and median values are very close to each other.The relatively large standard errors of laser ablation analyses are representative of true inter-specimen variability and largely unrelated to analytical errors.Therefore, the relatively large standard errors do not present a limitation for how much can be interpreted from the data.The major difference between LA and "wet chemistry" data is that the latter method averages individual variability before analysis by measuring multiple dissolved shells in one go, while LA captures individual variability (which is large and real as argued above) and averages afterwards.
One could further argue that the uncertainty stemming from the analysis of culture water δ 11 B should also be propagated when plotting in "normal δ 11 B sw " space (Table S4 in the Supplement).The propagated error is, of course, large as it includes the individual δ 11 B variability of the foraminifers.It is important to acknowledge that this variability represents true data, which is largely unrelated to analytical uncertainty.We added a calcite vs. borate δ 11 B cross-plot (Fig. 4) to avoid the conversion into the seawater scale and making the error propagation obsolete.However, as not all studies report the parameters required for the calculation of δ 11 B of borate, we plotted for comparison in "normal δ 11 B sw " space (Fig. 3a) but did not propagate the error related to the analysis of culture water δ 11 B.

B / Ca ratios
The B / Ca ratio of O. universa shows a strong negative correlation (R 2 = 0.96) with C T , irrespective of the pH T of the culture media (Fig. 2a).It is also correlated to [CO 2 ] but to a lesser extent (R 2 = 0.64; Fig. 2c).B / Ca also decreases with increasing [CO 2− 3 ] in specimens grown under a pH T of 8.05, (Fig. 2e).However, the B / Ca ratio of specimens grown under lower pH T values (7.9 and 7.7) is negatively offset from the relationship found at pH T 8.05 and the overall correlation of B / Ca and [CO 2− 3 ] is very low (R 2 = 0.2; Fig. 2e).Of all the carbonate species, the B / Ca ratio exhibits the best negative relationship with increasing [HCO − 3 ], irrespective of the pH T of the culture medium (R 2 = 0.96; Fig. 2g).Plotted against the ratio of [B(OH) − 4 ], over each of the carbon species (Fig. 2b, d, f, h), the correlations are high for all combinations but highest for [B(OH . Based on first principles, we predict a positive correlation between www.biogeosciences.net/14/415/2017/Biogeosciences, 14, 415-430, 2017 , error bars represent standard error.
B / Ca and δ 11 B because at higher pH, not only does the isotopic composition of borate become heavier, but its concentration also increases.

Boron isotopic fractionation (δ 11 B)
Single, measured δ 11 B values of O. universa are given in Table S1 and errors are calculated according to Eq. ( 2).Median and converted values using Eqs.( 4) and ( 5) are shown in Fig. 3a and Table 3.The fractionation of boron isotopes in the shells of O. universa is dependent on the pH of the culture medium, increasing with pH T from 15 ‰ at pH T 7.7 to 18.8 ‰ at pH T 8.05.These values are close to the B fractionation curve of B(OH) − 4 obtained for artificial seawater by Klochko et al. (2006;Fig. 3a).δ 11 B increases with increasing [CO 2− 3 ] at constant pH T from 17.2 ‰ at 238 µmol kg −1 CO 2− 3 to 19.9 ‰ at 534 µmol kg −1 CO 2− 3 (Table 3; Fig. 3b).Applying ANOVA with a Bonferroni test, which is best suited for a limited number of pairs, the p value of the overall ANOVA is 0.00203, demonstrating a significant difference between two or more population means.The differences between the mean δ 11 B values of the [CO 2− 3 ] treatments 239 and 286 µmol kg −1 were close to significance but only between 239 and 534 µmol kg −1 ; the difference was significant (Table S3).Because this range in [CO 2− 3 ] is beyond that of the real ocean and because pH and [CO 2− 3 ] covary, we believe that this observation is only important for a better understanding of the δ 11 B controls and does not significantly impact existing calibrations., he also demonstrated that B / Ca in combination with δ 11 B can be used to fully constrain the carbonate system in downcore records.Nonetheless, he identified [CO 2− 3 ] as having a major (secondary) control on B / Ca in samples of foraminifera from down-core samples and core tops.A similar conclusion was reached by Allen et al. (2011) for O. universa.These authors demonstrated a trend of decreasing B / Ca with increasing pH and [CO 2− 3 ]; however, due to the co-variations of the carbonate system in natural seawater, it is difficult to identify the differential effects of the individual parameters.Allen and Hönisch (2012) conclude that the relationships between K D and seawater parameters can sometimes be driven by the denominator of the empirical boron partition coefficient ([B(OH) − 4 ] / [HCO − 3 ]), and not by B / Ca of seawater itself.Reconstructions based on such B / Ca-independent relationships are susceptible to being driven by other environmental parameters.They conclude that application of the empirical boron partition coefficient should be avoided until more is known about the relative influences of different chemical species on boron incorporation.
www.biogeosciences.net/14/415/2017/Biogeosciences, 14, 415-430, 2017 the δ 11 B values of borate ion in artificial seawater given by Klochko et al. (2006).This is probably caused by the suppression of the vital effects imposed by O. universa.Theoretical considerations demonstrate that at 10 × boron concentration compared to natural seawater, vital effects are suppressed and the isotopic value of biogenic calcite approaches the value of borate (Zeebe, 2003).This was confirmed by the comparison of the boron isotopic values of O. universa grown at low and high light (Hönisch et al., 2003) and supports the notion that borate is, indeed, the species being taken up.There is a trend of varying [CO 2− 3 ] on δ 11 B of samples grown at the same pH but, most importantly, in light of the results obtained for the B / Ca ratio, there is no effect of [HCO − 3 ] (Fig. 3c).

Proxy implications
A sound understanding of the effects of past carbon perturbations becomes increasingly urgent in an age where anthropogenic activities are producing such rapid changes in global climate (Bijma et al., 2013;Knoll and Fischer, 2011).The usefulness of biogeochemical proxies to reconstruct paleooceanographic conditions is well established for many environmental parameters (Wefer et al., 1999) but uncertainties remain for proxies related to pH and the carbonate system (Allen and Hönisch, 2012;Hönisch et al., 2007;Katz et al., 2010;Pagani et al., 2005).This study confirms the robustness of δ 11 B as an independent pH proxy and supports the growing body of evidence that B / Ca in planktonic foraminiferal calcite is mechanistically controlled by [HCO 2− 3 ] (Yu et al., 2007), thereby allowing researchers to fully constrain the carbonate system in combination with δ 11 B.
Based on our results and other culture studies, it becomes clear that despite strong biological effects on the ambient carbonate chemistry (Köhler- Rink andKühl, 2001, 2000;Rink et al., 1998;Wolf-Gladrow et al., 1999;Zeebe et al., 2008), the boron isotopic composition and the B / Ca are faithful predictors of seawater pH and bicarbonate ion concentration, respectively.Our results provide strong evidence that [HCO − 3 ] is the primary control of the B / Ca ratio.The correlation of the B / Ca ratio to [HCO − 3 ] rather than to [CO 2− 3 ] might have some implications for existing paleo-carbonate chemistry reconstructions based on this proxy such as the study by Foster (2008) and that of Yu et al. (2014), since it seems reasonable to assume that the same relationship probably holds for benthic foraminifers as for planktonics.
A wide range of [HCO − 3 ] was necessary to facilitate decoupling the carbonate system from pH T .The high [HCO − 3 ] in some of these treatments are unrealistic for natural seawater systems and more environmentally-relevant values should be used for future calibration experiments.The proxy should, therefore, be ground-truthed using water column and core top samples.
Recently, Henehan et al. (2015) showed that B / Ca in G. ruber collected with a plankton net was perfectly corre-  3).E. L. Howes et al.: Controls on the incorporation of boron into Orbulina universa boron partitioning in inorganic precipitation experiments increases with increasing growth rate, and early work by the pioneers of foraminiferal biology and calcification (e.g.Bé et al., 1982;Bé, 1965;Caron et al., 1982;Hemleben et al., 1987;Jørgensen et al., 1985;Spero and Parker, 1985) clearly demonstrated the huge impact of symbionts on foraminiferal shell growth.Interestingly, Babila et al. (2014) state: "The seasonal cycle of B / Ca in G. ruber white was more strongly correlated with light intensity than with temperature.Both observations suggest that the presence of symbionts in G. ruber and seasonal variability in their photosynthetic activity act to modify the internal pH during calcification, by up to 0.2 units relative to ambient seawater."This supports our line of argumentation above.
In another recent paper on B / Ca, Salmon et al. ( 2016) write: "We provide the first evidence for a strong positive relationship between area density (test thickness) and B / Ca, and reveal that this is consistent in all species studied, suggesting a likely role for calcification in controlling boron partitioning into foraminiferal calcite."Their conclusion also supports our reasoning, that, mechanistically, increased photosynthesis may lead to higher calcification rates.Remarkably, Salmon et al. (2016) show that B / Ca of the nonsymbiont-bearing species (Globigerina bulloides and Globigerina inflata) and even the symbiont-bearing species G. sacculifer are related to [CO 2− 3 ] and [B(OH) − 4 / HCO − 3 ].In our view, those results demonstrate the primary control by carbonate chemistry parameters not masked by symbiont photosynthesis.One could even argue that there is a positive trend for O. universa but that the natural range of [CO 2− 3 ] variability (or borate/bicarbonate) is small (ca.20 µmol kg −1 in the depth range 30 to 50 m) in comparison to the decoupling we carried out in controlled culture experiments.Interestingly, Henehan et al. (2016) propose a field calibration for O. universa that is very close to δ 11 B of borate, suggesting that their "vital effects" are muted in the real ocean, especially the symbiont impact of raising the calibration curve above δ 11 B of borate.This is supported by the observation of Hemleben and Bijma (1994) that O. universa occupies a subsurface maximum (in the Red Sea) between 20 and 60 m (Hemleben and Bijma, 1994, Fig. 5) and could explain why B / Ca in this species is not (completely) masked by symbiont photosynthesis (Salmon et al., 2016).
Our final conclusion is that, although controlled laboratory studies are the only means to clarify the mechanisms of proxy incorporation, field studies are required to determine to what extent vital effects determine species-specific offsets from the target parameters.For instance, the light level used in the culture experiments of Sanyal et al. (2001) was ∼ 380 µmol photons m −2 s −1 , providing a photon flux for maximum photosynthetic rates (P max ) of the symbionts (Spero and DeNiro, 1987;Spero and Parker, 1985;Spero and Williams, 1988).Consequently, the impact of photosynthesis on the G. sacculifer calibration of Sanyal et al. (2001) is fully expressed.However, in the real ocean this species may experience lower irradiance, shifting the calibration curve more towards the borate values.In our study, the average irradiance in the culture jars was about 290 µmol photons m −2 s −1 , which is well below P max of the symbionts and apparently closer to the irradiance conditions of their natural depth habitat.Therefore, the impact of photosynthesis is muted (Hönisch et al., 2003;Zeebe, 2003) and our laboratory calibration closer to the field calibration of (Henehan et al., 2016).
Laboratory experiments are usually carried out with foraminifera selected as model organisms for ease of availability and ability to be maintained in culture but, generally, state nothing about their suitability for paleo-studies.Field studies are much better to identify which species are best suited for down-core reconstructions.We agree with Henehan et al. (2015) that G. ruber is not a good choice for B / Ca as its primary relationship to carbonate chemistry parameters is, apparently, not very robust.However, other symbiontbearing species, non-symbiotic planktonic foraminifera and deep-sea benthic foraminifera, may still be a viable option to use B / Ca for carbonate chemistry reconstructions.

Data availability
Data are included in the Supplement and archived on the Pangaea Database (Howes et al., 2016).
The Supplement related to this article is available online at doi:10.5194/bg-14-415-2017-supplement.
Competing interests.The authors declare that they have no conflict of interest.
E. L. Howes et al.: Controls on the incorporation of boron into Orbulina universa 2 Methods
Figure S2 shows the individual B / Ca, δ 11 B pairs per treatment.As expected for individual LA shell analyses, the inter specimen variability is quite large.Individual B / Ca ratios vary by almost 50 % in each treatment and individual δ 11 B values vary by ca.4-6 ‰ per treatment (cf.Kaczmarek et al., 2015b).Although one could argue for a positive trend between B / Ca and δ 11 B in some of the treatments, we believe that the individual B / Ca, δ 11 B pairs within a treatment are uncorrelated.However, the average values for the four treatments with [CO 2− 3 ] between 238-297 µmol kg −1 do show a positive correlation between B / Ca and δ 11 B. The "outlier" (treatment at pH 8.05; [CO 2− 3 ] = 534) can be explained by the high [HCO − 3 ], relative to the other pH = 8.05 treatments.
showed that the partition coefficient K D for the B / Ca ratio is influenced by [CO 2− 3 ] (and temperature).Although complicating the application as a proxy related to [B(OH) − 4 ] / [HCO − 3 ]

Table 1 .
Average properties of the manipulated seawater culture medium from four samples (two from the start of the incubation and two from the end of the incubation).

Table 2 .
Instrumental operating conditions for the MC-ICP-MS and LA.