Search for: All records

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. 

In inertial-range turbulence, structure functions can diagnose transfer or dissipation rates of energy and enstrophy, which are difficult to calculate directly in flows with complex geometry or sparse sampling. However, existing relations between third-order structure functions and these rates only apply under isotropic conditions. We propose new relations to diagnose energy and enstrophy dissipation rates in anisotropic two-dimensional (2-D) turbulence. These relations use second-order advective structure functions that depend on spatial increments of vorticity, velocity, and their advection. Numerical simulations of forced-dissipative anisotropic 2-D turbulence are used to compare new and existing relations against model-diagnosed dissipation rates of energy and enstrophy. These simulations permit a dual cascade where forcing is applied at an intermediate scale, energy is dissipated at large scales, and enstrophy is dissipated at small scales. New relations to estimate energy and enstrophy dissipation rates show improvement over existing methods through increased accuracy, insensitivity to sampling direction, and lower temporal and spatial variability. These benefits of advective structure functions are present under weakly anisotropic conditions, and increase with the flow anisotropy as third-order structure functions become increasingly inappropriate. Several of the structure functions also show promise for diagnosing the forcing scale of 2-D turbulence. Velocity-based advective structure functions show particular promise as they can diagnose both enstrophy and energy cascade rates, and are robust to changes in the effective resolution of local derivatives. Some existing and future datasets that are amenable to advective structure function analysis are discussed.