# OVLI-TA: An Unmanned Aerial System for Measuring Profiles and Turbulence in the Atmospheric Boundary Layer

Sara Alaoui-Sosse, Pierre Durand, Patrice Medina, Philippe Pastor, Marie Lothon, Iuri Cernov

► To cite this version:

Sara Alaoui-Sosse, Pierre Durand, Patrice Medina, Philippe Pastor, Marie Lothon, et al.. OVLI-TA: An Unmanned Aerial System for Measuring Profiles and Turbulence in the Atmospheric Boundary Layer. Sensors, 2019, 19 (3), pp.581. 10.3390/s19030581. hal-02383020

HAL Id: hal-02383020

https://hal.science/hal-02383020

Submitted on 6 Jul 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# OVLI-TA: An Unmanned Aerial System for Measuring Profiles and Turbulence in the Atmospheric Boundary Layer

Sara Alaoui-Sosse<sup>1,2,*</sup>, Pierre Durand<sup>2</sup>, Patrice Medina<sup>2</sup>, Philippe Pastor<sup>1 ID</sup>, Marie Lothon<sup>2</sup>
and Iuri Cernov<sup>2</sup>

Received: 16 November 2018; Accepted: 28 January 2019; Published: 30 January 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** In recent years, we developed a small, unmanned aerial system (UAS) called OVLI-TA (Objet Volant Leger Instrumenté-Turbulence Atmosphérique) dedicated to atmospheric boundary layer research, in Toulouse (France). The device has a wingspan of 2.60 m and weighed 3.5 kg, including payload. It was essentially developed to investigate turbulence in a way complementary to other existing measurement systems, such as instrumented towers/masts. OVLI-TA's instrumental package includes a 5-hole probe on the nose of the airplane to measure attack and sideslip angles, a Pitot probe to measure static pressure, a fast inertial measurement unit, a GPS receiver, as well as temperature and moisture sensors in specific housings. In addition, the Pixhawk autopilot is used for autonomous flights. OVLI-TA is capable of profiling wind speed, wind direction, temperature, and humidity up to 1 km altitude, in addition to measuring turbulence. After wind tunnel calibrations, flight tests were conducted in March 2016 in Lannemezan (France), where there is a 60-m tower equipped with turbulence sensors. In July 2016, OVLI-TA participated in the international project DACCIWA (Dynamics-Aerosol-Chemistry-Clouds Interactions in West Africa), in Benin. Comparisons of the OVLI-TA observations with both the 60 m tower measurements and the radiosonde profiles showed good agreement for the mean values of wind, temperature, humidity, and turbulence parameters. Moreover, it validated the capacity of the drone to sample wind fluctuations up to a frequency of around 10 Hz, which corresponds to a spatial resolution of the order of 1 m.

**Keywords:** unmanned aerial vehicle (UAV); turbulence observations; atmospheric boundary layer (ABL); five-hole probe; OVLI-TA

## 1. Introduction

Technological advances related to the development of unmanned aircrafts for atmospheric observations make it possible to carry out measurements in areas at very low heights that are difficult to reach by piloted aircrafts, thus covering a wider airspace for collecting meteorological data. The unmanned aerial vehicle (UAV) offers an opportunity to investigate the atmospheric surface layer (ASL) in a way that is complementary to other platforms such as instrumented towers/masts, radiosondes, and piloted airplanes, thereby filling the gap between all available fixed and mobile platforms.

The ASL is the layer at the bottom of the atmospheric boundary layer (ABL), which directly interacts with the surface of the earth [1]. In general, it is some tens of meters thick [2]. It is characterized by a high variability, both in time (related to the diurnal cycle) and along the vertical; shearing stress is

<footnote n="1"><sup>1</sup> Institut Supérieur de l'Aéronautique et de l'Espace (ISAE-SUPAERO), Université de Toulouse, 31400 Toulouse, France; Philippe.PASTOR@isae.fr</footnote>

<footnote n="2"><sup>2</sup> Laboratoire d'Aérologie, Université de Toulouse, CNRS, UPS, 31400 Toulouse, France; pierre.durand@aeroobs-mip.fr (P.D.); patrice.medina@aeroobs-mip.fr (P.M.); marie.lothon@aeroobs-mip.fr (M.L.); iuri.cernov@gmail.com (I.C.)</footnote>

<footnote>* Correspondence: sara.alaoui-sosse@isae-supaero.fr; Tel.: +33-769-896-099</footnote>
constant in the vertical and the flow is insensitive to the earth's rotation. Profiles of meteorological parameters, such as temperature, moisture, wind, and turbulence, vary substantially from the surface to the top of the ASL. These follow for example, logarithmic laws or even more complex shapes according to the surface characteristics and turbulence conditions. The two sources of ASL turbulence are friction and buoyancy. The former results from the momentum extracted from the flow, and the latter from the difference in temperature (and, at a lower level, of moisture) between the surface and the air above it. Turbulence is three-dimensional and isotropic, at least for the smaller eddies, which are well described by the “K-41 theory” (Kolmogorov, 1941) on the inertial sub-range and its famous “-5/3 slope” for the spectral energy.

Very high level of ASL knowledge is required for both academic research and practical applications. For example, improving turbulent exchange parameterization in a numerical weather prediction model, or estimating wind turbine wake characteristics require appropriate observations in the ASL. The most important parameters to measure are the temperature, moisture, and wind, at a rate fast enough to capture eddies that significantly contribute to the turbulence energy. Starting from the know-how built up in past decades regarding observations aboard piloted airplanes (see e.g., [3]), and thanks to the miniaturization of sensors and acquisition systems, numerous unmanned aerial systems (UASs) have been developed in recent years all over the world. Similarly to on-board piloted airplanes, many developments were made to characterize the surface overflown, through remote sensing instruments. However, in this paper, we will focus on in situ observations of meteorological parameters and turbulence.

The UAVs were classified according to three main categories, depending on their weight and their payload capacity [4]. The first category regroups all UAVs weighing 10-30 kg, such as for example the Manta, ScanEagle, Aerosonde, and RPMSS (Robotic Plane Meteorological Sounding System); their advantage is their endurance and their highest payload capacity, but they are very expensive and often difficult to operate. The Manta (27.7 kg) and ScanEagle (22 kg), in particular, were used for turbulent flux measurements within terrestrial and marine ABL; wind components and humidity were measured up to 25 Hz and temperature up to 5 Hz [5]. Similarly, ALADINA (Application of Light-weight Aircraft for Detecting IN situ Aerosol) (25 kg) was elaborated for measuring boundary-layer properties, atmospheric particles, and solar radiation, as well as turbulence up to a frequency of ~7 Hz [6].

Secondly, category II includes vehicles that weigh more than 1 kg and less than 10 kg such as the Tempest, the meteorological mini unmanned aerial vehicle (M²AV), the NexSTAR, the multipurpose automatic sensor carrier (MASC), and the small multifunction autonomous research and teaching sonde (SMARTSonde). Although these UAVs are smaller, and have less payload capacity and autonomy of flight than the vehicles in category I, they can carry many of the sensors used for wind measurements. In addition to this, their cost is moderate, which makes them more easily deployable. For instance, M²AV (6 kg) measures the meteorological wind up to 40 Hz, equivalent to a spatial resolution of 55 cm at an airspeed of 22 m/s, by coupling a GPS and an inertial measurement unit (IMU) with a Kalman filter, and combining data with a five-hole probe for which the calibration had been obtained during a wind tunnel test [7]. M²AV and MASC were also used to assess the accuracy and frequency response of a multi-hole probe [8], and finally they can measure fluctuations up to 20 Hz after re-evaluating the pneumatic tubing setup and data acquisition. MASC was also operated for wind energy research [9]. Moreover, the BLUECAT5 (5 kg), developed for turbulence measurements in the ABL, could measure turbulence with a sampling of 60 Hz, by following a profile pattern at loiter radius [10]. However, the circular flight pattern prevents the conversion of time series into spatial scales, except for high rates (small scales) where air sampling can be considered as straight.

Finally category III assembles UAVs that weigh less than 1 kg, for instance the small unmanned meteorological observer (SUMO) and DataHawk, which have limited payload capacity and endurance compared with the two other categories. Their cost and facility of deployment enable small experiments to be carried out from anywhere. SUMO efficiently captures meteorological profiles of wind speed,
wind direction, temperature, and humidity [11,12]. Moreover, capabilities and limitations for turbulence observations are described in Reference [13].

Furthermore, different field campaigns have been conducted using UAVs. For instance, SUMO, MASC, and two multicopters participated in the Hailuoto 2017 field campaign in the arctic, in order to increase the understanding of the stable boundary layer by complementing the existing observation systems—like ground-based eddy covariance and automatic weather stations. This large set of data brought insights into the nature of turbulent events that induce rapid warming of layers close to the ice surface [14]. During the CLOUD-MAP campaign, thanks to a great amount of work on sensor integration and calibration/validation, atmospheric sampling of thermodynamic parameters and boundary-layer profiling were done with fixed and rotary wings UASs [15]. Additionally, UAVs are used for specific atmospheric issues as pollution and trace gas monitoring: a six-rotor UAS was used for studying air pollution episode characteristics and influential mechanisms that occurred in Nanjing during the 3–4 December 2017 [16]. On the other hand, three types of UAV including micro aerial vehicles (MAVs), vertical take-off and landing (VTOL), and low-altitude short endurance (LASE) systems, were evaluated to operate three kinds of atmospheric trace gas sensors. The best compromise was given by UAVs, which have wingspans <3 m for payloads <5 kg [17].

The unmanned aerial vehicle OVLI-TA ("Objet Volant Leger Instrumenté-Turbulence Atmosphérique", that is instrumented light aerial vehicle-atmospheric turbulence) belongs to the second category of UAVs. OVLI-TA was developed at Laboratoire d'Aérologie in Toulouse, France. The purpose was to develop a low cost system able to fly at low altitudes and easy to deploy; that may measure temperature, humidity, wind vector, and observe turbulence within the ABL. One of the unique attributes of the OVLI-TA system with respect to comparable UAVs is that its nose serves as a five-hole probe, which is 3D printed, easy to manufacture, and allows us to install the inertial platform and aerodynamic sensors in the same location. Thus, when solving the 3-component wind equation (see Section 2), even at a rate fast enough to capture turbulent scales, all the terms with angular rotations can be neglected. Whereas, they must be taken into account when using systems based on a multi-hole probe placed at the extremity of a boom. Furthermore, tubing length effects on pressure measurements are eliminated in our system because the transducers can be installed as close as possible to the pressure ports. Regarding UAV performance, our choice was guided by a compromise between a system small enough to be easy to operate, and large enough to allow a payload compatible with a good-quality sensors package and an endurance compatible with atmospheric boundary layer probing. With our system, the autonomy of flight was at least one hour. Moreover, even if OVLI-TA is lighter (3.5 kg) than other category 2 UAVs (e.g., 6 kg for the M²AV and 5 to 7.5 kg for the MASC), it permits comparable horizontal resolution. In addition to the hemispherical 5-hole probe situated on the nose of the airplane to measure attack and sideslip angles, there is a Pitot probe, a fast IMU, a GPS receiver, as well as temperature and moisture sensors placed in specific housings. Its cost is affordable, and hence we can risk flying in areas with turbulent conditions. OVLI-TA is mainly a system for profiling the atmospheric boundary layer, it also has the capacity to measure turbulence properly. Several flights were conducted in Lannemezan (France), where there is an equipped 60 m tower that was a reference to our measurements. The drone then participated in the international project DACCIWA (Dynamics-Aerosol-Chemistry-Clouds Interactions in West Africa), in Benin.

This paper gives a technical description of OVLI-TA, and highlights the results of profiling the atmospheric boundary layer obtained during flights in Lannemezan and in DACCIWA. The performance for measuring the mean meteorological parameters was evaluated against the 60 m tower and radiosonde profiles. Furthermore, turbulence observations on the tower and with the UAV were compared through a time-space conversion based on Taylor's hypothesis applied to the two platforms.
## 2. OVLI-TA Unmanned Aerial Vehicle and Methods

### 2.1. OVLI-TA Airframe Characteristics

OVLI-TA (Figure 1) is based on the commercial airborne model Techpod, manufactured by HobbyUAV. It has a wingspan of 2.60 m and a fuselage length of 1.14 m. Its cruise speed is between 12 and 20 m/s and the maximum speed reaches 28.3 m/s. This UAV weighs 3.5 kg, including the payload and battery (1.25 kg without). The wings and fuselage were made from Expanded Polyolefin (EPO). The drone was powered with electrical propulsion: MT2216-V2 motor and two 5200 mAh 11.1 V LiPo batteries, controlled by an RCTimer 50 A electronic speed controller (ESC), that enabled two hours of endurance. No catapult was needed for take-off, and the plane was launched by hand. Characteristics of OVLI-TA are summarized in Table 1.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** OVLI-TA UAV (a,b) based on the commercial airborne model Techpod, manufactured by HobbyUAV.</figcaption></figure>

<table><caption>**Table 1.** Technical characteristics of OVLI-TA.</caption><thead><tr><th>Wingspan/Wing Area</th><th>Fuselage Length</th><th>Dry Weigh/Payload</th><th>Max Cruise Speed</th><th>Ascent Range</th><th>Endurance</th></tr></thead><tbody><tr><td>2.6 m/0.4 m²</td><td>1.14 m</td><td>1.25/2.25 kg</td><td>28.3 m/s</td><td>3–10 m/s</td><td>1–2 h</td></tr></tbody></table>

### 2.2. Wind Vector

Calculating wind vector is a complex task because it requires a large number of parameters, each of them needing to be measured very accurately. The velocity of the air with respect to earth $V_w$ (wind) is the sum of the velocity of the aircraft with respect to earth $V_g$ (ground speed), and the velocity of the air with respect to the aircraft $V_{air}$ (airspeed):

$$V_w = V_g + V_{air} \qquad (1)$$

The three components of ground speed ($U_p$, $V_p$, and $W_p$ toward the east, north, and upward vertical directions, respectively) were provided by the GPS-IMU coupled system. The airspeed vector was first computed in the aircraft-related frame (x', y', and z', see Figure 2), and then projected into the earth-related frame with the Euler angles: true heading $\psi$, roll angle $\phi$, and pitch angle $\theta$. In the aircraft-related frame, the three components of the airspeed $V_{air}$ were calculated from the following three parameters: the velocity of the airflow $V_a$, the attack ($\alpha$), and sideslip ($\beta$) angles, with the following Equations [3]:

$$V_{air} = \begin{cases} -V_a D^{-1} \\ -V_a D^{-1} \tan \beta \\ -V_a D^{-1} \tan \alpha \end{cases} \qquad (2)$$

where D is expressed as:

$$D = \sqrt{1 + \tan \alpha^2 + \tan \beta^2} \qquad (3)$$
<figure><img src="image_5.png" /><figcaption>**Figure 2.** The earth and airplane coordinate systems were used for calculating the air velocity with respect to the earth referential. (x, y, z): earth coordinate system (i.e., east, north, upward); (x', y', z'): UAV coordinate system. The curved arrows indicate positive angles (see the text for angle definitions). Adapted from Reference [3].</figcaption></figure>

$V_a$ is computed from the dynamic pressure ($P_{dyn}$) measured with the Pitot tube, according to Bernoulli's equation:

$$V_a = \sqrt{\frac{2P_{dyn}}{\rho}} \qquad (4)$$

where $\rho$ is the air density. This equation is based on incompressibility assumption, which is still valid given the airspeed of the UAS (~13 m/s). In the above equation, $\rho$ is computed with perfect gas law, with the static pressure measured by the Pitot tube installed on the left side of the fuselage, and temperature measured in the hood as described below. The influence of humidity is not taken into account here.

The wind vector with respect to the earth reference was thus given as [3]:

$$U = -V_a D^{-1}[\sin\psi\cos\theta + \tan\beta(\cos\psi\cos\phi + \sin\psi\sin\theta\sin\phi) + \tan\alpha(\sin\psi\sin\theta\cos\phi - \cos\psi\sin\phi)] + U_p \quad (5)$$

$$V = -V_a D^{-1}[\cos\psi\cos\theta - \tan\beta(\sin\psi\cos\phi - \cos\psi\sin\theta\sin\phi) + \tan\alpha(\cos\psi\sin\theta\cos\phi + \sin\psi\sin\phi)] + V_p \quad (6)$$

$$W = -V_a D^{-1}[\sin\theta - \tan\beta\cos\theta\sin\phi - \tan\alpha\cos\theta\cos\phi] + W_p \quad (7)$$

where U, V, and W (Up, Vp, and Wp, respectively) are the wind components (UAV groundspeed components) pointing toward east, north, and up, respectively. In the above equations, it is assumed that the airflow and inertial measurements were made at the same place, which is a reasonable assumption since the IMU is located inside the nose of the aircraft. For attack, sideslip, pitch, and roll angles, small angle approximations can often be made (in particular during straight and level portions of flight patterns), and hence the above equations can be simplified as [3]:

$$U = -V_a \sin(\psi + \beta) + U_p \quad (8)$$

$$V = -V_a \cos(\psi + \beta) + V_p \quad (9)$$

$$W = -V_a \sin(\theta - \alpha) + W_p \quad (10)$$
## 2.3. Autopilot and Acquisition System

OVLI-TA flies by the 3DR Pixhawk autopilot that ensures stable navigation. The autopilot was composed of an independent data acquisition system with a GPS for measuring the absolute position. This was coupled with the ADIS16448 IMU including a triaxial digital gyroscope, a triaxial digital accelerometer, and a triaxial digital magnetometer, in order to measure the three components of the angular speed, the linear acceleration, and the magnetic field, respectively. In addition, it contained a digital barometer from which the fast altitude variations can be computed, and an embedded temperature sensor. All these parameters were delivered to the acquisition system at a rate of 100 Hz. The GPS data was updated with a frequency of 5 Hz. The Pixhawk telemetry system sent in-flight data to the ground station, and also recorded them on an SD card. The data acquisition system is described in Figure 3.

<figure><img src="image_6.png" /><figcaption>Figure 3. Diagram of connections between different devices on OVLI-TA.</figcaption></figure>

## 2.4. Meteorological Instrumentation of OVLI-TA

OVLI-TA is equipped with two temperature and humidity sensors SHT75 (Sensirion), recorded at 2.5 Hz. They are placed inside dedicated hoods (Figure 4). These 3D printed hoods were designed in such a way that droplets (if any) go straightforward by inertia through the inlet without any contact with the sensor, whereas the measured air parcels were deflected inside the body of the hood with a stable velocity relative to the sensor.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** The hoods where the temperature and humidity sensors are placed: (a) hoods installed on the top of OVLI-TA fuselage; and (b) 3D design of the hoods.</figcaption></figure>

The attack and sideslip angles were computed from differential pressures according to the system proposed in 1983 by Brown et al. [18] on a piloted airplane. A 3D printed five-hole hemispherical probe was situated on the nose of the drone (see Figure 4), the holes were connected to three differential pressure transducers (HCEM010, sampled at 100 Hz) by silicone tubes. There were five pressure ports: the center hole was in the longitudinal axis of the airplane, the four others were placed above and below, and to the right and left of the central hole, making an angle of 45° to the center of the sphere. The diameter of the hemisphere and holes were 10 cm and 1 mm, respectively. The first pressure sensor made it possible to measure the difference between the total pressure (from the center hole of the five-hole probe) and the static pressure (from the Pitot tube) in order to calculate the dynamic pressure and thus the airspeed, as a redundant system to the Pitot tube. The other sensors measured the differential pressure between the up and down (resp. left and right) holes in order to calculate the angle of attack $α$ (resp. angle of sideslip $β$). These sensors need a calibration procedure (see below). The acquisition rate of pressure data was 100 Hz. The payload was installed at the front of OVLI-TA. The overall equipment is presented in Figure 5.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Exploded diagram of the electronic devices on OVLI-TA.</figcaption></figure>

## 2.5. *The Five-Hole Probe Calibrations*

The wind tunnel test for the five-hole probe calibration was conducted at the IMFT (Institut de Mécanique des Fluides de Toulouse, Toulouse, France) in 2014. The range of wind speed produced by the wind tunnel was from 1 to 30 ms⁻¹, and the section of the tunnel was 2.4 m × 2.0 m.

The data acquisition system was in the same configuration as during flights. Figure 6 shows how the UAV was installed in the wind tunnel.
<figure><img src="image_9.png" /><figcaption>Figure 6. OVLI-TA during the wind tunnel test: (a) angle of attack testing; and (b) angle of sideslip testing [19].</figcaption></figure>

Theoretically, and according to Brown et al. [18], the equation for the attack angle $\alpha$ is:

$$ \Delta \hat{P}_{\alpha} = \hat{P}_{1} - \hat{P}_{3} = -\frac{9}{4} * \sin 2\omega * \sin 2\alpha, \qquad (11) $$

where $\hat{P}_1$ and $\hat{P}_3$ are the pressures at port 1 and 3 of the five-hole probe corresponding to the holes above and below the central hole, respectively, normalized by the dynamic pressure; and $\omega$ is the angle between the longitudinal axis of the 5-hole probe and the axis passing through one port and center of the sphere (Figure 7). The same equation stands for sideslip angle, replacing $\alpha$ with $\beta$ and ports 1 and 3 with ports 4 and 2, respectively. Therefore, by assuming low values for the attack and sideslip angles, and given that $\omega = \pi/4$, Equation (11) can be rewritten as:

$$ \begin{aligned} \alpha &= k_{\alpha}^{-1} * \Delta \hat{P}_{\alpha} \\ \beta &= k_{\beta}^{-1} * \Delta \hat{P}_{\beta}, \end{aligned} \qquad (12) $$

with a theoretical value of $k_{\alpha} = k_{\beta} = 4.5 \text{ rad}^{-1} (0.0785 \text{ deg}^{-1})$

<figure><img src="image_10.png" /><figcaption>Figure 7. The five-hole probe; (a) front view; and (b) right view.</figcaption></figure>

In fact, $k_\alpha$ and $k_\beta$ can depart from the theoretical value for several reasons; like a non-perfect hemispherical shape of the nose, a shift in the position of the holes with respect to the theoretical position, or a pressure perturbation related to the fine shape of the pressure ports. The exact value of the coefficients was therefore determined from wind tunnel tests. The airplane was fixed on a structure, the pitch of which was varied and measured. The principle of the calibration was to consider that in the wind tunnel flow, when the roll is zero, the pitch and attack angle must be identical. Similarly, when the airplane was rolled by 90°, the pitch of the structure was equal to the sideslip angle.

The angle of attack varied between -5° and 5°, with a zero angle of sideslip, so as to determine the sensitivity factor $k_\alpha$. The relationship between the differential normalized pressure and IMU's angle of attack is shown in Figure 8a. From this plot we deduced that the sensitivity factor of the attack angle was $k_\alpha = 0.0662 \text{ deg}^{-1}$. Similarly, the angle of sideslip varied between -5° and 5° by maintaining the angle of attack equal to zero, in order to determine the sensitivity factor $k_\beta$. The relationship between
the differentially normalized pressure and the angle of sideslip was $k_β = 0.0658 \text{ deg}^{-1}$ as shown in Figure 8b.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** (a) Left plot: normalized differential pressure as a function of attack angle; and (b) right plot: normalized differential pressure as a function of sideslip angle.</figcaption></figure>

## 2.6. Flight Sites and Patterns

OVLI-TA flew on two different sites. Firstly, flight tests were conducted on March 24, 2016 in Lannemezan, in the south of France. This site was the instrumental platform of the "Laboratoire d'Aérologie", belonging to the Pyrenean Platform for the Observation of the Atmosphere (P2OA, p2oa.aero.obs-mip.fr). A 60 m tower was equipped at three levels (30 m, 45 m, and 60 m) for turbulence measurements with sonic anemometers that gave us the three wind components and the so-called sonic temperature at a frequency of 10 Hz. In addition to that, relative humidity and temperature were measured at these levels, and recorded at 1 Hz. The flight pattern chosen was a succession of straight and level paths back and forth, which was the most appropriate type of flight plan for atmospheric turbulence measurements (Figure 9a). Each straight line of flight was 300–400 m long and lasted around 30 s; hence 61 sequences of flight in a straight line were selected and used for turbulence calculation. OVLI-TA flew for almost 1 h, between 60 m and 130 m but not exactly in the footprint of the tower.

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Flight patterns: (a) Flight conducted in Lannemezan site. The 60 m tower is situated around 200 m to the south of the lower left corner of the picture (length of straight flight: around 500 m); (b) flight conducted in DACCIWA site on June 29, 2016 (Radius circle: 100 m); and (c) flight conducted in DACCIWA site on July 15, 2016 (radius circle: 100 m; and length of straight flight: 700 m).</figcaption></figure>

Afterwards, the drone participated in the DACCIWA campaign in West Africa (Savè, Benin). The site was equipped with radio soundings, remotely piloted airplane system, cloud camera, cloud radars, Doppler lidar, radiometer, sodar, UHF wind profiler, and meteorological/turbulence/chemistry ground stations. The ground-based field campaign took place from June 14 to July 30, 2016 [20]. We
chose to present results from the flights on June 29, 2016 and July 15, 2016. On June 29, 2016, the drone
flew between 10 m and 700 m agl, with ascending and helical patterns (Figure 9b) in order to compare
OVLI-TA profiles of wind, temperature, and humidity with those obtained with radiosondes regarded
as a reference. On July 15, 2016, the OVLI-TA flight pattern consisted of straight and level runs of
similar lengths as shown in Figure 9c. The drone flew between 10 m and 1000 m agl, with ascending
helices and descending straight lines.

3. Results

3.1. *Mean Wind and ABL Profiles*

The 3D representation (Figure 10) of humidity, temperature, and wind vector provides an overview of the full flight pattern in Lannemezan. Temperature and humidity appeared consistent during the flight, with little variation along each individual straight leg. However, the altitude range analyzed was not wide enough to extract the variations, which should be observed on vertical profiles. That was the reason why we focused the analysis of this flight on wind and turbulence parameters, the temperature and humidity profiles being analyzed from the DACCIWA flights (see later).

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Relative humidity (a), temperature (b), and wind speed (c) measured by OVLI-TA during a flight in Lannemezan on 24 March 2016. The time window is from 12:35 to 13:27 (UTC).</figcaption></figure>

As mentioned above, wind (speed and direction) was among the most critical parameters extracted
from airplane measurements, because it required a lot of variables to be measured (see previous
Equations (5)–(7)). Furthermore, the parameters to be computed (wind components) were generally one order of magnitude lower than some variables in the equations, like the groundspeed or airspeed components. The resulting errors could thus be considerable, and evolve with the orientation of the airplane with respect to the wind’s direction. For example, overestimating true airspeed leads to an underestimation (resp. overestimate) of wind speed when the airplane is facing the wind (resp. has the wind at its back), whereas the impact on the wind direction estimates is much lower. On the other hand, when the airplane is flying crosswind, an error on the true airspeed impacts (at the first order) the estimate of the wind direction, whereas the wind speed estimate is biased through an error on the sideslip angle. If one assumes that the mean wind does not vary much during a back and forth sequence, one can thus extract information about the bias of those variables from the difference in the wind computed on the two paths. These biases can be computed from appropriate flight sequences (like back and forth runs), and the quality of the wind estimates thus improved.

In the flight presented on Figure 9, OVLI-TA flew 61 straight sequences over the same 500 m long ground track, with true headings around 210° and 360° (hereafter called “SSW runs” and “N runs”, respectively), at heights between 60 and 130 m. The mean wind speed and direction was computed on each of these sequences, and the global scatter of the results was analyzed. (Thereafter, the terms “mean wind direction”, that will be used several times for simplicity, represents the direction calculated from the mean wind components in an earth-related frame). Figure 11 presents the mean wind speed and direction for the 61 sequences according to the heading of the drone. Without any adjustment, it can be remarked that the cluster of SSW runs showed wind speed values higher on average than the cluster of N runs. The difference was less significant on the wind direction. Taking into consideration this flight configuration and once a correction (an offset equal to 0.055 rad) had been applied to the sideslip angle, the two clusters of wind direction get comparable shapes. The impact on the wind direction was weaker, with a reduction in the overall range of the computed values (from [263°–337°] without correction to [269°–327°]). The origin of this offset probably lies in a misalignment between the IMU and the five-hole probe, due to the successive mounting and dismounting of the system on the airplane after the wind tunnel calibration.

<figure><img src="image_14.png" /><figcaption>**Figure 11.** Wind speed (a) and wind direction (b) as a function of heading angle, estimated on the 61 flight sequences before (in blue) and after (in red) the bias correction.</figcaption></figure>

The improvement brought by the correction was evaluated by comparing the standard deviation of the wind speed and direction before and after correction (Table 2). Wind direction is a circular variable that requires a specific method for estimating its standard deviation, because we have to cope
with the passage at 0/360° (or +/- 180°). We used the method described by Yamartino [21], who proposed the following relations:

$$ \sigma_P = \sin^{-1}(\varepsilon)^*[1 + 0.1547*\varepsilon^3] \quad (13) $$

$$ \varepsilon = 1 - (S_a^2 + C_a^2) \quad (14) $$

$$ S_a = n^{-1}\Sigma_i \sin P_i \quad (15) $$

$$ C_a = n^{-1}\Sigma_i \cos P_i \quad (16) $$

where σ is the standard deviation of wind direction, and Pᵢ is the wind direction on the ith of the 61 legs.

<table><caption>**Table 2.** Standard deviation of corrected and non-corrected values of wind speed and wind direction computed on the 61 flight sequences by OVLI-TA.</caption><thead><tr><th colspan="2">Wind Speed (m/s)</th><th colspan="2">Wind Direction (Degrees)</th></tr><tr><th>Non-Corrected</th><th>Corrected</th><th>Non-Corrected</th><th>Corrected</th></tr></thead><tbody><tr><td>1.31</td><td>0.98</td><td>15.0</td><td>13.8</td></tr></tbody></table>

The values in Table 2 confirm that the correction applied to the sideslip angle had a major impact on the wind speed estimate, with a standard deviation reduced by 25% (from 1.31 m/s to 0.98 m/s) and a minor impact on the wind direction estimate by reducing the standard deviation by only 8% (from 15.0° to 13.8°). The flight configuration was close to the “cross-wind” one, and the adjustment on the sideslip angle had a direct impact on the wind speed estimate. The overall performance was regarded as very satisfactory for the wind estimate with a mobile platform.

Once optimized, thanks to the flight pattern, the north and east wind speeds measured by OVLI-TA were compared to those measured by the sonic anemometer mounted at the 60 m level on the tower (Figure 12). For each platform, we present two time series: the instantaneous values, which contained the turbulent fluctuations; and values averaged over the straight and level runs of OVLI-TA, and averaged over the same corresponding periods on the tower. On the instantaneous time series, we observed that the two measurements were consistent, even though the OVLI-TA signal presents a higher variability, related in part to poorer wind estimates during half turns between two straight flights; in part due to the altitude range flown by the drone (between 60 and 130 m, whereas the tower observations are made at 60 m only), and in part to the difference between the areas sensed by the two platforms (given the wind direction, the footprint of the tower lies to its WNW, whereas the flight track is to the north of the tower). Furthermore, if we look at the values averaged over the periods of time of the straight sequences (green and blue dots on Figure 12), we find a very satisfactory agreement, given the scatter induced by the difference between the two platforms regarding their footprint, and the altitude of air masses sensed.

The absolute accuracy of the sensor used in the tower was lower than 0.1 m/s (https://www.campbellsci.com/csat3), i.e., an order of magnitude lower than what can be expected with a mobile platform like a drone. For each of the 61 samples (blue and green dots in Figure 12), we therefore compared the east and north mean components between the drone and tower. For each component, the difference between the two platforms was quantified through the mean absolute error (MAE) and the standard deviation of the difference computed on the 61 values. MAE is defined as:

$$ \text{MAE} = \sum_{i=1}^{n} \frac{|U_{\text{Mast},i} - U_{\text{OVLI},i}|}{n} $$
<figure><img src="image_15.png" /><figcaption>**Figure 12.** Continuous lines: north and east wind speeds measured by OVLI-TA (in red) and on the tower at 60 m (in black), during the whole flight. Dots: mean of north and east wind speeds for each straight and level sequence of the flight measured by OVLI-TA (in green), and on the tower at 60 m for the same periods of time (in blue).</figcaption></figure>

Table 3 summarizes the obtained results. The order of magnitude of these errors was 1 m/s. Though the standard deviation of the difference between the two platforms was identical for the two components, we observed a higher MAE on the north component (1.2 m/s) compared to that on the east component (0.8 m/s). We propose two hypotheses to explain such a difference: first this could be due to the distinct local areas of where the 60 m tower is situated and where the drone was flying, because OVLI-TA was not exactly flying in the footprint of the 60 m tower (the drone tracks are some hundreds of meters away from the tower) and the levels of flight were between 60 m and 130 m (whereas the tower observations are at 60 m only); second, a bias might remain in the drone north wind estimate resulting from a non-perfect calibration of a parameter (in the present case, given that the aircraft headings are close to the N or S orientations, the true airspeed might be suspected). However, we considered that these errors belong to the acceptable range of the errors targeted for the wind estimates with our mobile platform, and must be regarded as the upper bound of the drone error.
<table><caption>**Table 3.** Mean absolute error of east and north wind speed. Standard deviation of the difference between east (or north) wind component measured respectively by the 60 m mast and OVLI-TA.</caption><thead><tr><th></th><th>East Wind Component (m/s)</th><th>North Wind Component (m/s)</th></tr></thead><tbody><tr><td>MAE</td><td>0.8</td><td>1.2</td></tr><tr><td>σ (U<sub>Mast</sub>-U<sub>OVLI</sub>)</td><td>1.0</td><td>1.0</td></tr></tbody></table>

After flights in Lannemezan, OVLI-TA was operated during the DACCIWA campaign for profiling temperature, humidity, wind speed, and direction in the ABL between 10 and 1000 m agl. Figure 13 illustrates the observations of temperature and wind during the flight on July 15, 2016 that started at 16:45 UTC and lasted for 40 min. Stacked back and forth sequences were made on a racetrack pattern, and a continuous profile was obtained from a helical pattern.

<figure><img src="image_16.png" /><figcaption>**Figure 13.** Air temperature (a) and wind vector (b) measured on July 15, 2016 at 16:45 (UTC).</figcaption></figure>

The quality of the measurements was checked by comparing profiles during the ascending and descending parts of the flight. This is shown in Figure 14, for the wind speed and direction, relative humidity, and temperature. The differences in relative humidity below 600 m were within the expected accuracy of the sensor. The large fluctuations above 600 m were related to the so-called entrainment zone, which lies in the upper part of the ABL, and contains a mixture of air parcels coming from the well-mixed layer below, and from the upper layer above, the latter being characterized by much warmer and drier conditions. The difference in the two temperature profiles was also the highest in this area, for the same reasons. In the mixed layer below, the temperature difference between the two profiles was of the order of 0.2–0.3 °C in the 100–600 m layer, and could reach 0.5–0.6 °C in the 50–100 m layer, which was also within the absolute accuracy of the sensor. Note that on this flight, the ascent speed of OVLI-TA was 1.5 m/s and the descent speed was 0.7 m/s.
<figure><img src="image_17.png" /><figcaption>Figure 14. Profiles of temperature (a), relative humidity (b), wind speed (c), and wind direction (d) measured during ascending (in green) and descending (in blue) flight sequences on July 15, 2016.</figcaption></figure>

During the campaign, radiosondes were launched under meteorological balloons. This system measured the profiles of temperature, relative humidity, and wind. The ascending speed of the radiosonde was about 5 m/s, the 1st kilometer of the atmosphere was therefore sensed in around 200 s, which can be considered as “instantaneous” with regard to atmospheric variability. Radiosondes are generally considered as a reference measurement system, though moisture observations were questioned for specific conditions. During the OVLI-TA flight of June 29, a radiosonde was launched at 16:59 UTC, whereas two profiles (ascending and descending) were made by the drone between 17:12 and 17:35 UTC. The comparison of the profiles measured by the two systems is shown in Figure 15. The radiosonde used (http://www.meteomodem.com/docs/fr/Brochure-m10.pdf) had a time resolution of 1 s, and the response time of the sensors was 1 s for temperature and 2 s for humidity; consequently we can expect a vertical delay in the order of 5 m for temperature and 10 m for humidity. Due to this delay, and given the vertical gradients observed in the ABL, the profiles were therefore biased by ~+0.04 °C for temperature and ~−0.3 % for relative humidity. Such values were low enough to consider the radiosonde profiles as the reference against which the drone profiles have to be compared. The agreement was good for the three parameters; the main discrepancy being observed was in the 500–700 m layer, humidity, and wind speed/direction. However, we observed a difference between ascending and descending profiles which could be due to sensor response time. This can be estimated from the value of difference and the ascending and descending speeds during the profiles. To do that, we focused on altitude ranges chosen in the so-called “mixed layer”, i.e., above the atmospheric surface layer where the vertical and horizontal homogeneity is large, and well below the top the ABL where strong vertical gradients exist. For the retained flight sequences, only 12 min elapsed between the start of the ascent and the end of the descent, so we can assume a steady-state situation. The response time $\tau_x$ is thus computed as:

$$ \tau_x = \frac{x_{up} - x_{down}}{\left(\frac{dx}{dz}\right)\left(v_{up} - v_{down}\right)}, \qquad (17) $$
where *x* being either the temperature T or relative humidity RH. The indexes up and down refer to the ascending and descending phases, respectively; *v* is the vertical speed of the drone; and *z* is the altitude. *v*<sub>up</sub> and *v*<sub>down</sub> were 1.4 and -5.0 m/s, respectively. The resulting response times were 8.9 s and 14.8 s for T and RH, respectively. The corresponding delays in the observed profiles during the ascent phases were therefore ~12 m and ~21 m for temperature and relative humidity, respectively. Such values were around twice those computed for the radiosonde, considered as the reference, and the corresponding biases were +0.10 °C for temperature and -0.6 % for relative humidity, which lay well within the absolute accuracy of the sensors. This means that one can easily cope with the biases resulting from the response time of the sensors provided that the ascending or descending speed of the drone is low enough.

<figure><img src="image_18.png" /><figcaption>**Figure 15.** Same as Figure 14, but for the flight on 29 June 2016 and together with the observations obtained from a radiosonde (RS). (a) Temperature (°C), (b) relative humidity (%), (c) wind speed (m.s⁻¹), (d) Wind direction (degrees). Lines: radiosonde measurements (in blue), measurements of OVLI-Ta during ascent (in pink) and descent (in green).</figcaption></figure>

## 3.2. Turbulence

The flight on the Lannemezan site was suitable to investigate the capability for turbulence observation. The three wind components were computed at a rate of 100 Hz. In this section, we
investigate the maximum frequency at which the signals can be used. Figure 16a–c shows the energy spectra of the north (V), east (U), and vertical (W) wind components, for a straight and level sequence. Given the meteorological conditions, and that the flight was performed in a convective boundary layer, a 3D turbulence was expected, with a –5/3 power slope in the inertial sub-range. This behavior was well observed on the spectra, except for the highest frequencies where the signal was noisy. We cannot investigate in detail the possible (and multiple) origins of this noise. Such a shape on spectra is frequently observed on wind signal computed with aircraft platforms, because a flying system constitutes an adverse environment for fine and high-frequency measurements. In order to improve the statistics, we averaged the individual energy spectra computed on 61 straight and level sequences. An identical period of time was retained on each sequence, in order to have the same frequencies on the individual spectra. The resulting spectrum for the vertical wind was thus presented (up to 10 Hz) in Figure 16d. We highlight the result on the vertical wind, because this component was involved in the estimate of the vertical turbulent fluxes through the eddy-covariance technique, as the horizontal turbulent fluxes were much lower in the ABL. The –5/3 slope, characteristic of the inertial sub-range, was observed on this spectrum up to 10 Hz.

<figure><img src="image_19.png" /><figcaption>**Figure 16.** (a) Spectra of vertical wind speed, (b) north wind speed, (c) east wind speed, and (d) mean spectrum of vertical wind speed for all sequences. Note that the frequency scale is different on (d).</figcaption></figure>

In order to go further into the qualification of turbulence measurements, we made a comparison with the 60 m tower observations. As mentioned above, there was a sonic anemometer installed at the top of the tower, measuring the three wind components at a rate of 10 Hz. We focused on the vertical wind, for the reasons mentioned above. In order to compare between tower and drone spectra, we must go back to spatial scales as observed by the two platforms. In other words, we must convert a frequency $f$ into the corresponding wavenumber $k$ assuming Taylor's hypothesis $k = 2\pi f/u$, where $u$ is the mean speed of the flow with respect to the sensor. Hence, we considered the mean airspeed for OVLI-TA, and the mean wind speed for the tower. Figure 17a shows the comparison for the first
sequence of straight flight using the Welch periodogram. The period of time T<sub>t</sub>, over which the tower observations were considered is defined according to the distance in the air mass sampled by the drone during this sequence: T<sub>t</sub> = T<sub>d</sub> × V<sub>a</sub> × V<sub>w</sub><sup>-1</sup>, where V<sub>a</sub> is the mean true airspeed measured by the drone, T<sub>d</sub> is the duration of the drone path, and V<sub>w</sub> is the mean wind speed measured on the tower. Given the differences between the two air masses sampled by the two platforms, the spectra cannot be identical. However, the general shape and the characteristic scales should be similar. This is confirmed in Figure 17a,b where it can be seen that the two platforms capture the same level of energy, even in the inertial sub-range, up to wavenumbers of 4–5 rad/m. This means that the drone was able to correctly sample the ABL down to spatial scales of the order of 1 m, i.e., was able to capture the scales, which significantly contribute to the turbulent kinetic energy.

<figure><img src="image_20.png" /><figcaption>**Figure 17.** (a) Vertical wind spectrum as a function of wavenumber computed during the first straight sequence. OVLI-TA measurements are in red and the corresponding 60 m tower in blue. The spectra are computed using Welch's method. (b) The scatter plot of the power density of the mast against OVLI-TA.</figcaption></figure>

A scatter-plot provides a more in-depth comparison of the two spectra. However, since the Taylor's hypothesis is used to convert the frequency spectra into wavenumber spectra, and given that flow speeds used for this conversion are different for the two platforms (the true airspeed for the aircraft and the wind speed for the tower), the wavenumbers at which the spectral energy is computed are different. Therefore, the energies cannot be directly compared. To cope with this difficulty, we computed for each platform the average energy (on the Welch spectra) in non-overlapping wavenumber ranges identical for the two platforms. Furthermore, the scatter was quantified with the standard deviation in each wavenumber bin. The result is presented in Figure 17b. The agreement was quantified through the determination coefficient (0.76). It can be seen that the energy computed agreed
well between the two platforms, except at the lowest energies for the aircraft where the values were higher. This is related to the noisy tail of the aircraft spectrum, as mentioned previously.

## 4. Discussion

The development of the small, unmanned aerial system OVLI-TA for probing the ASL/ABL makes it possible to explore atmospheric layers as a complementary method to towers and other platforms. OVLI-TA is easy to operate, in a category of UASs for which the air traffic rules remain acceptable, avoiding the constraints we have to deal with for heavier airplanes. Nevertheless, the payload capacity (~2 kg) allows us to embark high-performance instrumentation, composed of atmospheric sensors, an inertial platform, and an efficient autopilot. This instrument package has revealed its capacity to measure the temperature, humidity, and wind components all along the aircraft's trajectory. In particular, the three wind components, although hard to compute on mobile platforms, can be estimated not only for their mean value, but also for the turbulent fluctuations. In particular, the instantaneous values of the airspeed vector, required for these calculations, are computed from a multi-hole probe installed on the nose of the airplane, and calibrated in a wind-tunnel for the estimates of attack and sideslip angles.

The quality of the measurements was determined from self-consistency methods and also by comparison with other measurements on platforms considered as a reference. For temperature and humidity profiles, we checked the similarity between the observations collected during ascending and descending flight profiles, and compared the observations to those obtained with radiosonde profiles. From the difference between up- and down-ward profiles performed in a short time range, we estimated the response times of the temperature and relative humidity probes, as around nine and 15 s, respectively. With a profiling vertical speed around 1 m/s, this allows us to lower the bias of the observed profiles at a level twice that of the radiosonde under meteorological balloons. For wind speed and direction, the variability of the estimates during back and forth runs, allowed us to adjust calibration coefficients and reduce the scatter of the estimates. Comparison with the 60 m tower observations revealed an excellent agreement, even though the air mass sensed by the two platforms was not identical. The absolute accuracy of the wind components is better than 1 m/s, because this value, deduced from the difference between drone and tower estimates, involves at least in part the atmospheric variability. The turbulence spectra, computed on straight and level flight sequences, revealed the expected slope in the inertial sub-range, up to a frequency of about 10 Hz. We compared these spectra with those computed from the tower observations, once the transformation of frequency into wavenumber had been done according to the sample velocity of each platform. The spectra agree as well as it can be expected, given the difference between the footprints of the two platforms. In addition, a scatter plot of power density of the mast against OVLI-TA confirms that the agreement is good up to 10 Hz.

The OVLI-TA platform is able to sample ABL turbulence with a spatial resolution of around 1 m. We can question whether this allows us to capture most of the energetic turbulence. In fact, 3D turbulence spectra in the ASL, when represented in the unit of variance, are characterized by a peak occurring at a certain frequency, with the energy decreasing on each side of the peak [22]. We can consider that the turbulence observation is appropriate when the sampling frequency of the measures is at least one order of magnitude higher than that of the peak. Starting from the surface, this peak shifts continuously towards lower frequencies (larger wavelengths) when the observation height increases, because the size of eddies grows with the distance from the surface. For convective, daytime conditions, we can expect a peak wavelength of the order of 3–5 z, where z is the height above the ground [22,23]. That means that at z = 10 m, which is the lowest height at which the drone could be reasonably operated, we are able to resolve eddies of a size 30–50 times smaller than the most energetic ones, and therefore most of the significant turbulence energy can be captured.
## 5. Conclusions

OVLI-TA is a small affordable UAV, developed in France. It was instrumented for profiling the ABL and measuring turbulence. It flew in two different sites, Lannemezan (France) and Savé (Bénin). From the comparison between OVLI-TA and respectively the 60 m tower and "radiosondes" measurements, we demonstrated its efficiency in measuring the turbulence up to 10 Hz, which is equivalent to a horizontal resolution of around 1 m. In spite of the quite slow temperature and moisture sensors, the UAV is able to give reliable profiles in the ABL provided that the vertical speed of the aircraft is kept slow (~1 m/s).

In the near future, we plan to extend OVLI-TA's capabilities with fast response temperature and moisture sensors. This will give us the possibility to compute kinematic heat and moisture fluxes, which are essential terms in ABL dynamics. Furthermore, we are developing a new instrumentation on a heavier drone, to extend the payload to 5 kg and the endurance to 5-6 h. This will allow us to embark high-performance instruments, in particular for turbulence measurements, and to fly in more adverse conditions, such as strong winds and those with high turbulence. However, the operation of this new system would be more restrictive, especially with regard to flight authorization and operation. We will therefore benefit from the complementarity of the two systems.

**Author Contributions:** Conceptualization, P.D. and P.M.; formal analysis, S.A-S.; investigation, P.M., I.C., M.L., and P.D.; methodology, P.D.; software, S.A-S. and P.M.; supervision, P.D., P.P., and P.M.; writing—original draft preparation, S.A-S.; writing—review and editing, S.A-S., P.D., P.P., M.L, and P.M.

**Funding:** The DACCIWA project received funding from the European Union Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 603502. OVLI-TA development was supported by the CNRS-INSU LEFE program and Observatoire Midi-Pyrénées. P2OA facilities and staff are funded and supported by the University Toulouse III Paul Sabatier, France, and CNRS (Centre National de la Recherche Scientifique).

**Acknowledgments:** We are indebted to V. Gervaise and many other trainees who contributed to OVLI-TA's development. We also express our gratitude to J.M. Thouery and the "Agence nationale pour la formation professionnelle des adultes" who offered us several internships. Finally, we thank the two reviewers for their constructive comments.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## References

1. Kaimal, J.C.; Finnigan, J.J. Flow over flat uniform terrain. In *Atmospheric Boundary Layers Flows: Their Structure and Measurement*; Oxford University Press: New York, NY, USA, 1994; pp. 3–26.

2. Stull, R.B. Mean boundary layer characteristics. In *An Introduction to Boundary Layer Meteorology*; Kluwer Academic Publishes: Dordrecht, The Netherlands, 1988; pp. 1–26, ISBN 90-277-2768-6.

3. Lenschow, D.H. Aircraft Measurements in the boundary Layer. In *Probing the Atmospheric Boundary Layer*, 1st ed.; Lenschow, D.H., Ed.; American Meteorological Society: Boston, MA, USA, 1986; pp. 39–55.

4. Elston, J.; Argrow, B.; Stachura, M.; Weibel, D.; Lawrence, D.; Pope, D. Overview of small fixed-wing unmanned aircraft for meteorological sampling. *J. Atmos. Ocean. Technol.* **2015**, *32*, 97–115. [CrossRef]

5. Reineman, B.D.; Lenain, L.; Statom, N.M.; Melville, K.W. Development and Testing of Instrumentation for UAV-Based Flux Measurements within Terrestrial and Marine Atmospheric Boundary Layers. *J. Atmos. Ocean. Technol.* **2013**, *30*, 97–115. [CrossRef]

6. Bärfuss, K.; Pätzold, F.; Altstädter, B.; Kathe, E.; Nowak, S.; Bretschneider, L.; Bestmann, U.; Lampert, A. New Setup of the UAS ALADINA for Measuring Boundary Layer Properties, Atmospheric Particles and Solar Radiation. *J. Atmos.*** **2018**, *9*, 28. [CrossRef]

7. Van Den Kroonenberg, A.; Martin, T. Measuring the Wind Vector Using the Autonomous Mini Aerial Vehicle M²AV. *J. Atmos. Ocean. Technol.* **2008**, *25*, 1969–1981. [CrossRef]

8. Wildmann, N.; Ravi, S.; Bange, J. Towards higher accuracy and better frequency response with standard multi-hole probes in turbulence measurement with remotely piloted aircraft (RPA). *J. Atmos. Meas. Tech.* **2014**, *7*, 1027–1041. [CrossRef]

9. Wildmann, N.; Hofsäß, M.; Weimer, F.; Joos, A.; Bange, J. MASC—A small Remotely Piloted Aircraft (RPA) for wind energy research. *J. Adv. Sci. Res.* **2014**, *11*, 55–61. [CrossRef]
10. Witte, B.M.; Singler, R.F.; Bailey, S.C.C. Development of an Unmanned Aerial Vehicle for the Measurement of Turbulence in the Atmospheric Boundary Layer. *J. Atmos. 2017*, *8*, 195. [CrossRef]

11. Reuder, J.; Brisset, P.; Jonassen, M.; Müller, M.; Mayer, S. The Small Unmanned Meteorological Observer SUMO: A new tool for atmospheric boundary layer research. *J. Meteorol. 2009*, *18*, 141–147. [CrossRef]

12. Reuder, J.; Baserud, L.; Jonassen, M.O.; Kral, S.T.; Müller, M. Exploring the potential of the RPA system SUMO for multipurpose boundary-layer missions during the BLLAST campaign. *J. Atmos. Meas. Tech.* **2016**, *9*, 2675–2688. [CrossRef]

13. Båserud, L.; Reuder, J.; Jonassen, M.O.; Kral, S.T.; Paskyabi, M.B.; Lothon, M. Proof of concept for turbulence measurements with the RPAS SUMO during the BLLAST campaign. *J. Atmos. Meas. Tech.* **2016**, *9*, 4901–4913. [CrossRef]

14. Kral, S.T.; Reuder, J.; Vihma, T.; Suomi, I.; O'Connor, E.; Kouznetsov, R.; Wrenger, B.; Rautenberg, A.; Urbancic, G.; Jonassen, M.O.; et al. Innovative Strategies for Observations in the Arctic Atmospheric Boundary Layer (ISOBAR)—The Hailuoto 2017 Campaign. *J. Atmos.* **2018**, *9*, 268. [CrossRef]

15. Jacob, J.D.; Chilson, P.B.; Houston, A.L.; Weaver Smith, S. Considerations for Atmospheric Measurements with Small Unmanned Aircraft Systems. *J. Atmos. 2018*, *9*, 252. [CrossRef]

16. Zhou, S.; Peng, S.; Wang, M.; Shen, A.; Liu, Z. The Characteristics and Contributing Factors of Air Pollution in Nanjing: A Case Study Based on an Unmanned Aerial Vehicle Experiment and Multiple Datasets. *J. Atmos.* **2018**, *9*, 343. [CrossRef]

17. Schuyler, T.J.; Guzman, M.J. Unmanned Aerial Systems for Monitoring Trace Tropospheric Gases. *J. Atmos.* **2017**, *8*, 206. [CrossRef]

18. Brown, E.N.; Friehe, C.A.; Lenschow, D.H. The use of pressure fluctuations on the nose of an aircraft for measuring air motion. *J. Clim. Appl. Meteorol.* **1980**, *22*, 171–180. [CrossRef]

19. Cernov, I. Study of the Atmospheric Turbulence Measurement by a Probe on a Light Drone. Master's Thesis, ISAE-ENSICA, Toulouse, France, 2014.

20. Kalthoff, N.; Lohou, F.; Brooks, B.; Jegede, G.; Adler, B.; Babić, K.; Dione, C.; Ajao, A.; Amekudzi, L.K.; Aryee, J.N.A.; et al. An overview of the diurnal cycle of the atmospheric boundary layer during the West African monsoon season: Results from the 2016 observational campaign. *J. Atmos. Chem. Phys.* **2018**, *18*, 2913–2928. [CrossRef]

21. Yamartino, R.J. A comparison of several "Single Pass" estimators of the standard deviation of wind direction. *J. Clim. Appl. Meteorol.* **1984**, *23*, 1362–1366. [CrossRef]

22. Stull, R.B. Some mathematical and conceptual tools: Part2. Time series. In *An Introduction to Boundary Layer Meteorology*; Kluwer Academic Publishes: Dordrecht, The Netherlands, 1988; pp. 295–344.

23. Kaimal, J.C.; Finnigan, J.J. Spectra and cospectra over flat uniform terrain. In *Atmospheric Boundary Layers Flows: Their Structure and Measurement*; Oxford university press: New York, NY, USA, 1994; pp. 33–63, ISBN 0-19-506239-6.

<figure><img src="image_21.png" /></figure>