Description: Warm periods in Earth's history tend to cool more slowly than cool periods warm. Carbon cycle feedbacks play a major role in these dynamics, from the slower rate of recovery of ocean carbon export production, to the slower re- establishment of geosphere carbon reservoirs, relative to rates of loss. Here we explore one- differences in how the global ocean takes up and gives up heat and carbon in forced rapid warming and cooling climate scenarios. We force an intermediate- complexity earth system model using two atmospheric CO2 scenarios. A ramp-up (1% per year increase in atmospheric CO2 for 150 years) starts from an average global CO2 concentration of 285 ppm to represent warming of an icehouse climate. A ramp- down (1% per year decrease in atmospheric CO2 for 150 years) starts from an average global CO2 concentration of 1257 ppm to represent cooling of a greenhouse climate. Atmospheric CO2 is then held constant in each simulation and the model is integrated an additional 350 years. The ramp-down simulation shows a weaker response of surface air temperature to changes in radiative forcing relative to the ramp-up scenario. This weaker response is due to a relatively large and fast release of heat from the ocean to the atmosphere. This asymmetry in heat exchange in cooling and warming scenarios exists mainly because of differences in the response of the ocean circulation to forcing. In the ramp-up, increasing stratification and weakening of meridional overturning circulation slows ocean carbon and heat uptake. In the ramp-down, cooling accelerates meridional overturning and deepens vertical mixing, accelerating the release of carbon and heat stored at depth. Though idealized, our experiments offer insight into differences in ocean dynamics in icehouse and greenhouse climate transitions.

Description: Marine calcifiers as a plankton functional type (PFT) are a crucial part of the global carbon cycle, being responsible for much of the carbon export to the deep ocean entering via biological pathways. Deep ocean carbon export through calcifiers is controlled by physiological, ecological, and biogeochemical factors. This paper describes the implementation of a calcifying phytoplankton PFT in the University of Victoria Earth System Climate Model, version 2.9 (UVic ESCM), and mechanistic improvements to the representation of model carbon export (a full calcite tracer, carbonate chemistry dependent calcite dissolution rates, and a ballasting scheme). An iterative method for stabilizing and tuning the biogeochemistry is furthermore described. The UVic ESCM now fills a niche in Earth system modelling that was previously unoccupied in that it is relatively inexpensive to run, yet resolves the complete Earth system carbon cycle including prognostic calcium carbonate and a separate phytoplankton calcifier PFT. The model is now well suited to testing feedbacks between the carbonate and carbon cycles and the climate system as transient simulations. The modifications described improve the UVic ESCM's mechanistic realism without compromising performance with respect to observed carbon and nutrient fluxes. Primary production, export production, particulate organic carbon, and calcite fluxes all fall within independently observed estimates.

Description: Treatment of the underwater light field in ocean biogeochemical models has been attracting increasing interest, with some models moving towards more complex parameterisations. We conduct a simple sensitivity study of a typical, highly simplified parameterisation. In our study, we vary the phytoplankton light attenuation parameter over a range constrained by data during both pre-industrial equilibrated and future climate scenario RCP8.5. In equilibrium, lower light attenuation parameters (weaker self-shading) shift net primary production (NPP) towards the high latitudes, while higher values of light attenuation (stronger shelf-shading) shift NPP towards the low latitudes. Climate forcing magnifies this relationship through changes in the distribution of nutrients both within and between ocean regions. Where and how NPP responds to climate forcing can determine the magnitude and sign of global NPP trends in this high CO2 future scenario. Ocean oxygen is particularly sensitive to parameter choice. Under higher CO2 concentrations, two simulations establish a strong biogeochemical feedback between the Southern Ocean and low-latitude Pacific that highlights the potential for regional teleconnection. Our simulations serve as a reminder that shifts in fundamental properties (e.g. light attenuation by phytoplankton) over deep time have the potential to alter global biogeochemistry.

Description: Autotrophy is largely resource-limited in the modern ocean. Paleo evidence indicates this was not necessarily the case in warmer climates, and modern observations as well as standard metabolic theory suggest continued ocean warming could shift global ecology towards heterotrophy, thereby reducing autotrophic nutrient limitation. Such a shift would entail strong nutrient recycling in the upper ocean and high rates of net primary production (NPP), yet low carbon export to the deep ocean and sediments. We demonstrate transition towards such a state in the early 22nd century as a response to business-as-usual representative concentration pathway forcing (RCP8.5) in an intermediate complexity Earth system model in three configurations; with and without an explicit calcifier phytoplankton class and calcite ballast model. In all models nutrient regeneration in the near-surface becomes an increasingly important driver of primary production. The near-linear relationship between changes in NPP and global sea surface temperature (SST) found over the 21st century becomes exponential above a 2–4${\;}^{\circ }{\rm{C}}$ global mean SST change. This transition to a more heterotrophic ocean agrees roughly with metabolic theory.

Description: Conventional integration of Earth system and ocean models can accrue considerable computational expenses, particularly for marine biogeochemical applications. "Offline" numerical schemes in which only the biogeochemical tracers are time stepped and transported using a pre-computed circulation field can substantially reduce the burden and are thus an attractive alternative. One such scheme is the "transport matrix method" (TMM), which represents tracer transport as a sequence of sparse matrix–vector products that can be performed efficiently on distributed-memory computers. While the TMM has been used for a variety of geochemical and biogeochemical studies, to date the resulting solutions have not been comprehensively assessed against their "online" counterparts. Here, we present a detailed comparison of the two. It is based on simulations of the state-of-the-art biogeochemical sub-model embedded within the widely used coarse-resolution University of Victoria Earth System Climate Model (UVic ESCM). The default, non-linear advection scheme was first replaced with a linear, third-order upwind-biased advection scheme to satisfy the linearity requirement of the TMM. Transport matrices were extracted from an equilibrium run of the physical model and subsequently used to integrate the biogeochemical model offline to equilibrium. The identical biogeochemical model was also run online. Our simulations show that offline integration introduces some bias to biogeochemical quantities through the omission of the polar filtering used in UVic ESCM and in the offline application of time-dependent forcing fields, with high latitudes showing the largest differences with respect to the online model. Differences in other regions and in the seasonality of nutrients and phytoplankton distributions are found to be relatively minor, giving confidence that the TMM is a reliable tool for offline integration of complex biogeochemical models. Moreover, while UVic ESCM is a serial code, the TMM can be run on a parallel machine with no change to the underlying biogeochemical code, thus providing orders of magnitude speed-up over the online model.

Description: Phytoplankton calcifiers contribute to global carbon cycling through their dual formation of calcium carbonate and particulate organic carbon (POC). The carbonate might provide an efficient export pathway for the associated POC to the deep ocean, reducing the particles' exposure to biological degradation in the upper ocean and increasing the particle settling rate. Previous work has suggested ballasting of POC by carbonate might increase in a warming climate, in spite of increasing carbonate dissolution rates, because calcifiers benefit from the widespread nutrient limitation arising from stratification. We compare the biogeochemical responses of three models containing (1) a single mixed phytoplankton class, (2) additional explicit small phytoplankton and calcifiers, and (3) additional explicit small phytoplankton and calcifiers with a prognostic carbonate ballast model, to two rapid changes in atmospheric CO2. The first CO2 scenario represents a rapid (151-year) transition from a stable icehouse climate (283.9 ppm) into a greenhouse climate (1263 ppm); the second represents a symmetric rapid transition from a stable greenhouse climate into an icehouse climate. We identify a slope change in the global net primary production response with a transition point at about 3.5 ∘C global mean sea surface temperature change in all models, driven by a combination of physical and biological changes. We also find that in both warming and cooling scenarios, the application of a prognostic carbonate ballast model moderates changes in carbon export production, suboxic volume, and nitrate sources and sinks, reducing the long-term model response to about one-third that of the calcifier model without ballast. Explicit small phytoplankton and calcifiers, and carbonate ballasting, increase the physical separation of nitrate sources and sinks through a combination of phytoplankton competition and lengthened remineralization profile, resulting in a significantly higher global nitrate inventory in this model compared to the single phytoplankton type model (15 % and 32 % higher for icehouse and greenhouse climates). Higher nitrate inventory alleviates nitrate limitation, increasing phytoplankton sensitivity to changes in physical limitation factors (light and temperature). This larger sensitivity to physical forcing produces stronger shifts in ocean phosphate storage between icehouse and greenhouse climates. The greenhouse climate is found to hold phosphate and nitrate deeper in the ocean, despite a shorter remineralization length scale than the icehouse climate, because of the longer residence times of the deep water masses. We conclude the global biogeochemical impact of calcifiers extends beyond their role in global carbon cycling, and that the ecological composition of the global ocean can affect how ocean biogeochemistry responds to climate forcing.