Abstract. The consideration of marine biogeochemistry is essential for simulating the carbon cycle in an Earth system model. Here we present the implementation and evaluation of a marine biogeochemical model, Model of Oceanic Pelagic Stoichiometry (MOPS) in the Flexible Ocean and Climate Infrastructure (FOCI) climate model. FOCI-MOPS enables the simulation of marine biological processes, the marine carbon, nitrogen and oxygen cycles, air-sea gas exchange of CO2 and O2, and simulations with prescribed atmospheric CO2 or CO2 emissions. A series of experiments covering the historical period (1850–2014) were performed following the DECK (Diagnostic, Evaluation and Characterization of Klima) and CMIP6 (Coupled Model Intercomparison Project 6) protocols. Overall, modelled biogeochemical tracer distributions and fluxes, as well as transient evolution in surface air temperature, air-sea CO2 fluxes, and changes of ocean carbon and heat, are in good agreement with observations. Modelled inorganic and organic tracer distributions are quantitatively evaluated by statistically-derived metrics. Results of the FOCI-MOPS model, also including sea surface temperature, surface pH, oxygen (100–600 m), nitrate (0–100 m), and primary production, are within the range of other CMIP6 model results. Overall, the evaluation of FOCI-MOPS indicates its suitability for Earth climate system simulations.