Articles | Volume 16, issue 2
Research article
20 Jan 2023
Research article |  | 20 Jan 2023

Detection of turbulence occurrences from temperature, pressure, and position measurements under superpressure balloons

Richard Wilson, Clara Pitois, Aurélien Podglajen, Albert Hertzog, Milena Corcos, and Riwal Plougonven

This article deals with the detection of small-scale turbulence from in situ meteorological measurements performed under superpressure balloons (SPBs). These balloons allow long-duration flights (several months) at a prerequisite height level. The data set is gathered from the Strateole-2 probationary campaign during which eights SPBs flew in the tropical tropopause layer at altitudes of around 19 and 20.5 km from November 2019 to March 2020.

Turbulence is not directly measured by the instrument set onboard the SPBs. Nonetheless, there is the potential to derive information about the occurrence of turbulence from the temporally well-resolved measurements of pressure, temperature, and position. It constitutes a challenge to extract the aforementioned information from a measurement set that was not designed for quantifying turbulence, and the paper explains the methodology developed to overcome this difficulty.

It is observed that SPBs oscillate quasi-periodically around their equilibrium positions. The oscillation periods, which are 220 s on average with a range of 130 to 500 s, are close to but noticeably smaller than the Brunt–Väisälä period (∼300 s). The amplitude of these vertical motions is ±15 m, inducing large fluctuations in all quantities, whether measured (e.g., pressure, temperature and position) or inferred (e.g., density and potential temperature). The relationships between the changes in these quantities and the vertical displacements of the balloons are used to infer properties of the flow in which the SPBs drift.

In the case of active turbulence, the vertical stratification as well as the wind shear are likely to be reduced by mixing. Hence, the increments of potential temperature, δθ, and of the vertical displacements of the balloon, δzB, are expected to be uncorrelated because θ/z0. Moreover, the local vertical gradients of measured quantities, temperature (T) and horizontal velocities (u and v), are estimated from the covariance of the increments of the considered quantity with δzB. The Richardson number of the flow is deduced.

Several binary indexes (true or false) to describe the state of the flow, laminar or turbulent, are evaluated. These turbulence indexes, based either on correlations between δθ and δzB or on estimates of the local Richardson number, are found to be consistent, as they differ in less than 3 % of cases. The flow is observed to be turbulent for about 5 % of the time, with strong inhomogeneities along the longitude.

1 Introduction

The vertical transport of heat, momentum, and minor constituents in the tropical upper troposphere–lower stratosphere (UTLS) is an important issue, as this region is recognized as the gateway for tropospheric air into the stratosphere (Fueglistaler et al.2009). Above 15 km altitude, this vertical transport is believed to mainly result either from the mean tropical upwelling associated with the Brewer–Dobson circulation or from small-scale turbulence. The relative contribution of turbulent mixing to this vertical transport is highly uncertain, partly owing to the lack of observations. Turbulence observations in the tropical UTLS are indeed sparse: they mostly come from two large very high frequency (VHF) radars, from relatively few radiosondes, or from research aircraft.

From measurements of the Equatorial Atmospheric Radar (EAR), located in West Sumatra (0.20 S, 100.32 E), Indonesia, Fujiwara et al. (2003) observed intermittent turbulence near the tropical tropopause with significant enhancements (a factor of 5) in the turbulent kinetic energy (TKE) lasting several days. Such enhancements in the TKE are believed to result from the breaking of large-scale Kelvin waves. Using the same EAR data set, Yamamoto et al. (2003) showed that eastward vertical wind shear around the equatorial tropopause frequently generates turbulence through Kelvin–Helmholtz instabilities (KHIs). Mega et al. (2010) presented detailed structures of KHIs in the equatorial UTLS from both high-resolution EAR measurements (using an interferometric imaging method) and radiosondes. From the VHF radar located in Gadanki (13.5 N, 79.2 E), India, Satheesan and Murthy (2002) and Satheesan and Krishna Murthy (2004) described turbulence characteristics in the tropical UTLS. These authors estimated TKE and TKE dissipation rates using several methods. Interestingly, they did not observe a clear variability in the turbulence intensity within the UTLS with altitude.

Sunilkumar et al. (2015) and Muhsin et al. (2016) presented the characteristics of turbulence in the tropical UTLS from GPS radiosonde observations obtained during more than 3 years at two stations located on the Indian Peninsula: Trivandrum (8.5 N, 76.9 E) and Gadanki (13.5 N, 79.2 E). The turbulent layers were detected using Thorpe's analysis (Thorpe1977) and following the procedure proposed by Wilson et al. (2010, 2011, 2013). The statistics of various turbulence parameters as well as the Brunt–Väisälä frequencies and the vertical shears are described for the convective troposphere and for the UTLS. The parameters describing the turbulence are either directly measured, such as the Thorpe lengths (an outer scale of turbulence) and the frequency of appearance of unstable layers, or are inferred on the basis of physical assumptions, like the TKE dissipation rates and the eddy diffusivity. Muhsin et al. (2020) extended the two previous studies by analyzing the soundings of six stations in South India, adding the Cochin (10 N, 76.3 E), Coimbatore (10.9 N, 76.9 E), Goa (15.5 N, 73.8 E), and Hyderabad (17.5 N, 78.6 E) stations to the two previously mentioned. These data were acquired during 4 years, from August 2013 to December 2017. All of these studies, based on radiosondes measurements, consistently show that the probability of occurrence of instability, either estimated from the gradient Richardson number (Ri) or from the squared Brunt–Väisälä frequency (N2), decreases with altitude above an altitude of 15 km, i.e., in the UTLS. They found that the probability of occurrence for Ri to be less than 0.25 is between 0 % and 5 % in the height range from 18 to 25 km and that the probability of occurrence of unstable regions (N2<0) is quasi-null. Interestingly, Muhsin et al. (2016) did not observe a clear diurnal cycle in the tropical UTLS.

Many studies of turbulence in the free atmosphere have been based on measurements from research aircraft (see Dörnbrack et al.2022, and references therein); however, to our knowledge, only one study has involved the tropical UTLS. Podglajen et al. (2017) used high-resolution (20 Hz) airborne measurements to study the occurrence and properties of small-scale (≲100 m) wind and temperature fluctuations in the tropical UTLS over the Pacific Ocean. They showed that wind fluctuations are very intermittent and appear to be localized in shallow layers, typically with ∼100 m thickness. Furthermore, active turbulent events appear to be more frequent at relatively low altitude and near deep convection. They observed that the motions are close to 3D isotropic and that the power spectra follow a -5/3 power law scaling. The diffusivity induced by turbulent bursts is estimated to be in the order of 10−1 m2 s−1 and decreases from the bottom to the top of the tropical UTLS.

Apart from convective regions, many observations of the free atmosphere show successions of strata, i.e., alternating layers in which the flow is turbulent, whose depth varies from a few tens to a few hundreds of meters, separated by stable static regions, i.e., regions where the flow is laminar (e.g., Luce et al.2015; Podglajen et al.2017; Wilson et al.2018). The numerical simulations of Fritts et al. (2003) support such observations. Turbulent and laminar strata are expected to exist in the tropical UTLS. As turbulence, by nature, is a dissipative process, the turbulent layers have a finite lifetime, depending on the mesoscale processes yielding energy (shear instability or wave breaking). Thus, by observing the flow at a given altitude level, we expect to find an alternation of turbulent and laminar episodes.

The Strateole-2 project has been set up in recent years in order to better understand the dynamics, transport, microphysics, and dehydration of the tropical UTLS (Haase et al.2018). It is an international project involving several research groups, mostly in France and the USA, led by the Laboratoire de Météorologie Dynamique (LMD) and Centre National d'Études Spatiales (the French space agency). The uniqueness and strength of the Strateole-2 project come from the fact that the atmospheric measurements are obtained under superpressure balloons (SPBs) that can fly for several months (typically 3 months) at an approximately constant density level. Standard measurements performed under all of the Strateole-2 SPBs, including temperature, pressure, and GPS position, allow one to describe meso- and small-scale dynamic processes at the flight level, in a quasi-Lagrangian way (Hertzog and Vial2001; Hertzog et al.2012; Podglajen et al.2016; Corcos et al.2021). The SPBs usually carry a second gondola with several configurations of instruments, allowing measurements of constituents, aerosols, cloud heights, and radiation, among others.

The main purpose of this paper is to present the methods allowing one to determine the dynamical state of the flow, laminar or turbulent, in which the SPBs are flowing. It is based on the estimation of the local stratification of the flow using measurements of temperature, pressure, and GPS position. The rest of the paper is organized as follows: the data set used is presented in Sect. 2, the methods for detecting turbulence and for estimating various indicators of the flow stability are described in Sect. 3, some results are presented (Sect. 4) and discussed (Sect. 5), and a concluding section (Sect. 6) summarizes the main findings of the paper.

2 The data set

2.1 The Strateole-2 C0 campaign

The probatory Strateole-2 campaign, called Strateole-2 C0, was held from November 2019 to March 2020. During this campaign, eight SPBs were launched from Mahé (437 S, 5527 E), Seychelles, and flew eastward following the quasi-biennial oscillation (QBO) phase that prevailed in the lower stratosphere at this time (Fig. 1). The main characteristics of the flights are provided in Table 1. Four of these flights carried only the TSEN (Thermodynamics SENsors) instrument and its associated GPS receiver, whereas the other four carried additional scientific instruments. At first order, SPBs fly at a quasi-constant density level (Vincent and Hertzog2014). TTL and STR flights were associated with different carried masses and/or balloon sizes, and they drifted at altitudes of ∼19 and ∼20.5 km, respectively.

Figure 1Trajectories of the eight SPBs.

Table 1Main characteristics of Strateole-2 C0 flights. All superpressure balloons are equipped with TSEN (Thermodynamics SENsors) instruments.

Download Print Version | Download XLSX

2.2 The TSEN measurements

In this study, we use the TSEN measurements that are performed at the balloon flight level for each flight. They consist of two respective temperature measurements, one with a thermistor sensor (TS) and one with a thermocouple sensor (TC), and a pressure measurement (P). They are closely associated with GPS observations that provide the position (long; lat; and altitude above sea level, a.s.l.) of the balloon as well as with the solar zenith angle (SZA). All of these measurements are acquired every 30 s, except for the pressure measurements that are acquired every 1 s.

These raw measurements enable us to compute two estimates of the balloon potential temperature (θS and θC) and density (ρS and ρC), respectively, obtained with the thermistor and thermocouple temperatures. The horizontal velocity components (u and v) are evaluated from successive GPS positions.

In the following, we will make use of the increments of measured quantities, namely the differences between consecutive measurements: δTS will, for instance, stand for the increments of the thermistor temperatures. Note that two independent estimates of the balloon vertical displacements are available: (1) δzGPS obtained with the raw GPS altitudes and (2) δzP computed with the pressure measurements and further assuming hydrostaticity.

TSEN measurements were also performed on the profiling unit of the Reeldown Aerosol Cloud Humidity and Temperature Sensors (RACHuTS) instrument (Kalnajs et al.2021). RACHuTS notably allows one to obtain vertical temperature profiles of 2 km length below the balloon by reeling down (and then up) the profiling unit at a vertical velocity of about 1 m s−1. Only one TSEN temperature sensor (a thermistor) is implemented in the profiling unit. Measurements are performed with a nominal sampling rate of 1 Hz, i.e., the vertical resolution of the RACHuTS temperature profiles is ∼1 m. We degraded its vertical resolution to about 30 m in order to (1) improve upon the 1 m raw precision of altitude measurements and (2) obtain a similar vertical resolution from TSEN and RACHuTS measurements (as the amplitude of the balloon oscillations is typically ±15 m).

Instrumental noise

The instrumental noise of the TSEN measurements is an important characteristic that has to be taken into account to assess whether the flow is turbulent. It is assumed that this noise is a zero-mean, uncorrelated process, i.e., a white noise, contributing to the measured signal. A way to estimate the white noise level is to compute the variance of the nth-order increments of the time series (n⪆6), such as high-order differentiation performing a high-pass filtering of the time series. As shown in Appendix A, the variance of the nth-order increments tends to the variance of the uncorrelated signal weighted by the sum of the squared binomial coefficients of order n−1. Here, noise levels are estimated on time segments of 21 samples (10 min); thus, they are estimated on more than 10 000 segments for a flight lasting ∼100 d. This method works well if the signal spectra exhibit a quasi-constant floor, i.e., if white noise clearly contributes to the measured signal for frequencies smaller than the Nyquist frequency, as is frequently the case for lidar or radar signals. However, we do not observe any white noise level for the TSEN measurements acquired with a sampling period of 30 s (TS, z, u) (see, for instance, the power spectra of the vertical displacements δz shown in Fig. 6). Therefore, the variance of the nth-order increments cannot be interpreted as being solely due to uncorrelated noise, even if some noise is expected to contribute to this high-frequency variance. However, a rough estimate of the noise level has been evaluated from the data segments showing the smallest variances, i.e., segments for which a comparatively smaller contribution of the atmospheric signal is expected. Tables 2 and 3 display estimates of the noise levels for quantities measured or inferred, obtained with the average of the smallest 10 % of variance values of the sixth-order increments. Such estimates, although not fully satisfying, nevertheless provide valuable information on the relative quality of the data, allowing one to compare noise levels between sensors, between flights, or between night and day measurements.

Table 2Estimated noise (standard deviation) of temperature, pressure, and dry potential temperature. For each flight, the first row indicates the flight-mean values (bold), and the second row indicates the respective nighttime and daytime average values.

Download Print Version | Download XLSX

Table 3Same as Table 2 but for altitudes derived from the pressure and altitude observations as well as zonal and meridional velocities.

Download Print Version | Download XLSX

The noise levels of the two temperature sensors, the thermocouple and thermistor, are quasi-identical. We note that the noise level of temperature measurements is a factor of 3–5 larger during daytime than during nighttime for all flights. Such an increase in temperature noise during the day very likely results from the random passage of the sensors in the wake of their electrical wires or mechanical support. These devices, which are significantly thicker than the sensors themselves, are heated by solar radiation during daytime and are consequently warmer than the air or even the sensor temperature. We also observe that the flight 03_TTL3 temperature measurements are the noisiest (14 mK vs. 2–9 mK). Unlike the other flights, during which temperature sensors were hanging at the very bottom of the flight train as far as possible from large elements (e.g., the gondolas), sensors on flight 03_TTL3 were located within the flight train. Thus, they were more prone to be affected by the warm wake of other devices in the flight train during daytime.

The noise estimates with respect to the pressure measurements and to the vertical position measurements (from both pressure differences and GPS) are slightly larger during daytime than during nighttime. The reason for this slight increase in the noise level is probably not the consequence of an increase in instrumental noise. It may rather be the consequence of the greater amplitudes of the balloon oscillations during daytime, which would impact the energy density of P and z at high frequency. The noise levels of the horizontal wind components estimated from GPS positions are very small, less than 1 cm s−1. They do not show any night–day difference.

2.3 Oscillating movements of the balloons around their equilibrium positions

Superpressure balloons drift with the winds, following a quasi-Lagrangian behavior, which is useful to document the detailed evolution of a given air mass. Their displacements, however, differ from those of an air parcel in at least two ways. First and most importantly, the balloons follow isopycnic trajectories. The balloon envelope is almost inextensible, with the balloon diameter varying by less than 1 %; hence, the volume is fixed as long as a superpressure is present. As the total mass of the flight train is fixed, the density of the balloon remains constant, implying an isopycnic trajectory. In contrast, air parcels follow isentropic trajectories in the absence of diabatic forcing. The second difference comes from the existence of natural oscillations of the balloon around its equilibrium density surface (EDS), where it achieves neutral buoyancy. These oscillations are of minor importance for the study of phenomena on timescales larger than a few tens of minutes (Vincent and Hertzog2014). However, they will be central to the current study, as an oscillating balloon samples short vertical profiles of key meteorological variables. The methodology developed in the present study precisely aims at exploiting these short “unintended” vertical profiles to diagnose the occurrence of turbulence.

The atmospheric density (ρ) at the flight level is estimated from the measurements of temperature (T) and pressure (P) using the ideal gas law: ρ=P/RaT, where Ra is the ideal gas constant for dry air per unit mass. As two independent measurements of T exist, two estimates of ρ can be calculated, denoted as ρS and ρC, using the thermistor or thermocouple temperatures, respectively. Figure 2 displays the systematic balloon fluctuations in density about their EDS (Fig. 2a), as well as the related variations in balloon altitudes (Fig. 2b) during 4 h of flight 01_STR1. The fluctuations in density are about ±5×10-4 kg m−3, corresponding to relative fluctuations of ±0.5 %. The corresponding amplitude of the balloon vertical displacements is ±20 m in this case. Thus, SPBs move around their EDS, exploring the atmosphere a few tens of meters above and below their position. These oscillatory motions induce fluctuations in the measured (P, T, zGPS) and inferred quantities (ρ, θ, zP). We shall use these variations to describe some properties of the flow by hypothesizing that the observed variability in a quantity depends on both the balloon vertical displacement and the local vertical gradient in the quantity.

Figure 2(a) Density fluctuations observed from 06:00 to 10:00 UTC on 19 November. (b) Fluctuations in the height of the balloon for the same time interval.


Periods and amplitudes of the vertical oscillations

Three angular frequencies, which are close to each other, first need to be distinguished: the Brunt–Väisälä frequency (N), the neutral buoyant oscillation frequency (ωNBO), and the observed frequency of SPB oscillations (ωB). The frequency at which a spherical balloon with constant volume oscillates about its EDS, the neutral buoyant oscillation (NBO) frequency, reads as follows (Hanna and Hoecker1971; Nastrom1980; Vincent and Hertzog2014):

(1) ω NBO 2 = 4 π 2 t NBO 2 = 2 g 3 T T z + g R a ,

where tNBO is the NBO period, and g is the acceleration of gravity. The squared Brunt–Väisälä (BV) frequency, N2=g/T(T/z+g/cp), can be expressed as a function of ωNBO2:

(2) N 2 = 3 2 ω NBO 2 - g γ H = 3 2 ω NBO 2 - 5 7 g H ,

where cp is the air specific heat capacity at constant pressure, γ is the heat capacity ratio, and H=RaT/g is the atmospheric scale height.

From Eq. (2), it can be shown that NωNBO as long as ωNBO210g/7H, i.e., tNBO2.8π2H/g120 s, or T/z8g/7Ra39 K km−1. Such conditions are met very frequently in the lower stratosphere, if not always. Note that both ωNBO and N increase with the vertical gradient of temperature T/z.

Figure 3a displays the observed periods of the balloon vertical oscillations (tB) for the whole 01_STR1 flight. The gray and white stripes correspond to nights and days, respectively, and the orange curve shows the daytime and nighttime averages. The histogram of the oscillation periods is shown in Fig. 3b, and the periods are observed to range from ∼135 to 800 s. No clear day–night variation is visible. Figure 3c shows the cumulative distribution function (CDF) of the periods: the median tB value is close to 220 s, and the mean is ∼250 s. Hence, ωB is larger than N as long as N is smaller than 2.5×10-2 rad s−1, which is a typical value in the tropical lower stratosphere. In the following, we shall see that ωB is systematically larger than N.

Figure 3(a) Estimates of the SPB oscillation periods during flight 01_STR1; the orange staircase curve shows daytime and nighttime averages. (b) Histogram of the oscillation periods. (c) Cumulative distribution function of the oscillation periods.


Figure 4 shows the corresponding amplitudes of the observed balloon oscillations as well as their probability and cumulative distribution function. The observed amplitudes range from ∼0 to 150 m, with an average (median) of 14 m (15 m). Some large amplitude oscillations (>100 m) are observed. We found that they are most often associated with depressurization events. A weak but clear day–night variability is found: the oscillation amplitudes are about 20 % larger during daytime than during nighttime. Moreover, greater variability in altitude is observed during daytime at high frequencies, as reported in Table 3.

Figure 4(a) Estimates of SPB oscillation amplitudes during flight 01_STR1; the orange staircase curve shows daytime and nighttime averages. (b) Histogram of the oscillation amplitudes. (c) CDF of the oscillation amplitudes.


3 Data processing

3.1 Methods for detecting the occurrence of turbulence from the TSEN measurements

The timescales of the turbulent fluctuations are expected to be smaller than the Brunt–Väisälä period, tN. (In the following, “high frequencies” refer to frequencies larger than N/(2π).) As tB is close to or even smaller than tN, the high-frequency variability up to the Nyquist frequency (1/60 Hz) will be dramatically affected by the balloon oscillations regardless of the state of the flow (laminar or turbulent). Hence, the detection of turbulence from either the variance or the spectral characteristics of the TSEN measurements at high frequencies appears difficult, if not impossible.

Diagnostics on the dynamical state of the flow can, therefore, only be based on the local properties of atmospheric stratification, either stable or neutral/unstable. A consequence of turbulence is to restore stability from a preceding unstable state of the flow, which is achieved by locally mixing the fluid. In this case, the conservative quantities (potential temperature, specific humidity, and momentum) are expected to exhibit weak horizontal and vertical variability within the turbulent layer. Note that neutral or even unstable stratification conditions may precede turbulence, but such conditions cannot persist and will necessarily cause turbulence. Conversely, laminar flow is expected to exhibit significant variability in the same quantities along the vertical direction.

The first way to diagnose that the flow is turbulent is, thus, to detect a null gradient in potential temperature. In such a case, the correlation between the vertical displacements (δzB) and the corresponding potential temperature increments (δθ) is expected to be null. Therefore, the implementation consists of testing the H0 hypothesis of a null correlation between δzB and δθ. This method is hereafter called the correlation method. The second way to diagnose a turbulent flow is based on the estimation of the local Richardson number Ri=N2/S2, where S2=(u/z)2+(v/z)2 is the squared shear. The flow is assumed to be turbulent when the Richardson number becomes less that 0.25. This method is hereafter called the Ri method.

3.2 Implementation of the two methods of turbulence detection

3.2.1 Relationship between the increments of measured quantities and the increments of vertical displacements

The two proposed methods of turbulence detection are based on either estimates of correlations (for the correlation method) or linear regressions (for the Richardson method) between increments. The increments are simply defined as the differences between consecutive measures, separated by δt=30 s. Such a differentiation realizes a high-pass filtering of the time series. Both the correlations and the linear regressions are computed over time periods of 1 h, i.e., with 120 observations that are 30 s long. The choice of this number of observations enables one to obtain relatively small uncertainties on the estimates of correlation coefficients and slopes, at the price of being unable to detect turbulence layers with timescales significantly shorter than 1 h (in the frame moving with the wind).

The oscillatory motions of SPBs at frequency ωB are not expected to occur on isentropic surfaces nor on isopycnic surfaces. If the vertical displacement of the balloon between two different times separated by δt=30 s is δzB=z2-z1, the change in the measured potential temperature δθ depends on both δzB and the vertical displacement of the atmosphere δζθ during δt:

(3) δ θ = θ z ( z 2 - z 1 ) = θ z ( δ z B - δ ζ θ ) ,

with δζθ=z1-z1 corresponding to the change in the height of the isentropic surface during δt (see Fig 5). Equation (3) must be modified for nonconservative quantities, such as temperature (T) or horizontal velocities (u and v), in order to take into account the change in the considered quantity during δt on the isentropic surface. For instance, the temperature increment reads as follows:

(4) δ T = T z δ z B - T z + g c p δ ζ θ .

Figure 6 shows the power spectral density (PSD) of the vertical displacements (δzB) for flight 02_STR2, either estimated from the GPS altitudes (δzGPS) or from the pressure observations (δzP). Both PSDs look very similar over the whole frequency range. A large spectral peak, centered at 4×10-3 Hz, is observed on these PSDs. This spectral peak corresponds to the SPB oscillatory motions about their EDS with a period of ∼220 s, as described in Sect. 2.3. Apart from long-duration balloons, radar is one of the few techniques that is able to infer the air vertical velocity in the free atmosphere. Few radar-borne PSDs of the vertical velocity have been published (Ecklund et al.1986; VanZandt et al.1991; Satheesan and Murthy2002). They either show a weak enhancement at frequencies close to N with respect to the spectral level at ωN or no enhancement at all. Thus, it is believed that the high-frequency vertical displacements in the balloon observations, 15 m on average, mostly result from the balloons' oscillating motions about their EDS, rather than from the isentropic vertical displacements of air parcels. In other words, δzBδζθ in Eqs. (3) and (4).

Figure 5Schematic of the vertical displacement of an SPB. Between the measurement times t1 and t1+30 s, the vertical displacement of the SPB is δzB, and the vertical displacement of the isentropic surface is δζθ.


Figure 6Power spectra of vertical increments of SPB heights, i.e., of vertical displacements between consecutive measurements (30 s), for flight 02_STR2.


3.2.2 The correlation method

The covariance of two quantities X and Y, is defined as follows:


Here, E is the mathematical expectation that is estimated by

(5) Cov ̃ ( X , Y ) = 1 N - 1 i = 1 N ( X - X ) ( Y - Y ) ,

where N is the sample size. The Pearson correlation coefficient ρP(X,Y)=Cov(X,Y)σXσY is estimated as follows:

(6) r P ( X , Y ) = Cov ̃ ( X , Y ) s X s Y ,

where sX is the estimate of the standard deviation σX of quantity X, i.e., sX=Coṽ(X,X). A second correlation coefficient, the Spearman rank correlation coefficient, rS(X,Y), has also been used in the present study. It is a measure of the statistical dependence between the rank of two variables X and Y. It is computed as the Pearson correlation between the rank values of those two variables:

(7) r S ( X , Y ) = r P ( rk X , rk Y ) ,

where rkX is the rank of variable X, i.e., sample X is replaced by the rank of X in the expression of rP (Eq. 6). The use of the nonparametric Spearman correlation makes it possible to get rid of outliers in the time series (Spearman1904).

If the flow is stably stratified, the increments of potential temperature δθ are related to the vertical displacements of the SPB δzB (Eq. 3). The covariance Cov(δθ,δzB) reads as follows:

(8) Cov δ θ , δ z B = E δ θ × δ z B = E θ z δ z B × δ z B - δ z B × δ ζ θ .

Noting that (i) δζθδzB and (ii) the vertical oscillations of the balloon and the displacements of the isentropic surfaces are expected to be nonsynchronous because ωB does not correspond to any atmospheric frequency (ωB>N), we assume that the covariance E[δzB×δζθ] is negligible compared with the variance of δzB, i.e.,

(9) Cov δ θ , δ z B = E θ z δ z B × δ z B .

Notice that, under the above hypothesis, the covariance of a nonconservative measured quantity X (T, u, v) reads as follows:

(10) Cov δ X , δ z B = E X z δ z B × δ z B .

If θz is positive in the time interval during which the covariance is estimated, Cov(δθ,δzB) is expected to be positive. Moreover, if θz is strictly constant during the time interval,

(11) θ z = Cov ( δ θ , δ z B ) Var δ z B .

If θz is not strictly constant, as is very likely the case, the Cov(δθ,δzB)/Var[δzB] ratio can be interpreted as an estimate of the mean gradient θz during the considered time interval at the flight level of the SPB. For a nonconservative quantities X,

(12) X z = Cov ( δ X , δ z B ) Var [ δ z B ] .

In the ideal case where θz is constant and there is no instrumental noise, σδθ=θzσδzB and the correlation coefficient rP(δθ,δzB)=1, for all θz>0. This conclusion also holds for the Spearman correlation rS if a linear relation is assumed between δθ and δzB. If θz>0 but is no longer constant during the considered time interval, the correlation coefficient is smaller than 1 but is still positive, as the product δθ×δzB is always positive. The correlation is further reduced, although remains positive, in the presence of uncorrelated noise (see Appendix B).

Figure 7(a) Time series of correlation coefficients rP and rS applied to the increments' potential temperature derived from TC and vertical displacements from zP. The correlations are evaluated on 1 h time segments (120 samples). The time series runs for 7 d, from 8 to 15 December 2019. The gray vertical stripes indicate nighttime. The light blue shaded areas show the ±1 standard deviation interval for the Pearson correlations. (b) Histograms of the two correlation coefficients obtained during the week. The blue shading and thick red line correspond to all of the 1 h time segments, and the thin green and orange lines correspond to the calm (nonturbulent) and turbulent time intervals, respectively.


The first method to detect turbulent mixing is to check for a nonpositive correlation, i.e., a null or negative correlation, between δθ and δzB. In order to infer if the null hypothesis H0 has to be rejected – that is, if a positive correlation exists – a standard hypothesis test is performed. A confidence interval for the correlation coefficient is based on a Fisher transformation of the correlation estimate (Hotelling1953). Here, the H0 hypothesis is rejected if the correlation estimate exceeds half of its standard deviation (unilateral test), with the corresponding confidence level being 69.15 %. If the H0 hypothesis cannot be rejected, turbulent mixing is diagnosed, and the turbulence flag is set to “true”. Two turbulence indexes are built, either from the Pearson or Spearman correlations.

Figure 7 illustrates the time series (Fig. 7a) and histograms (Fig. 7b) of rP and rS during 1 week of flight 02_STR2. Both correlations are estimated over time segments of 1 h. The light blue shaded areas show the ±1 standard deviation interval for the Pearson correlations only. The thick green lines at the top of Fig. 7a display the time interval during which turbulence is detected, i.e., when the turbulence indexes are set to true.

As expected, the correlation coefficients are found to be positive most of the time. On a few occurrences during the time period shown in Fig. 7, turbulence is detected as the correlation coefficients approach zero. A good agreement in the time periods identified as turbulent is found regardless of the correlation method used.

3.2.3 The impact of instrumental noise on the correlation coefficients

The impact of measurement noise can be evaluated using the averages of the largest correlation coefficients, all corresponding to θz>0 a priori. The reader should recall that, without instrumental noise, the correlation (δθ,δzB) is ideally 1 if θz>0 and remains constant over 1 h. As shown by Eq. (B1) in Appendix B, it is the ratio of the instrumental noise to the variance of the “geophysical” signal that appears in the correlation coefficient. In other words, the larger the amplitude of the balloon oscillations about their EDS, the less the measurement noise reduces the correlation coefficient.

The Spearman and Pearson correlation coefficients have been estimated for both potential temperature increments, δθC and δθS, and for both altitude increments, δzP=-δP/ρC and δzGPS. Thus, four turbulence indexes coefficients are estimated. We have then averaged the 50 % of correlation coefficients with the largest values, i.e., the quantiles 50–100 of the coefficients, for each of the eight flights (Table 4). These averaged coefficients give information, at least relatively, on the quality of the data and on the performance of the turbulence estimators: the larger the correlations, the smaller the impact of instrumental noise.

Table 4Averages of the 50 % of correlation coefficients with the largest values between δθ and δz. The Spearman and Pearson coefficients are displayed for the four combinations between δθS or δθC and δzP or δzGPS. About 1200 coefficients are averaged here for each flight. The largest correlation coefficient is displayed in bold for each flight. The larger the correlations, the smaller the impact of instrumental noise (see text).

Download Print Version | Download XLSX

Table 4 again illustrates that flight 03_TTL3 is abnormally noisy because of the poorer (potential) temperature measurements (see Table 2). Systematically, we observe the following:

  1. the correlations are larger when using θC rather than θS,

  2. there is a slight amelioration in using δzP rather than δzGPS,

  3. the Spearman correlations (rS) are slightly larger than the Pearson correlations (rP).

We shall, therefore, prefer the Spearman estimator, using δTC (and derived δθC) combined with δzP, in evaluating both the correlations and the vertical gradients from which the turbulent indexes are estimated, as the impact of instrumental noise appears to be lower in this estimator.

3.2.4 The Richardson method

The second detection method is based on an evaluation of the local Richardson number (Ri). Notice that θz (Eq. 11) is the slope of a least squares linear regression between the 120 potential temperature increments (δθ) and vertical displacement increments (δzB). This method of estimating the vertical gradient of the potential temperature is labeled the least squares fitting (LSF) method in the following. An alternative fitting method, the so-called Theil–Sen fitting (TSF) method, has also been used (Sen1968). The Theil–Sen estimator is defined as the median of the slopes of all lines through pair of points of the sample. It is a nonparametric and robust method that is almost insensitive to outliers. Both fitting methods have been applied to estimate the vertical gradients of several measured quantities, such as T and P or the horizontal wind velocities u and v (Eq. 12).

Two estimates of the mean Brunt–Väisälä frequency are thus computed with either the LSF or the TSF, applied to δTC and δzP, which were identified as the less noisy estimates of the air temperature and vertical displacements, respectively (see Tables 2 and 3). The squared shear, evaluated from the vertical gradients of the horizontal velocities, is obtained similarly. Hence, two turbulence indexes are built using LSF and TSF estimates of the Richardson number and are set to true if Ri1/4.

Figure 8a displays the vertical gradients of temperature for the same 1-week time series as that shown in Fig. 7. Estimates are shown for only the thermocouple temperatures and vertical displacements obtained with the pressure increments. Time series (not shown) obtained from other pairs of variables (using TS or TC or using zGPS or zP) are very similar to the one presented, regardless of the fitting method used. The uncertainty (±1 standard deviation) is shown in light red for the LSF estimator. The turbulence detection criteria based on the Richardson number are also displayed (thick green lines at the top of the figure). Figure 8b shows the probability distribution functions (PDFs) of the vertical gradients of temperature shown above (blue shading and thick curve). PDFs corresponding to turbulent and calm (i.e., nonturbulent) periods according to the Richardson number criterion are also displayed using thin lines.

First, it should be noticed that the choice of the LSF or TSF estimates for the vertical gradient of temperature or for the turbulence index has only a minor impact. Similarly, the choice of parameters (zP or zGPS, TC or TS) and any combination of those has little influence on the detection of turbulence layers (not shown). Figure 8 also clearly illustrates that the time periods flagged as turbulent by the Richardson number index are associated with the small-value tail of the temperature gradient PDF (dT/dz-10 K km−1). Finally, those time periods are essentially the same as the one shown in Fig. 7, obtained using the correlation method.

Figure 8(a) Time series of T/z estimated by two respective fitting methods, least squares fitting (LSF) and Theil–Sen fitting (TSF), from TC and zP. Linear fits are performed on 1 h time segments (120 samples). The light red shaded areas show the ±1 standard deviation interval for the LSF estimates. (b) Histograms of the temperature gradients from the two respective methods. The histograms corresponding to the turbulent and laminar time segments are also plotted (light curves).


4 Results

4.1 Turbulence detection with the correlation method

Figure 9a displays the time series of both the Spearman and Pearson correlation coefficients between δθC and δzP for the whole flight 02_STR2. As done previously, the correlations are estimated on time segments of 1 h (120 samples). They range from -0.5 to 0.96. Both time series exhibit very similar variations, showing time intervals lasting several days with relatively large correlations (>0.7) and short bursts of time with low or even negative correlations. Two turbulence indexes inferred from the Spearman and Pearson correlation methods are shown as thick green lines at the top of Fig. 9a. The fraction of time during which the flow is detected as turbulent for flight 02_STR2 is 4.5 % (4.4 %) from the Spearman (Pearson) correlations (see Table 5). Time series and histograms of correlations obtained using measurements of δθS or δzGPS (not shown) provide time series that are almost indistinguishable from those shown in Fig. 9.

Figure 9(a) Time series of the Spearman and Pearson correlation coefficients between δθC and δzP during flight 02_STR2 (103 d). The turbulence indexes deduced from the two correlations are shown (thick green lines at the top of the plot). (b) Histograms of the Spearman (thick line) and Pearson (filled) coefficients for flight 02_STR2. Also shown are the histograms corresponding to the turbulent and laminar cases (thin lines). The time series and histograms of the two correlation coefficients are very similar, even though a few differences in turbulence detection are visible.


Figure 9b shows the PDFs of the two correlation coefficients over the whole flight (blue shading and thick red line). Both distributions are observed to be almost identical, although rS is slightly and systematically larger than rP for correlations larger than 0.5. The average values for both correlation coefficients are close: rS[δθC,δzP]=0.50 and rP[δθC,δzP]=0.48. The thin lines show the PDFs of rS for solely turbulent (orange) and solely laminar time periods (green). The PDFs associated with the laminar time periods are almost identical to those of the whole flight, as they correspond to about 95 % of the data. The turbulent time periods correspond only to the small-value tail of the overall distribution.

The fractions of time during which the flow is found to be turbulent using the correlation method are reported in Table 5 for the eight flights of the campaign. Overall, they look very consistent, typically ranging between 4.5 % and 6.5 %, with the exception of flight 03_TTL3, which is associated with noisier temperature measurements.

Table 5Fraction of time (percent) during which the flow is found to be turbulent. Turbulence is diagnosed from the correlation method based of δθC and δzP (second and third columns) and from the Ri(δTC,δzP) criterion (last column).

Download Print Version | Download XLSX

4.2 Estimations of the vertical temperature gradient with the fitting methods

Figure 10 shows estimates of T/z for flight 02_STR2 inferred from linear fitting of the increments of temperature and vertical displacements during hourly intervals. More precisely, Fig. 10a shows the vertical temperature gradient obtained with the TSF method applied to δTC and δzP. Time series obtained using other combinations of measured variables and fitting methods (not shown) are almost indistinguishable from the one shown in Fig. 10. The dry adiabatic lapse rate Γ=-g/cP is also indicated as a black line. The turbulence index based on Ri is shown as a thick green discontinuous line at the top of Fig. 10a. Figure 10b shows the PDFs of T/z for the whole flight and for the respective calm and turbulent time periods. The distribution of T/z is found to be very asymmetric, with the mode being close to −4 K km−1. Most of the turbulent cases are associated with T/z<Γ, with a few of them being associated with T/z>Γ. The cumulative distribution functions (CDFs) of the temperature gradients are displayed in Fig. 10c. About 80 % of the detected turbulent cases, associated with Ri<0.25, correspond to a super-adiabatic temperature gradient.

Figure 10(a) Two estimates of TC/zP for flight 02_STR2, from the linear fitting of TC and δzP at hourly intervals. The black line shows the adiabatic lapse rate -g/cP. The Ri turbulence index is drawn on the top of the figure as a thick green discontinuous line. (b) Histograms of TC/zP for all of the time intervals (blue shading) and for the respective laminar and turbulent cases (thin lines). (c) CDFs of the temperature gradients.


Based on T/z estimates, the Brunt–Väisälä and NBO frequencies can be evaluated, see Eqs. (1) and (2). Figure 11 displays the squared Brunt–Väisälä (N2), the theoretical NBO (ωNBO2), and the observed balloon (ωB2) frequencies for the whole 02_STR2 flight. The Brunt–Väisälä and NBO frequencies are deduced from the hourly estimates of T/z. The observed balloon frequencies are directly estimated from observations of the balloon oscillations (Fig. 3). A 60 min running average is then applied to the raw ωB time series. The three frequencies are close to but distinct from one another, with the observed balloon frequency ranging between N and ωNBO. The CDF of frequency ratios (Fig. 11c) reveals that ∼95 % of the measured ωB are smaller than the ωNBO estimates and that 19 % of ωB are smaller than N. This property is consistent with numerical simulations of the motion of a spherical SPB assuming that atmospheric forcing occurs at frequencies lower than the Brunt–Väisälä frequency (Podglajen et al.2016, in particular see their Supplement).

Figure 11Times series (a) and histograms (b) of the Brunt–Väisälä, NBO, and observed balloon frequencies. The Brunt–Väisälä and NBO frequencies are deduced from the estimates of T/z from Cov(δTC,δzP). (c) CDF of the ωB2/N2 and ωB2/ωNBO2 ratios.


Figure 12 shows the histograms of Tz=TC/zP for flight 3 obtained from the TSEN measurements (TSF method) and from the RACHuTS temperature profiles – down to 2 km below the balloon. The vertical gradients of the RACHuTS temperature are estimated using 30 m vertical segments. The RACHuTS histogram has common features with the TSEN histograms: an asymmetric distribution, the same negative modes, and a sharp transition around −10C km−1. These distributions are considered to be consistent, despite the fact that they are not obtained in the same altitude domain nor simultaneously.

Figure 12Histograms of T/zP from the covariance of the TSEN measurements (blue shading) and from the RACHuTS temperature profile (black line). The two histograms have common characteristics, despite the fact that they are not obtained at the same altitude levels.


4.3 Turbulence detection with the Ri method

Figure 13 shows estimates of the Richardson number along flight 02_STR2, obtained using either the LSF or TSF estimates, with TC, zP, u, and v observations. The corresponding two turbulence indexes (Ri1/4) are also shown in Fig. 13a (green thick line). Figure 13b displays the histograms of the two Ri estimates. The time series and histograms of the two Ri estimates look very similar, even though a few differences in the turbulence detection are visible.

Figure 13(a) Richardson number estimates from TC, u, v, and zP using the respective least squares fitting (LSF) and Theil–Sen fitting (TSF) methods. The thick black line shows the threshold Ri=0.25. The two inferred turbulence indexes are shown (thick green lines at the top of the plot). (b) Histograms of the two Richardson number estimates.


As shown in Table 5, the fraction of time during which the flow is found to be turbulent according to the Ri criterion is 3.7 % for flight 02_STR2. More generally, the detection of turbulent flow during the Strateole-2 C0 campaign ranges from 3.3 % to 5.7 % of the time, with the exception of flight 03_TTL3.

5 Discussion

5.1 Possible impact of wake on the temperature measurements

Due to the vertical oscillations of the balloon, the T sensors were possibly in the wake of the balloon or of the flight chain. Note that we expect the balloon wake to be warm during daytime and cold during nighttime, with the balloons being cooler than the ambient air during nighttime.

The temperature sensors were located 27 m below the balloon base (except for during the 03_TTL3 flight) and 15 m below the EUROS gondola. The diameter of the balloons was either 11 m (TTL) or 13 m (STR). On all but the TTL3 flights, the T sensors were located 7 m below the last gondola in the flight chain. Flight 03_TTL3, carrying the RACHuTS system, is an exception, as the temperature sensors were located 30 cm from the EUROS gondola.

The question of the possible impact of the wake on temperature measurements has been taken into account. Indeed, we have calculated the statistics, vertical gradients, and correlations, considering only the downward phases of the oscillations, i.e., when the temperature sensors (which are located at the lower end of the flight chain) sample “fresh” air if a minimum shear exists. The resulting time series (not shown) are noisier because we only consider about half of the samples. However, both time series of correlations and temperature gradients have similar characteristics to those calculated when considering all samples, showing the same succession of stable and unstable periods. Therefore, we conclude that the impact of the wake does not significantly affect the estimated statistics, correlations, nor covariances, and, due to the increase in noise, we choose to consider all of the samples.

5.2 Comparison of the turbulence indexes

Table 6 shows the percentage of identical detections (laminar or turbulent) of the four turbulence indexes (rP, rS, RiLSF, and RiTSF) for the eight balloon flights. The percentage of similar detections ranges from 97.9 % (rS vs. RiLSF) to 99.07 %. However, as pointed out by an anonymous reviewer, such overall agreement does not imply equally good agreement between estimators when turbulence is detected. Choosing RiTSF as a reference, we compared the diagnoses with the other three estimators (RiLSF, rP, and rS). Figure 14 shows the percentages of similar and different detections when the flow is diagnosed as turbulent (T) or laminar (L). It reveals that the detections are similar for more than 99 % of cases if the flow is diagnosed as laminar. If the flow is diagnosed as turbulent, the rate of identical detections drops to about 80 %, with the comparison being the worst for rP (76 %) and comparable for rS and RiLSF (83.2 % and 85.4 %, respectively). We believe that these differences mainly result from the fact that the threshold values, zero correlation or Ri=0.25, correspond to the tails of the distributions of these estimates (see the histograms in Figs. 9 and 13). When the atmosphere is weakly stratified, threshold effects are likely to be important, leading to some differences in the diagnoses of flow conditions.

We also found that the time fraction of turbulent episodes obtained by the Ri criterion is almost always smaller than that obtained by the correlation methods (Table 5). This can be partly due to the threshold values of the hypothesis tests of a null correlation (i.e., to the choice of a confidence interval).

Table 6Percentage of identical turbulence detections from four turbulence index using (TC, zP) for flight 02_STR2.

Download Print Version | Download XLSX

Figure 14Percentages of true (similar) and false (different) detections for turbulent (T) and laminar (L) episodes diagnosed by the RiTSF criterion. The three estimators, RiLSF, rP, and rS, are compared to RiTSF. The rate of agreement is larger than 99 % in the case of laminar flow and about 80 % in the case of turbulent flow.


5.3 Occurrences of negative values for Ri and N2

The time series of N2 (Fig. 11) and Ri (Fig. 13) show some negative values. For the considered flight (02_STR2), the occurrence frequency of negative N2 is 3.4 % (from the Theil–Sen regression performed on TC and zP). Such a negative Ri (N2) can result from both the dispersion of the temperature gradient estimates and the occurrence of episodes of unstable stratification.

Negative estimates of N2 (or Ri) could be due to the precision of the temperature gradient estimates, which are expected to be scattered around a value close to −10C km−1 in the case of neutral stratification. Temperature gradients are estimated from the covariance of temperature increments and displacements, with these covariances scattered around their mean values. As a result, N2 estimates can be negative even if the stratification is neutral or nearly neutral.

However, unstable stratifications (N2<0) seem to occur in the lower stratosphere, as they have been reported in the literature. For instance, the detection of turbulence by the Thorpe method from in situ measurements is based on observations of θ/z<0, i.e., N2<0 (Thorpe1977). The probability of occurrence of such unstable layers likely depends on the vertical resolution of the profiles (e.g., Wilson et al.2011) but is not zero. In the lower stratosphere, KHIs are expected to be the main source of instability. For KHIs, turbulence is expected to be triggered for 0<Ri<1/4, i.e., for N2>0; however, once it is developed, the stratification can become almost neutral (N2≈0) or even unstable (N2<0) as a result of stirring and mixing. Therefore, it is plausible that the occurrence of unstable episodes may also contribute to negative values of N2 or Ri estimated from covariances.

Numerical simulations indicate that large temperature gradients are expected at the edges of turbulent layers (Fritts et al.2003; Werne and Fritts1999). These strong temperature gradients are also commonly observed from radiosonde profiles when turbulence detection is performed by the Thorpe method (e.g., Luce et al.2002). It is clear that sampling by balloons drifting within an air mass, and not cutting vertically through it as done by a radiosonde, will not allow one to identify turbulence from these sharp temperature gradients at the layer edges. Only the central part of the turbulent region, in which stratification is almost zero, can be identified as turbulent.

5.4 Spatial inhomogeneity of turbulence detections

The detection of turbulent episodes is far from uniform around the globe. Figure 15 shows the positions of turbulence detections (green dots) for the eight flights of the C0 campaign. Turbulence detections seem to be very rare over some regions (e.g., the southern Atlantic Ocean) and quite frequent over other regions (e.g., the Maritime Continent or the western Pacific). At the present stage, this result is only qualitative; however, interestingly, these inhomogeneities are likely related to the causes of the turbulence occurrences (see for instance Fritts and Alexander2003). A study examining processes that cause turbulence in the UTLS is in progress.

Figure 15The positions of turbulent patches are shown as green dots for the eight flights of the C0 campaign. Each dot corresponds to a 1 h time interval. The detections are based on the RiTSF criterion. The detected turbulence episodes are far from uniform around the globe.

6 Summary and concluding remarks

The present paper dealt with the detection of turbulence on superpressure balloon flights that drift for several months in the UTLS. During the Strateole-2 C0 campaign, eight SPBs were launched. Statistical methods to infer the flow regime, either laminar or turbulent, in which the SPBs were drifting are described. Some properties of the local stratification of the flow are also inferred. These methods are based on in situ GPS altitude, pressure, and temperature measurements, which were performed on all of the SPBs of the Strateole-2 C0 campaign with a time resolution of 30 s.

We make use of the SPB oscillations about their equilibrium density surface, which typically enable SPBs to vertically explore the atmosphere over 30 m (peak-to-peak displacements). The observed periods of oscillation, ∼220 s on average, are significantly smaller than the Brunt–Väisälä period. The large amplitude of these balloon motions at frequencies higher than the Brunt–Väisälä frequency (where turbulence is expected to occur) makes it very difficult, if not impossible, to detect turbulence from the direct characterization (i.e., from the variance or power spectra) of high-frequency fluctuations. On an other hand, thanks to the vertical motions of the SPBs around their EDS, the vertical gradients of any measured quantity can be estimated, either via the covariance between the increments of this quantity and those of the vertical displacements or, equivalently, via the linear fit between these two increments. For the present study, covariances and linear fits have been estimated on data segments of 1 h, i.e., 120 data samples.

Several turbulence indexes (true or false) are defined and compared. A first index is based on an inference test on the correlation between potential temperature and altitude increments. A null correlation is expected in the case of turbulent mixing, as the vertical gradient of potential temperature tends to zero. Two correlation coefficients, the Pearson and Spearman coefficients, have been tested and compared for the eight balloon flights. Alternatively, based on a linear fit between increments, vertical temperature gradients and horizontal wind can be evaluated at 1 h intervals, from which the local Richardson number can be deduced. A second turbulence index is then based on the criterion Ri<0.25. These different turbulence indexes compare well, as they coincide for more than 97 % of cases.

SPBs sample the atmosphere without any spatial or temporal sampling bias, drifting nights and days, over oceans and continents, and above convective and non-convective regions. The fraction of time during which the flow is found to be turbulent in the lower stratosphere appears to be quite small, ranging from 3.3 % to 6.6 %. Because of the lack of sampling bias, such a fraction of time can be interpreted as a fraction of space. The probability of occurrence of instabilities is far from being uniform in time (or space). Periods of several days are frequently observed during which the atmosphere is stable, i.e., without any instability. On the other hand, during certain periods, the frequency of instabilities appears to be quite high. At first sight, these differences could be attributable to the underlying deep convection, as observed, for example, when the balloons are flying over the Maritime Continent. However, this remains a very preliminary conclusion that needs to be substantiated in future work.

Appendix A: Estimation of the uncorrelated noise

The measured signal is assumed to contain an uncorrelated and centered noise contribution. This noise level is estimated on short data segments (∼20 data samples), with the useful signal being described by a polynomial fit of degree n. The time series of quantity X reads

(A1) X i ¯ = X i + ξ i 1 i n ,

where Xi¯ is the measured signal, and ξi is an uncorrelated noise of variance σξ. The measured first increment reads

(A2) δ X i ¯ = X i + 1 - X i + ξ i + 1 - ξ i .

If Xi is constant, i.e., Xi+1-Xi=0, the variance of δXi¯ reduces to

(A3) Var δ X i ¯ = 2 Var [ ξ ] = 2 σ ξ 2 .

The measured second increment reads

(A4) δ 2 X i ¯ = δ X i + 1 ¯ - δ X i ¯ = X i + 2 - 2 X i + 1 + X i + ξ i + 2 - 2 ξ i + 1 + ξ i .

If Xi varies according a linear trend, i.e., Xi+2-2Xi+1+Xi=0, the variance of δ2Xi¯ reduces to

(A5) Var δ 2 X i ¯ = Var [ ξ ] + 4 Var [ ξ ] + Var [ ξ ] = k = 0 2 2 k 2 σ ξ 2 = 6 σ ξ 2 .

The measured nth increment reads

(A6) δ n X i ¯ = δ n - 1 X i + 1 ¯ - δ n - 1 X i ¯ = k = 0 n ( - 1 ) n - k n k X i + k + k = 0 n ( - 1 ) n - k n k ξ i + k .

If Xi is described by a polynomial of degree n−1, the first term on the right-hand side of Eq. (A6) cancels, and the variance of δnXi¯ reduces to

(A7) Var [ δ n X i ¯ ] = k = 0 n n k 2 σ ξ 2 .

Increasing the order of differentiation enhances the relative contribution of the uncorrelated signal in the time series. After several differentiations, the variance of the differentiated time series is expected to converge to the weighted variance of the uncorrelated noise.

Appendix B: Impact of instrumental noise on the correlation coefficients

The Pearson or Spearman correlation coefficients (ρP and ρS, respectively) will be reduced because of uncorrelated noise in the time series of δzB and δθ. A simplistic model may help to illustrate this assertion. Let us note θ¯, the measured potential temperature, and zB¯, the measured altitude of the balloon. Assume that θ¯ can be analyzed as θ¯=θ+η, where θ is the real value and η is a centered random noise. Similarly, zB¯=zB+ζ. The variance of the measured increments δX¯ (X=θ or zB) reads


where σX is the standard deviation of the random noise on X. Also, Cov[δθ¯,δzB¯]=Cov[δθ,δzB]. The expectation of the Pearson correlation coefficient ρP[δθ¯,δzB¯]=Cov[δθ¯,δzB¯]/(Var[δθ¯]Var[δzB¯])1/2 reads

(B1) ρ P [ δ θ ¯ , δ z B ¯ ] = ρ P [ δ θ , δ z B ] 1 + 2 σ ζ 2 Var [ δ z B ] + 2 σ θ 2 Var [ δ θ ] + 4 σ ζ 2 σ θ 2 Var [ δ θ ] Var [ δ z B ] 1 / 2 ,

where σθ2 and σζ2 are the variances of the noises on θ and zB, respectively. Thus, the correlation coefficients are expected to be significantly reduced due to the instrumental noise. However, they will retain the same sign as the correlation coefficients without instrumental noise.

Data availability

The balloon-borne TSEN and RACHuTS data used in this study were collected as part of Strateole-2, which is sponsored by CNES, CNRS/INSU, NSF, and ESA. The Strateole-2 data set is available at (Hertzog2021).

Author contributions

RW developed the data processing method and wrote the majority of the manuscript. CP participated in conditioning and processing the data. AP processed the RACHuTS data. AH is PI of Strateole-2 and preprocessed the raw data, making them available. AP and AH made important contributions with respect to extensive exchanges about the physics of superpressure balloon-borne measurements. All authors (RW, CP, AP, AH, MC, and RP) participated in numerous discussions to elaborate and improve the overall understanding of the physics of measurements, and all co-authors contributed significantly to writing the paper.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


The authors warmly thank the three reviewers for their thorough review of this article. Their remarks and suggestions have helped to improve the paper in a very significant way. The authors acknowledge support from the Agence Nationale de la Recherche (TuRTLES project). Clara Pitois is funded via a Sorbonne Université doctoral fellowship.

Financial support

This research has been supported by the Agence Nationale de la Recherche (grant no. ANR-21-CE01-0016-01) and the Centre National d'Études Spatiales (grant no. TOSCA).

Review statement

This paper was edited by Markus Rapp and reviewed by Ulrich Schumann and two anonymous referees.


Corcos, M., Hertzog, A., Plougonven, R., and Podglajen, A.: Observation of Gravity Waves at the Tropical Tropopause Using Superpressure Balloons, J. Geophys. Res.-Atmos., 126, e2021JD035165,, 2021. a

Dörnbrack, A., Bechtold, P., and Schumann, U.: High-Resolution Aircraft Observations of Turbulence and Waves in the Free Atmosphere and Comparison With Global Model Predictions, J. Geophys. Res.-Atmos., 127, e2022JD036654,, 2022. a

Ecklund, W. L., Gage, K. S., Nastrom, G. D., and Balsley, B. B.: A Preliminary Climatology of the Spectrum of Vertical Velocity Observed by Clear-Air Doppler Radar, J. Appl. Meteorol. Clim., 25, 885–892, 1986. a

Fritts, D. C. and Alexander, M. J.: Gravity wave dynamics and effects in the middle atmosphere, Rev. Geophys., 41, 1003,, 2003. a

Fritts, D. C., Bizon, C., Werne, J. A., and Meyer, C. K.: Layering accompanying turbulence generation due to shear instability and gravity-wave breaking, J. Geophys. Res., 108, 8452,, 2003. a, b

Fueglistaler, S., Dessler, A. E., Dunkerton, T. J., Folkins, I., Fu, Q., and Mote, P. W.: Tropical tropopause layer, Rev. Geophys., 47, RG1004,, 2009. a

Fujiwara, M., Yamamoto, M., Hashiguchi, H., Horinouchi, T., and Fukao, S.: Turbulence at the tropopause due to breaking Kelvin waves observed by the Equatorial Atmosphere Radar, Geophys. Res. Lett., 30, 1171,, 2003. a

Haase, J., Alexander, M., Hertzog, A., Kalnajs, L., Davis, S., Ploufonven, R., Cocquerez, P., and Venel, S.: Around the world in 84 days, EOS, 99,, 2018. a

Hanna, S. R. and Hoecker, W. H.: The Response of Constant-Density Balloons to Sinusoidal Variations of Vertical Wind Speeds, J. Appl. Meteorol. Clim., 10, 601–604, 1971. a

Hertzog, A.: STRATEOLE2 - C0: TSEN - Thermodynamics SENsor, LMD/IPSL [data set],, 2021. a

Hertzog, A. and Vial, F.: A study of the dynamics of the equatorial lower stratosphere by use of ultra-long-duration balloons 2. Gravity waves, J. Geophys. Res., 106, 22745–22761, 2001. a

Hertzog, A., Alexander, M. J., and Plougonven, R.: On the Intermittency of Gravity Wave Momentum Flux in the Stratosphere, J. Atmos. Sci., 69, 3433–3448, 2012. a

Hotelling, H.: New Light on the Correlation Coefficient and its Transforms, J. Roy. Stat. Soc. B Met., 15, 193–225, 1953. a

Kalnajs, L. E., Davis, S. M., Goetz, J. D., Deshler, T., Khaykin, S., St. Clair, A., Hertzog, A., Bordereau, J., and Lykov, A.: A reel-down instrument system for profile measurements of water vapor, temperature, clouds, and aerosol beneath constant-altitude scientific balloons, Atmos. Meas. Tech., 14, 2635–2648,, 2021. a

Luce, H., Fukao, S., Dalaudier, F., and Crochet, M.: Strong mixing events observed near the tropopause with the MU radar and hight-resolution Balloon techniques, J. Atmos. Sci., 59, 2885–2896, 2002.  a

Luce, H., Wilson, R., Dalaudier, F., Nishi, N., Shibagaki, Y., and Nakajo, T.: Simultaneous observations or tropospheric turbulence from radiosondes using Thorpe analysis and the VHF MU radar, Radio Sci., 49, 1106–1123, 2015. a

Mega, T., Yamamoto, M., Luce, H., Tabata, Y. ndd Hashiguchi, H., Yamamoto, M., Yamanaka, M., and Fukao, S.: Turbulence generation by Kelvin‐Helmholtz instability in the tropical tropopause layer observed with a 47 MHz range imaging radar, J. Geophys. Res., 115, D18115,, 2010. a

Muhsin, M., Sunilkumar, S., Venkat Ratnam, M., Parameswaran, K., Krishna Murthy, B., Ramkuma, G., and Rajeev, K.: Diurnal variation of atmospheric stability and turbulence during different seasons in the troposphere and lower stratosphere derived from simultaneous radiosonde observations at two tropical stations in the Indian Peninsula, Atmos. Res., 180, 12–23, 2016. a, b

Muhsin, M., Sunilkumar, S., Ratnam, M. V., Parameswaran, K., Mohankumar, K., Mahadevan, S., Murugadass, K., Muraleedharan, P., Kumar, B. S., Nagendra, N., Emmanuel, M., Chandran, P., Koushik, N., Ramkumar, G., and Murthy, B.: “Contrasting features of tropospheric turbulence over the Indian peninsula”, J. Atmos. Sol.-Terr. Phys., 197, 105179,, 2020. a

Nastrom, G. D.: The response of superpressure balloons to gravity waves, J. Appl. Meteorol., 19, 1013–1019, 1980. a

Podglajen, A., Hertzog, A., Plougonven, R., and Legras, B.: Lagrangian temperature and vertical velocity fluctuations due to gravity waves in the lower stratosphere, Geophys. Res. Lett., 43, 3543–3553, 2016. a, b

Podglajen, A., Bui, T., Dean-Day, J., Pfister, L., Jensen, E., Alexander, M., Hertzog, A. Kärcher, B., Plougonven, R., and Randel, W.: Small-Scale Wind Fluctuations in the Tropical Tropopause Layer from Aircraft Measurements: Occurrence, Nature, and Impact on Vertical Mixing, J. Atmos. Sci., 74, 3847–3869,, 2017. a, b

Satheesan, K. and Krishna Murthy, R.: Estimation of turbulence parameters in the lower atmosphere from MST radar observations, Q. J. Roy. Meteor. Soc., 130, 1235–1249, 2004. a

Satheesan, K. and Murthy, R. F. K.: Turbulence parameters in the tropical troposphere and lower stratosphere, J. Geophys. Res., 107, ACL 2-1–ACl 2-13,, 2002. a, b

Sen, P.: Estimates of the Regression Coefficient Based on Kendall's Tau, J. Am. Stat. Assoc., 63, 1379–1389,, 1968. a

Spearman, C.: The Proof and Measurement of Association between Two Things, Am. J. Psychol., 15, 72–101,, 1904. a

Sunilkumar, S., Muhsin, M., Parameswaran, K., Venkat Ratnam, M., Ramkumar, G., Rajeev, K., Krishna Murthy, B., Sambhu Namboodiri, K., Subrahmanyam, K., Kishore Kumar, K., and Shankar Das, S.: Characteristics of turbulence in the troposphere and lower stratosphere over the Indian Peninsula, J. Atmos. Sol.-Terr. Phys., 133, 36–53, 2015. a

Thorpe, S. A.: Turbulence and mixing in a Scottish Lock, Philos. T. Roy. Soc. A, 286, 125–181, 1977. a, b

VanZandt, T. E., Nastrom, G. D., and Green, J. L.: Frequency spectra of vertical velocity from Flatland VHF radar data, J. Geophys. Res.-Atmos., 96, 2845–2855, 1991. a

Vincent, R. A. and Hertzog, A.: The response of superpressure balloons to gravity wave motions, Atmos. Meas. Tech., 7, 1043–1055,, 2014. a, b, c

Werne, J. and Fritts, D. C.: Stratified shear turbulence: Evolution and statistics, Geophys. Res. Lett., 26, 439–442,, 1999. a

Wilson, R., Luce, H., Dalaudier, F., and Lefrère, J.: Turbulent Patch Identification in Potential Density/Temperature Profiles, J. Atmos. Ocean. Tech., 26, 977–993, 2010. a

Wilson, R., Dalaudier, F., and Luce, H.: Can one detect small-scale turbulence from standard meteorological radiosondes?, Atmos. Meas. Tech., 4, 795–804,, 2011.  a, b

Wilson, R., Luce, H., Hashiguchi, H., Shiotani, M., and Dalaudier, F.: On the effect of moisture on the detection of tropospheric turbulence from in situ measurements, Atmos. Meas. Tech., 6, 697–702,, 2013. a

Wilson, R., Hashiguchi, H., and Yabuki, Y.: Vertical Spectra of Temperature in the Free Troposphere at Meso-and-Small Scales According to the Flow Regime: Observations and Interpretation, Atmosphere, 9, 415,, 2018. a

Yamamoto, M. K., Fujiwara, M., Horinouchi, T., Hashiguchi, H., and Fukao, S.: Kelvin-Helmholtz instability around the tropical tropopause observed with the Equatorial Atmosphere Radar, Geophys. Res. Lett., 30, 1476,, 2003. a

Short summary
Strateole-2 is an French–US initiative designed to study atmospheric events in the tropical upper troposphere–lower stratosphere. In this work, data from several superpressure balloons, capable of staying aloft at an altitude of 18–20 km for over 3 months, were used. The present article describes methods to detect the occurrence of atmospheric turbulence – one efficient process impacting the properties of the atmosphere composition via stirring and mixing.