ALBERT

Description: The boron isotope composition (d11B) of planktic foraminiferal calcite, which reflects seawater pH, is a well-established proxy for reconstructing palaeo-atmospheric CO2 and seawater carbonate chemistry. However, to translate d11B measurements determined in calcareous fossils into pH we need to know the boron isotope composition of the parent seawater (d11Bsw). While a number of d11Bsw reconstructions exist, the discrepancies between them reveals uncertainties and deficiencies that need to be addressed. Here we present a new d11Bsw record based on the d11B difference between planktic and benthic foraminifera and an estimate of the pH gradient between surface and deep water. We then calculate d11Bsw two different ways. One variant of our method assumes that the pH gradient between surface and deep has remained the same as today over the past 23 Ma; the other uses the d13C gradient between surface and deep to represent change in the pH gradient through time. The results of these two methods of calculating d11Bsw are broadly consistency with each other, however, based on extensive carbon cycle modelling using CYCLOPS and GENIE we favour the d13C gradient method. In our favoured d11Bsw reconstruction, d11Bsw is around 2 per mil lower than today at ~37.5 per mil during the early and middle Miocene and increases to the modern value (39.61 per mil) by ~5 Ma. A similar pattern of change is evident in the seawater composition of three other stable isotope systems, Mg, Li and Ca. Concurrent shifts in the seawater isotopic composition of all four of these elements during the late Miocene, suggest a common forcing mechanism. We hypothesise the most likely cause of these shifts is a change in the isotopic composition of the riverine input, potentially driven by an increase in secondary mineral formation since ~15 Ma.

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Climate of the Past 13 (2017): 149-170, doi:10.5194/cp-13-149-2017.

Description: The boron isotope composition (δ11B) of foraminiferal calcite reflects the pH and the boron isotope composition of the seawater the foraminifer grew in. For pH reconstructions, the δ11B of seawater must therefore be known, but information on this parameter is limited. Here we reconstruct Neogene seawater δ11B based on the δ11B difference between paired measurements of planktic and benthic foraminifera and an estimate of the coeval water column pH gradient from their δ13C values. Carbon cycle model simulations underscore that the ΔpH–Δδ13C relationship is relatively insensitive to ocean and carbon cycle changes, validating our approach. Our reconstructions suggest that δ11Bsw was  ∼  37.5 ‰ during the early and middle Miocene (roughly 23–12 Ma) and rapidly increased during the late Miocene (between 12 and 5 Ma) towards the modern value of 39.61 ‰. Strikingly, this pattern is similar to the evolution of the seawater isotope composition of Mg, Li and Ca, suggesting a common forcing mechanism. Based on the observed direction of change, we hypothesize that an increase in secondary mineral formation during continental weathering affected the isotope composition of riverine input to the ocean since 14 Ma.

Description: The work was supported by NERC grants NE/I006176/1 (Gavin L. Foster and Caroline H. Lear), NE/H006273/1 (Gavin L. Foster), NE/I006168/1 and NE/K014137/1 and a Royal Society Research Merit Award (Paul A. Wilson), a NERC Independent Research Fellowship NE/K00901X/1 (Mathis P. Hain) and a NERC studentship (Rosanna Greenop).

Description: During the four most recent glacial maxima, atmospheric CO2 has been lowered by about 90–100 ppm with respect to interglacial concentrations. It is likely that most of the atmospheric CO2 deficit was stored in the ocean. Changes in the biological pump, which are related to the efficiency of the biological carbon uptake in the surface ocean and/or of the export of organic carbon to the deep ocean, have been proposed as a key mechanism for the increased glacial oceanic CO2 storage. The biological pump is strongly constrained by the amount of available surface nutrients. In models, it is generally assumed that the ratio between elemental nutrients, such as phosphorus, and carbon (C∕P ratio) in organic material is fixed according to the classical Redfield ratio. The constant Redfield ratio appears to approximately hold when averaged over basin scales, but observations document highly variable C∕P ratios on regional scales and between species. If the C∕P ratio increases when phosphate availability is scarce, as observations suggest, this has the potential to further increase glacial oceanic CO2 storage in response to changes in surface nutrient distributions. In the present study, we perform a sensitivity study to test how a phosphate-concentration-dependent C∕P ratio influences the oceanic CO2 storage in an Earth system model of intermediate complexity (cGENIE). We carry out simulations of glacial-like changes in albedo, radiative forcing, wind-forced circulation, remineralization depth of organic matter, and mineral dust deposition. Specifically, we compare model versions with the classical constant Redfield ratio and an observationally motivated variable C∕P ratio, in which the carbon uptake increases with decreasing phosphate concentration. While a flexible C∕P ratio does not impact the model's ability to simulate benthic δ13C patterns seen in observational data, our results indicate that, in production of organic matter, flexible C∕P can further increase the oceanic storage of CO2 in glacial model simulations. Past and future changes in the C∕P ratio thus have implications for correctly projecting changes in oceanic carbon storage in glacial-to-interglacial transitions as well as in the present context of increasing atmospheric CO2 concentrations.

Description: During the four most recent glacial cycles, atmospheric CO2 during glacial maxima has been lowered by about 90–100 ppm with respect to interglacials. There is widespread consensus that most of this carbon was partitioned in the ocean. It is however still debated which processes were dominant in achieving this increased carbon storage. In this paper, we use an Earth system model of intermediate complexity to constrain the range in ocean carbon storage for an ensemble of ocean circulation equilibrium states. We do a set of simulations where we run the model to pre-industrial equilibrium, but where we achieve different ocean circulation by changing forcing parameters such as wind stress, ocean diffusivity and atmospheric heat diffusivity. As a consequence, the ensemble members also have different ocean carbon reservoirs, global ocean average temperatures, biological pump efficiencies and conditions for air-sea CO2 disequilibrium. We analyse changes in total ocean carbon storage and separate it into contributions by the solubility pump, the biological pump and the CO2 disequilibrium component. We also relate these contributions to differences in strength of ocean overturning circulation. In cases with weaker circulation, we see that the ocean's capacity for carbon storage is larger. Depending on which ocean forcing parameter that is tuned, the origin of the change in carbon storage is different. When wind stress or ocean vertical diffusivity is changed, the response of the biological pump gives the most important effect on ocean carbon storage, whereas when atmospheric heat diffusivity or ocean horizontal diffusivity is changed, the solubility pump and the disequilibrium component are also important and sometimes dominant. Finally, we do a drawdown experiment, where we investigate the capacity for increased carbon storage by maximising the efficiency of the biological pump in our ensemble members. We conclude that different initial states for an ocean model result in different capacities for ocean carbon storage, due to differences in the ocean circulation state. This could explain why it is difficult to achieve comparable responses of the ocean carbon pumps in model intercomparison studies, where the initial states vary between models. The drawdown experiment highlights the importance of the strength of the biological pump in the control state for model studies of increased biological efficiency.

Description: The study of the mechanical energy budget of the oceans using the Lorenz available potential energy (APE) theory is based on knowledge of the adiabatically rearranged Lorenz reference state of minimum potential energy. The compressible and nonlinear character of the equation of state for seawater has been thought to cause the reference state to be ill defined, casting doubt on the usefulness of APE theory for investigating ocean energetics under realistic conditions. Using a method based on the volume frequency distribution of parcels as a function of temperature and salinity in the context of the seawater Boussinesq approximation, which is illustrated using climatological data, the authors show that compressibility effects are in fact minor. The reference state can be regarded as a well-defined one-dimensional function of depth, which forms a surface in temperature, salinity, and density space between the surface and the bottom of the ocean. For a very small proportion of water masses, this surface can be multivalued and water parcels can have up to two statically stable levels in the reference density profile, of which the shallowest is energetically more accessible. Classifying parcels from the surface to the bottom gives a different reference density profile than classifying in the opposite direction. However, this difference is negligible. This study shows that the reference state obtained by standard sorting methods is equivalent to, though computationally more expensive than, the volume frequency distribution approach. The approach that is presented can be applied systematically and in a computationally efficient manner to investigate the APE budget of the ocean circulation using models or climatological data.

Description: A decomposition of meridional overturning circulation (MOC) cells into geostrophic vertical shears, Ekman, and bottom pressure–dependent (or external mode) circulation components is presented. The decomposition requires the following information: 1) a density profile wherever bathymetry changes to construct the vertical shears component, 2) the zonal-mean zonal wind stress for the Ekman component, and 3) the mean depth-independent velocity information over each isobath to construct the external mode. The decomposition is applied to the third-generation Hadley Centre Coupled Ocean–Atmosphere General Circulation Model (HadCM3) to determine the meridional variability of these individual components within the Atlantic Ocean. The external mode component is shown to be extremely important where western boundary currents impinge on topography, and also in the area of the overflows. The Sverdrup balance explains the shape of the external mode MOC component to first order, but the time variability of the external mode exhibits only a very weak dependence on the wind stress curl. Thus, the Sverdrup balance cannot be used to determine the external mode changes when examining temporal change in the MOC. The vertical shears component allows the time-mean and the time-variable upper North Atlantic MOC cell to be deduced at 25°S and 50°N. A stronger dependency on the external mode and Ekman components between 8° and 35°N and in the regions of the overflows means that hydrographic sections need to be supplemented by bottom pressure and wind stress information at these latitudes. At the decadal time scale, variability in Ekman transport is less important than that in geostrophic shears. In the Southern Hemisphere the vertical shears component is dominant at all time scales, suggesting that hydrographic sections alone may be suitable for deducing change in the MOC at these latitudes.