The GFDL Global Ocean and Sea Ice Model OM4.0: Model Description and Simulation Features

• Published in: Journal of Advances in Modeling Earth Systems
• Main Authors: Adcroft, Alistair, Anderson, Whit, Balaji, V., Blanton, Chris, Bushuk, Mitchell, Dufour, Carolina O., Dunne, John P., Griffies, Stephen M., Hallberg, Robert, Harrison, Matthew J., Held, Isaac M., Jansen, Malte F., John, Jasmin G., Krasting, John P., Langenhorst, Amy R., Legg, Sonya, Liang, Zhi, McHugh, Colleen, Radhakrishnan, Aparna, Reichl, Brandon G., Rosati, Tony, Samuels, Bonita L., Shao, Andrew, Stouffer, Ronald, Winton, Michael, Wittenberg, Andrew T., Xiang, Baoqiang, Zadeh, Niki, Zhang, Rong
• Format: Article in Journal/Newspaper
• Language: English
• Published: Journal of Advances in Modeling Earth Systems 2019
• Subjects: ocean circulation model; CORE; hybrid coordinates; Sea ice
• Online Access: http://hdl.handle.net/1828/11899; https://doi.org/10.1029/2019MS001726

We document the configuration and emergent simulation features from the Geophysical Fluid Dynamics Laboratory (GFDL) OM4.0 ocean/sea ice model. OM4 serves as the ocean/sea ice component for the GFDL climate and Earth system models. It is also used for climate science research and is contributing to the Coupled Model Intercomparison Project version 6 Ocean Model Intercomparison Project. The ocean component of OM4 uses version 6 of the Modular Ocean Model and the sea ice component uses version 2 of the Sea Ice Simulator, which have identical horizontal grid layouts (Arakawa C‐grid). We follow the Coordinated Ocean‐sea ice Reference Experiments protocol to assess simulation quality across a broad suite of climate‐relevant features. We present results from two versions differing by horizontal grid spacing and physical parameterizations: OM4p5 has nominal 0.5° spacing and includes mesoscale eddy parameterizations and OM4p25 has nominal 0.25° spacing with no mesoscale eddy parameterization. Modular Ocean Model version 6 makes use of a vertical Lagrangian‐remap algorithm that enables general vertical coordinates. We show that use of a hybrid depth‐isopycnal coordinate reduces the middepth ocean warming drift commonly found in pure z * vertical coordinate ocean models. To test the need for the mesoscale eddy parameterization used in OM4p5, we examine the results from a simulation that removes the eddy parameterization. The water mass structure and model drift are physically degraded relative to OM4p5, thus supporting the key role for a mesoscale closure at this resolution. We thank the leadership of NOAA/GFDL (Director V. Ramaswamy and Associate Director W. Anderson) for their ongoing support during this model development project. We also thank the many GFDL scientists and engineers who provided valuable resources for algorithm development and code optimization, thus allowing OM4 to run efficiently and robustly across a suite of computational platforms. A. A., S. L. and B. G. R. were supported by Award NA14OAR4320106 ...