Comparing ocean surface boundary vertical mixing schemes including langmuir turbulence

Six recent Langmuir turbulence parameterization schemes and five traditional schemes are implemented in a common single‐column modeling framework and consistently compared. These schemes are tested in scenarios versus matched large eddy simulations, across the globe with realistic forcing (JRA55‐do, WAVEWATCH‐III simulated waves) and initial conditions (Argo), and under realistic conditions as observed at ocean moorings. Traditional non‐Langmuir schemes systematically underpredict large eddy simulation vertical mixing under weak convective forcing, while Langmuir schemes vary in accuracy. Under global, realistic forcing Langmuir schemes produce 6% (−1% to 14% for 90% confidence) or 5.2 m (−0.2 m to 17.4 m for 90% confidence) deeper monthly mean mixed layer depths than their non‐Langmuir counterparts, with the greatest differences in extratropical regions, especially the Southern Ocean in austral summer. Discrepancies among Langmuir schemes are large (15% in mixed layer depth standard deviation over the mean): largest under wave‐driven turbulence with stabilizing buoyancy forcing, next largest under strongly wave‐driven conditions with weak buoyancy forcing, and agreeing during strong convective forcing. Non‐Langmuir schemes disagree with each other to a lesser extent, with a similar ordering. Langmuir discrepancies obscure a cross‐scheme estimate of the Langmuir effect magnitude under realistic forcing, highlighting limited understanding and numerical deficiencies. Maps of the regions and seasons where the greatest discrepancies occur are provided to guide further studies and observations

Challenges and Prospects in Ocean Circulation Models

We revisit the challenges and prospects for ocean circulation models following Griffies et al. (2010). Over the past decade, ocean circulation models evolved through improved understanding, numerics, spatial discretization, grid configurations, parameterizations, data assimilation, environmental monitoring, and process-level observations and modeling. Important large scale applications over the last decade are simulations of the Southern Ocean, the Meridional Overturning Circulation and its variability, and regional sea level change. Submesoscale variability is now routinely resolved in process models and permitted in a few global models, and submesoscale effects are parameterized in most global models. The scales where nonhydrostatic effects become important are beginning to be resolved in regional and process models. Coupling to sea ice, ice shelves, and high-resolution atmospheric models has stimulated new ideas and driven improvements in numerics. Observations have provided insight into turbulence and mixing around the globe and its consequences are assessed through perturbed physics models. Relatedly, parameterizations of the mixing and overturning processes in boundary layers and the ocean interior have improved. New diagnostics being used for evaluating models alongside present and novel observations are briefly referenced. The overall goal is summarizing new developments in ocean modeling, including: how new and existing observations can be used, what modeling challenges remain, and how simulations can be used to support observations

Challenges and Prospects in Ocean Circulation Models

We revisit the challenges and prospects for ocean circulation models following Griffies et al. (2010). Over the past decade, ocean circulation models evolved through improved understanding, numerics, spatial discretization, grid configurations, parameterizations, data assimilation, environmental monitoring, and process-level observations and modeling. Important large scale applications over the last decade are simulations of the Southern Ocean, the Meridional Overturning Circulation and its variability, and regional sea level change. Submesoscale variability is now routinely resolved in process models and permitted in a few global models, and submesoscale effects are parameterized in most global models. The scales where nonhydrostatic effects become important are beginning to be resolved in regional and process models. Coupling to sea ice, ice shelves, and high-resolution atmospheric models has stimulated new ideas and driven improvements in numerics. Observations have provided insight into turbulence and mixing around the globe and its consequences are assessed through perturbed physics models. Relatedly, parameterizations of the mixing and overturning processes in boundary layers and the ocean interior have improved. New diagnostics being used for evaluating models alongside present and novel observations are briefly referenced. The overall goal is summarizing new developments in ocean modeling, including: how new and existing observations can be used, what modeling challenges remain, and how simulations can be used to support observations.Peer reviewe

Evaluation of global ocean–sea-ice model simulations based on the experimental protocols of the Ocean Model Intercomparison Project phase 2 (OMIP-2)

We present a new framework for global ocean- sea-ice model simulations based on phase 2 of the Ocean Model Intercomparison Project (OMIP-2), making use of the surface dataset based on the Japanese 55-year atmospheric reanalysis for driving ocean-sea-ice models (JRA55-do).We motivate the use of OMIP-2 over the framework for the first phase of OMIP (OMIP-1), previously referred to as the Coordinated Ocean-ice Reference Experiments (COREs), via the evaluation of OMIP-1 and OMIP-2 simulations from 11 state-of-the-science global ocean-sea-ice models. In the present evaluation, multi-model ensemble means and spreads are calculated separately for the OMIP-1 and OMIP-2 simulations and overall performance is assessed considering metrics commonly used by ocean modelers. Both OMIP-1 and OMIP-2 multi-model ensemble ranges capture observations in more than 80% of the time and region for most metrics, with the multi-model ensemble spread greatly exceeding the difference between the means of the two datasets. Many features, including some climatologically relevant ocean circulation indices, are very similar between OMIP-1 and OMIP- 2 simulations, and yet we could also identify key qualitative improvements in transitioning from OMIP-1 to OMIP- 2. For example, the sea surface temperatures of the OMIP- 2 simulations reproduce the observed global warming during the 1980s and 1990s, as well as the warming slowdown in the 2000s and the more recent accelerated warming, which were absent in OMIP-1, noting that the last feature is part of the design of OMIP-2 because OMIP-1 forcing stopped in 2009. A negative bias in the sea-ice concentration in summer of both hemispheres in OMIP-1 is significantly reduced in OMIP-2. The overall reproducibility of both seasonal and interannual variations in sea surface temperature and sea surface height (dynamic sea level) is improved in OMIP-2. These improvements represent a new capability of the OMIP-2 framework for evaluating processlevel responses using simulation results. Regarding the sensitivity of individual models to the change in forcing, the models show well-ordered responses for the metrics that are directly forced, while they show less organized responses for those that require complex model adjustments. Many of the remaining common model biases may be attributed either to errors in representing important processes in ocean-sea-ice models, some of which are expected to be reduced by using finer horizontal and/or vertical resolutions, or to shared biases and limitations in the atmospheric forcing. In particular, further efforts are warranted to resolve remaining issues in OMIP-2 such as the warm bias in the upper layer, the mismatch between the observed and simulated variability of heat content and thermosteric sea level before 1990s, and the erroneous representation of deep and bottom water formations and circulations. We suggest that such problems can be resolved through collaboration between those developing models (including parameterizations) and forcing datasets. Overall, the present assessment justifies our recommendation that future model development and analysis studies use the OMIP-2 framework

Evaluation of global ocean–sea-ice model simulations based on the experimental protocols of the Ocean Model Intercomparison Project phase 2 (OMIP-2)

We present a new framework for global ocean–sea-ice model simulations based on phase 2 of the Ocean Model Intercomparison Project (OMIP-2), making use of the surface dataset based on the Japanese 55-year atmospheric reanalysis for driving ocean–sea-ice models (JRA55-do). We motivate the use of OMIP-2 over the framework for the first phase of OMIP (OMIP-1), previously referred to as the Coordinated Ocean–ice Reference Experiments (COREs), via the evaluation of OMIP-1 and OMIP-2 simulations from 11 state-of-the-science global ocean–sea-ice models. In the present evaluation, multi-model ensemble means and spreads are calculated separately for the OMIP-1 and OMIP-2 simulations and overall performance is assessed considering metrics commonly used by ocean modelers. Both OMIP-1 and OMIP-2 multi-model ensemble ranges capture observations in more than 80 % of the time and region for most metrics, with the multi-model ensemble spread greatly exceeding the difference between the means of the two datasets. Many features, including some climatologically relevant ocean circulation indices, are very similar between OMIP-1 and OMIP-2 simulations, and yet we could also identify key qualitative improvements in transitioning from OMIP-1 to OMIP-2. For example, the sea surface temperatures of the OMIP-2 simulations reproduce the observed global warming during the 1980s and 1990s, as well as the warming slowdown in the 2000s and the more recent accelerated warming, which were absent in OMIP-1, noting that the last feature is part of the design of OMIP-2 because OMIP-1 forcing stopped in 2009. A negative bias in the sea-ice concentration in summer of both hemispheres in OMIP-1 is significantly reduced in OMIP-2. The overall reproducibility of both seasonal and interannual variations in sea surface temperature and sea surface height (dynamic sea level) is improved in OMIP-2. These improvements represent a new capability of the OMIP-2 framework for evaluating process-level responses using simulation results. Regarding the sensitivity of individual models to the change in forcing, the models show well-ordered responses for the metrics that are directly forced, while they show less organized responses for those that require complex model adjustments. Many of the remaining common model biases may be attributed either to errors in representing important processes in ocean–sea-ice models, some of which are expected to be reduced by using finer horizontal and/or vertical resolutions, or to shared biases and limitations in the atmospheric forcing. In particular, further efforts are warranted to resolve remaining issues in OMIP-2 such as the warm bias in the upper layer, the mismatch between the observed and simulated variability of heat content and thermosteric sea level before 1990s, and the erroneous representation of deep and bottom water formations and circulations. We suggest that such problems can be resolved through collaboration between those developing models (including parameterizations) and forcing datasets. Overall, the present assessment justifies our recommendation that future model development and analysis studies use the OMIP-2 framework.This research has been supported by the Integrated Research Program for Advancing Climate Models (TOUGOU) of the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan (grant nos. JPMXD0717935457 and JPMXD0717935561), the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) (grant no. 274762653), the Helmholtz Climate Initiative REKLIM (Regional Climate Change) and European Union's Horizon 2020 Research & Innovation program (grant nos. 727862 and 800154), the Research Council of Norway (EVA (grant no. 229771) and INES (grant no. 270061)), the US National Science Foundation (NSF) (grant no. 1852977), the National Natural Science Foundation of China (grant nos. 41931183 and 41976026), NOAA's Science Collaboration Program and administered by UCAR's Cooperative Programs for the Advancement of Earth System Science (CPAESS) (grant nos. NA16NWS4620043 and NA18NWS4620043B), and NOAA (grant no. NA18OAR4320123).Peer ReviewedPostprint (published version

Parameterization of ocean eddies: potential vorticity mixing, energetics and Arnold's first stability theorem

A family of eddy closures is studied that flux potential vorticity down-gradient and solve an explicit budget for the eddy energy, following the approach developed by Eden and Greatbatch (2008, Ocean Modelling). The aim of this manuscript is to demonstrate that when energy relates naturally to the instability or stability of the flow as described by Arnold's first stability theorem. The resultant family of eddy closures therefore possesses some of the ingredients necessary to parameterize the gross effects of eddies in both forced-dissipative and freely-decaying turbulence. These ideas are illustrated through their application to idealized, barotropic wind-driven gyres in which the maximum eddy energy occurs within the viscous boundary layers and separated western boundary currents, and to freely-decaying turbulence in a closed barotropic basin in which inertial Fofonoff gyres emerge as the long-time solution. The result that these eddy closures preserve the relation between the growth or decay of eddy energy and the instability or stability of the flow provides further support for their use in ocean general circulation models

Simulating Water Residence Time in the Coastal Ocean: A Global Perspective

Exchanges between coastal and oceanic waters shape both coastal ecosystem processes and signatures that they impart on global biogeochemical cycles. The timescales of these exchanges, however, are poorly represented in current‐generation, coarse‐grid climate models. Here we provide a novel global perspective on coastal residence time (CRT) and its spatio‐temporal variability using a new age tracer implemented in global ocean models. Simulated CRTs range widely from several days in narrow boundary currents to multiple years on broader shelves and in semi‐enclosed seas, in agreement with available observations. Overall, CRT is better characterized in high‐resolution models (1/8° and 1/4°) than in the coarser (1° and 1/2°) versions. This is in large part because coastal and open ocean grid cells are more directly connected in coarse models, prone to erroneous coastal flushing and an underestimated CRT. Additionally, we find that geometric enclosure of a coastal system places an important constraint on CRT