Abstract. We present a newly developed reduced-order 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 with limited additional computational cost. The reduced-order 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 open-ocean 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 providing a detailed description of BFM17, we couple it with the one-dimensional Princeton Ocean Model (POM) for validation using observational data from the Sargasso Sea. Results show good agreement with the observational data and with corresponding results from BFM56, including the ability to capture the subsurface chlorophyll maximum and bloom intensity. In comparison to previous reduced-order models of similar size, BFM17 provides improved correlations between model output and field data, indicating that significant improvements in the reproduction of in situ data can be achieved with a low number of variables, while maintaining the functional group approach.

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.