Search for: All records

There are two distinct parameterizations for the restratification effect of mesoscale eddies: the Greatbatch and Lamb (1990, GL90,https://journals.ametsoc.org/view/journals/phoc/20/10/1520-0485_1990_020_1634_opvmom_2_0_co_2.xml?tab_body=abstract-display) parameterization, which mixes horizontal momentum in the vertical, and the Gent and McWilliams (1990, GM90,https://journals.ametsoc.org/view/journals/phoc/20/1/1520-0485_1990_020_0150_imiocm_2_0_co_2.xml) parameterization, which flattens isopycnals adiabatically. Even though these two parameterizations are effectively equivalent under the assumption of quasi‐geostrophy, GL90 has been used much less than GM90, and exclusively inz‐coordinate models. In this paper, we compare the GL90 and GM90 parameterizations in an idealized isopycnal coordinate model, both from a theoretical and practical perspective. From a theoretical perspective, GL90 is more attractive than GM90 for isopycnal coordinate models because GL90 provides an interpretation that is fully consistent with thickness‐weighted isopycnal averaging, while GM90 cannot be entirely reconciled with any fully isopycnal averaging framework. From a practical perspective, the GL90 and GM90 parameterizations lead to extremely similar energy levels, flow and vertical structure, even though their energetic pathways are very different. The striking resemblance between the GL90 and GM90 simulations persists from non‐eddying through eddy‐permitting resolution. We conclude that GL90 is a promising alternative to GM90 for isopycnal coordinate models, where it is more consistent with theory, computationally more efficient, easier to implement, and numerically more stable. Assessing the applicability of GL90 in realistic global ocean simulations with hybrid coordinate schemes should be a priority for future work.

 

Abstract. We describe an idealized primitive-equation model for studying mesoscale turbulence and leverage a hierarchy of grid resolutions to make eddy-resolving calculations on the finest grids more affordable.The model has intermediate complexity, incorporating basin-scale geometry with idealized Atlantic and Southern oceans and with non-uniform ocean depth to allow for mesoscale eddy interactions with topography.The model is perfectly adiabatic and spans the Equator and thus fills a gap between quasi-geostrophic models, which cannot span two hemispheres, and idealized general circulation models, which generally include diabatic processes and buoyancy forcing.We show that the model solution is approaching convergence in mean kinetic energy for the ocean mesoscale processes of interest and has a rich range of dynamics with circulation features that emerge only due to resolving mesoscale turbulence. 

We develop a parameterization for representing the effects of submesoscale symmetric instability (SI) in the ocean interior. SI may contribute to water mass modification and mesoscale energy dissipation in flow systems throughout the World Ocean. Dense gravity currents forced by surface buoyancy loss over shallow shelves are a particularly compelling test case, as they are characterized by density fronts and shears susceptible to a wide range of submesoscale instabilities. We present idealized experiments of Arctic shelf overflows employing the GFDL‐MOM6 inz* and isopycnal coordinates. At the highest resolutions, the dense flow undergoes geostrophic adjustment and forms bottom‐ and surface‐intensified jets. The density front along the topography combined with geostrophic shear initiates SI, leading to onset of secondary shear instability, dissipation of geostrophic energy, and turbulent mixing. We explore the impact of vertical coordinate, resolution, and parameterization of shear‐driven mixing on the representation of water mass transformation. We find that in isopycnal and low‐resolutionz* simulations, limited vertical resolution leads to inadequate representation of diapycnal mixing. This motivates our development of a parameterization for SI‐driven turbulence. The parameterization is based on identifying unstable regions through a balanced Richardson number criterion and slumping isopycnals toward a balanced state. The potential energy extracted from the large‐scale flow is assumed to correspond to the kinetic energy of SI which is dissipated through shear mixing. Parameterizing submesoscale instabilities by combining isopycnal slumping with diapycnal mixing becomes crucial as ocean models move toward resolving mesoscale eddies and fronts but not the submesoscale phenomena they host.

 

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 purez^*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.