skip to main content


Title: BFM17 v1.0: a reduced biogeochemical flux model for upper-ocean biophysical simulations
Abstract. We present a newly developed upper-thermocline, open-ocean biogeochemical flux model that is complex and flexible enough to capture open-ocean ecosystem dynamics but reduced enough to incorporate into highly resolved numerical simulations and parameter optimization studies with limited additional computational cost. The model, which is derived from the full 56-state-variable Biogeochemical Flux Model (BFM56; Vichi et al., 2007), follows a biological and chemical functional group approach and allows for the development of critical non-Redfield nutrient ratios. Matter is expressed in units of carbon, nitrogen, and phosphate, following techniques used in more complex models. To reduce the overall computational cost and to focus on upper-thermocline, open-ocean, and non-iron-limited or non-silicate-limited conditions, the reduced model eliminates certain processes, such as benthic, silicate, and iron influences, and parameterizes others, such as the bacterial loop. The model explicitly tracks 17 state variables, divided into phytoplankton, zooplankton, dissolved organic matter, particulate organic matter, and nutrient groups. It is correspondingly called the Biogeochemical Flux Model 17 (BFM17). After describing BFM17, we couple it with the one-dimensional Princeton Ocean Model for validation using observational data from the Sargasso Sea. The results agree closely with observational data, giving correlations above 0.85, except for chlorophyll (0.63) and oxygen (0.37), as well as with corresponding results from BFM56, with correlations above 0.85, except for oxygen (0.56), including the ability to capture the subsurface chlorophyll maximum and bloom intensity. In comparison to previous models of similar size, BFM17 provides improved correlations between several model output fields and observational data, indicating that reproduction of in situ data can be achieved with a low number of variables, while maintaining the functional group approach. Notable additions to BFM17 over similar complexity models are the explicit tracking of dissolved oxygen, allowance for non-Redfield nutrient ratios, and both dissolved and particulate organic matter, all within the functional group framework.  more » « less
Award ID(s):
1924658 1924636 1258995
NSF-PAR ID:
10295655
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Geoscientific Model Development
Volume:
14
Issue:
5
ISSN:
1991-9603
Page Range / eLocation ID:
2419 to 2442
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Surface ocean marine dissolved organic matter (DOM) serves as an important reservoir of carbon (C), nitrogen (N), and phosphorus (P) in the global ocean, and is produced and consumed by both autotrophic and heterotrophic communities. While prior work has described distributions of dissolved organic carbon (DOC) and nitrogen (DON) concentrations, our understanding of DOC:DON:DOP stoichiometry in the global surface ocean has been limited by the availability of DOP concentration measurements. Here, we estimate mean surface ocean bulk and semi‐labile DOC:DON:DOP stoichiometry in biogeochemically and geographically defined regions using newly available marine DOM concentration databases. Global mean surface ocean bulk (C:N:P = 387:26:1) and semi‐labile (C:N:P = 179:20:1) DOM stoichiometries are higher than Redfield stoichiometry, with semi‐labile DOM stoichiometry similar to that of global mean surface ocean particulate organic matter (C:N:P = 160:21:1) reported in a recent compilation. DOM stoichiometry varies across ocean basins, ranging from 251:17:1 to 638:43:1 for bulk and 83:15:1 to 414:49:1 for semi‐labile DOM C:N:P, respectively. Surface ocean DOP concentration exhibits larger relative changes than DOC and DON, driving surface ocean gradients in DOC:DON:DOP stoichiometry. Inferred autotrophic consumption of DOP helps explain intra‐ and inter‐basin patterns of marine DOM C:N:P stoichiometry, with regional patterns of water column denitrification and iron supply influencing the biogeochemical conditions favoring DOP use as an organic nutrient. Specifically, surface ocean marine DOM exhibits increasingly P‐depleted stoichiometries from east to west in the Pacific and from south to north in the Atlantic, consistent with patterns of increasing P stress and alleviated iron stress.

     
    more » « less
  2. Abstract

    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
  3. Abstract. Coupled physical–biogeochemical models can fill thespatial and temporal gap in ocean carbon observations. Challenges ofapplying a coupled physical–biogeochemical model in the regional oceaninclude the reasonable prescription of carbon model boundary conditions,lack of in situ observations, and the oversimplification of certainbiogeochemical processes. In this study, we applied a coupledphysical–biogeochemical model (Regional Ocean Modelling System, ROMS) to theGulf of Mexico (GoM) and achieved an unprecedented 20-year high-resolution(5 km, 1/22∘) hindcast covering the period of 2000 to 2019. Thebiogeochemical model incorporated the dynamics of dissolved organic carbon(DOC) pools and the formation and dissolution of carbonate minerals. Thebiogeochemical boundaries were interpolated from NCAR's CESM2-WACCM-FV2solution after evaluating the performance of 17 GCMs in the GoM waters. Modeloutputs included carbon system variables of wide interest, such aspCO2, pH, aragonite saturation state (ΩArag), calcitesaturation state (ΩCalc), CO2 air–sea flux, and carbon burialrate. The model's robustness is evaluated via extensive model–datacomparison against buoys, remote-sensing-based machine learning (ML)products, and ship-based measurements. A reassessment of air–sea CO2flux with previous modeling and observational studies gives us confidencethat our model provides a robust and updated CO2 flux estimation, andNGoM is a stronger carbon sink than previously reported. Model resultsreveal that the GoM water has been experiencing a ∼ 0.0016 yr−1 decrease in surface pH over the past 2 decades, accompanied by a∼ 1.66 µatm yr−1 increase in sea surfacepCO2. The air–sea CO2 exchange estimation confirms in accordance with severalprevious models and ocean surface pCO2 observations that theriver-dominated northern GoM (NGoM) is a substantial carbon sink, and theopen GoM is a carbon source during summer and a carbon sink for the rest ofthe year. Sensitivity experiments are conducted to evaluate the impacts ofriver inputs and the global ocean via model boundaries. The NGoM carbonsystem is directly modified by the enormous carbon inputs (∼ 15.5 Tg C yr−1 DIC and ∼ 2.3 Tg C yr−1 DOC) from theMississippi–Atchafalaya River System (MARS). Additionally,nutrient-stimulated biological activities create a ∼ 105 timeshigher particulate organic matter burial rate in NGoM sediment than in thecase without river-delivered nutrients. The carbon system condition of theopen ocean is driven by inputs from the Caribbean Sea via the Yucatan Channeland is affected more by thermal effects than biological factors. 
    more » « less
  4. Abstract

    Dissolved iron (dFe) plays an important role in regulating marine productivity. In high nutrient, low chlorophyll regions (>33% of the global ocean), iron is the primary growth limiting nutrient, and elsewhere iron can regulate nitrogen fixation by diazotrophs. The link between iron availability and carbon export is strongly dependent on the phytoplankton iron quotas or cellular Fe:C ratios. This ratio varies by more than an order of magnitude in the open ocean and is positively correlated with ambient dFe concentrations in field observations. Representing Fe:C ratios within models is necessary to investigate how ocean carbon cycling will interact with perturbations to iron cycling in a changing climate. The Community Earth System Model ocean component was modified to simulate dynamic, group‐specific, phytoplankton Fe:C that varies as a function of ambient iron concentration. The simulated Fe:C ratios improve the representation of the spatial trends in the observed Fe:C ratios. The acclimation of phytoplankton Fe:C ratios dampens the biogeochemical response to varying atmospheric deposition of soluble iron, compared to a fixed Fe:C ratio. However, varying atmospheric soluble iron supply has first order impacts on global carbon and nitrogen fluxes and on nutrient limitation spatial patterns. Our results suggest that pyrogenic Fe is a significant dFe source that rivals mineral dust inputs in some regions. Changes in dust flux and iron combustion sources (anthropogenic and wildfires) will modify atmospheric Fe inputs in the future. Accounting for dynamic phytoplankton iron quotas is critical for understanding ocean biogeochemistry and projecting its response to variations in atmospheric deposition.

     
    more » « less
  5. Abstract

    In contrast to its productive coastal margins, the open-ocean Gulf of Mexico (GoM) is notable for highly stratified surface waters with extremely low nutrient and chlorophyll concentrations. Field campaigns in 2017 and 2018 identified low rates of turbulent mixing, which combined with oligotrophic nutrient conditions, give very low estimates for diffusive flux of nitrate into the euphotic zone (< 1 µmol N m−2d−1). Estimates of local N2-fixation are similarly low. In comparison, measured export rates of sinking particulate organic nitrogen (PON) from the euphotic zone are 2 – 3 orders of magnitude higher (i.e. 462 – 1144 µmol N m−2d−1). We reconcile these disparate findings with regional scale dynamics inferred independently from remote-sensing products and a regional biogeochemical model and find that laterally-sourced organic matter is sufficient to support >90% of open-ocean nitrogen export in the GoM. Results show that lateral transport needs to be closely considered in studies of biogeochemical balances, particularly for basins enclosed by productive coasts.

     
    more » « less