Processing your ADCP data using structure function techniques

From Atomix
Revision as of 09:49, 23 May 2022 by Yuengdjern (talk | contribs)

To calculate the dissipation rate at a specific range bin and a specific time ensemble:

Schematic showing along-beam distance r and radial velocities.
  1. Extract or compute the along-beam bin center separation [δr0] based on the instrument geometry
  2. Calculate the along-beam velocity fluctuation time-series in each bin n, where [Failed to parse (syntax error): {\displaystyle b’(n, t)} ] from the along-beam velocity data that has met the QC criteria (i.e. the data in Level 2 of the netcdf file)
  3. Select the maximum distance (rmax) over which to compute the structure function based on conditions of the flow (e.g., expected max overturn). The corresponding number of bins is [nrmax=rmax/r0]
  4. Calculate the structure function Dll for all possible bin separations δ within rmax using either a bin-centred difference scheme or a forward-difference scheme.
  5. Perform a regression of Dll(n,δ) against (δr0)2/3 for the appropriate range of bins and δr0 separation distances. Be aware of special considerations for forward-difference, center-difference schemes in setting up the regression calculation. The regression is typically done as a least-squares fit, either as:

    Dll=a0+a1(δr0)2/3; or as
    Dll=a0+a1(δr0)2/3+a3((δr0)2/3)3

    the former being the canonical method that excludes non-turbulent velocity differences between bins, whereas the latter is a modified method that includes non-turbulent velocity differences between bins due to any oscillatory signal (e.g. surface waves, motion of the ADCP on a mooring).
  6. Use the coefficient a1 to calculate ε as

    ε=(a1C2)2/3

    where C2 is an empirical constant, typically taken as 2.0 or 2.1.




Next step: Apply quality-control on dissipation rates (QA2)

Previous step: Apply quality-control on velocity time series data (QA1)

Return to ADCP Flow Chart front page