skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: Simulations With the Marine Biogeochemistry Library (MARBL)

The Marine Biogeochemistry Library (MARBL) is a prognostic ocean biogeochemistry model that simulates marine ecosystem dynamics and the coupled cycles of carbon, nitrogen, phosphorus, iron, silicon, and oxygen. MARBL is a component of the Community Earth System Model (CESM); it supports flexible ecosystem configuration of multiple phytoplankton and zooplankton functional types; it is also portable, designed to interface with multiple ocean circulation models. Here, we present scientific documentation of MARBL, describe its configuration in CESM2 experiments included in the Coupled Model Intercomparison Project version 6 (CMIP6), and evaluate its performance against a number of observational data sets. The model simulates present‐day air‐sea CO2flux and many aspects of the carbon cycle in good agreement with observations. However, the simulated integrated uptake of anthropogenic CO2is weak, which we link to poor thermocline ventilation, a feature evident in simulated chlorofluorocarbon distributions. This also contributes to larger‐than‐observed oxygen minimum zones. Moreover, radiocarbon distributions show that the simulated circulation in the deep North Pacific is extremely sluggish, yielding extensive oxygen depletion and nutrient trapping at depth. Surface macronutrient biases are generally positive at low latitudes and negative at high latitudes. CESM2 simulates globally integrated net primary production (NPP) of 48 Pg C yr−1and particulate export flux at 100 m of 7.1 Pg C yr−1. The impacts of climate change include an increase in globally integrated NPP, but substantial declines in the North Atlantic. Particulate export is projected to decline globally, attributable to decreasing export efficiency associated with changes in phytoplankton community composition.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Marine dissolved organic phosphorus (DOP) serves as an organic nutrient to marine autotrophs, sustaining a portion of annual net community production (ANCP). Numerical models of ocean circulation and biogeochemistry have diagnosed the magnitude of this process at regional to global scales but have thus far been validated against DOP observations concentrated within the Atlantic basin. Here we assimilate a new marine DOP data set with global coverage to optimize an inverse model of the ocean phosphorus cycle to investigate the regionally variable role of marine DOP utilization by autotrophs contributing to ANCP. We find ∼25% of ANCP accumulates as DOP with a regionally variable pattern ranging from 8% to 50% across nine biomes investigated. Estimated mean surface ocean DOP lifetimes of ∼0.5–2 years allow for transport of DOP from regions of net production to net consumption in subtropical gyres. Globally, DOP utilization by autotrophs sustains ∼14% (0.9 Pg C yr−1) of ANCP with regional contributions as large as ∼75% within the oligotrophic North Atlantic and North Pacific. Shallow export and remineralization of DOP within the ocean subtropics contributes ∼30%–80% of phosphate regeneration within the upper thermocline (<300 m). These shallow isopycnals beneath the subtropical gyres harboring the preponderance of remineralized DOP outcrop near the poleward edge of each gyre, which when combined with subsequent lateral transport equatorward by Ekman convergence, provide a shallow overturning loop retaining phosphorus within the subtropical biome, likely helping to sustain gyre ANCP over multiannual to decadal timescales.

    more » « less
  2. Abstract. Recent earth system models predict a 10 %–20 % decrease in particulate organic carbon export from the surface ocean by the end of the21st century due to global climate change. This decline is mainly caused by increased stratification of the upper ocean, resulting in reducedshallow subsurface nutrient concentrations and a slower supply of nutrients to the surface euphotic zone in low latitudes. These predictions,however, do not typically account for associated changes in remineralization depths driven by sinking-particle size. Here we combinesatellite-derived export and particle size maps with a simple 3-D global biogeochemical model that resolves dynamic particle size distributions toinvestigate how shifts in particle size may buffer or amplify predicted changes in surface nutrient supply and therefore export production. We showthat higher export rates are empirically correlated with larger sinking particles and presumably larger phytoplankton, particularly in tropical andsubtropical regions. Incorporating these empirical relationships into our global model shows that as circulation slows, a decrease in export isassociated with a shift towards smaller particles, which sink more slowly and are thus remineralized shallower. This shift towards shallowerremineralization in turn leads to greater recycling of nutrients in the upper water column and thus faster nutrient recirculation into the euphoticzone. The end result is a boost in productivity and export that counteracts the initial circulation-driven decreases. This negative feedbackmechanism (termed the particle-size–remineralization feedback) slows export decline over the next century by ∼ 14 % globally (from −0.29to −0.25 GtC yr−1) and by ∼ 20 % in the tropical and subtropical oceans, where export decreases are currently predicted tobe greatest. Our findings suggest that to more accurately predict changes in biological pump strength under a warming climate, earth system modelsshould include dynamic particle-size-dependent remineralization depths. 
    more » « less
  3. Abstract

    Among marine organisms, gelatinous zooplankton (GZ; cnidarians, ctenophores, and pelagic tunicates) are unique in their energetic efficiency, as the gelatinous body plan allows them to process and assimilate high proportions of oceanic carbon. Upon death, their body shape facilitates rapid sinking through the water column, resulting in carcass depositions on the seafloor (“jelly‐falls”). GZ are thought to be important components of the biological pump, but their overall contribution to global carbon fluxes remains unknown. Using a data‐driven, three‐dimensional, carbon cycle model resolved to a 1° global grid, with a Monte Carlo uncertainty analysis, we estimate that GZ consumed 7.9–13 Pg C y−1in phytoplankton and zooplankton, resulting in a net production of 3.9–5.8 Pg C y−1in the upper ocean (top 200 m), with the largest fluxes from pelagic tunicates. Non‐predation mortality (carcasses) comprised 25% of GZ production, and combined with the much greater fecal matter flux, total GZ particulate organic carbon (POC) export at 100 m was 1.6–5.2 Pg C y−1, equivalent to 32–40% of the global POC export. The fast sinking GZ export resulted in a high transfer efficiency (Teff) of 38–62% to 1,000 m and 25–40% to the seafloor. Finally, jelly‐falls at depths >50 m are likely unaccounted for in current POC flux estimates and could increase benthic POC flux by 8–35%. The significant magnitude of and distinct sinking properties of GZ fluxes support a critical yet underrecognized role of GZ carcasses and fecal matter to the biological pump and air‐sea carbon balance.

    more » « less
  4. Abstract

    Phytoplankton exhibit diverse physiological responses to temperature which influence their fitness in the environment and consequently alter their community structure. Here, we explored the sensitivity of phytoplankton community structure to thermal response parameterization in a modelled marine phytoplankton community. Using published empirical data, we evaluated the maximum thermal growth rates (μmax) and temperature coefficients (Q10; the rate at which growth scales with temperature) of six key Phytoplankton Functional Types (PFTs): coccolithophores, cyanobacteria, diatoms, diazotrophs, dinoflagellates, and green algae. Following three well‐documented methods, PFTs were either assumed to have (1) the sameμmaxand the sameQ10(as in to Eppley, 1972), (2) a uniqueμmaxbut the sameQ10(similar to Kremer et al., 2017), or (3) a uniqueμmaxand a uniqueQ10(following Anderson et al., 2021). These trait values were then implemented within the Massachusetts Institute of Technology biogeochemistry and ecosystem model (called Darwin) for each PFT under a control and climate change scenario. Our results suggest that applying aμmaxandQ10universally across PFTs (as in Eppley, 1972) leads to unrealistic phytoplankton communities, which lack diatoms globally. Additionally, we find that accounting for differences in theQ10between PFTs can significantly impact each PFT's competitive ability, especially at high latitudes, leading to altered modeled phytoplankton community structures in our control and climate change simulations. This then impacts estimates of biogeochemical processes, with, for example, estimates of export production varying by ~10% in the Southern Ocean depending on the parameterization. Our results indicate that the diversity of thermal response traits in phytoplankton not only shape community composition in the historical and future, warmer ocean, but that these traits have significant feedbacks on global biogeochemical cycles.

    more » « less
  5. Abstract

    The elemental ratios of carbon, nitrogen, and phosphorus (C:N:P) within organic matter play a key role in coupling biogeochemical cycles in the global ocean. At the cellular level, these ratios are controlled by physiological responses to the environment. But linking these cellular‐level processes to global biogeochemical cycles remains challenging. We present a novel model framework that combines knowledge of phytoplankton cellular functioning with global scale hydrographic data, to assess the role of variable carbon‐to‐phosphorus ratios (RC:P) on the distribution of export production. We implement a trait‐based mechanistic model of phytoplankton growth into a global biogeochemical inverse model to predict global patterns of phytoplankton physiology and stoichiometry that are consistent with both biological growth mechanisms and hydrographic carbon and nutrient observations. We compare this model to empirical parameterizations relatingRC:Pto temperature or phosphate concentration. We find that the way the model represents variable stoichiometry affects the magnitude and spatial pattern of carbon export, with globally integrated fluxes varying by up to 10% (1.3 Pg C yr−1) across models. Despite these differences, all models exhibit strong consistency with observed dissolved inorganic carbon and phosphate concentrations (R2 > 0.9), underscoring the challenge of selecting the most accurate model structure. We also find that the choice of parameterization impacts the capacity of changingRC:Pto buffer predicted export declines. Our novel framework offers a pathway by which additional biological information might be used to reduce the structural uncertainty in model representations of phytoplankton stoichiometry, potentially improving our capacity to project future changes.

    more » « less