ALBERT

Description: This study focuses on an improved representation of the biological soft tissue pump in the global three-dimensional biogeochemical ocean model PISCES. We compare three parameterizations of particle dynamics: (1) the model standard version including two particle size classes, aggregation-disaggregation and prescribed sinking speed; (2) an aggregation-disaggregation model with a particle size spectrum and prognostic sinking speed; (3) a mineral ballast parameterization with no size classes, but prognostic sinking speed. In addition, the model includes a description of surface sediments and organic carbon early diagenesis. Model output is compared to data or data based estimates of ocean productivity, pe-ratios, particle fluxes, surface sediment bulk composition and benthic O2 fluxes. Model results suggest that different processes control POC fluxes at different depths. In the wind mixed layer turbulent particle coagulation appears as key process in controlling pe-ratios. Parameterization (2) yields simulated pe-ratios that compare well to observations. Below the wind mixed layer, POC fluxes are most sensitive to the intensity of zooplankton flux feeding, indicating the importance of zooplankton community composition. All model parameters being kept constant, the capability of the model to reproduce yearly mean POC fluxes below 2000 m and benthic oxygen demand does at first order not dependent on the resolution of the particle size spectrum. Aggregate formation appears essential to initiate an intense biological pump. At great depth the reported close to constant particle fluxes are most likely the result of the combined effect of aggregate formation and mineral ballasting.

Description: PISCES-v2 is a biogeochemical model which simulates the lower trophic levels of marine ecosystem (phytoplankton, microzooplankton and mesozooplankton) and the biogeochemical cycles of carbon and of the main nutrients (P, N, Fe, and Si). The model is intended to be used for both regional and global configurations at high or low spatial resolutions as well as for short-term (seasonal, interannual) and long-term (climate change, paleoceanography) analyses. There are twenty-four prognostic variables (tracers) including two phytoplankton compartments (diatoms and nanophytoplankton), two zooplankton size-classes (microzooplankton and mesozooplankton) and a description of the carbonate chemistry. Formulations in PISCES-v2 are based on a mixed Monod–Quota formalism: on one hand, stoichiometry of C/N/P is fixed and growth rate of phytoplankton is limited by the external availability in N, P and Si. On the other hand, the iron and silicium quotas are variable and growth rate of phytoplankton is limited by the internal availability in Fe. Various parameterizations can be activated in PISCES-v2, setting for instance the complexity of iron chemistry or the description of particulate organic materials. So far, PISCES-v2 has been coupled to the NEMO and ROMS systems. A full description of PISCES-v2 and of its optional functionalities is provided here. The results of a quasi-steady state simulation are presented and evaluated against diverse observational and satellite-derived data. Finally, some of the new functionalities of PISCES-v2 are tested in a series of sensitivity experiments.

Description: This study focuses on an improved representation of the biological soft tissue pump in the global three-dimensional biogeochemical ocean model PISCES. We compare three parameterizations of particle dynamics: (1) the model standard version including two particle size classes, aggregation-disaggregation and prescribed sinking speed; (2) an aggregation-disaggregation model with a particle size spectrum and prognostic sinking speed; (3) a mineral ballast parameterization with no size classes, but prognostic sinking speed. In addition, the model includes a description of surface sediments and organic carbon early diagenesis. The integrated representation of material fluxes from the productive surface ocean down to the sediment-water interface allows taking advantage of surface ocean observations, sediment trap data and exchange fluxes at the sediment-water interface. The capability of the model to reproduce yearly averaged particulate organic carbon fluxes and benthic oxygen demand does at first order not dependent on the resolution of the particle size spectrum. Model results obtained with the standard version and with the one including a particle size spectrum and prognostic sinking speed are not significantly different. Both model versions overestimate particulate organic carbon between 1000 and 2000 m, while deep fluxes are of the correct order of magnitude. Predicted benthic oxygen fluxes correspond with respect to their large scale distribution and magnitude to data based estimates. Modeled particulate organic C fluxes across the mesopelagos are most sensitive to the intensity of zooplankton flux feeding. An increase of the intensity of flux feeding in the standard version results in lower mid- and deep-water particulate organic carbon fluxes, shifting model results to an underestimation of particulate organic carbon fluxes in the deep. The corresponding benthic oxygen fluxes are too low. The model version including the mineral ballast parameterization yields an improved fit between modeled and observed particulate organic carbon fluxes below 2000 m and down to the sediment-water interface. Our results suggest that aggregate formation alone might not be sufficient to drive an intense biological pump. The later is most likely driven by the combined effect of aggregate formation and mineral ballasting.

Description: Planktonic Foraminifera are a major contributor to the deep carbonate flux and their microfossil deposits form one of the richest databases for reconstructing paleoenvironments, particularly through changes in their taxonomic and shell composition. Using an empirically based planktonic foraminifer model that incorporates three known major physiological drivers of their biogeography – temperature, food and light – we investigate (i) the global redistribution of planktonic Foraminifera under anthropogenic climate change and (ii) the alteration of the carbonate chemistry of foraminiferal habitat with ocean acidification. The present-day and future (2090–2100) 3-D distributions of Foraminifera are simulated using temperature, plankton biomass and light from an Earth system model forced with a historical and a future (IPCC A2) high CO2 emission scenario. Foraminiferal abundance and diversity are projected to decrease in the tropics and subpolar regions and increase in the subtropics and around the poles. Temperature is the dominant control on the future change in the biogeography of Foraminifera. Yet food availability acts to either reinforce or counteract the temperature-driven changes. In the tropics and subtropics the largely temperature-driven shift to depth is enhanced by the increased concentration of phytoplankton at depth. In the higher latitudes the food-driven response partly offsets the temperature-driven reduction both in the subsurface and across large geographical regions. The large-scale rearrangements in foraminiferal abundance and the reduction in the carbonate ion concentrations in the habitat range of planktonic foraminifers – from 10–30 μmol kg−1 in their polar and subpolar habitats to 30–70 μmol kg−1 in their subtropical and tropical habitats – would be expected to lead to changes in the marine carbonate flux. High-latitude species are most vulnerable to anthropogenic change: their abundance and available habitat decrease and up to 10% of the volume of their habitat drops below the calcite saturation horizon.

Description: The ocean is a substantial source of nitrous oxide (N2O) to the atmosphere, but little is known about how this flux might change in the future. Here, we investigate the potential evolution of marine N2O emissions in the 21st century in response to anthropogenic climate change using the global ocean biogeochemical model NEMO-PISCES. Assuming nitrification as the dominant N2O formation pathway, we implemented two different parameterizations of N2O production which differ primarily under low-oxygen (O2) conditions. When forced with output from a climate model simulation run under the business-as-usual high-CO2 concentration scenario (RCP8.5), our simulations suggest a decrease of 4 to 12 % in N2O emissions from 2005 to 2100, i.e., a reduction from 4.03/3.71 to 3.54/3.56 TgN yr−1 depending on the parameterization. The emissions decrease strongly in the western basins of the Pacific and Atlantic oceans, while they tend to increase above the oxygen minimum zones (OMZs), i.e., in the eastern tropical Pacific and in the northern Indian Ocean. The reduction in N2O emissions is caused on the one hand by weakened nitrification as a consequence of reduced primary and export production, and on the other hand by stronger vertical stratification, which reduces the transport of N2O from the ocean interior to the ocean surface. The higher emissions over the OMZ are linked to an expansion of these zones under global warming, which leads to increased N2O production, associated primarily with denitrification. While there are many uncertainties in the relative contribution and changes in the N2O production pathways, the increasing storage seems unequivocal and determines largely the decrease in N2O emissions in the future. From the perspective of a global climate system, the averaged feedback strength associated with the projected decrease in oceanic N2O emissions amounts to around −0.009 W m−2 K−1, which is comparable to the potential increase from terrestrial N2O sources. However, the assessment for a potential balance between the terrestrial and marine feedbacks calls for an improved representation of N2O production terms in fully coupled next-generation Earth system models.

Description: This study compares spatial and temporal variability in net primary productivity (PP) and particulate organic carbon (POC) export production (EP) from three different coupled climate carbon cycle models (IPSL, MPIM, NCAR) with observation-based estimates derived from satellite measurements of ocean colour and inverse modelling. Satellite observations of ocean colour have shown that temporal variability of PP on the global scale is largely dominated by the permanently stratified, low-latitude ocean (Behrenfeld et al., 2006) ocite{Behrenfeld06} with stronger stratification (higher SSTs) leading to negative PP anomalies and vice versa. Results from all three coupled models confirm the role of the low-latitude, permanently stratified ocean for global PP anomalies. Two of the models also reproduce the inverse relationship between stratification (SST) and PP, especially in the equatorial Pacific. With the help of the model results we are able to explain the chain of cause and effect leading from stratification (SST) through nutrient concentrations to PP and finally to EP. There are significant uncertainties in observational PP and especially EP. Our finding of a good agreement between independent estimates from coupled models and satellite observations provides increased confidence that such models can be used as a first basis to estimate the impact of future climate change on marine productivity and carbon export.

Description: This model study addresses the change in pelagic calcium carbonate production (CaCO3, as calcite in the model) and dissolution in response to rising atmospheric CO2. The parameterization of CaCO3 production includes a dependency on the saturation state of seawater with respect to calcite. It was derived from laboratory and mesocosm studies on particulate organic and inorganic carbon production in Emiliania huxleyi as a function of pCO2. The model predicts values of CaCO3 production and dissolution in line with recent estimates. The effect of rising pCO2 on CaCO3 production and dissolution was quantified by means of model simulations forced with atmospheric CO2 increasing at a rate of 1% per year from 286 ppm to 1144 ppm. The simulation predicts a decrease of CaCO3 production by 27%. The combined change in production and dissolution of CaCO3 yields an excess uptake of CO2 from the atmosphere by the ocean of 5.9 GtC.

Description: The marine aragonite cycle has been included in the global biogeochemical model PISCES to study the role of aragonite in shallow water CaCO3 dissolution. Aragonite production is parameterized as a function of mesozooplankton biomass and aragonite saturation state of ambient waters. Observation-based estimates of marine carbonate production and dissolution are well reproduced by the model and about 60% of the combined CaCO3 water column dissolution from aragonite and calcite is simulated above 2000 m. In contrast, a calcite-only version yields a much smaller fraction. This suggests that the aragonite cycle should be included in models for a realistic representation of CaCO3 dissolution and alkalinity. For the SRES A2 CO2 scenario, production rates of aragonite are projected to notably decrease after 2050. By the end of this century, global aragonite production is reduced by almost one third and total CaCO3 production by 19% relative to pre-industrial. Geographically, the effect from increasing atmospheric CO2, and the subsequent reduction in saturation state, is largest in the subpolar and polar areas where the modeled aragonite production is projected to decrease by 65% until 2100.