An official website of the United States government
Abstract 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
Approximating Isoneutral Ocean Transport via the Temporal Residual Mean
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
- PAR ID:
- 10177876
- Date Published:
- Journal Name:
- Fluids
- Volume:
- 4
- Issue:
- 4
- ISSN:
- 2311-5521
- Page Range / eLocation ID:
- 179
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Flow past disperse solid particles or bubbles induces fluctuations in carrier fluid velocity, which correlate with temperature fluctuations in non-isothermal flows resulting in the pseudo-turbulent heat flux (PTHF). In the Eulerian-Eulerian (EE) two-fluid (TF) model, the transport of PTHF is shown to be an important contributor to the overall energy budget, and is modeled using a pseudo-turbulent thermal diffusivity (PTTD). The PTHF and PTTD were originally quantified using particle-resolved direct numerical simulation (PR-DNS) data, and correlations were developed over a range of solid volume fraction (0.1 ≤ 𝜀𝑠 ≤ 0.5) and mean slip Reynolds number (1 ≤ 𝑅𝑒𝑚 ≤ 100) for a Prandtl number of 0.7. However, the original PTTD correlation diverges to infinity as the solid volume fraction goes to zero, which is physically unrealistic. This singular behavior is problematic for EE TF simulations at particle material fronts where solid volume fraction values can fall below the lower limit of existing data (𝜀𝑠 =0.1) to zero in the pure carrier phase. In this work, additional PR-DNS data are reported for 𝜀𝑠 < 0.1, and improved correlations are developed for the PTHF and PTTD. The new PTTD correlation is non- singular, and both the PTHF and PTTD decay exponentially to zero as the solid volume fraction approaches zero, which is physically reasonable. This improves prediction of PTHF transport in dilute flow using EE TF heat transfer simulations.more » « less
-
Abstract Coastal physical processes are essential for the cross‐shore transport of meroplanktonic larvae to their benthic adult habitats. To investigate these processes, we released a swarm of novel, trackable, subsurface vehicles, the Mini‐Autonomous Underwater Explorers (M‐AUEs), which we programmed to mimic larval depth‐keeping behavior. The M‐AUE swarm measured a sudden net onshore transport of 30–70 m over 15–20 min, which we investigated in detail. Here, we describe a novel transport mechanism of depth‐keeping plankton revealed by these observations. In situ measurements and models showed that, as a weakly nonlinear internal wave propagated through the swarm, it deformed surface‐intensified, along‐isopycnal background velocities downward, accelerating depth‐keeping organisms onshore. These higher velocities increased both the depth‐keepers' residence time in the wave and total cross‐shore displacement, leading to wave‐induced transports twice those of fully Lagrangian organisms and four times those associated with the unperturbed background currents. Our analyses also show that integrating velocity time series from virtual larvae or mimics moving with the flow yields both larger and more accurate transport estimates than integrating velocity time series obtained at a point (Eulerian). The increased cross‐shore transport of organisms capable of vertical swimming in this wave/background‐current system is mathematically analogous to the increase in onshore transport associated with horizontal swimming in highly nonlinear internal waves. However, the mechanism described here requires much weaker swimming speeds (mm s−1vs. cm s−1) to achieve significant onshore transports, and meroplanktonic larvae only need to orient themselves vertically, not horizontally.more » « less
-
Abstract Realistic computational simulations in different oceanic basins reveal prevalent prograde mean flows (in the direction of topographic Rossby wave propagation along isobaths; aka topostrophy) on topographic slopes in the deep ocean, consistent with the barotropic theory of eddy-driven mean flows. Attention is focused on the western Mediterranean Sea with strong currents and steep topography. These prograde mean currents induce an opposing bottom drag stress and thus a turbulent boundary layer mean flow in the downhill direction, evidenced by a near-bottom negative mean vertical velocity. The slope-normal profile of diapycnal buoyancy mixing results in downslope mean advection near the bottom (a tendency to locally increase the mean buoyancy) and upslope buoyancy mixing (a tendency to decrease buoyancy) with associated buoyancy fluxes across the mean isopycnal surfaces (diapycnal downwelling). In the upper part of the boundary layer and nearby interior, the diapycnal turbulent buoyancy flux divergence reverses sign (diapycnal upwelling), with upward Eulerian mean buoyancy advection across isopycnal surfaces. These near-slope tendencies abate with further distance from the boundary. An along-isobath mean momentum balance shows an advective acceleration and a bottom-drag retardation of the prograde flow. The eddy buoyancy advection is significant near the slope, and the associated eddy potential energy conversion is negative, consistent with mean vertical shear flow generation for the eddies. This cross-isobath flow structure differs from previous proposals, and a new one-dimensional model is constructed for a topostrophic, stratified, slope bottom boundary layer. The broader issue of the return pathways of the global thermohaline circulation remains open, but the abyssal slope region is likely to play a dominant role.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/fluids4040179
