Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: The Averaged Hydrostatic Boussinesq Ocean Equations in Generalized Vertical Coordinates
Due to their limited resolution, numerical ocean models need to be interpreted as representing filtered or averaged equations. How to interpret models in terms of formally averaged equations, however, is not always clear, particularly in the case of hybrid or generalized vertical coordinate models, which limits our ability to interpret the model results and to develop parameterizations for the unresolved eddy contributions. We here derive the averaged hydrostatic Boussinesq equations in generalized vertical coordinates for an arbitrary thickness‐weighted average. We then consider various special cases and discuss the extent to which the averaged equations are consistent with existing ocean model formulations. As previously discussed, the momentum equations in existing depth‐coordinate models are best interpreted as representing Eulerian averages (i.e., averages taken at fixed depth), while the tracer equations can be interpreted as either Eulerian or thickness‐weighted isopycnal averages. Instead we find that no averaging is fully consistent with existing formulations of the parameterizations in semi‐Lagrangian discretizations of generalized vertical coordinate ocean models such as MOM6. A coordinate‐following average would require “coordinate‐aware” parameterizations that can account for the changing nature of the eddy terms as the coordinate changes. Alternatively, the model variables can be interpreted as representing either Eulerian or (thickness‐weighted) isopycnal averages, independent of the model coordinate that is being used for the numerical discretization. Existing parameterizations in generalized vertical coordinate models, however, are not always consistent with either of these interpretations, which, respectively, would require a three‐dimensional divergence‐free eddy tracer advection or a form‐stress parameterization in the momentum equations.  more » « less
Award ID(s):
1912163 1912332
PAR ID:
10563672
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Volume:
16
Issue:
12
ISSN:
1942-2466
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; (, Journal of Advances in Modeling Earth Systems)
    Abstract 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. 
    more » « less
  2. ; ; (, Journal of Physical Oceanography)
    null (Ed.)
    Abstract Theories of the Beaufort Gyre (BG) dynamics commonly represent the halocline as a single layer with a thickness depending on the Eulerian-mean and eddy-induced overturning. However, observations suggest that the isopycnal slope increases with depth, and a theory to explain this profile remains outstanding. Here we develop a multilayer model of the BG, including the Eulerian-mean velocity, mesoscale eddy activity, diapycnal mixing, and lateral boundary fluxes, and use it to investigate the dynamics within the Pacific Winter Water (PWW) layer. Using theoretical considerations, observational data, and idealized simulations, we demonstrate that the eddy overturning is critical in explaining the observed vertical structure. In the absence of the eddy overturning, the Ekman pumping and the relatively weak vertical mixing would displace isopycnals in a nearly parallel fashion, contrary to observations. This study finds that the observed increase of the isopycnal slope with depth in the climatological state of the gyre is consistent with a Gent–McWilliams eddy diffusivity coefficient that decreases by at least 10%–40% over the PWW layer. We further show that the depth-dependent eddy diffusivity profile can explain the relative magnitude of the correlated isopycnal depth and layer thickness fluctuations on interannual time scales. Our inference that the eddy overturning generates the isopycnal layer thickness gradients is consistent with the parameterization of eddies via a Gent–McWilliams scheme but not potential vorticity diffusion. This study implies that using a depth-independent eddy diffusivity, as is commonly done in low-resolution ocean models, may contribute to misrepresentation of the interior BG dynamics. 
    more » « less
  3. ; ; ; (, Geoscientific Model Development)
    null (Ed.)
    Abstract. Eastern boundary upwelling systems (EBUSs) are physically and biologically active regions of the ocean with substantial impacts on ocean biogeochemistry, ecology, and global fish catch. Previous studies have used models of varying complexity to study EBUS dynamics, ranging from minimal two-dimensional (2-D) models to comprehensive regional and global models. An advantage of 2-D models is that they are more computationally efficient and easier to interpret than comprehensive regional models, but their key drawback is the lack of explicit representations of important three-dimensional processes that control biology in upwelling systems. These processes include eddy quenching of nutrients and meridional transport of nutrients and heat. The authors present the Meridionally Averaged Model of Eastern Boundary Upwelling Systems (MAMEBUS) that aims at combining the benefits of 2-D and 3-D approaches to modeling EBUSs by parameterizing the key 3-D processes in a 2-D framework. MAMEBUS couples the primitive equations for the physical state of the ocean with a nutrient–phytoplankton–zooplankton–detritus model of the ecosystem, solved in terrain-following coordinates. This article defines the equations that describe the tracer, momentum, and biological evolution, along with physical parameterizations of eddy advection, isopycnal mixing, and boundary layer mixing. It describes the details of the numerical schemes and their implementation in the model code, and provides a reference solution validated against observations from the California Current. The goal of MAMEBUS is to facilitate future studies to efficiently explore the wide space of physical and biogeochemical parameters that control the zonal variations in EBUSs. 
    more » « less
  4. (, Fluids)
    Ocean volume and tracer transports are commonly computed on density surfaces because doing so approximates the semi-Lagrangian mean advective transport. The resulting density-averaged transport can be related approximately to Eulerian-averaged quantities via the Temporal Residual Mean (TRM), valid in the limit of small isopycnal height fluctuations. This article builds on a formulation of the TRM for volume fluxes within Neutral Density surfaces, (the “NDTRM”), selected because Neutral Density surfaces are constructed to be as neutral as possible while still forming well-defined surfaces. This article derives a TRM, referred to as the “Neutral TRM” (NTRM), that approximates volume fluxes within surfaces whose vertical fluctuations are defined directly by the neutral relation. The purpose of the NTRM is to more closely approximate the semi-Lagrangian mean transport than the NDTRM, because the latter introduces errors associated with differences between the instantaneous state of the modeled/observed ocean and the reference climatology used to assign the Neutral Density variable. It is shown that the NDTRM collapses to the NTRM in the limiting case of a Neutral Density variable defined with reference to the Eulerian-mean salinity, potential temperature and pressure, rather than an external reference climatology, and therefore that the NTRM approximately advects this density variable. This prediction is verified directly using output from an idealized eddy-resolving numerical model. The NTRM therefore offers an efficient and accurate estimate of modeled semi-Lagrangian mean transports without reference to an external reference climatology, but requires that a Neutral Density variable be computed once from the model’s time-mean state in order to estimate isopycnal and diapycnal components of the transport. 
    more » « less
  5. ; ; ; ; ; (, Journal of Advances in Modeling Earth Systems)
    Abstract Ocean circulations are strongly constrained by rotation and stratification, making potential vorticity (PV) a key dynamical property. Specifically, PV is conserved following water parcels in the absence of diabatic forcing and friction—a desirable feature for isopycnal ocean models. Here, we devise a finite‐volume discretization for PV fluxes in the vector‐invariant formulation of shallow‐water equations. The scheme adopts a seventh‐order weighted essentially non‐oscillatory reconstruction of PV on the Arakawa C‐grid while ensuring PV conservation. To prevent PV singularities arising from vanishing layer thickness—typical in isopycnal ocean models—we propose a thickness weighting for the PV reconstruction. The scheme is implemented in the Modular Ocean Model 6 (MOM6) dynamical core and tested through idealized simulations of baroclinic geostrophic turbulence and wind‐driven double gyres. Compared with existing discretizations in MOM6 that conserve energy and/or potential enstrophy, the new scheme not only substantially reduces dispersive errors but also yields a weak potential enstrophy dissipation near grid scales, and thus permits stable time integration without superfluous horizontal viscosity. Additionally, the scheme supports a strong inverse energy cascade and exhibits an elevated eddy kinetic energy. These properties in turn lead to more realistic ocean energetics in eddy‐permitting simulations for which the energy input by surface wind stress is dissipated through bottom drag. In terms of computational efficiency, the scheme leads to notable improvement in energy spectra in eddy‐permitting simulations with a minor increase in total cost by . Results here are promising for improving the simulation of ocean eddies and large‐scale circulations for climate applications. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Text Encoded Conversion
  • XML
  • Save / Share this Record
  • Send to Email