TY - JOUR
T1 - Comparing Two Parameterizations for the Restratification Effect of Mesoscale Eddies in an Isopycnal Ocean Model
AU - Loose, Nora
AU - Marques, Gustavo M.
AU - Adcroft, Alistair
AU - Bachman, Scott
AU - Griffies, Stephen M.
AU - Grooms, Ian
AU - Hallberg, Robert W.
AU - Jansen, Malte F.
N1 - Publisher Copyright:
© 2023 The Authors. Journal of Advances in Modeling Earth Systems published by Wiley Periodicals LLC on behalf of American Geophysical Union.
PY - 2023/12
Y1 - 2023/12
N2 - 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 in z-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.
AB - 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 in z-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.
KW - isopycnal coordinates
KW - mesoscale eddies
KW - ocean model
KW - parameterization
KW - thickness-weighted average
KW - vertical viscosity
UR - http://www.scopus.com/inward/record.url?scp=85179183859&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85179183859&partnerID=8YFLogxK
U2 - 10.1029/2022MS003518
DO - 10.1029/2022MS003518
M3 - Article
AN - SCOPUS:85179183859
SN - 1942-2466
VL - 15
JO - Journal of Advances in Modeling Earth Systems
JF - Journal of Advances in Modeling Earth Systems
IS - 12
M1 - e2022MS003518
ER -