## 1. Introduction

Atmospheric turbulence is a complex phenomenon, characterized by the presence of a plethora of scales (eddies) ranging from small viscosity-dominated to large, energy-containing structures. Turbulence may undergo space and time variations due to rapidly changing external conditions, it may be locally suppressed or enhanced. To describe characteristic features of turbulence, statistical theories are sought for. In this context, a number of recent research works (Vassilicos 2015; Mazellier and Vassilicos 2008; Bos and Rubinstein 2017; Bos et al. 2007; McComb et al. 2010) address the problem of the equilibrium Taylor law (Taylor 1935) and its failure in the presence of rapid changes of the system. A new, nonclassical, although universal scaling is introduced to describe the latter.

*ϵ*with the statistics determined mostly by large turbulence eddies: the root-mean-square of the turbulent velocity fluctuations

*s*= (∂

_{ij}*u*∂

_{i}/*x*+ ∂

_{j}*u*∂

_{j}/*x*)

_{i}*/*2,

*u*is the velocity fluctuation component,

_{i}*k*is the turbulence kinetic energy, and

*ν*is the kinematic viscosity.

*C*= const underlies the well-known Kolmogorov–Richardson equilibrium cascade picture. It is also a foundation of many eddy-viscosity turbulence models, such as

_{ϵ}*k*–

*ϵ*(Jones and Launder 1972), where the characteristic length scale of turbulence eddies is estimated from the turbulence kinetic energy

*k*and the dissipation rate

*ϵ*. Another conclusion that follows from Eq. (1) is that

The general validity of the “equilibrium” equations, Eqs. (1) and (2), was coming into question, as laboratory studies by Antonia and Pearson (2000) and Burattini et al. (2005) revealed that *C _{ϵ}* varies considerably. This behavior was first explained by the dependence on the inflow and boundary conditions. In another work by Bos et al. (2007) it was shown based on theoretical arguments that

*C*takes different values in forced and decaying turbulence, even though the classical, Kolmogorov −5

_{ϵ}*/*3 form of the spectrum was assumed.

*C*is not constant, but depends on the inlet conditions and local Reynolds number

_{ϵ}_{0}=

*U*

_{∞}

*L*/

_{b}*ν*,

*U*

_{∞}is the flow speed at the inlet,

*L*is the length scale defined by the grid spacing,

_{b}*m*and

*n*are power coefficients and it was found experimentally that both are close to unity. Moreover, instead of (2), another, striking relation was observed:

*increased or decreased.*

_{λ}Rubinstein and Clark (2017) argued that turbulence characterized by the classical laws (1) and (2) should be called “equilibrium” turbulence. In such type of flow the energy spectrum takes a self-similar form and the classical −5*/*3 scaling holds (Kolmogorov 1941). Nonequilibrium states of a flow field appear after a sudden change of external conditions (e.g., forcing) when the system evolves toward another equilibrium (Mahrt and Bou-Zeid 2020; Wacławczyk 2021, 2022). In the nonequilibrium turbulence, the turbulence dissipation coefficient *C _{ϵ}* and the integral-to-Taylor-scale ratios are described by relations (3) and (4) and a nonequilibrium correction should be added to the spectrum. Bos and Rubinstein (2017) derived the form of this correction, as well as the nonequilibrium scaling laws for

*C*and

_{ϵ}The nonequilibrium relations (3) and (4) are substantially different from their equilibrium counterparts in (1) and (2). Hence, classification of turbulence in a given area can be made on this basis. Previous research works which focus on such analysis concern laboratory or numerical experiments (e.g., decaying grid turbulence, turbulent boundary layer) under controlled conditions (Valente and Vassilicos 2015; Obligado et al. 2016; Cafiero and Vassilicos 2019; Obligado et al. 2022). The main idea of this work is to use this procedure to analyze data from in situ airborne measurements of atmospheric turbulence, where the flow conditions in the investigated region are unknown a priori and only limited information, that is 1D intersections of the flow field along the flight track are available. In principle, such analysis could be performed on any high-resolution observations of the atmospheric boundary layer. In this study we emphasize theoretical aspects of the method to be used; hence, we decided to analyze data already known from previous campaigns, but in a way that adds one more important aspect of turbulence. Our choice is the stratocumulus-topped boundary layer (STBL). Several research campaigns and subsequent data analyses were devoted to the study of the turbulence in STBL (see, e.g., Stevens et al. 2003; Malinowski et al. 2013; Siebert et al. 2010; Jen-La Plante et al. 2016). We focus on the data from the Azores Stratocumulus Measurements of Radiation, Turbulence and Aerosols (ACORES) campaign collected in the eastern North Atlantic (Siebert et al. 2021). Nowak et al. (2021) compared the properties of turbulence between the cases of coupled and decoupled STBLs [see Wood et al. (2015) for the definition of decoupling]. They suggested that the observed differences in inertial range scaling might be related to different stages of turbulence lifetime (e.g., development and decay) but did not study these mechanisms in detail. In this work, from the same dataset, we derive all turbulence statistics necessary to calculate *C _{ϵ}*,

*. The aim is to divide the flow area into regions where turbulence decays, develops or is in a stationary state, by investigating values of*

_{λ}*C*. Moreover, by plotting

_{ϵ}*C*and

_{ϵ}*we want to identify whether turbulence is in or out of equilibrium. Additionally, by the curve fitting of the data, estimation of the initial (or equilibrium) Reynolds number becomes possible. We also propose a procedure to roughly estimate the characteristic time scale of the turbulence decay/development. To the best of the authors’ knowledge, this is the first application of such analysis in the study of atmospheric turbulence.*

_{λ}We perform a sensitivity study to show that all statistics, including *ϵ*, can be estimated from the wind velocity records, even in the case of nonequilibrium, provided that the resolution of the data is good enough and the flight segments are sufficiently long. These conditions are fulfilled by the ACORES data. Moreover, we extend analysis of Bos (2020) and estimate the upper bound of *C _{ϵ}* for flows with nonzero buoyancy.

The present paper is structured as follows. In the following section 2 formal definitions of the “equilibrium” and “nonequilibrium” turbulence are presented and the theoretical derivations of Bos and Rubinstein (2017) are recalled. In section 3 the investigated data are briefly presented. This is followed by the detailed description of the methods applied in the data analysis. Results are presented in section 4, and finally, conclusions and perspectives are discussed in section 5.

## 2. Theory of nonequilibrium turbulence

### a. Nonequilibrium scaling laws

This section is devoted to the theory presented previously in the works of Bos and Rubinstein (2017) and Rubinstein and Clark (2017). Their derivations are briefly recalled here for the sake of clarity.

*C*= const refers to the case where turbulence energy spectrum takes a (possibly time-dependent) self-similar form:

_{ϵ}*κ*is the wavenumber,

*L*(

*t*) is a characteristic length scale (e.g., integral), and Φ is a function of a “fixed” form. All considerations are performed under the assumption of isotropy of the flow field. With the use of a spectral model those authors conclude the Kolmogorov K41 −5

*/*3 scaling is the solution for a stationary flow. In the case of nonstationarity, a small, subdominant correction appears; however, the K41 inertial range can still be treated as an approximate solution. Equation (5) hence describes turbulence in a state of a possibly unsteady equilibrium. As shown by Bos et al. (2007) in such state

*C*does not depend on Re

_{ϵ}*; however, it can still take different values depending on whether turbulence is stationary or decaying. In fact, in the latter case the observed*

_{λ}*C*was ca. twice higher than in the forced case. This observation concerned final stages of decay in the grid turbulence, where the self-similar state was reached. Another conclusion follows for the region closer to the grid, where

_{ϵ}*C*≠ const and turbulence is out of equilibrium. According to Rubinstein and Clark (2017) such a state occurs when turbulence evolves rapidly from one self-similar state to another. Those authors propose to express the energy spectrum as a sum of equilibrium and nonequilibrium parts

_{ϵ}*E*

_{0}(

*κ*,

*t*) satisfies the equilibrium relation Eq. (5). Bos and Rubinstein (2017) derived the Kolmogorov’s relation for

*E*

_{0}in the inertial range

*C*≈ 1.5 is the Kolmogorov’s constant, and the following form of the nonequilibrium part

_{K}*κ*<

_{L}*κ*<

*κ*, see Fig. 1. Outside this range

_{η}*E*(

*κ*,

*t*) was assumed to be zero. Explicit form of

*E*(

*κ*,

*t*) was also derived in Yang et al. (2018), using somewhat different approach, which includes dissipative part of the spectrum.

*C*:

_{ϵ}*n*and

*m*in Eq. (3) are close to 1. Similarly, the ratio

All considerations of Bos and Rubinstein (2017) were performed under the assumptions of homogeneity and local isotropy. These assumptions are usually employed to calculate turbulence dissipation rate from in situ measurements performed by aircrafts, where it is also assumed that the energy spectrum follows the −5/3 law. We also use the assumption of local isotropy in our analyses; however, following Bos and Rubinstein (2017) we account for possible deviations from the Kolmogorov’s scaling due to nonstationarity. The next step, left for further work, would be to account for inhomogeneity and anisotropy. This could be done along the line of the studies of Chen and Vassilicos (2022), where transport equations for the second-order structure functions are considered. We note that the nonequilibrium scaling relations (11) and (12) are also present in the nonhomogeneous flows, like, e.g., turbulent boundary layers or wakes, cf. Obligado et al. (2016, 2022), Nedić et al. (2017); hence, we could assume that the dependence of *C _{ϵ}* and

*will not be much affected by anisotropy in the atmospheric boundary layer flows, although the proportionality constants may be somewhat different.*

_{λ}### b. Estimation of C_{ϵ} from turbulence models

*C*, based on common turbulence models. For this, Eq. (9) was used:

_{ϵ}*ϵ*was replaced by

*ϵ*

_{0}, as the nonequilibrium correction

*k*–

*ϵ*closure equations [cf. Pope (2000) and appendix B], which model the time evolution of the kinetic energy and

*ϵ*. It is to note, however, that the transport equation for

*ϵ*is purely empirical and its form cannot be derived from first principles; therefore, it should be used only due to a lack of better options.

*k*–

*ϵ*model for buoyant flows, as suggested by Rodi (1987) and assume that transport terms of

*ϵ*are negligible in comparison to other components in the

*ϵ*equation. With this we have

*P*is the shear production of the turbulence kinetic energy,

*G*is the buoyancy production (see appendix B), and the model constants take the following values:

*A*

_{1}= 1.44,

*A*

_{2}= 1.92, and

*A*

_{3}= 1.14. It should be noted that the latter appears to be highly case dependent and different values and modifications of the model were reported in the literature. The one used here was suggested by Baum and Caponi (1992).

*exceeds a critical value 0.2–0.25. However, some small-scale “supercritical” turbulent motions can still survive up to Ri*

_{f}*= 0.5. We can now consider different values of*

_{f}*P*/

*ϵ*

_{0}and Ri

*and estimate corresponding*

_{f}*C*/

_{ϵ}*C*

_{ϵ}_{0}from Eq. (16). Ri

*= 0.5 will be considered as the upper limit of the flux Richardson number. As reported by Bos (2020) for the case*

_{f}*G*= 0, when the production is equal to the dissipation

*P*=

*ϵ*

_{0}

*k*–

*ϵ*model considered predicts stationary state for production slightly exceeding dissipation. For the case without buoyancy,

*C*is bounded from above by the value (Bos 2020)

_{ϵ}*P*/

*ϵ*= 0. This follows exactly from Eq. (15), with

*G*/

*ϵ*= 0. We estimated the ratio

*C*/

_{ϵ}*C*

_{ϵ}_{0}for nonzero buoyancy from Eq. (16). Results are presented in Fig. 2.

Values of *C _{ϵ}*/

*C*

_{ϵ}_{0}estimated from Eq. (16).

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

Values of *C _{ϵ}*/

*C*

_{ϵ}_{0}estimated from Eq. (16).

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

Values of *C _{ϵ}*/

*C*

_{ϵ}_{0}estimated from Eq. (16).

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

Generally when turbulence becomes locally stronger, *C _{ϵ}* is smaller than its equilibrium counterpart. The opposite is true when turbulence becomes locally weaker. It is to note that the value

*C*= 1.8

_{ϵ}*C*

_{ϵ}_{0}remains an upper bound also in the presence of negative buoyancy flux (positive Ri

*), which follows from the fact that (*

_{f}*A*

_{1}−

*A*

_{3}Ri

*) in Eq. (16) remains positive for the maximal value of Ri*

_{f}*= 0.5. It could, however, be possible that some limited regions of zero turbulence production*

_{f}*P*and negative buoyancy

*G*< 0 would be present. We can expect turbulence would be in a strong nonequilibrium, decaying fast in such zones.

Larger values of *C _{ϵ}* were obtained, e.g., in the study of Nedić et al. (2017) from direct numerical simulations data of boundary layer flows at relatively low Reynolds number. However, the derivations of Bos and Rubinstein (2017) were performed for simplified high-Re form of the spectrum, cf. Fig. 1. Hence, the results may not be valid at low Re, where the dissipative part of the spectrum is nonnegligible. Low values of Re

*are unlikely to be found in the atmospheric boundary layer turbulence; hence, we would rather expect*

_{λ}*C*to remain bounded from above by 1.8

_{ϵ}*C*

_{ϵ}_{0}.

## 3. Data and methods

To verify, whether the presented above analysis can be applied to atmospheric turbulence, we used high-resolution data collected in July 2017 during the ACORES campaign in the eastern North Atlantic around the island of Graciosa (Siebert et al. 2021). Measurements were performed with the Airborne Cloud Turbulence Observation System (ACTOS) (Siebert et al. 2006) mounted 170 m below the BO-105 helicopter. We analyze the vertical wind velocity *w* corrected for platform motion and attitude (cf. Nowak et al. 2021) provided by the ultrasonic anemometer–thermometer (Siebert and Muschinski 2001). The sensor was mounted on ACTOS and has the sampling frequency of 100 Hz. The standard deviations due to uncorrelated noise for velocity measurements are 0.002 m s^{−1}. Relatively low velocity of the helicopter (∼20 m s^{−1}) and high frequency of the sensor allows us to reliably resolve small scales of turbulence down to ∼0.5 m. The further analysis is based on the vertical velocity component only, as the horizontal components still showed some influence of platform attitude and motion which could not be completely removed with standard procedures. These issues caused problems with Reynolds decomposition and subsequent analysis. We note that Obligado et al. (2022) compared estimates of *C _{ϵ}* with

Helicopter flights during ACORES were performed over the ocean inside the 10 km × 10 km square adjacent to the Graciosa island. The typical flight duration was 2 h. The trajectory included vertical profiles up to 2000 m and horizontal legs of the length of several kilometers. In our analysis we use data from the horizontal segments for the same flights as selected in Nowak et al. (2021): flight 5 on 8 July 2017 and flight 14 on 18 July 2017, distinctive by the stratocumulus presence and STBL stratification (considerably well mixed in flight 5, considerably decoupled in flight 14). Detailed analysis of turbulence statistics based on these data was performed in Nowak et al. (2021). The study of nondimensional dissipation coefficient and

### a. Modified spectra

Due to finite frequency of ACTOS turbulence sensor, the wind velocity is measured with a spectral cutoff. As mentioned, ACORES data have the spatial resolution of ∼0.5 m, which is still larger than the typical size of the smallest Kolmogorov’s eddies (∼0.001 m to ∼0.01 m). With this, part of the energy spectrum remains unresolved and to recover value of *ϵ* indirect methods, which rely on the Kolmogorov’s hypotheses, should be used.

In this subsection we show that in spite of the presence of the nonequilibrium correction, the turbulence kinetic energy dissipation rate can still be calculated from the one-dimensional spectra by fitting to the −5/3 slope, provided that the fitting range is moved toward larger wavenumbers. This is possible, because the nonequilibrium correction affects the one-dimensional spectra only slightly and its influence is negligible for large wavenumbers. Hence, we can still estimate *ϵ* from the measured velocity signals even in the presence of the spectral cutoff at the Nyquist frequency *f _{s}*/2 = 50 Hz.

*E*(

*κ*,

*t*) is described in Eq. (6) and its equilibrium and nonequilibrium parts are given in Eqs. (7) and (8):

*P*/

*ϵ*and

*G*/

*ϵ*and the energy spectrum can be determined if

*k*

_{0}and

*ϵ*

_{0}are known. Plots of the one-dimensional spectra

*E*

_{11}, calculated numerically from Eq. (17) after substituting Eq. (18), corresponding to the case where

*P*/

*ϵ*= 2,

*G*/

*ϵ*= 0 (developing turbulence with the production twice larger than the dissipation) and

*P*/

*ϵ*= 0,

*G*/

*ϵ*= 0 (decaying turbulence, no turbulence production) are presented in Fig. 3. We used

*k*

_{0}= 0.06 m

^{2}s

^{−2}and

*ϵ*

_{0}= 5 × 10

^{−4}m

^{2}s

^{−3}, which are sample values (to the order of magnitude) measured in the STBL.

One-dimensional spectra calculated for the case *P*/*ϵ* = 2, *G*/*ϵ* = 0 (developing turbulence with the production twice larger than the dissipation): solid red line; *P*/*ϵ* = 0, *G*/*ϵ* = 0 (decaying turbulence, no turbulence production): dot–dashed blue line; compared with K41 spectrum: dashed black line.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

One-dimensional spectra calculated for the case *P*/*ϵ* = 2, *G*/*ϵ* = 0 (developing turbulence with the production twice larger than the dissipation): solid red line; *P*/*ϵ* = 0, *G*/*ϵ* = 0 (decaying turbulence, no turbulence production): dot–dashed blue line; compared with K41 spectrum: dashed black line.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

One-dimensional spectra calculated for the case *P*/*ϵ* = 2, *G*/*ϵ* = 0 (developing turbulence with the production twice larger than the dissipation): solid red line; *P*/*ϵ* = 0, *G*/*ϵ* = 0 (decaying turbulence, no turbulence production): dot–dashed blue line; compared with K41 spectrum: dashed black line.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

*C*

_{11}≈ 0.49. Linear least squares fit procedure in the selected range of wavenumbers is performed to estimate the coefficient and calculate

*ϵ*. The above formula is no longer true in the case of modified spectra; however, if the fitting range is moved toward larger wavenumbers, the influence of

*κ*than

*E*

_{0}. We used here the range of wavenumbers [1 − 10] m

^{−1}. Our reference case will be

*P*=

*ϵ*,

*G*= 0, for which we assume that

*E*

_{11}follows Eq. (19). For this case the linear least squares fit procedure provides

*ϵ*

_{ref}= 4.55 × 10

^{−4}m

^{2}s

^{−3}, which is somewhat smaller than the value prescribed a priori (equal to

*ϵ*

_{0}= 5 × 10

^{−4}m

^{2}s

^{−3}). This follows from the fact that Eq. (14) predicts

*P*=

*ϵ*. We next modified

*P*/

*ϵ*and

*G*/

*ϵ*and calculated the corresponding

*ϵ*in the same range. Results of

*ϵ*/

*ϵ*

_{ref}and

*C*/

_{ϵ}*C*

_{ϵ}_{0}determined from Eq. (15) are given in Table 1. As it is seen, even in the case of production being locally 3 times larger than the dissipation, the overprediction of

*ϵ*is around 10%. Analogously, in case of no production

*ϵ*is underpredicted by 7%. Most of the values of

*C*/

_{ϵ}*C*

_{ϵ}_{0}, calculated for ACORES data, are within the range 0.6

*C*

_{ϵ}_{0}–1.8

*C*

_{ϵ}_{0}(see section 4). Hence, as it follows from the above estimation, we expect the over/underprediction of

*ϵ*due to modification of the spectra should be less than 10% for the chosen fitting range. To minimize this error it is desirable to move the fitting range toward possibly largest resolved wavenumbers, but such that the effect of aliasing due to finite sampling frequency is still negligible. Wacławczyk et al. (2020) argued that

*ϵ*estimates based on the structure function approach are better in such case. This is also confirmed in the experience of Siebert et al. (2006) and Nowak et al. (2021). For this reason, in the subsequent analysis we estimate

*ϵ*using the transverse second-order structure function approach, by fitting the data to the formula

*D*

_{2}≈

*C*

_{2}(

*ϵr*)

^{2/3}, where

*C*

_{2}= 2.8. Figure 4 presents structure functions estimated from airborne measurements in the decoupled STBL. The wind velocity was recorded at height 230 m in a region of large turbulence production and at 990 m (in a region inside the cloud with locally weak turbulence production). As it is seen

*S*

_{2}profiles clearly deviate from the Kolmogorov’s scaling at high

*r*.

One-dimensional transverse structure functions estimated from airborne measurements, from the vertical wind velocity component recorded at height 280 m (region of large turbulence production): solid red line; at height 990 m (region inside the cloud with locally weak turbulence production): dot–dashed blue line; compared with K41 spectrum: dashed black line. Vertical lines mark the fitting range.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

One-dimensional transverse structure functions estimated from airborne measurements, from the vertical wind velocity component recorded at height 280 m (region of large turbulence production): solid red line; at height 990 m (region inside the cloud with locally weak turbulence production): dot–dashed blue line; compared with K41 spectrum: dashed black line. Vertical lines mark the fitting range.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

One-dimensional transverse structure functions estimated from airborne measurements, from the vertical wind velocity component recorded at height 280 m (region of large turbulence production): solid red line; at height 990 m (region inside the cloud with locally weak turbulence production): dot–dashed blue line; compared with K41 spectrum: dashed black line. Vertical lines mark the fitting range.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

Parameters estimate for different values of *P*/*ϵ* and *G*/*ϵ*.

### b. Characteristic time scales

*τ*defined above is a measure of nonequilibrium. The equilibrium time scale of turbulence will be defined by

*P*/

*ϵ*and

*G*/

*ϵ*as prescribed in Table 1. Results are presented in Table 1, last column. As it is seen, the further the system is from its stationary state, the smaller the time scale

*τ*becomes. Hence,

*τ*can be interpreted as the time scale of turbulence “adjustment” toward equilibrium.

### c. Effects of Reynolds decomposition

To calculate turbulence statistics the instantaneous velocity should first be decomposed into the mean and fluctuating parts. Here, the ensemble average is approximated by the space average along the flight track. The size of averaging window used for detrending, AW* _{D}*, should be much larger than the characteristic turbulence length scale but much smaller than the length associated with the mean changes. In practice, in the atmospheric turbulence, the presence of large-scale convective motions, internal waves, and changes of atmospheric conditions along the flight track makes the choice of the averaging window difficult. In this work we perform the analysis for vertical velocity components, for which the changes of the mean component were of lesser magnitude. In Nowak et al. (2021), AW

*= 50 s, which correspond to, approximately 1 km length, were used for the horizontal segments. Such window was approximately 10 times larger than the estimated integral turbulence length scales which were of order 100 m.*

_{D}*affects mostly the estimation of the integral length scale and the turbulence kinetic energy. We investigated how the values of dimensionless dissipation*

_{D}*C*determined from the formula

_{ϵ}*. The integral length scale*

_{D}*ϵ*was estimated from the second-order structure function in the fitting range [0.6–6] m, which corresponds to the range of wavenumbers considered in section 3a. Finally

*w*.

The analysis is performed for the subcloud segment F05LEG307 from flight 5, characterized by low values of *C _{ϵ}*, which suggests nonequilibrium, developing turbulence, a part of subcloud segment F14LEG287 from flight 14, where turbulence is close to equilibrium and cloud segment F14LEG992 from flight 14 with

*C*>

_{ϵ}*C*

_{ϵ}_{0}. For the first two segments the conditions remain approximately constant along the flight track. In the third segment, larger variations of turbulence properties along the flight track are observed; however, the mean value of

*C*still exceeds the equilibrium value, see section 6. After detrending another averaging window AW

_{ϵ}*should be chosen for the calculation of statistics. This window should be sufficiently long to reduce the bias and the random errors to acceptable levels (Lenschow et al. 1994). In the following subsection, we estimate the size of AW*

_{S}*sufficient to estimate*

_{S}*C*with a good accuracy. Here, we perform averaging over the whole considered parts of the segments, that is, AW

_{ϵ}*= 270 s for flight 5 and AW*

_{S}*= 300 s for both segments of flight 14. These lengths correspond to at least*

_{S}Results are presented in Table 2. As can be seen, the larger AW* _{D}* is chosen, the larger resulting

*even further we would get slower growth and eventually convergence of*

_{D}*, changes of*

_{D}*C*are much smaller. E.g., once AW

_{ϵ}*is reduced from 75 to 50 s,*

_{D}*C*changes by 0.01 or 0.02. If AW

_{ϵ}*is further reduced, the value of*

_{D}*C*becomes close to equilibrium. This would be expected, as both

_{ϵ}*is too much reduced. We decided to use AW*

_{D}*= 50 s, the same as in Nowak et al. (2021), for further analysis. Such value still allows us to detect differences in*

_{D}*C*levels, and, with moderate

_{ϵ}*, could be taken small enough to observe variations of turbulence properties along flight tracks.*

_{S}Estimates of turbulence statistics for different AW* _{D}*.

### d. Effect of averaging windows

Keeping AW* _{D}* = 50s constant we next investigated influence of the results to the choice of AW

*, that is, the window used to calculate statistics from the previously detrended signal. The analysis is performed for the first half of F14LEG287 from flight 14 in the decoupled STBL, where*

_{S}*C*remains approximately constant. The largest AW

_{ϵ}*= 200 s corresponds to, approximately*

_{S}*, as otherwise changes of*

_{S}*C*are averaged out. On the other hand, AW

_{ϵ}*should be large enough to calculate statistics with a good accuracy. We gradually decreased AW*

_{S}*to 150, 100, and 50 s which corresponds to, approximately*

_{S}*C*) for the investigated segment increase with the decreasing AW

_{ϵ}*and equal, respectively, std(*

_{S}*C*) = 0.036, 0.045, 0.064, and 0.21. The last number is clearly unacceptable (which is also visible in Fig. 5), as it should be possible to detect changes of

_{ϵ}*C*on the order of 0.1 to draw conclusions about the state of turbulence in the investigated region. We decided to further use such averaging window AW

_{ϵ}*which correspond to the length larger than*

_{S}*C _{ϵ}*/

*C*

_{ϵ}_{0}calculated for different AW

*from part of signal from flight 14, F14LEG287. Horizontal coordinate*

_{S}*x*denotes the position of the center of the averaging window.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

*C _{ϵ}*/

*C*

_{ϵ}_{0}calculated for different AW

*from part of signal from flight 14, F14LEG287. Horizontal coordinate*

_{S}*x*denotes the position of the center of the averaging window.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

*C _{ϵ}*/

*C*

_{ϵ}_{0}calculated for different AW

*from part of signal from flight 14, F14LEG287. Horizontal coordinate*

_{S}*x*denotes the position of the center of the averaging window.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

## 4. Results: Analysis of turbulence states in STBL

We first address and compare two horizontal flight segments of comparable, low altitudes of 287 m (decoupled STBL, F14LEG287) and 307 m (coupled STBL, F05LEG307). In the former case the averaging window AW* _{S}* = 150 s was used, in the latter, the window was increased to AW

*= 170 s due to larger length scales. The window was moved every 5 s (≈100 m) along the flight track and each time turbulence statistics were calculated. Figure 6a presents values of*

_{S}*C*as a function of the position of the center of the averaging window. In the figure two parallel lines of a constant

_{ϵ}*C*

_{ϵ}_{0}= 0.45 and

*C*= 0.8 are plotted. The first one would correspond to the equilibrium, stationary case. We note that this equilibrium value may in fact differ slightly for different flows, as it is weakly influenced by the large turbulence structures and their scaling (Valente and Vassilicos 2012; Bos et al. 2007). This influence was not taken into account in the theoretical considerations in section 2, as it was assumed that the energy density has sharp spectral cutoffs at

_{ϵ}*κ*and

_{L}*κ*. Moreover, possible flow anisotropy may also play a role and modify the proportionality constant. We choose

_{η}*C*

_{ϵ}_{0}= 0.45 as it is the value obtained for F14LEG287 at small AW

*= 15 s, see Table 2, for which nonequilibrium effects were mostly filtered out, as it was discussed in section 5d. The second value,*

_{D}*C*= 0.8 corresponds to a freely decaying turbulence with

_{ϵ}*P*= 0, where the equilibrium, self-similar state was reached (see Table 1).

*C _{ϵ}* as a function of the position of the center of the averaging window. Red rectangles mark nonequilibrium region. Data for decoupled STBL: (a) F14LEG287, (b) F14LEG448, (c) F14LEG992.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

*C _{ϵ}* as a function of the position of the center of the averaging window. Red rectangles mark nonequilibrium region. Data for decoupled STBL: (a) F14LEG287, (b) F14LEG448, (c) F14LEG992.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

*C _{ϵ}* as a function of the position of the center of the averaging window. Red rectangles mark nonequilibrium region. Data for decoupled STBL: (a) F14LEG287, (b) F14LEG448, (c) F14LEG992.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As it is seen in Fig. 6a, initially, *C _{ϵ}* ≈

*C*

_{ϵ}_{0}. As the averaging window is moved further along the flight track, a moderate increase of

*C*is first observed, followed by a sharp decrease toward smaller values. It would suggest that a weak turbulence decay, followed by a region of strong turbulence production was detected. In the analogous figure for the coupled STBL, Fig. 7a, only

_{ϵ}*C*<

_{ϵ}*C*

_{ϵ}_{0}are detected, which can be again interpreted as a region of strong turbulence production.

As in Fig. 6, but for coupled STBL: (a) F05LEG307, (b) F05LEG553, (c) F05LEG819.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 6, but for coupled STBL: (a) F05LEG307, (b) F05LEG553, (c) F05LEG819.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 6, but for coupled STBL: (a) F05LEG307, (b) F05LEG553, (c) F05LEG819.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

We next compared two subcloud flight segments: F14LEG448 in the decoupled STBL of height 448 m and F05LEG553 in the coupled STBL of height 553 m. In Nowak et al. (2021) nearly zero buoyancy production at those levels was reported. The weak turbulence production is manifested through larger average values of *C _{ϵ}*, both in decoupled and coupled case, cf. Figs. 6b and 7b, respectively. Insufficient turbulence production at this altitude may lead to decoupling, when the eddies fail to mix air over the entire depth of the STBL, and consequently STBL separates into two parts: cloud driven at the top and surface driven at the bottom, cf. Nowak et al. (2021).

Finally, we consider in-cloud segments F14LEG992 in the decoupled STBL and F05LEG819 in the coupled STBL. The signals were recorded at heights 992 and 819 m. In both cases values of *C _{ϵ}* vary considerably, cf. Figs. 6c and 7c; however, they are somewhat larger in the decoupled case which indicates weaker turbulence production and is again in line with the observations of Nowak et al. (2021). The results of

*C*can be compared with experimental estimations of Zheng et al. (2021) in high-Re turbulence behind a grid, where

_{ϵ}*C*first decreased in the production region and next increased in the nonequilibrium decay region toward value

_{ϵ}*C*≈ 0.8. We also mention here the direct numerical simulation study of Gallana et al. (2022), who calculated

_{ϵ}*C*in stably and unstably stratified flows. In the former case

_{ϵ}*C*increased toward

_{ϵ}*C*≈ 0.8, indicating turbulence suppression and in the unstable conditions decreased toward

_{ϵ}*C*≈ 0.4, below the equilibrium value. In STBL turbulence is manly driven by radiative and evaporative cooling at the cloud top, leading to instability. However, boundary layer is also capped from above by a stably stratified layer. This may explain large variations of

_{ϵ}*C*along the in-cloud flight tracks.

_{ϵ}To classify whether turbulence is in its equilibrium or nonequilibrium state, we calculated the length scales and the local Reynolds numbers and plotted *C _{ϵ}* as a function of Re

*. The data points were compared with the theoretical predictions (lines). For*

_{λ}*C*, in the equilibrium case in high-Re atmospheric flows, we should obtain

_{ϵ}*C*= const. In the nonequilibrium, dependence of

_{ϵ}*C*on Re

_{ϵ}*should be described by Eq. (11). As we do not have the knowledge about the equilibrium value Re*

_{λ}

_{λ}_{0}, the coefficient in Eq. (11) was calculated from the linear least squares algorithm. In the lowest flight segment F14LEG287 in decoupled STBL statistics calculated at

*x*< 4500 m follow, approximately the equilibrium line

*C*= const (black dots in Fig. 8a). However, when the position of the center of the averaging window exceeds 4500 m the nonequilibrium states can be detected (red stars in Fig. 8a). The values of

_{ϵ}*C*calculated for coupled STBL for F05LEG307 decrease slightly with Re

_{ϵ}*, which suggests the presence of nonequilibrium scaling (cf. Fig. 9a).*

_{λ}*C _{ϵ}* as a function of Re

*. Solid black lines: equilibrium scalings*

_{λ}*C*

_{ϵ}_{0}= 0.45 and

*C*

_{ϵ}_{1}= 1.8

*C*

_{ϵ}_{0}; dashed red line: nonequilibrium scaling, Eq. (11); calculated statistics: symbols. Decoupled STBL: (a) F14LEG287, (b) F14LEG448, (c) F14LEG992.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

*C _{ϵ}* as a function of Re

*. Solid black lines: equilibrium scalings*

_{λ}*C*

_{ϵ}_{0}= 0.45 and

*C*

_{ϵ}_{1}= 1.8

*C*

_{ϵ}_{0}; dashed red line: nonequilibrium scaling, Eq. (11); calculated statistics: symbols. Decoupled STBL: (a) F14LEG287, (b) F14LEG448, (c) F14LEG992.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

*C _{ϵ}* as a function of Re

*. Solid black lines: equilibrium scalings*

_{λ}*C*

_{ϵ}_{0}= 0.45 and

*C*

_{ϵ}_{1}= 1.8

*C*

_{ϵ}_{0}; dashed red line: nonequilibrium scaling, Eq. (11); calculated statistics: symbols. Decoupled STBL: (a) F14LEG287, (b) F14LEG448, (c) F14LEG992.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 8, but for coupled STBL: (a) F05LEG307, (b) F05LEG553, (c) F05LEG819. In (c) color coded are different local equilibrium states.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 8, but for coupled STBL: (a) F05LEG307, (b) F05LEG553, (c) F05LEG819. In (c) color coded are different local equilibrium states.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 8, but for coupled STBL: (a) F05LEG307, (b) F05LEG553, (c) F05LEG819. In (c) color coded are different local equilibrium states.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

Figures 8b and 9b refer to the flight segments F14LEG448 and F05LEG553 below the cloud. Figure 8b suggests that part of the signal in the decoupled STBL is a record of equilibrium turbulence, while statistics of another part fit the nonequilibrium formulas. For the coupled STBL the statistics follow rather the equilibrium predictions, cf. Fig. 9b. A large scatter of results is visible in Fig. 9b such that it is not easy to draw conclusions. Note that the flight segment F05LEG553 is relatively short and so is the range of the detected Re* _{λ}* values. Equilibrium scaling may suggest that in spite of the weak turbulence production, the changes of atmospheric conditions are relatively slow, such that the turbulence spectra already relaxed to the self-similar form (5). Similar, nearly constant levels of

*C*were reported, e.g., by Neunaber et al. (2022) in the study of turbulence in wakes behind turbines.

_{ϵ}Particularly interesting is the study of the in-cloud segments F14LEG992 and F05LEG819. In the decoupled STBL most of the points follow the nonequilibrium scaling relations, cf. Fig. 8c. Exceptions are several points corresponding to large *C _{ϵ}* ≈ 1.8

*C*

_{ϵ}_{0}, which would indicate regions of no turbulence production and/or stable stratification. In the coupled STBL on the other hand, in spite of large variations of

*C*, the points seem to follow three different lines of constant

_{ϵ}*C*, namely,

_{ϵ}*C*≈ 0.8,

_{ϵ}*C*≈ 0.6, and

_{ϵ}*C*≈ 0.4, cf. Fig. 9c. We may call such state of the system as a quasi-equilibrium state, which means that the changes of external conditions are slow enough, such that the system undergoes through several local equilibria.

_{ϵ}We note here that for certain points the assignment to equilibrium or nonequilibrium region may be somewhat arbitrary. This concerns especially those points which represent statistics calculated at the boundary between the equilibrium and nonequilibrium region. Sometimes this assignment becomes more definite when the dependence of the length scale ratio * _{λ}* is studied, see Figs. 10 and 11. In the equilibrium case, the statistics should follow Eq. (2), while in the nonequilibrium

*are inversely proportional and follow Eq. (12). Both functional dependencies are plotted in Figs. 10 and 11 as the black and red lines, respectively. We note that we did not perform another curve fitting here, the constants are those estimated for data from Figs. 8 and 9 for respective flight segments and divided by 15, as follows from the Eqs. (11) and (12). The good agreement confirms validity of the theoretical derivations of Bos and Rubinstein (2017). After the classification of statistics in Figs. 8–11 the nonequilibrium parts were marked in Figs. 6 and 7 (as red rectangles).*

_{λ}* _{λ}*. Solid black line: equilibrium scaling, Eq. (2); dashed red line: nonequilibrium scaling, Eq. (12); calculated statistics: symbols. Decoupled STBL: (a) F14LEG287, (b) F14LEG448, (c) F14LEG992.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

* _{λ}*. Solid black line: equilibrium scaling, Eq. (2); dashed red line: nonequilibrium scaling, Eq. (12); calculated statistics: symbols. Decoupled STBL: (a) F14LEG287, (b) F14LEG448, (c) F14LEG992.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

* _{λ}*. Solid black line: equilibrium scaling, Eq. (2); dashed red line: nonequilibrium scaling, Eq. (12); calculated statistics: symbols. Decoupled STBL: (a) F14LEG287, (b) F14LEG448, (c) F14LEG992.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 10, but for coupled STBL: (a) F05LEG307, (b) F05LEG553, (c) F05LEG819. In (c) color codes are different local equilibrium states.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 10, but for coupled STBL: (a) F05LEG307, (b) F05LEG553, (c) F05LEG819. In (c) color codes are different local equilibrium states.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 10, but for coupled STBL: (a) F05LEG307, (b) F05LEG553, (c) F05LEG819. In (c) color codes are different local equilibrium states.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

Nonequilibrium scaling relations were found in many laboratory studies, as examples we mention here the works of Zheng et al. (2021), which concerns high-Re flows behind a grid, and of Nedić et al. (2017), which concerns boundary layer flows.

*C*≠

_{ϵ}*C*

_{ϵ}_{0}can be distinguished. The part on the left, where

*C*>

_{ϵ}*C*

_{ϵ}_{0}suggests turbulence decay. Statistics calculated at the final part of the signal, from F14LEG143, with

*C*<

_{ϵ}*C*

_{ϵ}_{0}, imply on the other hand the dominant role of the production. On the middle part of the curve in Fig. 12,

*C*oscillates around the equilibrium value of 0.45. These equilibrium values are marked as black dots in Fig. 13a. Apparently, they do not depend on Re

_{ϵ}*. Moreover, the ratio*

_{λ}*, as predicted by the equilibrium relation (2), see Fig. 13b. However, the remaining points, marked as red and blue stars in Fig. 13, do not follow one nonequilibrium line. We can argue that they correspond to two distinct regions, with different initial conditions and different Re*

_{λ}

_{λ}_{0}. By using the least squares algorithm we calculated two coefficients of the nonequilibrium relation (11),

*A*refers to the points with

*C*>

_{ϵ}*C*

_{ϵ}_{0}(red stars in Fig. 13) and

*B*refers to the points with

*C*<

_{ϵ}*C*

_{ϵ}_{0}(blue stars in Fig. 13). We found Re

_{λ}_{0}

*≈ 1.2Re*

_{B}

_{λ}_{0}

*. Assuming*

_{A}*C*

_{ϵ}_{0}≈ 0.45 it is possible to recover the initial-state Re

_{λ}_{0}from the calculated

*A*and

*B*coefficients. We found Re

_{λ}_{0}

*≈ 7500 and Re*

_{A}

_{λ}_{0}

*≈ 9000 as the equilibrium (initial) Reynolds numbers in the respective regions. We next plotted the red points for rescaled coordinate*

_{B}As in Fig. 6, but for the decoupled STBL, F14LEG143.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 6, but for the decoupled STBL, F14LEG143.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 6, but for the decoupled STBL, F14LEG143.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

(a) *C _{ϵ}* as a function of Re

*, (b)*

_{λ}*. Decoupled STBL, F14LEG143. Solid black lines: equilibrium scaling; dashed red line and dot–dashed blue lines: nonequilibrium scaling, Eq. (11); calculated statistics: symbols.*

_{λ}Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

(a) *C _{ϵ}* as a function of Re

*, (b)*

_{λ}*. Decoupled STBL, F14LEG143. Solid black lines: equilibrium scaling; dashed red line and dot–dashed blue lines: nonequilibrium scaling, Eq. (11); calculated statistics: symbols.*

_{λ}Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

(a) *C _{ϵ}* as a function of Re

*, (b)*

_{λ}*. Decoupled STBL, F14LEG143. Solid black lines: equilibrium scaling; dashed red line and dot–dashed blue lines: nonequilibrium scaling, Eq. (11); calculated statistics: symbols.*

_{λ}Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 13, but for rescaled Re* _{λ}*.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 13, but for rescaled Re* _{λ}*.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

As in Fig. 13, but for rescaled Re* _{λ}*.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0028.1

## 5. Conclusions

The main purpose of this work was to show that certain information about the temporal variations of turbulence and the budget of the turbulence kinetic energy in the STBL can be acquired by analyzing statistics of wind velocity based on in situ measurement data. The key indicators are the nondimensional dissipation coefficient *C _{ϵ}* and the integral-to-Taylor-scale ratio

*C*takes largest values, up to around

_{ϵ}*C*≈ 0.8 when turbulence decays, reaches its equilibrium

_{ϵ}*C*≈

_{ϵ}*C*

_{ϵ}_{0}≈ 0.45 in the stationary state and becomes smaller

*C*<

_{ϵ}*C*

_{ϵ}_{0}when turbulence becomes stronger. We showed here that the upper limit

*C*≈ 1.8

_{ϵ}*C*

_{ϵ}_{0}should not be exceeded even under a stable stratification, due to the presence of the critical Richardson number, above which turbulence is suppressed.

Moreover, as it was discussed in the previous works (Bos and Rubinstein 2017) two different states of turbulence system can be identified: equilibrium, where the turbulence kinetic energy spectrum has a self-similar shape, (5), and a nonequilibrium where a correction to the spectrum should be considered; see Eq. (6). The equilibrium state can possibly be nonstationary; however, the changes of external conditions (forcing) are slow enough such that the system has enough time to relax and the energy spectrum takes the self-similar form (5). To detect the nonequilibrium states we studied dependence of *C _{ϵ}* and

*. In the equilibrium*

_{λ}*C*≈ const and

_{ϵ}*. In the nonequilibrium turbulence a different, although universal scaling is present, namely,*

_{λ}*C*/

_{ϵ}*C*

_{ϵ}_{0}.

In this work we focus on the analysis of airborne measurement data, which often suffer from various deficiencies, related to finite frequencies of the sensors or short available averaging windows (Wacławczyk et al. 2017; Akinlabi et al. 2019). We showed that despite these deficiencies both *C _{ϵ}* and

*ϵ*it is crucial to move the fitting range toward smaller scales (larger wavenumbers) to minimize the influence of nonequilibrium correction to the spectrum. This −7/3 correction is generally subdominant and influences larger scales. We also showed that the indicator functions can be estimated with a sufficient accuracy if the size of the averaging window used to calculate statistics is larger than

*C*and

_{ϵ}*as well as on the initial conditions Re*

_{λ}

_{λ}_{0}, which is shown by a proper rescaling of the data. Moreover, with a curve fitting of the data we were able to recover values of Re

_{λ}_{0}.

The analysis performed here allows for additional comparison of turbulence in the coupled and decoupled STBL’s, previously studied by Nowak et al. (2021). In both cases, *C _{ϵ}* is smaller at lower altitudes due to strong shear and/or buoyancy production. On the other hand, in the subcloud parts

*C*increases above

_{ϵ}*C*

_{ϵ}_{0}which suggests that in these regions turbulence decays. Large variations of

*C*were detected inside the cloud in both cases. Apparently, along the flight track, regions where turbulence production surpass or is in balance with the dissipation alternate with areas of turbulence decay. Here, the main difference between the coupled and decoupled STBL was observed. Namely, in the former, turbulence seems to undergo a series of equilibrium states, where

_{ϵ}*C*≈ const and

_{ϵ}*C*was higher in the decoupled than in the coupled case which indicates the weaker turbulence production in the decoupled STBL.

_{ϵ}In this work only limited amount of data from two ACTOS flights were investigated. In future, more thorough analysis to better characterize turbulence in the coupled and decoupled STBL is foreseen. The same analysis can be repeated for various measurement data in the atmospheric boundary layers where the impact of transient states is nonnegligible (cf. Mahrt and Bou-Zeid 2020), e.g., to study the collapse of turbulent convective daytime boundary layer (El Guernaoui et al. 2019; Wang et al. 2020; Lothon et al. 2014), in and between various types of the clouds in order to better understand their dynamics and even in clear air turbulence to characterize its dynamic properties.

Another interesting aspect is the relation of the present study to the K62 theory (Kolmogorov 1962), which takes into account temporal and spatial variations of the instantaneous dissipation rate. This phenomenon is known as the “internal intermittency” and affects the scaling of higher-order structure functions. A careful analysis of the fine-scale structure of turbulence performed by Siebert et al. (2010) based on ACTOS data revealed the intermittency coefficient in cloud turbulence is consistent with values obtained in laboratory experiments. A question to be asked is how the nonequilibrium structure functions are affected by the small-scale intermittency and whether temporal variations of statistics considered here can be linked to the premises of the K62 theory.

The presence of nonequilibrium in the stratocumulus clouds indicate that the common turbulence closures like the Smagorinsky LES model may fail to predict the dynamics of such clouds correctly. In particular, this may influence the prediction of boundary layer decoupling or cloud dissipation. A nonequilibrium turbulence model may largely improve these predictions. This is another direction of future studies.

## Acknowledgments.

This research was supported by the National Science Centre, Poland, Project 2020/37/B/ST10/03695 (theory and methodology), and by the European Union’s Horizon 2020 research and innovation programme under Grant Agreement 101003470, project Next Generation Earth Modelling Systems (NextGEMS) (analyses of measurement data of turbulence in the stratocumulus-topped boundary layer). The field campaign was supported by the Deutsche Forschungsgesellschaft (DFG; Grant SI 1543/4-1).

## Data availability statement.

All data from the ACORES study are available from the authors on request.

## APPENDIX A

### Derivation of Nonequilibrium Scaling Laws

*C*≈ 1.5. As a consequence, turbulence properties can also be decomposed into the equilibrium and nonequilibrium parts, that is,

_{K}*E*

_{0}by

*κ*≪

_{L}*κ*. Further, as

_{η}*C*and Re

_{ϵ}*in the following form:*

_{λ}*ϵ*≈

*ϵ*

_{0}due to the inequality

_{λ}_{0}above denotes the equilibrium Reynolds number. Next, after substituting (A10) and assuming that

*I*

_{0}is a first-order term of the Taylor-series expansion, the final formula for the nonequilibrium

*C*was obtained as follows:

_{ϵ}## APPENDIX B

### Transport Equations for the Turbulence Kinetic Energy and Dissipation Rate

*b*stands for the fluctuations of the buoyancy. The vertical component of the fluctuating velocity

*u*

_{3}will also be denoted by

*w*; hence, the alternative notation for the buoyancy flux is

*G*= ⟨

*wb*⟩.

*k*–

*ϵ*closures a gradient diffusion hypothesis is used to model the turbulent transport and the production terms, while another equation is solved for the turbulence kinetic energy dissipation rate (Jones and Launder 1972). This equation, adopted in Rodi (1987) for the study of buoyant flows reads

*ν*reads

_{T}*C*= 0.09,

_{μ}*A*

_{1}= 1.44,

*A*

_{2}= 1.92, and

*σ*= 1 are standard, the value of

_{ϵ}*A*

_{3}= 1.14 used in this work was proposed by Baum and Caponi (1992).

## APPENDIX C

### Estimation of the Integral Length Scale from Measurement Data

*L*denotes longitudinal vector component (i.e., along the flight track) and the subscript

*N*the transverse component (perpendicular to the flight track). In the isotropic turbulence

*f*(

*r*,

*t*) and

*g*(

*r*,

*t*) holds (Pope 2000):

*E*(

*κ*,

*t*) [see also Eq. (A5)]

*g*(

*r*,

*t*) described by the above formula is integrated over

*r*from 0 to

*g*(

*r*,

*t*) is calculated by averaging product of velocity components recorded at two points placed at a distance

*r*from each other. We gradually increase

*r*and obtain an approximate form of the function

*g*. Results are deteriorated for large

*r*, due to small size of the sample; however, we expect

*g*(

*r*,

*t*) is calculated with a good accuracy for small

*r*. Next we find the distance

*r*=

*r*

_{0}at which

*g*(

*r*,

*t*) crosses 0 for the first time. We integrate

*g*(

*r*,

*t*) from 0 to

*r*

_{0}. The calculated value is next divided by 0.57 and the result should be an estimation of the integral length scale

## REFERENCES

Akinlabi, E. O., M. Wacławczyk, J. P. Mellado, and S. P. Malinowski, 2019: Estimating turbulence kinetic energy dissipation rates in the numerically simulated stratocumulus cloud-top mixing layer: Evaluation of different methods.

,*J. Atmos. Sci.***76**, 1471–1488, https://doi.org/10.1175/JAS-D-18-0146.1.Antonia, R. A., and B. R. Pearson, 2000: Effect of initial conditions on the mean energy dissipation rate and the scaling exponent.

,*Phys. Rev. E***62**, 8086, https://doi.org/10.1103/PhysRevE.62.8086.Baum, E., and E. A. Caponi, 1992: Modeling the effects of buoyancy on the evolution of geophysical boundary layers.

,*J. Geophys. Res.***97**, 15 513–15 527, https://doi.org/10.1029/92JC01715.Bos, W. J. T., 2020: Production and dissipation of kinetic energy in grid turbulence.

,*Phys. Rev. Fluids***5**, 104607, https://doi.org/10.1103/PhysRevFluids.5.104607.Bos, W. J. T., and R. Rubinstein, 2017: Dissipation in unsteady turbulence.

,*Phys. Rev. Fluids***2**, 022601, https://doi.org/10.1103/PhysRevFluids.2.022601.Bos, W. J. T., L. Shao, and J.-P. Bertoglio, 2007: Spectral imbalance and the normalized dissipation rate of turbulence.

,*Phys. Fluids***19**, 045101, https://doi.org/10.1063/1.2714079.Burattini, P., P. Lavoie, and R. A. Antonia, 2005: On the normalized turbulent energy dissipation rate.

,*Phys. Fluids***17**, 098103, https://doi.org/10.1063/1.2055529.Cafiero, G., and J. C. Vassilicos, 2019: Non-equilibrium turbulence scalings and self-similarity in turbulent planar jet.

*Proc. Roy. Soc.**London*,**A475**, 20190038, https://doi.org/10.1098/rspa.2019.0038.Chen, J. G., and J. C. Vassilicos, 2022: Scalings of scale-by-scale turbulence energy in non-homogeneous turbulence.

,*J. Fluid Mech.***938**, A7, https://doi.org/10.1017/jfm.2022.153.El Guernaoui, O., J. Reuder, I. Esau, T. Wolf, and B. Maronga, 2019: Scaling the decay of turbulence kinetic energy in the free-convective boundary layer.

,*Bound.-Layer Meteor.***173**, 79–97, https://doi.org/10.1007/s10546-019-00458-z.Gallana, L., S. Abdunabiev, M. Golshan, and D. Tordella, 2022: Diffusion of turbulence following both stable and unstable step stratification perturbations.

,*Phys. Fluids***34**, 065122, https://doi.org/10.1063/5.0090042.Galperin, B., S. Sukoriansky, and P. S. Anderson, 2007: On the critical Richardson number in stably stratified turbulence.

,*Atmos. Sci. Lett.***8**, 65–69, https://doi.org/10.1002/asl.153.Grachev, A. A., E. L. Andreas, C. W. Fairall, P. S. Guest, and P. O. G. Persson, 2013: The critical Richardson number and limits of applicability of local similarity theory in the stable boundary layer.

,*Bound.-Layer Meteor.***147**, 51–82, https://doi.org/10.1007/s10546-012-9771-0.Jen-La Plante, I., and Coauthors, 2016: Physics of Stratocumulus Top (POST): Turbulence characteristics.

,*Atmos. Chem. Phys.***16**, 9711–9725, https://doi.org/10.5194/acp-16-9711-2016.Jones, W. P., and B. E. Launder, 1972: The prediction of laminarization with a two-equation model of turbulence.

,*Int. J. Heat Mass Transfer***15**, 301–314, https://doi.org/10.1016/0017-9310(72)90076-2.Kolmogorov, A. N., 1941: Dissipation of energy in the locally isotropic turbulence.

*Proc. Roy. Soc.**London*,**A434**, 15–17, https://doi.org/10.1098/rspa.1991.0076.Kolmogorov, A. N., 1962: A refinement of previous hypotheses concerning the local structure of turbulence in a viscous incompressible fluid at high Reynolds number.

,*J. Fluid Mech.***13**, 82–85, https://doi.org/10.1017/S0022112062000518.Lenschow, D. H., J. Mann, and L. Kristensen, 1994: How long is long enough when measuring fluxes and other turbulence statistics?

,*J. Atmos. Oceanic Technol.***11**, 661–673, https://doi.org/10.1175/1520-0426(1994)011<0661:HLILEW>2.0.CO;2.Lothon, M., and Coauthors, 2014: The BLLAST field experiment: Boundary-layer late afternoon and sunset turbulence.

,*Atmos. Chem. Phys.***14**, 10 931–10 960, https://doi.org/10.5194/acp-14-10931-2014.Mahrt, L., and E. Bou-Zeid, 2020: Non-stationary boundary layers.

,*Bound.-Layer Meteor.***177**, 189–204, https://doi.org/10.1007/s10546-020-00533-w.Malinowski, S. P., and Coauthors, 2013: Physics of Stratocumulus Top (POST): Turbulent mixing across capping inversion.

,*Atmos. Chem. Phys.***13**, 12 171–12 186, https://doi.org/10.5194/acp-13-12171-2013.Mazellier, N., and J. C. Vassilicos, 2008: The turbulence dissipation constant is not universal because of its universal dependence on large-scale flow topology.

,*Phys. Fluids***20**, 015101, https://doi.org/10.1063/1.2832778.McComb, W. D., A. Berera, M. Salewski, and S. Yoffe, 2010: Taylor’s (1935) dissipation surrogate reinterpreted.

,*Phys. Fluids***22**, 061704, https://doi.org/10.1063/1.3450299.Nedić, J., S. Tavoularis, and I. Marusic, 2017: Dissipation scaling in constant-pressure turbulent boundary layers.

,*Phys. Rev. Fluids***2**, 032601, https://doi.org/10.1103/PhysRevFluids.2.032601.Neunaber, I., J. Peinke, and M. Obligado, 2022: Application of the Townsend–George theory for free shear flows to single and double wind turbine wakes—A wind tunnel study.

,*Wind Energy Sci.***7**, 201–219, https://doi.org/10.5194/wes-7-201-2022.Nowak, J. L., H. Siebert, K.-E. Szodry, and S. P. Malinowski, 2021: Coupled and decoupled stratocumulus-topped boundary layers: Turbulence properties.

,*Atmos. Chem. Phys.***21**, 10 965–10 991, https://doi.org/10.5194/acp-21-10965-2021.Obligado, M., T. Dairay, and J. C. Vassilicos, 2016: Nonequilibrium scalings of turbulent wakes.

,*Phys. Rev. Fluids***1**, 044409, https://doi.org/10.1103/PhysRevFluids.1.044409.Obligado, M., C. Brun, J. H. Silvestrini, and E. B. C. Schettini, 2022: Dissipation scalings in the turbulent boundary layer at moderate

*Reθ.*,*Flow Turbul. Combust.***108**, 105–122, https://doi.org/10.1007/s10494-021-00270-1.Pope, S. B., 2000:

1st ed. Cambridge University Press, 802 pp.*Turbulent Flows.*Rodi, W., 1987:

IAHR, 104 pp.*Turbulence Models and Their Application in Hydraulics: A State of the Art Review.*Rubinstein, R., and T. T. Clark, 2017: “Equilibrium” and “non-equilibrium” turbulence.

,*Theor. Appl. Mech. Lett.***7**, 301–305, https://doi.org/10.1016/j.taml.2017.09.010.Seoud, R. E., and J. C. Vassilicos, 2007: Dissipation and decay of fractal-generated turbulence.

,*Phys. Fluids***19**, 105108, https://doi.org/10.1063/1.2795211.Siebert, H., and A. Muschinski, 2001: Relevance of a tuning-fork effect for temperature measurements with the Gill Solent HS ultrasonic anemometer–thermometer.

,*J. Atmos. Oceanic Technol.***18**, 1367–1376, https://doi.org/10.1175/1520-0426(2001)018<1367:ROATFE>2.0.CO;2.Siebert, H., K. Lehmann, and M. Wendisch, 2006: Observations of small-scale turbulence and energy dissipation rates in the cloudy boundary layer.

,*J. Atmos. Sci.***63**, 1451–1466, https://doi.org/10.1175/JAS3687.1.Siebert, H., R. A. Shaw, and Z. Warhaft, 2010: Statistics of small-scale velocity fluctuations and internal intermittency in marine stratocumulus clouds.

,*J. Atmos. Sci.***67**, 262–273, https://doi.org/10.1175/2009JAS3200.1.Siebert, H., and Coauthors, 2021: Observations of aerosol, cloud, turbulence, and radiation properties at the top of the marine boundary layer over the eastern North Atlantic Ocean: The ACORES campaign.

,*Bull. Amer. Meteor. Soc.***102**, E123–E147, https://doi.org/10.1175/BAMS-D-19-0191.1.Stevens, B., and Coauthors, 2003: Dynamics and Chemistry of Marine Stratocumulus—DYCOMS-II.

,*Bull. Amer. Meteor. Soc.***84**, 579–594, https://doi.org/10.1175/BAMS-84-5-579.Taylor, G. I., 1935: Statistical theory of turbulence.

,*Proc. Roy. Soc. London***A151**, 421–444, https://doi.org/10.1098/rspa.1935.0158.Valente, P. C., and J. C. Vassilicos, 2011: The decay of turbulence generated by a class of multiscale grids.

,*J. Fluid Mech.***687**, 300–340, https://doi.org/10.1017/jfm.2011.353.Valente, P. C., and J. C. Vassilicos, 2012: Universal dissipation scaling for nonequilibrium turbulence.

,*Phys. Rev. Lett.***108**, 214503, https://doi.org/10.1103/PhysRevLett.108.214503.Valente, P. C., and J. C. Vassilicos, 2015: The energy cascade in grid-generated non-equilibrium decaying turbulence.

,*Phys. Fluids***27**, 045103, https://doi.org/10.1063/1.4916628.Vassilicos, J. C., 2015: Dissipation in turbulent flows.

,*Annu. Rev. Fluid Mech.***47**, 95–114, https://doi.org/10.1146/annurev-fluid-010814-014637.Wacławczyk, M., Y.-F. Ma, J. M. Kopeć, and S. P. Malinowski, 2017: Novel approaches to estimating the turbulent kinetic energy dissipation rate from low- and moderate-resolution velocity fluctuation time series.

,*Atmos. Meas. Tech.***10**, 4573–4585, https://doi.org/10.5194/amt-10-4573-2017.Wacławczyk, M., A. S. Gozingan, J. Nzotungishaka, M. Mohammadi, and S. P. Malinowski, 2020: Comparison of different techniques to calculate properties of atmospheric turbulence from low-resolution data.

,*Atmosphere***11**, 199, https://doi.org/10.3390/atmos11020199.Wacławczyk, T., 2021: Modeling of non-equilibrium effects in intermittency region between two phases.

,*Int. J. Multiphase Flow***134**, 103459, https://doi.org/10.1016/j.ijmultiphaseflow.2020.103459.Wacławczyk, T., 2022: On differences between deterministic and statistical models of the interphase region.

,*Acta Mech. Sin.***38**, 722045, https://doi.org/10.1007/s10409-022-22045-w.Wang, D., I. S. Stachlewska, X. Song, B. Heese, and A. Nemuc, 2020: Variability of the boundary layer over an urban continental site based on 10 years of active remote sensing observations in Warsaw.

,*Remote Sens.***12**, 340, https://doi.org/10.3390/rs12020340.Wood, R., and Coauthors, 2015: Clouds, aerosols, and precipitation in the marine boundary layer: An ARM Mobile Facility deployment.

,*Bull. Amer. Meteor. Soc.***96**, 419–440, https://doi.org/10.1175/BAMS-D-13-00180.1.Yang, P.-F., A. Pumir, and H. Xu, 2018: Generalized self-similar spectrum and the effect of large-scale in decaying homogeneous isotropic turbulence.

,*New J. Phys.***20**, 103035, https://doi.org/10.1088/1367-2630/aae72d.Zheng, S., P. J. K. Bruce, C. Cuvier, J.-M. Foucaut, J. M. R. Graham, and J. C. Vassilicos, 2021: Nonequilibrium dissipation scaling in large Reynolds number turbulence generated by rectangular fractal grids.

,*Phys. Rev. Fluids***6**, 054613, https://doi.org/10.1103/PhysRevFluids.6.054613.Zilitinkevich, S., and A. Baklanov, 2002: Calculation of the height of the stable boundary layer in practical applications.

,*Bound.-Layer Meteor.***105**, 389–409, https://doi.org/10.1023/A:1020376832738.