Model output and code from a simulation of global change to 2100 for a marine ecosystem model
This version of the model was runs with the gud package (https://gitlab.com/jahn/gud).
The model includes the same ecosystem model and parameters as Dutkiewicz et al (2021), but here is used within a coarser physical ocean configuration that simulates a future warmer ocean.
The ecosystem component resolves the cycling of carbon, phosphorus, nitrogen silica, iron, and oxygen through inorganic, living, dissolved and particulate organic phases (including CDOM). The biogeochemical and biological tracers are transported and mixed by the MIT general circulation model (MITgcm, Marshall et al., 1997) driven by an earth system of intermediate complexity (MESM, Dutkiewicz et al, 2005; Moneir et al 2013; Dutkiewicz et al., 2019) . The three dimensional ocean has coarse resolution (2o×2.5o horizontally) and 22 levels ranging from 10m in the surface to 500m at depth.
We resolve 35 phytoplankton types, covering several functional groups (picophytoplankton, coccolithophores, diazotrophs, diatoms, mixotrophic dinoflagellates), and 16 size classes (from 0.6 to 228um equivalent spherical diameter). There are 16 grazer size classes from 6.6 to 2280um equivalent spherical diameter. The phytoplankton types differ in the types of nutrients they require (e.g. diatoms require silica), maximum growth rate, nutrient half saturation constants, sinking rates, and palatability to grazers.

Note that the atmospheric-ocean component of the earth system model used to produce the forcing fields has its own internal interannual variability. Though the frequency of phenomena such as ENSO are captured, they do not occur in same years as the real ocean

This particular version of the repository has model output used in Henson et al (Nature Communications, 2021), Cael et (Science Advances, 2021), and Bardon et al (Geophysical Research Letters, in press)

References:

Bardon, L., B.A. Ward, S. Dutkiewicz, and B.B. Cael. Testing the skill of a species distribution model in a virtual environment. Geophysical Research Letters, in press.

Cael, B.B., S. Dutkiewicz and S.A. Henson, 2021. Abrupt shifts in 21st century plankton communities. Science Advances, 7, eabf8593.

Dutkiewicz, S., P. Boyd, U. Riebesell 2021. Exploring biogeochemical and ecological redundancy in phytoplankton communities in the global ocean. Global Change Biology, 27, 1196-1213, https://doi.org/10.1111/gcb.15493

Dutkiewicz, S., A.E. Hickman, O. Jahn, E. Moneir, S. Henson, and C. Beaulieu, 2019. Ocean colour signature of climate change. Nature Communications, 10, doi:10.1038/s41467-019-08457-x.

Dutkiewicz. S., Sokolov, A., Scott, J. & Stone, P. A Three-Dimensional Ocean-Seaice-Carbon Cycle Model and its Coupling to a Two-Dimensional Atmospheric Model: Uses in Climate Change Studies. Report 122, Joint Program of the Science and Policy of Global Change, M.I.T., Cambridge, MA (2005)

Henson, S.A., S.R. Allen, B.B. Cael, and S. Dutkiewicz (2021). Future phytoplankton diversity in a changing climate. Nature Communications,12, 5372, doi:10.1038/s41467-021-25699-w .

Marshall, J., Adcroft, A., Hill, C. N., Perelman, L., and Heisey, C: A finite-volume, incompressible Navier–Stokes model for studies of the ocean on parallel computers, J. Geophys. Res., 102, 5753–5766, 1997.

Monier E., Scott, J.R., Sokolov, A.P., Forest, C.E., Schlosser, C.A. An integrated assessment modeling framework for uncertainty studies in global and regional climate change: the MIT IGSM-CAM (version 1.0). Geosci. Model Dev., 6, 2063-2085, doi:10.5194/gmd-6-2063-2013 (2013).