skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: An Eulerian scheme for identifying fronts and vortices in quasi-balanced flows
Abstract The identification of vortices in a fluid flow is a dynamically interesting problem that has practical applications in oceanography due to the outsized role eddies play in water mass, heat, and tracer transport. Here a new Eulerian scheme is developed to detect both vortices and strongly strained fronts, which are both ubiquitous in the world ocean. The new scheme is conceptually linked to the well-known Okubo-Weiss parameter, but is extended to quasigeostrophic flows by recognizing the strong role played by vertical shear in ocean dynamics. Adapted from the λ 2 -criterion for vortex identification (Jeong and Hussain 1995), the scheme considers the curvature of the pressure field as the differentiator between vortical and strained flow structures, and it is shown that its underlying geometry also exhibits characteristics of quasigeostrophic flow. The uses and skill of the scheme are demonstrated using a high-resolution regional ocean simulation, and prospects for its use with observational products are discussed.  more » « less
Award ID(s):
1912420
PAR ID:
10352093
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of Physical Oceanography
ISSN:
0022-3670
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. An idealized storm scenario is examined in which a wind-generated inertial wave interacts with a turbulent baroclinic quasigeostrophic flow. The flow is initialized by spinning up an Eady model with a stratification profile based on observations. The storm is modeled as an initial value problem for a mixed layer confined, horizontally uniform inertial oscillation. The primordial inertial oscillation evolves according to the phase-averaged model of Young and Ben Jelloul. Waves feed back onto the flow by modifying the potential vorticity. In the first few days, refraction dominates and wave energy is attracted (repelled) by regions of negative (positive) vorticity. Wave energy is subsequently drained down into the interior ocean guided by anticyclonic vortices. This drainage halts as wave energy encounters weakening vorticity. After a week or two, wave energy accumulates at the bottom of negative vorticity features, that is, along filamentary structures at shallow depths and in larger anticyclonic vortices at greater depths. Wave feedback tends to weaken vortices and thus slows the penetration of waves into the ocean interior. This nonlinear effect, however, is weak even for vigorous storms. 
    more » « less
  2. Abstract Horizontal boundary layer roll vortices are a series of large-scale turbulent eddies that prevail in a hurricane’s boundary layer. In this paper, a one-way nested sub-kilometer-scale large-eddy simulation (LES) based on the Weather Research and Forecasting (WRF) Model was used to examine the impact of roll vortices on the evolution of Hurricane Harvey around its landfall from 0000 UTC 25 August to 1800 UTC 27 August 2017. The simulation results imply that the turbulence in the LES can be attributed mainly to roll vortices. With the representation of roll vortices, the LES provided a better simulation of hurricane wind vertical structure and precipitation. In contrast, the mesoscale simulation with the YSU PBL scheme overestimated the precipitation for the hurricane over the ocean. Further analysis indicates that the roll vortices introduced a positive vertical flux and thinner inflow layer, whereas a negative flux maintained the maximum tangential wind at around 400 m above ground. During hurricane landfall, the weak negative flux maintained the higher wind in the LES. The overestimated low-level vertical flux in the mesoscale simulation with the YSU scheme led to overestimated hurricane intensity over the ocean and accelerated the decay of the hurricane during landfall. Rainfall analysis reveals that the roll vortices led to a weak updraft and insufficient water vapor supply in the LES. For the simulation with the YSU scheme, the strong updraft combined with surplus water vapor eventually led to unrealistic heavy rainfall for the hurricane over the ocean. 
    more » « less
  3. We present an Eulerian vortex method based on the theory of flow maps to simulate the complex vortical motions of incompressible fluids. Central to our method is the novel incorporation of the flow-map transport equations forline elements, which, in combination with a bi-directional marching scheme for flow maps, enables the high-fidelity Eulerian advection of vorticity variables. The fundamental motivation is that, compared to impulsem, which has been recently bridged with flow maps to encouraging results, vorticityωpromises to be preferable for its numerical stability and physical interpretability. To realize the full potential of this novel formulation, we develop a new Poisson solving scheme for vorticity-to-velocity reconstruction that is both efficient and able to accurately handle the coupling near solid boundaries. We demonstrate the efficacy of our approach with a range of vortex simulation examples, including leapfrog vortices, vortex collisions, cavity flow, and the formation of complex vortical structures due to solid-fluid interactions. 
    more » « less
  4. Abstract Several new methods are proposed that can diagnose the interscale transfer (or spectral flux) of kinetic energy (KE) and other properties in oceanic and broader geophysical systems, using integrals of advective structure functions and Bessel functions (herein “Bessel methods”). The utility of the Bessel methods is evaluated using simulations of anisotropic flow within two-dimensional (2D), surface quasigeostrophic (SQG), and two-layer QG systems. The Bessel methods diagnose various spectral fluxes within all of these systems, even under strong anisotropy and complex dynamics (e.g., multiple cascaded variables, coincident and opposing spectral fluxes, and nonstationary systems). In 2D turbulence, the Bessel methods capture the inverse KE cascade at large scales and the downscale enstrophy cascade (and associated downscale energy flux) at small scales. In SQG turbulence, the Bessel methods capture the downscale buoyancy variance cascade and the coincident upscale wavenumber-dependent KE flux. In QG turbulence, the Bessel methods capture the upscale kinetic energy flux. It is shown that these Bessel methods can be applied to data with limited extent or resolution, provided the scales of interest are captured by the range of separation distances. The Bessel methods are shown to have several advantages over other flux-estimation methods, including the ability to diagnose downscale energy cascades and to identify sharp transition scales. Analogous Bessel methods are also discussed for third-order structure functions, along with some caveats due to boundary terms. Significance StatementBig ocean eddies play an important role in Earth’s energy cycle by moving energy to both larger and smaller scales, but it is difficult to measure these “eddy energy fluxes” from oceanic observations. We develop a new method to estimate eddy energy fluxes that utilizes spatial differences between pairs of points and can be applied to various ocean data. This new method accurately diagnoses key eddy energy flux properties, as we demonstrate using idealized numerical simulations of various large-scale ocean systems. 
    more » « less
  5. One-quarter of the world’s tropical cyclones (TCs) occur in the Indian Ocean (IO) basin.The mechanisms for TC initiation in the IO are varied, but one recently discovered process involves the flow around the steep topography of Sumatra.  When the low-level flow impinges on Sumatra, it is blocked and the flow splits under typical environmental stratification.  As a result, wake vortices commonly develop at northern and southern island tips of the island. For the case of easterly flow, these circulationssubsequently move downstream over the IO.  The wake vortices emanating from the island tips are counter-rotating, but since Sumatra straddles the equator, the circulations are cyclonic in both hemispheres and thus have the potential for TC development. Using data from2.5yearsof observations from DYNAMO and YOTC, it is found that approximately 25% of the TCsthat occurred overIO basin during that periodwere initiated by Sumatra-induced wake vortices.Additional analysis of vortex statistics for the period 2008-17 has found that vortex counts are highest near Madden-Julian Oscillation (MJO) phase 1 when low-level easterlies are strongest across southern Sumatra. A secondary peak in vortex formation occurs during MJO phase 4 when low-level westerlies exist near the equator west of Sumatra. The latter finding suggests that MJO-related, low-level westerly surges on the equator impinging on Sumatracontribute to an increase in wake vortex development.  Numerical simulations have shown that circulations farther upstream such aswestern Pacific remnant TCs and the Borneo vortex can influence the development of Sumatra wake vortices and their growth into TCs over the IO. 
    more » « less