Turbulent diapycnal fluxes as a pilot Essential Ocean Variable

1 Introduction

In a quiescent fluid, the transfer of a tracer (e.g., dissolved gas) concentration only occurs via the slow process of molecular diffusion. In the ocean, three-dimensional turbulence significantly accelerates molecular diffusion by increasing the surface area between density layers. Turbulence is a non-linear eddying state of motion that cascades energy towards progressively smaller scales (see Frisch, 1995). It drives subsurface fluxes across density layers (i.e., diapycnal fluxes), contributing to the spatial redistribution of water properties.

Three-dimensional turbulent events usually occur at small spatial (<10 m) and temporal (seconds to minutes) scales and promote the redistribution of different oceanic waters, leading to the irreversible transformation of water masses. Although the processes driving these fluxes only affect scales of a few meters or less, the irreversible transformation associated with these fluxes impacts the large-scale ocean circulation at time scales ranging from seconds to decades (Moum, 2021).

Turbulent mixing is a key mechanism in the global overturning circulation (Munk, 1966; Wunsch and Ferrari, 2004) because it provides the mechanical energy necessary to upwell deep waters across isopycnals – the densest water masses at the bottom of the ocean gain buoyancy by mixing with lighter water above, providing a pathway by which water can return to the ocean surface after sinking at the poles. Mixing also sets the distributions of dissolved gases, nutrients and pollutants, which impact the Earth’s climate system and the global carbon cycle (Ellison et al., 2023), as well as the productivity of ecosystems (Bindoff et al., 2019; Melet et al., 2022). Diapycnal turbulent fluxes play an important role in emerging ocean industries such as deep-sea mining and marine Carbon Dioxide Removal (mCDR), expanding the need for observations. For example, deep ocean turbulence controls the scale of environmental impact of sediment plume deposition in the wake of deep-sea mining (Peacock and Ouillon, 2023) as well as the rate and permanence of carbon sequestration (National Academies of Sciences, Engineering, and Medicine, 2021).

Ocean mixing occurs at scales much smaller than the discretized grid cells of ocean and climate models, so it must be parameterized. It is common for model developers to assess the sensitivity of ocean circulation model performance to a particular turbulence closure scheme (e.g. Adcroft et al., 2019); more rare, however, is the direct comparison of the parameterized turbulence quantities with observed ones (e.g. Luneva et al., 2019; de Lavergne et al., 2020; Savelyev et al., 2022; Trossman et al., 2022). Since available datasets are limited to certain regions and times of the year, the parameterizations have regional and temporal biases. Additionally, the assumptions underlying the parameterizations often do not hold where mixing is known to be especially active (e.g., near the seafloor and active sub-mesoscale areas, see Section 5). In some cases, models are specifically tuned to reflect the higher mixing in these regions (Heuzé et al., 2015). Thus, these tunings require measurements in the first place, and it is unclear how they will hold as the ocean and climate system change. The inaccuracies associated with mixing parameterizations make them one of the major uncertainties in climate models (Hazeleger and Haarsma, 2005; Melet and Meyssignac, 2015; Exarchou et al., 2018; Zhu and Zhang, 2019; Deppenmeier et al., 2020); improving climate predictions thus requires more routine measurements of turbulence fluxes for parameterization development and validation.

Given the importance of mixing on climate projections and the protection of ocean health (see Section 2), we propose that turbulent diapycnal fluxes (the drivers of ocean mixing) become a pilot Essential Ocean Variable (EOV), as defined by the Global Ocean Observing System (GOOS). Ocean mixing measurement technology has matured significantly over the last two decades. Many commercial offerings now exist to measure turbulence, and the measurement platforms have expanded beyond ship-based profilers and bottom-landers. Autonomous platforms such as gliders, self-propelled vehicles, moorings, drifting profilers, and Argo floats now routinely collect datasets spanning weeks to several months.

Data processing and standards for archiving of turbulence measurements have also matured. In late 2020, the Scientific Committee on Oceanic Research (SCOR) approved Working Group #160 on “Analysing ocean Turbulence Observations to quantify MIXing” (ATOMIX) to develop best practices and quality control procedures for mixing data. The group advocates for “Turbulent Diapycnal Fluxes” to be considered an Essential Ocean Variable. Below, in section 2, we provide more context about the rationale of the proposed subvariables and derived variables outlined in Table 1.

Table 1 Ocean Turbulent Mixing variable and its sub-variables.

Overall, a “Turbulent Diapycnal Fluxes” EOV would allow GOOS to stimulate and coordinate ocean mixing science and technology to work toward significant improvements in ocean forecasts, climate projections, and the protection of ocean health (see Section 2). This document aims to present subsurface turbulent fluxes as a natural addition to the current EOV list because they are both crucial and feasible to measure globally. The following sections describe: 1) the need for sustained ocean mixing observations; 2) the methods for obtaining turbulent fluxes measurements, and the many turbulence sensors and platforms that are currently operational or under development; 3) the current and future management practices of turbulence data; 4) the feasibility and cost of creating large networks of turbulence-sensing instruments; 5) the international coordination required for supporting these networks.

2 Impact of turbulent fluxes on global climate, ocean health, and operational ocean services

By increasing the effective diffusivity across water mass boundaries, turbulent fluxes play a direct role in warming the surface ocean, producing the steric effects that drive sea level rise, and impacting coastal (within 100km) communities where ∼40% of the world population lives. The deep ocean is also strongly influenced by subsurface turbulent fluxes. Of the excess heat stored in the climate system (ocean and atmosphere), the ocean takes up 90%, with about 29% and 9% stored in the 700–2000m and below-2000m layers, respectively (Cheng et al., 2021).

At small scales, turbulent mixing accelerates the heat and momentum transfers between neighboring water masses. It disperses any particles or solutes they may contain, affecting local chemical and biological processes. At the largest scales, ocean mixing is a mechanical driver of the global ocean overturning circulation, upwelling the dense water that sinks at the poles (Figure 1). Added together, mixing events flux buoyancy downwards into abyssal water masses, allowing them to upwell across deep isopycnals and close the global overturning circulation (Talley, 2013). In this way, the upwelling of deep water is balanced by the downward mixing of buoyancy, which raises the potential energy stored in the water column (Munk, 1966; Wunsch and Ferrari, 2004; Rahmstorf, 2006). In an energetic sense, mixing is essential to the meridional overturning circulation. Through the downward mixing of heat, ocean mixing also increases the heat stored in the deep ocean, which has important impacts on Earth’s climate and its associated ecosystems and economies. The following sections provide brief examples of the impact of subsurface turbulent fluxes on a diverse range of processes.

Figure 1 Turbulent fluxes play a first-order role in the global overturning circulation system. The widespread zones of mixing-driven upwelling (red circles) balance the formation of North Atlantic Deep Water (formed at L and N locations) and Antarctic Bottom Water (formed at R and W locations). Figure copied from Kuhlbrodt et al. (2007).

2.1 Turbulent fluxes control Earth’s mean climate

Global climate models are too coarsely-spaced to resolve turbulent mixing processes, such that turbulent buoyancy fluxes and corresponding diffusivity coefficients must be parameterized. The substantial effect of imposing imprecise diffusivities on climate predictions was illustrated in a recent coupled global climate model analysis in which global average diapycnal diffusivity was varied between two realistic values (i.e., 0.9×10⁻⁴ m² s⁻¹ and 1.7×10⁻⁴ m² s⁻¹, Hieronymus et al., 2019). The difference between the least and most diffusive runs leads to a 3.6°C difference in volume-mean ocean temperature, a 2.4°C difference in sea surface temperature, and a 3°C difference in global-mean air temperature 2m above the ocean surface. In models too coarse to resolve the oceanic mesoscale eddy field (such as comprehensive Earth System Models), isopycnal diffusivities must also be parameterized (Redi, 1982) and these have a leading-order effect on the ventilation and structure of oceanic water masses (Jones and Abernathey, 2019).

2.2 Turbulent fluxes draw down anthropogenic heat and carbon

Ocean mixing influences the efficiencies of heat and carbon uptake from the atmosphere to the ocean through its impact on the stratification (e.g., Tatebe et al., 2018). Enhanced diapycnal mixing erodes the stratification in the upper ocean leading to a weaker and more diffuse thermocline, which enables more downward heat transport (e.g., Kuhlbrodt and Gregory, 2012; Melet et al., 2016) and more carbon uptake (Schmittner et al., 2009; Ehlert et al., 2017).

Using an eddy-resolving ocean model that assimilates ocean observations, Ellison et al. (2023) compared the effects of imposing two different background diffusivities. They showed that by altering the prescribed mixing rates from 10⁻⁴ m² s⁻¹ to 10⁻⁵ m² s⁻¹ (a common range of diapycnal diffusivities used in models), they could create a 40% change in Southern Ocean air–sea fluxes in only a few years. The different diffusivities led to an altered distribution of dissolved inorganic carbon, alkalinity, temperature and salinity, all of which affect the surface flux of CO₂ (Figure 2). Furthermore, direct turbulence observations – in regions where outgassing occurs in the subpolar Southern Ocean – showed strong episodic CO₂ outgassing events driven by storms (Nicholson et al., 2022).

Figure 2 Results from Ellison et al. (2023) showing how applying a fixed versus variable vertical mixing diffusivity to the Southern Ocean State Estimate impacts air-sea carbon fluxes on 1 to 6 years timescales. The Biogeochemical Southern Ocean State Estimate (BSOSE) (Verdy and Mazloff, 2017) highlights the importance of vertical mixing for air-sea fluxes. BSOSE includes carbon, oxygen, and nutrient cycles and is constrained with observations from profiling floats, shipboard data, underway measurements, and satellites, while maintaining closed budgets and obeying dynamical and thermodynamic balances. The ‘standard’ BSOSE employs a constant background diffusivity of 10⁻⁴ m² s⁻¹ along with a parameterization for mixed-layer dynamics. BSOSE with variable mixing is called BSOSE_mix. (A) Annual mean carbon flux for BSOSE (positive = ocean outgassing, negative = carbon uptake). Magenta lines show annual maximum and minimum sea-ice extents. (B) Annual mean difference between BSOSE_mix and BSOSE; positive = reduced carbon uptake or increased outgassing. (C) Cumulative sum of carbon fluxes from 70°S northward to 30°S for BSOSE (red) and BSOSE_mix (blue).

A remarkable emergent property of Earth system models is that, due to compensation between the effects of anthropogenic heat and carbon uptake, the peak global-mean surface warming is proportional to cumulative anthropogenic carbon emissions (Matthews et al., 2009). While this result underlies global climate mitigation policy (Drake and Henderson, 2022), the all-important constant of proportionality— known as the Transient Climate Response (TCR)—remains frustratingly uncertain (Matthews et al., 2020), with about 50% of the geophysical uncertainty attributable to oceanic tracer uptake processes (Lutsko and Popp, 2019) and the remainder due to radiative feedbacks. Using a simple climate-economic model, Hope (2015) estimates that the benefit of halving uncertainty in the TCR is valued at $10 trillion (USD), suggesting that the potential societal value of ocean mixing research could easily be in the many billions of dollars due to climate considerations alone.

These results illustrate the fundamental role that ocean mixing has in driving ocean warming, sea level rise, and ocean acidification, directly impacting coastal communities where 40% of the world’s population lives and indirectly impacting the rest of the global population through changes in regional climate.

2.3 Turbulent fluxes maintain healthy ecosystems

Stratification inhibits the vertical exchanges of nutrients and dissolved gases. Turbulent mixing processes supply nutrients to the biologically productive upper ocean, which helps regulate the net primary productivity of the open ocean and coastal environments. Turbulent vertical exchanges may also be the main source of dissolved oxygen replenishment in the deeper layers of the water column. In some environments, the lack of mixing precludes dissolved oxygen replenishment at depth (Bourgault et al., 2012), leading to hypoxic conditions that are detrimental to the ecosystem. Turbulence measurements are required regularly to document the physical processes that enhance these exchanges before they can be adequately linked with biogeochemical processes over seasonal timescales (e.g., Rippeth et al., 2009; Bluteau et al., 2021) and over time as the climate changes.

Mixing has an especially important role in maintaining the concentrations of dissolved oxygen, nitrate, and other nutrients in coastal and estuarine environments. A lack of mixing in the deepest parts of inlets and estuaries causes them to become hypoxic (Bourgault et al., 2012), while elevated mixing at sills may generate biodiversity hotspots by re-oxygenating bottom waters and transporting nutrients upwards into the euphotic (well-lit) surface layers. Estuarine circulations and tidal flow over topography cause mixing that brings nitrate up into the surface layer. The relative strength of this type of tidally-induced mixing varies depending on other biogeochemical processes at play, like the seasonal strength of riverine nitrate influx to the system, but it can be as important as river discharge (Bluteau et al., 2021).

One of the robust oceanic signals of climate change is an increase in upper ocean stratification (Bindoff et al., 2019; Melet et al., 2022), both from increased thermal stratification and increased salinity stratification in high latitudes. Projections suggest that the net primary productivity of the ocean will very likely decrease by between 4% and 11% by 2100 under an unmitigated climate change scenario, mostly because the enhanced upper ocean stratification reduces the supply of nutrients via turbulent mixing across the pycnocline (Bindoff et al., 2019; Melet et al., 2022).

2.4 Turbulent fluxes govern air–sea interactions

The upper boundary layer of the ocean and the processes occurring in the surface mixed layer figure prominently in the global climate system because it mediates momentum, heat, and gas fluxes between the ocean and the atmosphere (Frankignoul and Hasselmann, 1977; Bopp et al., 2015; Pörtner et al., 2019). Ultimately, these fluxes are governed by the dynamics of turbulent vertical mixing in this boundary layer (D’Asaro, 2014).

As an example, the sea surface temperature (SST) is a critical control on the atmosphere (Xie, 2004) and impacts equatorial climate by controlling precipitation (Xie et al., 2015). Over the African Sahel, precipitation changes with the SST difference between the neighboring subtropical North Atlantic and global tropics (Giannini et al., 2013). In the tropics, the SST rises by 2°C during boreal spring, when heating of the upper ocean by the atmosphere exceeds cooling by mixing from below. In boreal summer, SST decreases because cooling from below exceeds heating from above (Moum et al., 2013). Such a quantitative assessment of how mixing (computed from χ measurements) varies on timescales longer than a few weeks clearly shows its controlling influence on the seasonal cooling of SST in a critical oceanic regime. This case study is but one example of the indirect impacts of ocean mixing through coupled Earth system interactions and teleconnections.

2.5 Turbulent fluxes factor into natural resource assessments

In the coming decades, there will be greatly increased attention on anthropogenic use of the ocean. One example is deep-sea mining of abundant critical mineral resources (e.g., nickel, cobalt) in the abyssal ocean (Peacock and Alford, 2018). Such activities will generate benthic plumes of sediment, and other biogeochemical factors (e.g., dissolved heavy metals), that have the potential to negatively impact the ocean environment (Peacock and Ouillon, 2023). Recent modeling studies identified turbulent vertical mixing as a critical parameter that influences model predictions of plume extent, which will be the basis of decision-making by regulators (Chen et al., 2023).

While deep-sea mining is projected to be a multi-billion dollar industry, marine Carbon Dioxide Removal (mCDR) is projected to become an even larger trillion-dollar industry (National Academies of Sciences, Engineering, and Medicine, 2021). Several different mCDR technologies are being proposed (e.g., alkalinity enhancement, artificial upwelling, seaweed cultivation). For all proposed technologies, turbulent vertical mixing in the vicinity of the upper mixed layer is the critical governing physical process by which enhanced uptake of atmospheric carbon dioxide occurs, setting the efficacy of any mCDR technology. Critically, this nascent global industry will rely on model predictions informed by field data to identify potentially promising mCDR locations and to create the so-called MRV (Measurement, Reporting & Verification) value chain that will be the basis for a global market of mCDR carbon credits. Without correct parameterization of turbulent vertical mixing in mCDR modeling systems, there will be no way to achieve confidence in the mCDR MRV process – or to reliably assess the potentially negative impacts of proposed mCDR approaches. Global standards for turbulent diapycnal fluxes as an EOV will be essential.

Another technology that utilizes the ocean and is growing rapidly is offshore wind farms. There is, for example, a major expansion of activities in the United States (Musial et al., 2021), that lags behind Europe in this space (Fernández, 2023), which itself is scaling up activities in the North Sea. The introduction of fixed installations in shallow waters and floating structures in offshore waters creates strong turbulence in the wake of the installations (Dorrell et al., 2022). These enhanced turbulence levels can impact the surrounding marine environment and this needs to be well characterized in order to inform environmental assessments (Dorrell et al., 2022).

3 Definition of the proposed essential ocean variable

3.1 Selection of Essential Ocean Variables

The Global Ocean Observing System (GOOS) – a program co-sponsored by four organizations: three United Nations agencies, Intergovernmental Oceanographic Commission of UNESCO, the World Meteorological Organization, United Nations Environment Programme, along with a non-governmental organization, the International Council of Science – has the mission to support an international community of ocean observing organizations and set global standards for maintaining sustained observations and outputting data. Since 2011, GOOS has implemented a Framework for Ocean Observing that serves as a road map for supporting the ocean observing system (UNESCO, 2012). The Framework includes a list of Essential Ocean Variables (EOVs), which are measurable parameters considered vital to inform the three GOOS delivery areas: climate, ocean health, and forecasts and warnings. The selection of EOVs relies on two criteria: the impact (e.g., scientific, ocean services and health) and the feasibility (e.g., technological, political, economical) of making sustained measurements.

Once an EOV is selected, GOOS facilitates and coordinates the sustained operation of observation programs of the EOVs at global scales. It achieves this by encouraging funding from international or national agencies for observation programs and technological developments that would improve the sustainability of EOV measurements. The current list of EOVs contains variables linked to ocean circulation and the distribution and transport of heat, salt, and other water properties (Table 2). These include temperature, salinity, ocean currents, and sea-surface fluxes–variables that are continuously sampled through the global ocean today.

Table 2 Current list of Essential Ocean Variables.

3.2 Turbulent diapycnal fluxes as an Essential Ocean Variable

Our proposed variables (see Table 1) can be derived from the turbulent dissipation rate of kinetic energy (ϵ) and the turbulent dissipation rate of temperature (or another scalar) variance (χ), which are more easily measurable than direct covariance estimates from co-located velocity and tracer variance measurements (i.e., “eddy correlation” methods). Turbulent dissipation rates are obtained from measuring perturbations – ranging from a few centimeters to a few meters – of velocities and scalars such as temperature, conductivity (salinity), or dissolved oxygen. These variables are already part of the EOV list (Table 2), and ϵ is mentioned as a derived variable of Subsurface Currents in its EOV Specification Sheet. We argue, however, that GOOS should consider establishing an EOV specific to turbulent diapycnal fluxes for the following reasons:

1. Although turbulent diapycnal flux variables are related to existing EOVs (e.g., subsurface temperature and currents), they require different sampling protocols and instrument setups than traditional observations. Namely, selecting the time and spatial resolution depends on the turbulence theories invoked when deriving ϵ and/or χ from the observations, in addition to anticipated turbulence “energy” levels. The ultimate use for the observations also influences the choice of technique and sensors.

2. Turbulence variables depend on quadratic moments of existing EOVs at small spatial scales and fast time scales, so their calculation does not commute with averaging operations; crucially, this means that turbulence variables cannot simply be derived from existing EOV measurements, which sample at much lower frequencies and/or wavenumbers than needed to resolve J, ϵ, or χ (see Figure 3).

3. The EOV list should reflect the purpose of flux (J) or dissipation rate measurements (ϵ and χ) since they enable describing ocean dynamics beyond what can be documented from non-turbulent velocity and scalars. For instance, nutrient availability in the euphotic zone may be sporadically intense or input at a slower but more steady rate from turbulent exchanges. The magnitude of these turbulent exchanges, relative to other sources and sinks, cannot be gleaned solely by regularly monitoring currents and nutrients. Turbulent vertical fluxes may be the primary source for providing nutrients in the surface layer in many regions and they must be measured purposefully.

4. Global and sustained turbulent mixing observations are required to improve ocean mixing parameterizations in climate models. Ocean mixing is one of the three major uncertainties in climate models for sea level rise, along with ice sheet and cloud feedbacks (Melet and Meyssignac, 2015). Global climate models have begun using flow-dependent parameterizations for diffusivity, which require more direct diffusivity measurements to constrain them and assess the validity of their underlying assumptions.

Figure 3 From Bluteau et al. (2016b). Example spectral observations of shear spanning from ε ∼ 10⁻¹⁰ to ε ∼ 10⁻⁴ W kg⁻¹. The colored lines correspond to a theoretical shear spectra (e.g., Nasmyth spectrum) matching with the observations. Colored squares are placed at k_10, i.e., where 10% of the variance in the viscous range is resolved, while colored circles are placed at k₂₅. The thick lines delineate the wavenumber range that needs to be resolved to apply the integration method. The variance missing from the observation is corrected using the theoretical spectra. The bracket below the x-axis shows that the range of wavenumbers typically resolved by CTD and ADCP measurements is lower than the range required to resolve turbulent dissipation rates. A similar approach is used for the computation of χ with observed wavenumber spectra of temperature gradient matching the Batchelor spectrum. The reader should refer to the original study for the other features of the figure.

4 Measuring turbulent fluxes ($\overline{{\mathrm{U}}^{\prime }{C}^{\prime }}$) and turbulent diffusivity (K)

For turbulent diapycnal fluxes or mixing to occur, a source of energy must overcome the stable background stratification (e.g., through shear or convective instability) and create a growing overturning cell. A forward energy cascade occurs as the instability becomes non-linear and transitions to turbulence. As the cascade reaches scales where molecular processes become leading order, momentum, heat, salt, and passive tracers are diffused and get mixed between neighboring water masses, progressively lowering their total variance. In statistical equilibrium, the mean tracer variance $\overline{{(C^{\prime })}^2}$ is (on average) produced by fluctuating tracer fluxes acting upon mean gradients $\nabla \overline{C}$ at the same rate that it is dissipated by molecular diffusion:

(1)

where κ_m is the tracer’s molecular diffusivity and overbars denote averages. For the sake of simplicity, and mirroring the vast majority of previous ocean mixing studies, we have assumed here that the divergent transport terms can be neglected.

Turbulence is inherently a 3D process, but, at larger scales, it is constrained along isopycnals (approximately horizontal) by the overwhelming influences of rotation and stratification; only the relatively small, fast, and isotropic motions contribute to diapycnal (approximately vertical) fluxes. Thus, it is quite common to focus on the vertical component of Equation 1 which simplifies even further to

Here, w^′ and C^′ represent the microscale perturbations of vertical velocity and tracer concentration. In this formulation, the balance is between the average advective flux (i.e., turbulent diapycnal flux) of a tracer $(J_c\equiv \overline{w^{\prime }{C}^{\prime }})$ and the molecular dissipation of that tracer.

Focusing on the left hand side of equation 2, a bulk estimate of the turbulent diapycnal fluxes can be determined by invoking a Fick’s law approach, in which the flux is proportional to the concentration gradient of the tracer. The turbulent fluxes producing a tracer’s variance can be parameterized by the introduction of an effective turbulent diapycnal diffusivity K_C that shortcuts molecular dissipation by acting directly on mean tracer gradients, i.e.

The units of diffusivity are m²s⁻¹; it can be thought of as the rate at which the size of the tracer patch expands due to diffusion (Taylor, 1922; Ruan and Ferrari, 2021; Drake et al., 2022). While turbulent diffusivities can be separately defined for individual tracers (e.g., K_θ for potential temperature, K_S for salinity), it is often assumed that they are approximately equal in sufficiently turbulent environments – consistent with a mixing length argument (Prandtl, 1925).

The Fickian parameterization (equation 3) can be applied to the vertical fluxes of heat $(\rho c_p\overline{w^{\prime }{\theta }^{\prime }})$, salt $(\overline{w^{\prime }{S}^{\prime }})$ and any scalar $(\overline{w^{\prime }{C}^{\prime }})$, e.g., nutrients or dissolved gases (Gregg, 1987):

where J_q, J_S and J_C are the vertical fluxes of heat, salinity and a passive tracer C, respectively. This strategy enables the estimation of vertical fluxes for any parameter for which an accurate vertical gradient can be computed, provided K is known. There is no requirement to co-locate a fast-response sensor for the targeted parameter (e.g., oxygen) with a point-velocity sensor, which is the basis of existing “eddy correlation” field techniques (Pond et al., 1971; Lorrai et al., 2010; Bluteau et al., 2018).

Now focusing on the right hand side of equation 2, the molecular dissipation of a tracer can also be described by χ_C, which is the rate of dissipation of tracer variance, or the rate at which fluctuations in the tracer are smoothed out:

Using equations 2 and 3, the eddy diffusivity terms, K, can then be calculated directly from observations of χ:

with χ_θ and χ_S the dissipation of thermal and salinity variance, respectively. χ_C is the turbulent dissipation of tracer variance and quantifies the destruction of a tracer’s gradient through turbulence.

Similar relations can be derived for momentum from the turbulent kinetic energy equation (Osborn, 1980). Following the same assumption of statistical equilibrium and non-divergence of turbulent kinetic energy, the turbulent flux of momentum is balanced by the turbulent dissipation of kinetic energy ϵ and the turbulent buoyancy flux J_b, leading to:

with $N^2=-\frac{ɡ}{{\rho }_0}\frac{\partial \overline{\rho }}{\partial z}$ K_ρ as the turbulent diffusivity of buoyancy, and Γ as the mixing coefficient defined as a function of J_b and the momentum fluxes $\overline{w^{\prime }{u}^{\prime }}\frac{\partial \overline{\mathbf{u}}}{\partial z}$. In situations where temperature is the main control of the density distribution, it is possible to combine equations 4 and 8 to find that $\mathrm{\Gamma }=\frac{N^2}{2{\partial }_z{\overline{T}}^2}\frac{X_{\theta }}{\varepsilon }$. Typically, the mixing coefficient Γ is assumed to be ≈ 0.2 in the ocean, but its variability and dependence on flow properties remain the subject of continued research (Gregg et al., 2018). However, its variations are typically of order one, much less than the several orders of magnitude over which other turbulence quantities, such as ϵ and χ, vary in the ocean.

Turbulent flux equations (4) and (7) are simplified with diffusivity terms (6) and (8) by assuming that the production of variance is instantaneously balanced by the dissipation of variance (Osborn and Cox, 1972; Osborn, 1980). Although this assumption provides an indirect method of measuring turbulent diffusivities from ϵ or χ, values of K obtained in this manner are consistent with bulk diffusivities inferred from dye release experiments within accepted uncertainty ranges (Ledwell et al., 1998). Computing diffusivity by tracking dye provides a ground-truth of diffusivity, but measuring ϵ and χ is much more pragmatic for large-scale implementation in ocean observation networks.

As mentioned in section 3.2, measuring ϵ and χ requires data collection at centimeter-scale resolution or smaller. The required resolution cannot be obtain with most classic temperature, salinity and ocean velocity measurements, but with data collected at the appropriate resolution, ϵ and χ are “simply” the result of integrating the velocity and scalar (e.g., temperature) variations over the turbulence subranges – namely the inertial and the viscous subrange (Figure 3). The inertial subrange is a wavenumber band where the shear or temperature gradient increases toward smaller scale without being influenced by viscosity, and the viscous subrange corresponds to a wavenumber range where the signal is dissipated due to viscosity, creating a spectral roll-off. The shape and amplitude of the spectral signal can be modeled with theoretical spectra: the Batchelor spectrum for temperature gradient (Batchelor, 1959) and, among others, the Nasmyth spectrum for the shear (Nasmyth, 1970).

Practically, the integration of ϵ and χ is done assuming local isotropy, and integrating the one-dimensional spectrum in one direction, e.g., vertical direction for a vertical profiler or horizontal direction for a horizontal towed or moored sensor:

where ${\varphi }_{u_z}$ and ${\varphi }_{{\theta }_z}$ are the vertical (or horizontal) wavenumber spectra of velocity and temperature gradient (e.g., Lueck, 2022a, Lueck, 2022b), respectively. The spectra are integrated between k₀, the instrumentation-dependent lowest wavenumber accessible through the measurement observations, and k_c, a frequency cut-off separating the instrument noise from the signal. κ_T and ν are the thermal diffusivity and the kinematic viscosity, respectively. In the case of vertical profilers, we use the factors of 7.5 and 6 on the right-hand sides of Equation 9 and 10 to extrapolate the single observed dimension to three dimensions (note that these constants assume isotropic 3D turbulence).

Because ϵ and χ are spatial variance measurements, the sensors must record small variations in velocity and temperature within the microstructure range as a function of time while moving through the water. The time series can then be converted into spatial gradients using the mean flow passing the sensors. On profiling platforms, this is often taken as the instrument fall speed using the pressure data or nearby current data. Hence, measuring χ requires high-frequency thermistors (FP07) moving steadily through the water, i.e., a moving platform or flow passing through the sensor. Measuring ϵ can be done using shear probe (piezo-electric beams embedded in silicon) moving through the water, but is also achievable with other types of sensors and methods (Lueck et al., 2002; Le Boyer et al., 2021). For instance, ϵ can be determined from commercially available acoustic Doppler current profilers (ADCPs) using structure-function methods (e.g., Guerra and Thomson, 2017; McMillan and Hay, 2017), provided their sampling programs are specifically configured for turbulence measurements. Alternatively, point velocimeters that rapidly sample velocities at a fixed point in space can be used to obtain ϵ (Bluteau et al., 2011). Velocity-based methods typically rely on larger scales of turbulence within the inertial subrange, which are less taxing to measure but require the underlying assumptions in (9) to hold at larger scales. Hence, these techniques are typically limited to higher energy environments, but have the advantage of simultaneously providing information about the subsurface currents.

Applying eddy-correlation techniques to measurements from fast-sampling velocity-based instruments is the most direct way to obtain turbulent fluxes of a scalar (e.g., dissolved gases and salt). When paired with fast-response sensors such as dissolved oxygen (Bluteau et al., 2018) or temperature (Polzin et al., 2021), velocity-based instruments can be used to obtain turbulence quantities beyond ϵ and χ, such as the velocity-fluctuations (i.e., Reynolds stresses) or turbulent fluxes (〈w^′C^′〉). The scalar C^′ and the vertical velocities w^′ must be sampled sufficiently fast enough to resolve their time-averaged covariance over the flux-contributing time and length scales. In practice, this requires sensors that are sufficiently small and have response times of less than a second to measure the smallest flux-contributing scales while determining the largest flux-contributing scales to ensure the non-turbulent motions are excluded from the calculations (McGinnis et al., 2008; Lorrai et al., 2010). The relatively slow response times of most scalars (e.g., dissolved oxygen) hinder the application of the eddy-correlation method in high energy flows (high ϵ) because the flux-contributing scales are smaller with increasing ϵ. Because of these theoretical and technological difficulties, eddy-correlation methods have been limited to measurements from bottom landers in relatively quiescent environments; in particular, to examine near-bottom vertical fluxes of dissolved oxygen (e.g., Lorrai et al., 2010) and heat (e.g., Davis and Monismith, 2011). These difficulties are the reason that the scientific community has focused on collecting bulk turbulent flux measurements through ϵ and/or χ, rather than direct turbulent flux measurements.

Other methods also exist to estimate ϵ, which rely on even larger scale measurements than the velocitybased techniques above. Although easier to measure, they come with larger biases than measuring ϵ from shear probes or fast-sampling velocities. For example, using lower-resolution density profiles, ϵ can be computed from Thorpe scale overturns at meters scales. Thorpe scales highlight gravitationally unstable regions of the water column, indicating turbulence (Mater et al., 2015) with the caveat that, to do so, the overturn is assumed to be fully turbulent, which is not always the case and can induce bias in the estimation of ϵ. Alternatively, finescale parameterizations can also be used from lower resolution observations of shear, strain or kinetic energy (Polzin et al., 1995; Polzin et al., 2014; Whalen et al., 2015). Parameterizations estimate ϵ by comparing the observed internal wave field characteristics (e.g., shear, strain) with the Garrett and Munk spectrum – a description of an average internal wave field (Garrett and Munk, 1972). These estimates are known to differ by order of magnitude or more at some locations in the deeper ocean (Figure 4, Klymak et al., 2008; Waterman et al., 2014). Parameterizations rely on the assumption that most mixing is driven by internal waves, and is only applicable in the open ocean (Polzin et al., 2014). Until parameterizations can be built for more specific cases, there will always be a discrepancy between parameterized and measured turbulence values at ocean boundaries.

Figure 4 Evidence of discrepancies between fine-scale parameterization and direct observation from Waterman et al. (2014). Left: A comparison of the station-averaged vertical profiles of ϵ_fine/ϵ_micro (finescale parameterization for ϵ/measured microstructure ϵ) ratio as a function of height above the bottom for various implementations of the finescale parameterizations. The strain parameterization (red lines) is the only one that can be used with Argo float data. Right: Scatterplot of log₁₀(ϵ_micro) versus log₁₀(ϵ_fine) colored by the local predicted lee-wave energy flux.

5 Instrumentation maturity

Turbulence measurements require a high sampling frequency compared to the other EOVs. In this section, we present existing instrumentation allowing for the computation of ϵ or χ (section 5.1), for the direct measurement of turbulent fluxes (section 5.2, or for the integration of both methods on a single platform 5.4).

5.1 Microstructure-based measurements

The first in-situ observations of ocean turbulence were made by researchers at the Pacific Naval Laboratory of Canada’s Defence Research Board using hot film anemometers (Grant et al., 1962). Since then, the ocean mixing community has continued to design and improve instruments that are increasingly resilient and easy to deploy. The rapid evolution of instrumentation has enabled more autonomous measurements and the mobile electronics industry has progressed tremendously over the past decade, reducing the size and power of electronic components while increasing data storage capacity. The types of platforms used to measure turbulence have expanded from profilers and towed vehicles to moorings, gliders, AUVs, autonomous profilers, and Argo floats.

The basic sensor technology available to make direct mixing observations in the microstructure range has been in place since the 1970s (see exhaustive review by Lueck et al., 2002). Airfoil shear probes and FP07 thermistors sample small-scale gradients of velocity and temperature, respectively, at frequencies of 100-512 Hz, which equates to a spatial resolution of centimeters. FP07 thermistors and shear probes each have advantages and disadvantages for measuring turbulence. Shear probes measure small-scale velocity fluctuations, directly measuring the loss of kinetic energy associated with turbulent motions. FP07 thermistors measure small-scale temperature fluctuations and thus provide an indirect measurement of how turbulent motions affect background temperature gradients. Shear probes are more sensitive to platform speeds and vibrations than FP07 thermistors. Conversely, thermistors have a slower response time, limiting their ability to resolve the smallest eddies in areas of strong turbulence at reasonable fall speeds. Both types of sensors are very delicate and can be damaged if they encounter large zooplankton in the water column. A new type of plastic membrane (e.g., polyvinylidene fluoride, PVDF) with similar piezo-electric properties to the currently used ceramic materials has the potential to greatly improve the resilience of shear probes over the next five years. These plastic membranes are flexible and, therefore, more resilient to shocks and pressure forces.

Since the early turbulence measurements, the sensors themselves have not changed significantly, but the electronics and data acquisition systems have improved resulting in higher signal-to-noise ratios and reduced power consumption. Onboard data storage has also increased significantly, allowing the on-board processing of the turbulent dissipation rates (ϵ,χ). Autonomous platforms (e.g., gliders and floats) already use on-board processing schemes (e.g., Hughes et al., 2023). Using vibration sensors and data processing algorithms has also reduced signal contamination from sources such as platform vibrations (Goodman et al., 2006). One can also rely on the larger scales of the microstructure subrange to estimate ϵ (Bluteau et al., 2016b) and χ (Bluteau et al., 2018) to avoid vibration issues by sampling less rapidly. Furthermore, progress has been made in estimating the statistical uncertainty of a dissipation estimate and the quality of a spectrum objectively (Lueck, 2022a,b), which allows automated and reliable quality control of dissipation estimates. These improvements, combined with the recent progress in real-time data processing, have 1) allowed more compact instruments to be built, 2) enabled the transmission of data over satellite connections, and 3) allowed for instruments to be deployed for longer periods (e.g., Rainville et al., 2017). They also enable more data to be recovered from autonomous platforms deployed in high-traffic areas where platforms are sometimes damaged or lost.

5.2 Velocity-based measurements and platforms

Acoustic-Doppler Current Profilers (ADCPs) and point-velocity measurements are instruments routinely deployed on bottom landers and moorings. Estimating ϵ from these sensors requires measuring velocities at a sufficiently high rate to sample within the inertial subrange (seconds and several cms). The point-velocity measurements are more mature than using ADCPs to obtain ϵ as it relies on fitting the inertial subrange of spectra derived from time-series over a small sampling volume of a few centimeters. The first estimates of ϵ were from electro-magnetic sensors (Bowden and Fairbairn, 1956) or drag spheres (Lumley and Terray, 1983) until acoustic-Doppler velocimeters with sufficiently low measurement noise became common in the late 1990s (e.g., Voulgaris and Trowbridge, 1998; Kim et al., 2000). Initially, turbulence estimates were only possible from fixed platforms since high-frequency platform motion (e.g., vibrations) could not be easily removed from the velocity signals. Over time, techniques to handle the changes in expected spectral forms from surface waves were generalized (Lumley and Terray, 1983; Feddersen et al., 2007), while reducing the impact of motion contamination on moored point-velocity measurements (e.g., Bluteau et al., 2016a).

The use of ADCPs to derive ϵ, can provide estimates over a wider range of depths than one point-velocity meters. The technique was first introduced by Wiles et al. (2006) and relies on the structure-function first employed in meteorology in the late 1960s (Sauvageot, 1992). Deriving ϵ requires differencing instantaneous turbulent velocities at different separation distances along an ADCP beam. As with the point-velocities, the sampling rate must be sufficiently fast to measure within the inertial subrange (1-2 Hz typically). The noise and measurement quality must be much higher than for measuring mean currents, and so the ADCPs must be programmed in pulse-to-pulse coherent mode. This mode results in higher resolution measurements but over a small spatial range of typically less than 10m. The method also requires that the bin separation is sufficiently small so that a few bins are fully contained within the inertial subrange. For relatively weak turbulence levels of ϵ ∼ 10⁻⁸ Wkg⁻¹, a bin size of 20cm is sufficiently small. In more energetic flows, the inertial subrange extends to smaller scales but also to larger scales such that the bin size does not need to be reduced at the expense of the ADCP’s vertical sampling range.

Several commercial offerings currently exist to measure ϵ over a wide range of turbulence levels and spatial scales with the structure-function. For example, the structure-function has been applied to ADCPs available from Nortek and RDI Teledyne (see Guerra and Thomson, 2017, for comparisons between instruments). The structure-function method was been applied to moored ADCPs using off-the-shelf instruments (Lucas et al., 2014). Velocity profiles of sufficient quality for turbulence analysis are also possible from off-the-shelf instruments mounted on Lagrangian (drifting) floats (Shcherbina et al., 2018).

Another type of acoustic measurement uses travel-time velocimeter technology (TTV, Polzin et al., 2021). Similar to the point velocities, TTV instruments measure a small volume of water using a pair of transducers “pinging” at each other. This technology allows for a very fast sampling in low-scattering environments and can be used for direct observation of turbulent fluxes (eddy-correlation methods). The eddy-correlation method is relatively new in the deep ocean science community and requires a time-scale separation that could lead to bias in the turbulent fluxes computation.

5.3 Platforms

Microstructure measurements have been made using several different instruments mounted on a variety of platforms (Figure 5). The first measurements were made using horizontal profilers that were towed by a ship (e.g., Grant et al., 1962). By the late 1970s, most microstructure measurements were made using vertical profilers developed by over 10 different research groups around the world (Lueck et al., 2002). Currently, microstructure instruments are made commercially by Rockland Scientific International (Canada) and Sea and Sun Technology (Germany), and by research groups at Oregon State University (OSU), Scripps Institution of Oceanography (SIO), and Woods Hole Oceanographic Institution (WHOI).

Figure 5 Examples of platforms measuring ϵ, χ or turbulent fluxes directly.

Most of these instruments are designed to make simultaneous measurements of both shear and temperature microstructure from vertical profilers (autonomous or tethered) for a variety of depth ratings (up to 11,000m) and modular instruments that can be mounted on autonomous platforms such as gliders, AUVs, and drifting or moored profilers (e.g., Argo floats or Wirewalkers). Vertical profilers remained the dominant method of obtaining microstructure measurements until the 2010s, when instruments began to be mounted on moorings, gliders, AUVs, and surface-following platforms (e.g., Moum and Nash, 2009; Fer et al., 2014; Hughes et al., 2020; Iyer et al., 2021; Zippel et al., 2021). These new platforms, combined with more versatile instruments, have led to a significant increase in the volume of data being collected and the duration of deployments. Deployments on autonomous platforms are power-limited and typically last up to 45 days. This equates to 300 glider profiles to 1000m.

Longer turbulence deployments are possible from moored platforms, although measurements are limited to discrete depths rather than the full water column. χpods, which contain FP07 sensors mounted horizontally, have provided near-continuous mixing measurements at the equator on several of the Tropical Atmosphere–Ocean (TAO) array moorings since late 2005 (Moum et al., 2013) in the Pacific Ocean, and later in the Atlantic Ocean (PIRATA array) and Indian Ocean (RAMA array). The TAO/TRITON (Triangle Trans-Ocean Buoy Network) program evolved from the original TAO project started in 1985 and was sponsored by CLIVAR, GOOS, and SCOR (https://www.pmel.noaa.gov/gtmba/, McPhaden et al. (2010)). PIRATA launched in 1997 (Bourlès et al., 2008), RAMA evolved from the Indian Ocean Observing System (IndOOS) with original deployment in 2006 (McPhaden et al., 2009; Beal et al., 2020). Only a few sensors allowing for direct measurements of turbulent fluxes (e.g., MAVS Figure 5 using TTV technology) can be mounted on mooring and record velocity-tracers fluctuations over a wide range of frequencies.

CTD-χpods – an adaptation of χ-pods suitable for mounting on a standard rosette (Figure 6) – have been deployed on cruises that are part of the international GO-SHIP Repeat Hydrography Program since 2014. Through this effort, CTD-χpods have provided χ and turbulence measurements across sections in the Pacific, Atlantic, Indian, and Arctic Oceans with the intent to obtain repeated full-depth measurements on a global scale.

Figure 6 Left: CTD-rosette equipped with CTD-χpods deployed during a recurrent GO-SHIP transect. Right: CTD-χpod transects as of 2022.

5.4 Argomix

One of the most relevant programs to our discussion is the international Argo program, which collects ocean data using a fleet of autonomous instruments that drift with the ocean currents and move up and down between the surface and a mid-water level. The Argo floats are part of a global coordinated effort of over 30 nations with more than 3500 autonomous ocean-observing profilers currently active in the network (http://www.jcommops.org/). The floats provide the scientific community with near-real-time ocean measurements including temperature and salinity profiles, so they are already an important platform for physics EOVs. The ongoing integration of turbulence sensors on Argo floats has the potential to provide turbulent dissipation rate measurements on a daily basis at a global scale. The present discussion recognizes the need for the ocean mixing community to further coordinate its efforts to measure and distribute data, taking advantage of the means provided by the GOOS.

Mixing measurements from Argo floats (Argo-mix) are one of the emerging branches of the Argo program (Roemmich et al., 2019). This new branch, will provide the scientific community repeated microstructure measurements at a global scale. These measurements would offer insights into the impact of ocean mixing on water mass transformations, air-sea interactions and a plethora of other processes (Naveira Garabato and Meredith, 2022). Technological improvement needs to provide resilient and low-power new sensors so their integration does not modify the life expectancy or the core mission of a float– namely long-term deployments profiling the upper 2000 m.

An early effort obtained more than 1000 profiles over two weeks in the Bay of Bengal from two χSOLO floats, which are Argo floats equipped with fast thermistors (Shroyer et al., 2016). This instrument motivated development of the Flippin’ χSOLO (FCS), which has shear probes and fast thermistors at one end of the float and a communication antenna at the other (Moum et al., 2023). By flipping around at the top and bottom of a cast, turbulence profiling of undisturbed fluid is achieved on both down- and up-casts all the way through the sea surface. Similar to glider performances, FCS missions can last 45 days when profiling to 120m at a rate of roughly 100 profiles per day. Missions are controlled from a shore station by remote communications.

While χSOLOs have reached a mature level of operation, additional microstructure floats are under development (Figure 7). The microstructure sensors developed by Rockland Scientific International and SIO are currently being integrated in the MRV ALTO and APEX floats, respectively. Both of these floats are part of the Argo fleet. The deployments of these floats began in 2022 during a few pilot experiments that profiled to 2000 m. The first results show a good sensitivity range despite the slow motion of the floats. The integration of these turbulence sensors on Deep-SOLO floats – Argo float profiling down to 6000 m depth – is also underway. The probes used on the sensors are expected to survive the challenging pressure at such a depth since they show normal sensitivity and behavior after multiple bench-tested pressure cycles up to 10000 PSI (>7000 m depth).

Figure 7 Examples of Argo floats equipped with microstructure packages: An Apex float from the University of Washington with an Epsilometer from SIO (Left), A MicroAlto float from MRV equipped with a Rockland MAPLe™ (middle) and a SOLO float equipped with turbulence sensors from Oregon State University.

A new project from the National Oceanographic Partnership Program (NOPP) aims to improve the global simulation of internal waves and their impacts on ocean mixing. As part of this effort, the SQUID (“Sampling QUantitative Internal-wave Distributions”) sub-project sponsored by NSF will deploy 50 EMAPEX microstructure floats measuring velocity, temperature, salinity, and χ from cruises of opportunity throughout the world in 2023–2025 (Figure 8). Sampling will differ from the Argo paradigm by making bursts of 4-6 profiles to 2000 meters over two days, allowing the estimation of a statistical average and high-frequency perturbations in shear, strain, and dissipation, as well as energy flux by dominant internal wave components through low-mode fitting and harmonic analysis. Bursts will be separated by multi-day drifting at 1000m, with a target duration of one year per float. The overall SQUID array, while not covering the globe by any means, will aim to span the range of major internal wave forcing parameters (including internal tide amplitude as estimated from altimetry, bathymetric roughness, high-frequency wind forcing, mesoscale eddies, and latitude).

Figure 8 Map of possible deployment sites for 50 EM-APEX floats from cruises of opportunity in 2023 2025. These floats make up the SQUID (Sampling QUantitative Internal-wave Distributions) contribution to the NOPP Global Internal Waves experiment, aimed at validating internal wave models and their implications for mixing by internal waves. Courtesy of James Girton, APL-UW.

5.5 Error bar and uncertainties

Practical measurements of ϵ (and fluxes derived from these) have unavoidable uncertainties. Even an instrument with multiple (redundant) sensors will not necessarily produce the same values. Further, turbulence has short spatial and temporal scales and therefore point measurements are not necessarily representative of that part of the ocean on longer (e.g., >weekly) or larger (e.g., >10km) scales and requires multiple realizations to address the time and space variability turbulent diapycnal fluxes.

Lueck (2022a) showed that ϵ values measured by individual shear sensors on the same instrument agree in an average sense but that there is typically variability up to a factor of two. Specifically, for 90-cm segments, agreement between sensors was worse than a factor of two in 10% of cases. Shorter segments have larger variability and longer segments have smaller variability. Similar levels of agreement are described by Oakey (1982) who compared ϵ derived from a shear probe to that derived from a thin film thermometer, and by Kolås et al. (2022) who compared outputs from two shear probes installed on an AUV. (Further examples are given in Section 3.1 of the turbulence methodology review by Burchard et al. (2008)).

Averaged over hours to days, turbulence measurements from independent instruments—but from the same nominal part of the ocean—should give the same results. Moum et al. (1995) tested this with two free-falling profilers from two ships spaced within 11km of each other on the equator (0°, 140°W). At most depths, the agreement between the full 3.5-day means were within a factor of two. Agreement between individual data pairs was typically within a factor of three but sometimes beyond a factor of 10. In other words, the natural variability of geophysical turbulence over kilometer scales is comparable to, or larger than, the uncertainty associated with each point measurement of ϵ.

Clearly, if turbulence measurements at a given spot are to be considered representative, averaging is needed and/or uncertainty levels need to be quantified. As with any measurement, but especially for turbulence, more samples are better for achieving accurate statistics. Observed distributions of ϵ are often approximately lognormal. Hence, distributions of log(ϵ) are approximately Gaussian, which makes standard statistical quantities such as mean and standard deviation easy to calculate (e.g., Baker and Gibson, 1987). However, one can avoid assuming a Gaussian distribution by using bootstrapping, a statistical method that is agnostic to distribution shape. Indeed, bootstrapping is a standard procedure in turbulence analysis (e.g., Shay and Gregg, 1986; Lu et al., 2000; Greenan et al., 2001; Nash and Moum, 2001; Inall and Rippeth, 2002; Klymak and Moum, 2007; St Laurent and Thurnherr, 2007; Perlin and Moum, 2012; Whalen et al., 2012; Sutherland et al., 2013; Waterhouse et al., 2014; Wenegrat and McPhaden, 2015).

This discussion of uncertainties, which so far only considers ϵ, can be extended to turbulent fluxes, albeit with caution. As Moum (1997) emphasizes, microstructure measurements are not flux measurements, and the assumptions linking the two are ‘fraught with uncertainty’. That said, direct numerical simulations show that the Osborn and Cox (1972) and Osborn (1980) relations outperform their underlying assumptions in predicting diffusivity and hence the irreversible buoyancy flux (Taylor et al., 2019). This is consistent with the widely held notion from observations that Γ = 0.2 (see Gregg et al., 2018, and elsewhere in this paper).

6 Data processing and international coordination

The ocean mixing community is increasingly committed to integrating observation efforts and standardizing methods for data quality control and distribution. An archive for turbulence data is maintained by the CLIVAR’s Carbon Hydrographic Data Office (microstructure.ucsd.edu). This database was created by the Climate Process Team (CPT), funded by NSF and NOAA between 2010 and 2015. The team’s mission was to develop and test mixing parameterizations for climate models. To this end, they created an archive for turbulence data that compiles processed datasets of dissipation rates from research programs in a single repository. The format of the data included in this archive is standardized so that observed turbulent dissipation rates can be easily incorporated into models or process studies.

Creating a common dataset is a complex undertaking due, in part, to the sensitivity of the observed dissipation rates to the different algorithms independently developed and customized for the different platforms (up to five in the CPT database). The ATOMIX working group was founded in 2020 to establish a consensus regarding best practices for the derivation of ϵ from observed shear microstructure and acoustic measurements. These efforts come on the heels of a nearly 20-year history of official coordinated efforts of the ocean mixing community; a prior SCOR working group (Working group 121) began in 2002 with a goal to summarize the state of knowledge about ocean mixing and to identify the steps that would be required to improve mixing parameterizations in global climate models. Going forward, ATOMIX will establish and communicate measurement standards through a wiki and propose a set of validated data to build calibrated global turbulence measurements.

7 Feasibility and cost effectiveness

Besides the scientific relevance of observing subsurface turbulent fluxes at a global scale, the feasibility and the cost-effectiveness of these measurements are two of the main criteria that must be met to show that turbulent fluxes belong in the EOV list. In the EOV context,

● Feasibility implies that observing or deriving turbulent fluxes on a global scale is technically feasible using proven, scientifically understood methods.

● Cost effectiveness means that generating and archiving turbulence data is affordable, mainly relying on coordinated observing systems using proven technology, taking advantage where possible of historical datasets.

7.1 Feasibility and pilot programs

The TAO and CTD-χpods projects mentioned earlier have demonstrated the feasibility of collecting basinwide datasets of microstructure observations (Figure 6). These projects integrated turbulence measurements into the structure and culture of repeat hydropgraphy cruises and developed a methodology for processing and interpreting the acquired ocean mixing data. CTD-χpods are also being tested in coastal environments to possibly widen the scope of their usage. Expanding to global observations is a matter of putting more instruments in the ocean.

Measurements of ϵ are reaching a similar level of maturity within the Argo-mix context and the ongoing integration of turbulence packages measuring ϵ and χ inside Argo floats. The float controls a unique low-power electronic board that can collect data from analog channels and has the required memory space to compute the dissipation rates and a few quality flags. Such packages, as well as the use of easy-toreplace shear and FP07 probes, simplify the integration of microstructure sensors greatly. Once they are an integrated part of the Argo fleet, microstructure measurements at a global scale will automatically be available to a wide community. Before their full integration, a few regional pilot experiments (e.g., the SQUID project mentioned in section 5.4) will be performed to validate and prove that measurements from autonomous Argo floats are fit-for-purpose.

The recent success of the biogeochemical community in incorporating sensors onto global observation platforms like Argo floats can be used as a roadmap for the ocean mixing community. Autonomous floats or vehicles able to sample the turbulent fluxes and diffusivities during the lifetime of coherent oceanic structures (e.g., mesoscale eddies or equatorial cold tongue) can provide insights about local energy transfers inside structures impacting the global circulation. The scales (time and space) of these structures makes them ideal candidates for pilot experiments that will demonstrate that Argomix floats are also fit-for-purpose of the core Argo mission.

Argomix floats, ϵ measurements from already existing platforms (e.g., gliders, profilers, high-resolution ADCPs or ADVs), and current efforts to develop a consistent methodology for processing and interpreting measurements (e.g., ATOMIX SCOR group) all converge toward integrating turbulent flux observations into the structure and culture of repeated hydrography cruises and observing systems.

7.2 Cost-effectiveness

The cost-effectiveness of repeated microstructure measurements represents a balance between 1) the price of the turbulence sensors, 2) the cost of the platforms used for these measurements and 3) the funds dedicated to coordinate observing systems and ship opportunities in charge of maintaining these measurements. Over the past 40 years, many regional research projects have been funded by national research agencies in attempts to quantify turbulent dissipation (Waterhouse et al., 2014). Regional scale experiments provide high-resolution time series of ocean turbulence that shed light on its variability and the multi-scale interaction between ocean energy sources and energy dissipation. Understanding these interactions is of primary importance for future human activities (e.g., deep sea mining, energy harvesting, sea level rise). As the connection between ocean mixing and human activities has become more apparent, private organizations have begun supporting studies that include microstructure measurements (e.g., Schmidt Family Foundation and the Benioff Foundation, Muñoz-Royo et al., 2021). In 2018, the physical oceanography program at NASA identified ocean mixing as one of its three priority areas and funded several projects, including efforts associated with the χ measurements made on the international GO-SHIP hydrography program, described above.

The cost of autonomous platforms like Argo floats spans from $20k for a core (Temperature-SalinityPressure) float to about $100k for a biogeochemical Argo float with six additional sensors. As of today, the microstructure packages considered for the Argo float integration range around $30k–$40k in their current form. Progress in the electronics industry is triggering the development of new generations of turbulence sensors and new fabrication techniques are lowering the fabrication cost of these packages by a factor of three (Le Boyer et al., 2021). This will bring the cost of an “Argomix” float well within the price range of a core Argo float and a biogeochemical float. The necessary simplification of turbulence packages to fit into the tight space of the float, as well as any production scaling, will further reduce the price of Argo-mix floats.

Consequently, the ocean mixing community will be able to design regional scale experiments with sampling strategies using a small fleet of ∼10 microstructure sensors measuring both ϵ and χ from a variety of easy-to-deploy platforms with a<$500k budget within the next five years (e.g., SQUIDD project mentioned above).

8 Conclusions

Currently, global climate models do not resolve turbulent mixing, and instead use turbulent diffusivity schemes to parameterize turbulent fluxes using the physics of processes (e.g., internal waves, mesoscale and submesoscale processes) that are thought to drive part of the ocean mixing. In regions where turbulent production is dominated by the internal wave field, the finescale parameterization can be used to predict the patterns of ocean mixing based on the resolved model state. However, strong ocean mixing signals occur near the boundaries (surface and bottom) and in regions where mesoscale and submesoscale processes are important. The dynamics of these processes modify the internal wave field and violate the assumptions required to use the parameterization. As a result, models poorly represent ocean mixing where it matters the most, and the turbulent diffusivity scheme needs to be improved by observational inputs.

Not only are parameterizations of ocean mixing failing to estimate dissipation rates in key regions, but evidence is mounting that the picture drawn in the 1960s, presenting turbulent diffusivities as a largely dynamically passive phenomenon, does not provide an adequate explanation of many climatically important aspects of the ocean’s behavior. Major elements of the oceanic circulation are governed by dynamic interactions between ocean mixing processes and the large-scale ocean state, and are thus not well captured by the current generation of state-of-the-art Earth system–class ocean models. The conclusion seems inevitable that to generate a step change in our understanding and modeling capability, sustainable and global turbulence observations are required to describe and quantify the circulation-shaping role of ocean mixing (Naveira Garabato and Meredith, 2022).

The instrumentation for such observations has been developed since the 1960s and recently benefited greatly from the mobile technology revolution and the chip-miniaturization that went with it. Storage capacity, processing speeds and reduced power consumption allow for cheaper and modular instruments that can be embedded in autonomous platforms like Argo floats, gliders or AUVs. With their extended range and life expectancy, these platforms can now be used for large-scale and long-term experiments. The obvious and on-going next step for the ocean mixing community is to develop methods and coordinate efforts between research groups to ensure the creation of a sustained global observing system for mixing (Naveira Garabato and Meredith, 2022). The instruments and platforms mentioned in this document can be sorted, using the UNESCO readiness level grid, in terms of maturity for being deployed on a global scale. All of them are past the proof of concept level and are deployed in the context of pilot experiments (Figure 9).

Figure 9 Levels of maturity of turbulent diapycnal flux observations mapped onto the EOV framework: Requirements (i.e., instrumentation, Methods), Coordination (i.e., strong community) and Data management (i.e., database and processing best practices). The Framework for Ocean Observing evaluates possible new components to be included in the Global Ocean Observing System (GOOS) based on their readiness level (UNESCO, 2012). A variable in levels 1-3 is identified as societally important and observations have been made at basin scale levels. A variable in levels 4-6 has a verified sampling strategy with autonomous deployment and international collaboration. A fully mature variable is validated through peer review and has sustained measurements at global scales. Data products are routinely publicly available. Turbulent diapycnal fluxes fall in the range of “Pilot” level maturity (Levels 4-6), during which instruments and coordination are tested and made ready for large-scale implementation. The components are well on their way to the “Mature” range of readiness (Levels 7-9).

A number of initiatives and experiments participate in that goal. The Scientific Committee on Oceanic Research is supporting ATOMIX, a working group creating a universal set of standards for measuring and distributing ocean mixing data with appropriate quality control flags. Similarly, the National Science Foundation and National Oceanic and Atmospheric Administration supported the climate process team on internal-wave driven ocean mixing, which improved diffusivity parameterizations in global climate models (MacKinnon et al., 2017). A consequence of these efforts to coordinate research groups was the creation of an ocean turbulence database (Waterhouse et al., 2014). Along with these community efforts to centralize and standardize data production, sustained large-scale observation programs, like GO-SHIP and the Argo program, have either already integrated or are currently integrating turbulent mixing instrumentation.

Several research groups and companies are currently working on the integration of turbulent sensors in the Argo program. Floats with integrated turbulence packages are already part of regional experiments (Shroyer et al., 2016). The engineering challenges of sensor resilience, energy consumption and costs are being overcome thanks to new technologies and materials. At the time of writing, the turbulence sensing community is about five years away from producing instruments that can be easily deployed and survive hundreds of profiles, collecting quality-controlled estimates of ocean mixing. Keeping the pace of current development, major programs (e.g., Argo, NASA) are acknowledging that ocean mixing is relevant to climate and ocean health, as discussed in Section 1. Consistent global scale measurements will have an positive impact on the reliability of global climate models, our understanding of global circulation patterns, and our ability to map and predict physical, chemical, and biological parameter distributions on regional and global scales. As feasibility challenges are overcome through emerging technologies and the cost per unit comes down, the inclusion of ocean turbulent mixing as an Essential Ocean Variable should be on the horizon.

Author contributions

Lead authors AB and NC, designated by the symbol “*” and listed first in the list of authors, had primary responsibility for writing all sections of the manuscript. The other co-authors contributed in the following: Introduction – All co-authors with strong participation from HD, KH, CB; section 2 – HD, KH, CB, AN, MM, IF; section 3 – HD, CB, AN, TP, EF, AMa, LC, AMe, IF; section 4 – CB, KH, MA, CS, MD, IF; section 5 – MA, CB, IF, MD, AMo; section 6 – CB, KH, AMo; section 7 – AN, CB, AMo; Conclusion: All co-authors. All authors contributed to the article and approved the submitted version.

Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

Acknowledgments

All co-authors, and especially the lead authors acknowledge the tremendous input from J. Moum who helped re-oriented the main topic of this white paper toward what the ocean community is now considering the best option for turbulence data to become an EOV. We are also in debt of all the pioneers who participated in the development of today’s turbulence technology.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References