Articles | Volume 14, issue 3
Research article
26 Mar 2021
Research article |  | 26 Mar 2021

Suitability of fibre-optic distributed temperature sensing for revealing mixing processes and higher-order moments at the forest–air interface

Olli Peltola, Karl Lapo, Ilkka Martinkauppi, Ewan O'Connor, Christoph K. Thomas, and Timo Vesala

The suitability of a fibre-optic distributed temperature sensing (DTS) technique for observing atmospheric mixing profiles within and above a forest was quantified, and these profiles were analysed. The spatially continuous observations were made at a 125 m tall mast in a boreal pine forest. Airflows near forest canopies diverge from typical boundary layer flows due to the influence of roughness elements (i.e. trees) on the flow. Ideally, these complex flows should be studied with spatially continuous measurements, yet such measurements are not feasible with conventional micrometeorological measurements with, for example, sonic anemometers. Hence, the suitability of DTS measurements for studying canopy flows was assessed.

The DTS measurements were able to discern continuous profiles of turbulent fluctuations and mean values of air temperature along the mast, providing information about mixing processes (e.g. canopy eddies and evolution of inversion layers at night) and up to third-order turbulence statistics across the forest–atmosphere interface. Turbulence measurements with 3D sonic anemometers and Doppler lidar at the site were also utilised in this analysis. The continuous profiles for turbulence statistics were in line with prior studies made at wind tunnels and large eddy simulations for canopy flows. The DTS measurements contained a significant noise component which was, however, quantified, and its effect on turbulence statistics was accounted for. Underestimation of air temperature fluctuations at high frequencies caused 20 %–30 % underestimation of temperature variance at typical flow conditions. Despite these limitations, the DTS measurements should prove useful also in other studies concentrating on flows near roughness elements and/or non-stationary periods, since the measurements revealed spatio-temporal patterns of the flow which were not possible to be discerned from single point measurements fixed in space.

1 Introduction

The majority of the interaction between the atmosphere and Earth's surface takes place in a shallow air layer termed the atmospheric boundary layer (ABL). Insights on the atmospheric mixing processes in this layer are required in order to gain a better understanding on ecosystem–atmosphere feedbacks, air quality and weather-forecasting-related issues. Studies near the surface typically rely on time series analysis, since spatial details of the mixing close to the ground are difficult to measure with conventional instrumentation. However, similarity theories underlying the analysis of observations and models are posed in length scales, calling for spatially explicit sampling. In addition, turbulence statistics or spatial details of different flow structures are typically derived from a time series of observations by assuming ergodic hypothesis or Taylor's frozen turbulence hypothesis (Taylor1938). Yet, it has been recognised that these hypotheses are not universally valid (Mahrt et al.2009; Thomas2011; Higgins et al.2012, 2013; Cheng et al.2017).

Atmospheric boundary layer flows feature other flow modes besides turbulence, and these flow modes exhibit spatial patterns which are not directly related to their temporal counterparts. For instance, transient submesoscale motions may occur in the weak-wind stable boundary layer (Mahrt2014), which can travel in the opposite direction of the mean flow (Zeeman et al.2015) and interact with turbulence (Kang et al.2015; Sun et al.2015; Mahrt and Thomas2016; Vercauteren et al.2016), inflicting intermittent mixing and transport of gases, heat and momentum. Mechanistic understanding of these non-turbulent motions and related processes are limited, partly due to the missing spatial information on the flow.

Besides submeso motions in stable conditions, the flow in the roughness sublayer (RSL), by definition, exhibits spatial patterns that call for spatial sampling. For instance, flows within and above forest canopies are dominated by large coherent structures (Finnigan2000; Thomas and Foken2007a, b; Thomas et al.2008; Finnigan et al.2009) that are in continuous interaction with the surface roughness elements; this interaction can cause persistent spatial variability in the flow and turbulent transport (Bohrer et al.2009; Schlegel et al.2015) raising questions about the experimental approach deployed in global measurement networks monitoring the turbulent transport of gases above ecosystems (Baldocchi2014). Surface spatial heterogeneities also dominate the flow properties in the urban RSL, with urban surfaces exhibiting large variability in the surface heat flux and, hence, the thermal production of turbulence (Barlow2014). In general, the analysis of the effect of abrupt edges, or irregular discontinuities in surface characteristics, on the flow requires spatially explicit sampling.

Remote sensing instrumentation is capable of resolving the spatial detail of the flow (e.g. Newsom et al.2008; Träumner et al.2012); yet, small-scale features close to the ground are typically missing, and measurements near obstacles such as trees and buildings are not feasible. On the other hand, conventional precise in situ measurements circumvent these limitations but are fixed in space (e.g. at measurement towers). Consequently, spatial details of the flow can be deduced only by assuming Taylor's frozen turbulence hypothesis. Hence, there is an evident observational gap between the conventional measurement techniques (remote sensing vs. in situ), which results in limited understanding of flow processes falling in this observational gap. In situ measurements on moving platforms, such as tethered balloons or unoccupied aerial vehicles (UAVs) (Poulos et al.2002; Frehlich et al.2008; Higgins et al.2018; Egerer et al.2019), may partly fill the gap, yet, with these techniques, only non-continuous campaign type measurements are possible, resulting in low temporal coverage and representativeness. Furthermore, measurements close to or in-between roughness elements (trees and buildings) are typically not possible with these mobile platforms.

Distributed temperature sensing (DTS) has been utilised in environmental research since the first studies of Selker et al. (2006) and Tyler et al. (2009). The measurement method provides spatially continuous profiles along a fibre-optic cable which can be freely distributed in the measurement domain. The DTS data are provided at a similar temporal and spatial resolution along the cable to that of conventional in situ measurements, lending direct comparison between the measurement techniques and complementary analyses of joint measurements with multiple techniques. The measurement technique relies on optical time domain reflectometry and measurement of Raman backscattering of a light pulse travelling in the fibre-optic cable (Dakin et al.1985; Selker et al.2006). Due to its ability to provide spatio-temporal information at scales (down to 1 s and 25 cm) that are commensurate with the scales prevalent near the surface, the measurement method shows promise in answering many persistent, unanswered questions related to near-surface flow.

A growing body of research has already utilised DTS measurements in atmospheric studies (Keller et al.2011; Thomas et al.2012; Euser et al.2014; de Jong et al.2015; Sayde et al.2015; Zeeman et al.2015; Pfister et al.2017; Higgins et al.2018; Schilperoort et al.2018; Higgins et al.2019; Izett et al.2019; Mahrt et al.2019; Pfister et al.2019). The bulk of the studies have concentrated on nocturnal flows near the surface (Keller et al.2011; Thomas et al.2012; Zeeman et al.2015; Pfister et al.2017; Izett et al.2019; Mahrt et al.2019; Pfister et al.2019), and a few have concentrated on transition periods (Higgins et al.2018, 2019). By utilising different DTS measurement configurations, some studies have done distributed wind speed (heated cables) (Sayde et al.2015; Pfister et al.2017) or humidity (wetted cables) measurements (Euser et al.2014; Schilperoort et al.2018), whereas others have examined the radiation error of the cables (de Jong et al.2015; Sigmund et al.2017). However, thorough and critical analysis on the feasibility of DTS systems for measuring atmospheric scalar mixing has not been done since the pioneering study of Thomas et al. (2012). We complement the analysis made by Thomas et al. (2012) by comparing the DTS measurements against conventional in situ analysers within and above an aerodynamically rough forest canopy from the ground up to 120 m above the ground. Furthermore, we evaluate whether third-order statistics could be discerned from the DTS data to reveal important transport process information of the sweep ejection cycles created by the coherent structures. Deviations from Gaussian distribution are typically related to large coherent eddies or submeso air motions, whereas isotropic turbulence follows Gaussian distribution more closely, and hence, higher-order statistics are needed for studying these more organised flow patterns. We also assess the random uncertainties in the DTS measurements and their effect on the derived statistics. Finally, we demonstrate how a combination of vertical DTS measurements, together with 3D sonic anemometers and upward pointing Doppler lidar measurements, can be used to obtain continuous turbulence profiles, starting from the canopy sublayer all the way up to boundary layer top. Consequently, the DTS system bridges the scales between individual in situ measurement devices and remote sensing with lidar.

2 Materials and methods

2.1 Measurement site and instrumentation

The measurement campaign took place from 3 June to 8 July 2019 at the Hyytiälä SMEAR II (Station for Measuring Forest Ecosystem–Atmosphere Relations) station (Hari and Kulmala2005) located in central Finland (61.845 N, 24.289 E). The station is a member of several measurement networks, including ICOS (Integrated Carbon Observing System; Franz et al.2018) and ACTRIS (Aerosols, Clouds, and Trace gases Research Infrastructure), and thus, a wide range of atmospheric and biospheric measurements is conducted continuously. The site is dominated by a Scots pine (Pinus Sylvestris) stand, with an average tree height (h) of approximately 17 m, resulting in zero-plane displacement height (d) of 14 m. The forest canopy is between 10 and 17 m, and the all-sided leaf area index (LAI) was approximately 8 m2 m−2 in 2015. Below 10 m height, there is relatively open trunk space. The terrain is undulating, with the main slope (2) directed roughly in north–south direction. More details on the in-canopy turbulence and canopy structure are given in Launiainen et al. (2007) and on the terrain undulation in Alekseychik et al. (2013).

The measurements utilised in this study were conducted on the SMEAR II station's 125 m tall main mast. There are 3D ultrasonic anemometers installed at five heights on the mast (5.5, 24.6, 27, 68.3 and 125 m above ground), providing turbulent fluctuations of temperature and three wind speed components at 10 Hz sampling frequency. Different anemometers are present at different heights (USA-1, manufactured by METEK Meteorologische Messtechnik GmbH, Germany, at 5.5, 68.3 and 125 m; Solent Research model 1012R2 by Gill Instruments Limited, UK, at 24.6 m; HS-50 by Gill Instruments Limited, UK, at 27 m). The anemometer at 27 m height is part of the eddy covariance (EC) flux measurement setup for ICOS measurements (Rebmann et al.2018). Measurements at 5.5 m height started on 18 June 2019; all other 3D sonic anemometers were measuring continuously throughout the study period. The air temperature vertical profile below the ICOS EC measurement level was measured following ICOS procedures (Montagnani et al.2018). Radiation-shielded and ventilated platinum wire thermistors (PT100) were located at 3.3, 5.8, 8.8, 12.5, 16.8, 21.6 and 27 m heights above the ground. These measurements provide suitable data for comparison against the temperature data measured with DTS. Radiation components were measured with a four‐component net radiometer (model CNR4; Kipp & Zonen, Delft, the Netherlands) mounted on the mast at 67 m height.

Besides conventional in situ measurements on the mast, the atmospheric boundary layer was profiled with a scanning Doppler wind lidar (HALO Photonics StreamLine; Pearson et al.2009) located on the roof of a building approximately 400 m southwest from the mast. The wind lidar operates at 1.5 µm wavelength and, when pointing vertically, is configured to provide a profile of radial velocities approximately every 14 s with 30 m spatial resolution at Hyytiälä (Hirsikko et al.2014). The Doppler lidar provides data from 75 to 9585 m distance, and the clearest signal typically originates from the boundary layer where the aerosol loading is the highest.

DTS measurement setup

The DTS instrument (ULTIMA-S, 5 km variant; Silixa Ltd, Hertfordshire, UK) was housed in a wooden cabin located approximately 30 m away from the measurement mast. A total of two 50 L calibration baths were filled with water, equipped with reference thermometers (PT100; supplied with the DTS instrument and logged with the DTS instrument at 0.01 C precision), and the fibre-optic cable was guided twice through both baths. The baths were also equipped with aquarium pumps to ensure continuous mixing and prevent stratification of the water. Water temperature in the baths was controlled with two thermostats (RC 6 CS, Lauda Dr. R. Wobser GmbH & Co. KG, Domicile Lauda–Königshofen, Germany) and the bath temperatures were set at 5 and 30 C so that the temperatures bracketed the ambient air temperatures during the campaign. The baths were located next to the wooden cabin which housed the DTS instrument.

The measurements were conducted in double-ended mode, with signals from both directions stored separately (see Sect. 2.2 for processing). The path of the fibre-optic cable began by running through both calibration baths, then up and down the 125 m tall mast, and returning through the calibration baths again to finish. The cable was fastened to approximately 0.5 m long horizontal booms located every 5 to 10 m along the mast by squeezing the cable gently between two metallic plates covered with rubber slabs and soft electrical tape. Despite careful mounting of the cable, the cable holders caused disturbances in the data, and consequently, data points close to the holders were removed from subsequent analysis. Horizontal separation between the cable and the 3D sonic anemometers was 3.5 m at all levels. The vertical DTS measurements extended from 2 to 120 m above the ground. This setup enabled reference measurements at both the beginning and end of the cable from the calibration baths and provided a double measurement of the DTS vertical air temperature profile at 12.7 cm spatial and 1 Hz temporal resolution. Note, however, that, in double-ended mode, each direction along the cable is measured sequentially, so that the temperature data were available at 0.5 Hz resolution. A white, thin (outer diameter 0.9 mm) aramid-reinforced 50 µm multimode fibre-optic cable (AFL Telecommunications LLC, Duncan, SC 29334, US) was utilised in this study in order to minimise the impact of solar heating on the measurements (de Jong et al.2015; Sigmund et al.2017) and to maximise the high-frequency response of the setup.

2.2 Data processing

All high-frequency time series (3D sonic anemometers and DTS) were separated into the slowly varying mean and fluctuations around the mean as follows:

(1) T = T + T ,

where T is the measured high-frequency time series, overbar denotes averaging in time, and prime () fluctuations around the mean are calculated as a residual (T=T-T). Processing of 3D sonic anemometer data followed commonly applied procedures. The data were despiked in order to remove spurious outliers from the time series (Mauder et al.2013), and the wind vectors were aligned with the long-term mean wind field following the planar fit coordinate rotation method (Wilczak et al.2001). Turbulent fluxes and statistics (variance and skewness) were calculated with 30 min resolution throughout this study, unless otherwise noted. Note that, as H2O fluctuations were not measured at most heights, it was not possible to convert the sonic temperature to actual temperature (Schotanus et al.1983; Foken et al.2012). This caused a slight H2O variance and H2O flux-dependent bias in temperature variance and wT covariance, respectively. The stability of the atmospheric surface layer was evaluated based on the stability parameter as follows:

(2) ζ = z - d L ,

where z denotes height, d is displacement height, and L is Obukhov length. Negative values for ζ denote unstable stratification (buoyancy increases turbulence), positive values denote stable stratification (buoyancy diminishes turbulence) and zero denotes neutral stratification. A scaling parameter for temperature was used to describe the temperature variability in the atmospheric surface layer and was defined as follows:

(3) T * = w T u * ,

where u* is friction velocity. T* was calculated using data from the 3D sonic anemometer at 27 m height throughout the study.

DTS measurements were post-field calibrated using the measurements in the calibration baths and the following equation (e.g. van de Giesen et al.2012):

(4) T ( x , t ) = γ ln P s ( x , t ) P as ( x , t ) + C ( t ) + 0 x Δ α ( x ) d x ,

where T is the calibrated temperature as a function of position along the fibre (x) and time (t), γ is a DTS system-specific constant, Ps and Pas are the measured Stokes and anti-Stokes signals, C is a time-dependent calibration parameter, and the integral in the denominator is related to the differential attenuation of the anti-Stokes and Stokes signals along the cable. The double-ended configuration allowed the determination of the differential attenuation at each location along the cable separately, which was calculated following van de Giesen et al. (2012) separately for each 30 min averaging period. However, measurements from only one direction were utilised to calculate the temperature along the cable, i.e. measurements from the two directions were not averaged. After determining the differential attenuation, the calibration parameters γ and C were first determined by fitting Eq. (4) to the reference measurements made in the calibration baths for each time step separately producing time series for both parameters. Then, during the second processing stage, the data were calibrated again by fixing the value of γ to 476 K (mean of values obtained during the first processing stage) and letting C vary between time steps. We opted to process the data with this three-step procedure since (1) it resulted in less noisy calibration parameters, (2) theoretically γ should be constant and (3) a reduction in fit parameters is desirable given the number of calibration baths. The calibration parameter C showed only a slight variation in time (1.460 ± 0.003, mean ± SD), which resulted from a parabolic dependence on the indoor temperature of the cabin housing the DTS instrument. After quality filtering, 1513 (89 % of the whole period) 30 min periods of DTS data were available for further analysis.

The high-frequency response of the DTS system was evaluated by transforming the temperature time series to frequency domain with Fourier transform and comparing the power spectra against reference power spectra measured with the sonic anemometers. Following EC data-processing procedures (e.g. Ibrom et al.2007), the noise in the power spectra was subtracted from the signal prior to the comparison. By assuming that the DTS system as a whole behaves as a first-order response sensor, then, after normalisation, the ratio between the power spectra from DTS and the reference can be approximated by the following:

(5) T H ( f ) = 1 1 + 2 π f τ 2 ,

where f is frequency and τ is a first-order response time describing the high-frequency response of the system. An estimate for the attenuation of the DTS-derived temperature variance due to imperfect high-frequency response and lower sampling frequency can be derived by the following:

(6) AF = f 1 f 2 T H S TT , low d f f 1 f 3 S TT d f ,

where AF is the attenuation factor, STT is the power spectrum measured with the 3D sonic anemometer, STT,low is the power spectrum calculated from sonic anemometer data that were averaged in the time domain to match the temporal resolution of DTS prior to Fourier transform, and TH was calculated using a value for τ obtained experimentally with Eq. (5) (see Sect. 3.2). The integration limits f1, f2 and f3 are set by the length of the averaging period and the Nyquist frequency for DTS and 3D sonic anemometer data, respectively. By definition, AF is equal to one for a perfect sensor (τ=0 s) with high sampling frequency (i.e. f2=f3), and AF is equal to zero for a fully damped signal (τ→∞ s).

Determination of instrument noise

Instrument noise can significantly hinder the analysis of measurements if the signal-to-noise ratio is small. Here, we utilised the methodology developed by Lenschow et al. (2000) to extract turbulent statistics from noisy data. The method originally developed for lidar data has also been used for EC data (Mauder et al.2013; Rannik et al.2016) and is part of ICOS EC-data-processing routines (Nemitz et al.2018). The method has been shown to reliably estimate the influence of noise on second- and third-order statistics in turbulence measurements in prior studies (e.g. Lenschow et al.2000; Rannik et al.2016; Nakai et al.2020). The fluctuating part of the time series T(t) can be thought to consist of signal Ts(t) and noise ϵ(t) (T=Ts+ϵ). All the time series T, Ts and ϵ have zero means. The second-order autocovariance (M11) can then be written as follows:

(7) M 11 ( t l ) = T s + ϵ T s l + ϵ l ,

where overbar denotes time averaging, and subscript l denotes that the time series has been lagged with time lag tl (Lenschow et al.2000). Assuming ϵ and Ts are uncorrelated with each other, it can be found that, in the following:

(8) M 11 ( 0 ) = T s 2 + ϵ 2 ,

meaning that the time series variance M11(0) is a sum of signal variance Ts2 and noise variance ϵ2. Since noise is assumed to be uncorrelated, noise only contributes to the autocovariance at lag zero; hence, the signal variance can be estimated by extrapolating the autocovariance values to zero lag, and the noise variance can be estimated from the residual as follows:


Similarly, the effect of noise on third-order moments can be estimated using third-order autocovariance M21 as follows:


where it was further assumed that the product of Ts to any power and ϵ to an odd power are zero (Lenschow et al.2000). The time series skewness (Sk) was then estimated using the second- and third-order moments calculated with the equations above Sk=Ts3Ts23/2. The signal-to-noise ratio was used to estimate the magnitude of noise relative to the signal as follows:

(13) SNR = T s 2 ϵ 2 .
3 Results and discussion

3.1 DTS instrument noise determination

DTS measurements at high spatial and temporal resolution contain a significant noise contribution due to the small number of backscattered Raman photons; hence, duplicate or triplicate measurements are often conducted. In this study, the instrument noise was estimated using Lenschow et al. (2000, see Sect. 2.2.1) for each location along the fibre for each 30 min averaging period. Figure 1 shows the average dependence of the noise standard deviation, σnoise, with length along the fibre, LAF, together with the variability observed in the calibration baths. The noise increases exponentially with LAF, as expected, since the number of photons available for backscattering decreases along the fibre. After taking the mean of the measurements up and down the mast and then averaging over four consecutive points, the average noise decreased from 0.62 to 0.26 K, and was, in practise, independent of the position along the fibre and, thus, height. This decrease is similar to what would be expected for fully independent measurements (i.e. 0.62K/8=0.22 K). This averaging procedure was applied to all further analysis in this study and provided a temperature profile with approximately 0.5 m resolution from 2 up to 120 m above the ground. No temporal averaging was applied in order to preserve the temporal resolution of the signal. The noise level after averaging was below the typical temperature variability measured above the forest canopy during the campaign.

Figure 1Average dependence of noise (1σ), estimated with Eq. (10), for length along the fibre (LAF) and height above ground (top x axis). Circles show T variability in the calibration baths (mean ± SD). Dashed grey line shows the mean T standard deviation measured with 3D sonic anemometer at 27 m height during the measurement campaign. The vertical dotted black lines show the limits (bottom and top) of the vertical DTS measurements on the mast. The black line is a fit (σnoise=0.57K+0.01e8.6×10-3LAF, R2=0.99) to the noise estimates calculated from the cable going up (Tup) and down the mast (Tdown). The noise level decreased after averaging the up and down portions of the measurements(Tave; cyan line) and fell below the typical sonic T variability level further if four spatially consecutive measurements were averaged.


Based on the fit shown in Fig. 1, the noise at zero LAF was approximately 0.57 K, close to the instrument specifications. However, σnoise increased faster as a function of LAF than expected, likely due to the experimental setup which may have caused mechanical strain on the fibre (cable holders and wind load) and, hence, loss of signal. There were a few step losses (i.e. sudden decreases in Ps and Pas; see also Eq. 4) along the cable caused by the support structures, and these step losses added to the increases in the noise. The latter was linearly proportional to the magnitude of the step loss (Fig. 2). In addition to the mechanical stresses on the light-conducting glass core, the increase in noise as a function of LAF depends on the characteristics of the glass core in the cable and how it is coupled to the outer cable (sheath and protection). These are specific to the batch of the light-conducting glass core used during manufacturing of the cable. Nevertheless, the estimates obtained for instrument noise along the fibre can be used as a first-order estimate when designing future DTS measurement campaigns with inevitable signal artefacts due to the fibre-optic cable.

Figure 2Increase in noise at step losses (i.e. sudden decreases in Ps and Pas) along the cable versus step loss magnitude.


The value for σnoise at zero LAF is an intrinsic property of the measurement device and describes the noise level independently of fibre type or the measurement conditions. The variability in the instrument noise at zero LAF was evaluated using a similar fit to that shown in Fig. 1 for each 30 min measurement period. The obtained zero LAF noise estimates varied between 0.50 and 0.66 K (fifth and 95th percentiles of the obtained values) during the measurement campaign; the variability was not fully random as a systematic temporal component was also evident. Considering Eq. (4), the zero LAF noise in T is likely related to the variability in laser intensity or detector sensitivity. The DTS instrument contains an internal reference coil of fibre-optic cable used for a first instrument-based calibration using the variability in the Stokes signal (Ps,in) from this reference coil. We used this as a joint proxy for both the variability in emitted light and detector sensitivity, since the signal had not yet experienced any attenuation. The zero LAF noise displayed a linear dependence on this proxy (y=(0.10±0.001KdB-1)x-(0.44±0.001) K, R2=0.87, where y equals σnoise at LAF = 0 m, and x is the decrease in Ps,in (in decibels) from its maximum value). Interestingly, the gradient is similar to that obtained for σnoise increases at step losses (Fig. 2), indicative of a more general dependence between signal intensity and σnoise. It is known that laser intensity and detector responsivity can be sensitive to changes in temperature; the zero LAF noise showed a non-monotonic dependence on instrument internal temperature (calculated from data originating from the internal coil) and displayed additional changes when the instrument was cooling down or warming up (not shown). The dependence was not fully explained by the instrument internal temperature alone, and further studies are necessary to investigate this more deeply.

3.2 High-frequency response of the DTS system

The high-frequency response of the DTS system was evaluated by comparing the power spectra of DTS measurements against those of the co-located sonic anemometers as a reference (Fig. 3). The DTS power spectrum was dominated by white noise in the high-frequency part of the spectrum, whereas the reference followed the canonical inertial subrange slope. Following the commonly applied procedures for EC data processing (e.g. Ibrom et al.2007), the noise was removed from the DTS power spectrum by assuming that it was not correlated with the signal (see Sect. 2.2.1) and that the high-frequency response could be estimated from the ratio between the noise-removed power spectrum and the 3D sonic anemometer reference (see Sect. 2.2). Comparing DTS measurements to sonic anemometers at different heights (Table 1) provided slightly different values for the response time τ (Eq. 5), likely due to the uncertainty of the noise removal procedure. For reference, the values obtained are in the same range as those found for water vapour EC flux measurements in high relative humidity conditions (Ibrom et al.2007; Mammarella et al.2009; Runkle et al.2012) and are significantly higher than those for carbon dioxide (Mammarella et al.2009) or methane (Peltola et al.2014) EC flux measurements. When comparing the DTS power spectra to 3D sonic anemometer power spectra, it is important to recognise that the instruments were sampling at different rates (DTS with 0.5 Hz and 3D sonic anemometer with 10 Hz), and hence, the power spectra cannot fully be compared across all frequencies. The estimated response times and Eq. (5) should only be considered as a function for matching DTS and 3D sonic anemometer spectra and not as a transfer function describing the functioning of the DTS system.

Figure 3Ensemble-averaged, frequency-weighted and normalised temperature power spectrum estimated from the DTS data and co-located 3D sonic anemometer (reference – REF) at 25 m height above ground (1.5hc, where hc is canopy height). Only periods with high DTS SNR and moderate wind speed (between 1.7 and 3.3 m s−1) were selected for the ensemble. The DTS power spectrum is shown before and after the noise removal procedure, and the high-frequency attenuation is demonstrated by multiplying the reference power spectrum with a transfer function (Eq. 5).


Table 1Response times (τ) describing the high-frequency response of the DTS system derived via comparison with 3D sonic anemometers at different heights above ground (see Eq. 5), and attenuation factors (AFs) describing the attenuation of temperature variance due to the high-frequency loss of the signal. A total of 200 estimates for τ were estimated by bootstrap sampling the available data for each height, and the reported values are the medians (interquartile range) of the 200 estimates. AF was estimated for both height-dependent and constant τ for each 30 min time period during the measurement campaign with Eq. (6), and the reported values are the medians (interquartile range) of the AF time series.

Download Print Version | Download XLSX

Theoretically, the high-frequency response of the DTS system depends on (1) heat conduction across the quasi-laminar boundary layer surrounding the fibre cable, (2) heat conduction within the cable and (3) the spatial and temporal sampling resolution of the DTS instrument (Thomas et al.2012). The latter two can be considered independent of environmental conditions, whereas the first one depends on the depth of the boundary layer surrounding the cable, which is inversely dependent on wind speed. The DTS system high-frequency response time τ was calculated for different wind speed conditions at different heights, and no clear dependence on wind speed was found. This indicates that it is likely that the sampling resolution of the DTS instrument dominates the high-frequency response, and hence, τ can be considered constant for this particular fibre cable and DTS instrument combination.

Figure 4Dependence of AF on U/hc at different mast heights. The 30 min values (points), bin medians (filled circles) and the bin interquartile range (error bars) are shown. Data were screened based on flux stationarity (Foken and Wichura1996) and values for temperature variance and wT before plotting (|σT2|>0.02 K2 and |wT|>0.015 Km s−1).


The attenuation of the temperature variance due to the limited high-frequency response of the DTS system was estimated using the obtained values for τ and Eq. (6). The median values of the calculated AF values for different heights ranged between 0.74 and 0.86, indicating that the DTS temperature variances were typically underestimated by 20 %–30 % (Table 1). The signal attenuation is, in effect, a combination of attenuation (τ) and turbulent timescales (Horst1997); hence, AF was not constant while τ was. Turbulent timescales were estimated using roughness sublayer scaling, i.e. U/hc (Thomas and Foken2007a). The turbulent timescales are modulated by atmospheric stability; hence, AF should also depend on stability. When calculated with a fixed value for τ for each height, AF showed both an approximately linear decrease with U/hc and a dependence on stability (Fig. 4). However, the estimates of AF at different heights did not follow the same linear dependence, with clear differences between the dependence at 68 and 126 m (Fig. 4d and e) compared to lower heights, which may indicate that hc was not the physically meaningful length scale at these heights. The dependence of signal attenuation on turbulent timescales caused the median AF to increase weakly with height above the canopy, indicating smaller signal attenuation well above the canopy (Table 1). However, AF values from different heights were typically within 5 % of each other, indicating that the whole profile was attenuated in a similar fashion, and the reported values for AF serve as first-order estimates for the signal attenuation throughout the profile. Values of AF were higher below the canopy, indicating a reduction in the high-frequency attenuation of the signal; this was expected since organised (slow) motions have been previously found to dominate the scalar variability (and hence variance) within forests (Thomas and Foken2007a, b).

3.3 Second- and third-order statistics

Figure 5 and Table 2 show the agreement between temperature variances derived from DTS and the co-located 3D sonic anemometers. Note that statistics derived from 3D sonic anemometers were calculated from sonic temperature and, hence, slightly biased by H2O fluctuations (see Sect. 2.2). The DTS temperature variance was estimated using Eq. (9), where the noise variance was compensated for prior to the comparison. The bulk of the time series variance was related to fluctuations close to the peak of the power spectra, and the DTS system was able to resolve the variability at these frequencies (Fig. 3). Note that the gradient of the linear fits between the temperature variance estimates (Table 2) were close to the values obtained for the attenuation factors describing the high-frequency loss (Sect. 3.2), suggesting that the systematic mismatch between DTS and 3D sonic anemometer temperature variances was solely related to the limited high-frequency response of the DTS system.

Figure 5Comparison of DTS and 3D sonic anemometer (REF) 30 min values for temperature variance (a, d), wT covariance (b, e) and skewness (c, f) at two different heights on the mast. Upper panels (a–c) show below-canopy heights (0.3hc) and lower panels (d–f) show above-canopy heights (1.6hc). For skewness, only periods with SNR > 0.5 (Eq. 13) were used, otherwise all available data were included. Coefficients related to the linear fits are given in Table 2. Colour denotes the density of the point cloud. Note that no corrections for sensor separation effects (Horst and Lenschow2009) or signal attenuation (Sect. 3.2) were made prior to comparison.


Table 2Linear regression statistics between DTS and 3D sonic anemometers at different heights (y=ax+b, where y equals DTS, and x is a 3D sonic anemometer). The fits between data sets were performed using robust regression in order to minimise the effect of outliers on the regression coefficients. Standard errors for the coefficients are given in parentheses, and N shows the number of 30 min data points used in the fit. For skewness, only periods when the signal-to-noise ratio for DTS data was above 0.5 were utilised; for other statistics, all available data were used. Note that no corrections for sensor separation effects (Horst and Lenschow2009) or signal attenuation (Sect. 3.2) were made prior to comparison.

Download Print Version | Download XLSX

Following Thomas et al. (2012), the kinematic heat fluxes (i.e. wT) were calculated using the vertical wind speed from a 3D sonic anemometer and temperature either from the same 3D sonic anemometer or co-located DTS signal. Above the canopy, heat fluxes calculated from both systems showed relatively good agreement (Fig. 5e), indicating that the DTS system was able to resolve most of the eddies contributing to the vertical turbulent flux. However, within the forest canopy, the agreement was worse (Table 2 and Fig. 5b). There are at least two reasons for this. (1) There was a 3.5 m horizontal sensor separation between the 3D sonic anemometers and DTS fibre cable, which significantly contributed to dampening the high-frequency response of the joint anemometer–DTS flux calculation, especially close to the ground (Horst and Lenschow2009). (2) The forest floor is close, which suppresses the size of eddies dominating heat transfer, which is also influenced by the canopy elements, breaking the coherency of large eddies (spectral short circuit; Finnigan2000; Launiainen et al.2007), and DTS cannot capture the small-scale turbulence. This is supported by the finding that the in-canopy temperature variance was captured accurately, while the heat flux was not. The temperature fluctuations related to the large eddies sweeping into the in-canopy were captured accurately, yet their signal was decorrelated with respect to the vertical wind speed (i.e. vertical turbulent flux) by the canopy elements and large horizontal sensor separation. This is evident, for instance, in the comparison between, within and above forest canopy power spectra in Launiainen et al. (2007) measured at the same site. In general, scalars (e.g. temperature) behave very differently to the vectors within the canopy (Vickers and Thomas2014).

The third-order statistics, such as temperature skewness, were also compared. Again, the effect of noise was removed following the method by Lenschow et al. (2000), and the temperature skewness values were calculated using Eq. (11). Additionally, all data with SNR < 0.5 was discarded for the comparison presented in Fig. 5 and Table 2. The reason for this additional filtering can be seen in Fig. 6, where the discrepancy between DTS and 3D sonic anemometer temperature skewness increases rapidly as DTS SNR decreases below 0.5. The dependence of DTS SNR on height meant that better estimates of temperature skewness were available closer to the canopy than at higher elevations above it. Imposing a DTS SNR threshold of 0.5 still left 54 % of all DTS data measured during the campaign and 63 % of DTS data measured below 50 m height available for generating third-order statistics. Hence, we conclude that this DTS system can still be used after SNR filtering to monitor higher-order statistics and the non-Gaussian character of the flow, despite the higher noise floor of this instrument relative to the older 2 km variant used in Thomas et al. (2012).

Figure 6Difference in temperature skewness calculated from 3D sonic anemometer (REF) and DTS as a function of DTS signal-to-noise ratio (Eq. 13). Individual 30 min values (grey dots), bin medians (circles) and the interquartile range (error bars) are shown.


3.4 Profiles of first- through third-order statistics

A comparison of the profiles of turbulence statistics in different stability regimes are shown in Fig. 7. DTS provides continuous profiles, whereas conventional measurements provided point estimates; this is the unique strength of the spatially continuous DTS measurements. The mean potential temperature gradients from DTS exhibited some bias when compared against reference instrumentation (Fig. 7a), especially in mildly stable and strongly unstable cases. The fibre-optic cable was not radiation-shielded, and hence, it was exposed to shortwave and longwave radiation transfer with its surroundings. The 30 min mean temperature values were found to be biased by solar heating during daytime and longwave radiative cooling at night, in accordance with prior studies (de Jong et al.2015; Sigmund et al.2017). Due to the forest canopy, these radiation-induced biases differ below and above the canopy, introducing biases into the temperature gradient. In general, during night, DTS showed a greater decrease in temperature with height within the canopy than the reference profile, and the bias in the gradient depended on the amount of longwave radiative cooling (not shown). During daytime, DTS showed a weaker decrease in temperature with height due to canopy shading; this was particularly evident during morning and evening but less so during the middle of the day. The biases in across-canopy temperature differences were 0.07 K (median bias in temperature difference between 27 and 5.8 m heights) at night and 0.04 K at daytime, indicating that the DTS profile overestimated the across-canopy temperature gradient at night and underestimated the gradient during daytime. The unbiased median temperature gradients during these periods calculated from the reference measurements, at the heights mentioned above, were 0.34 and 0.29 K at night and daytime, respectively. Radiation shields around the cable could have presumably decreased these biases in gradients (Schilperoort et al.2018); however, the usage of screens would have invalidated the estimation of turbulent fluctuations from the DTS data since they disturb the turbulent airflow.

Figure 7Profiles of potential temperature gradient (θref= potential temperature at 27 m), normalised temperature variability (σT/|T*|) and temperature skewness derived from DTS (continuous lines) and reference (dots) measurements. The data were binned based on the stability parameter ζ estimated from EC measurements at 27 m height. For unstable conditions, σT/|T*| was also calculated based on Monin–Obukhov similarity theory (dashed lines) using σT/|T*|=2.3(1+9.5|ζ|)-1/3 (Cava et al.2008). Note that local scaling was not utilised, i.e. |T*| was calculated using measurements at a fixed height (27 m). Only periods with high SNR were used, and n denotes the number of 30 min periods in each stability bin. The right-hand axis displays height as a fraction of canopy height (hc), and the canopy height is denoted with a horizontal dashed line.


As shown above, DTS and 3D sonic temperature measurements showed good agreement in different mixing conditions within the canopy sublayer, slightly above the canopy and well above the forest. As expected, the temperature variance peaked at the canopy top due to (1) strong turbulence production originating from wind shear and (2) canopy heat source and/or sink related to solar heating (daytime) or radiative cooling (nighttime) of the canopy. The scaled temperature variability (σT/|T*|, where T* was calculated based on Eq. (3) and data from 27 m height) followed Monin–Obukhov (M–O) similarity scaling in the unstable and near-neutral regime above about 50 m height, indicating that surface layer scaling was valid between this height and the mast top. As expected, at lower heights (between 50 m height and ground level), the measured values for σT/|T*| departed from the M–O predictions. This is typical for the roughness sublayer, where the turbulence resembles more mixing layer than boundary layer turbulence (Raupach et al.1996; Finnigan2000; Finnigan et al.2009), meaning that the turbulent mixing is more efficient than in the surface layer and, hence, a smaller σT will yield the same turbulent flux. In other words, the correlation between w and T is higher in the roughness sublayer than in the surface layer (Patton et al.2010), and hence, M–O scaling is no longer valid. These results suggest that the height of the roughness sublayer at this site is about 50 m, which corresponds to approximately 3 times the height of the roughness elements (i.e. trees). This is in line with a previous study at this site (Rannik1998) and others, where estimates for the roughness sublayer height typically range between 2hc and 5hc (3hc being the most common estimate) (Garratt1980; Coppin et al.1986; Mölder et al.1999; Poggi et al.2004; Thomas et al.2006).

In near-neutral situations, the scaled temperature variability (σT/|T*|) exceeded the predictions made with M–O scaling since the heat fluxes (and hence also |T*|) decreased with ζ, yet the temperature variability (σT) did not decrease at the same rate. In other words, heat transfer efficiency approached zero at the neutral limit (e.g. Rannik1998). Under stable stratification, the scaled temperature variability showed consistent z dependence between the displacement height and approximately 30 m height, whereas above 30 m it was relatively independent of height, following z-less stratification, and any variability above this height was most likely related to flow processes other than interaction with the surface. These observations are in line with previous findings (e.g. Pahlow et al.2001). Note that, here, |T*| was calculated using measurements at a fixed height (27 m), i.e. local scaling was not utilised. The profiles shown in Fig. 7b for unstable and neutral periods could be expected to show a similar pattern, even if local scaling were used since turbulent fluxes can be conjectured to be constant with height in the bottom part of the ABL. However, for stable situations, the σT/|T*| calculated based on local scaling might depart from what is shown here since |T*| varies with z.

Temperature skewness was almost independent of height (approximately 0.5–0.6) in unstable conditions above 30–50 m (2hc3hc), which corresponded to the roughness sublayer height found above. Skewness increased as the instability increased. Skewness decreased when moving to the canopy sublayer, but stayed positive throughout the profile, indicating non-Gaussian flow. Non-zero skewness is often connected to an imbalance between the transport of the scalar in question by ejective and sweeping air motions (Katul et al.2018). Hence, based on these results, the ejections and sweeps were not in balance, and the ejective motions dominate heat transport in the unstable conditions throughout the air column. Non-zero skewness of scalar time series has also been linked to the influence of ABL-scale large eddies entraining air from the free troposphere and transporting it in downdraughts close to the ground (Mahrt1991; Couvreux et al.2007; van de Boer et al.2014). In contrast to the profiles in unstable conditions, in stable conditions above the forest canopy the flow was Gaussian (zero skewness), whereas skewness was positive below the canopy due to the sweeping motions penetrating through the canopy and bringing pulses of warm air from aloft into the cold below-canopy air.

3.5 Examples of organised patterns observed with the DTS system across vertical coupling regimes

The spatially continuous measurements of the DTS system enabled the detection and analysis of spatial temperature patterns in the flow. During well-developed turbulence, these patterns are interlinked with large-scale, organised turbulent air motions (see, e.g., Fig. 3 in Gao et al.1989). Figure 8 shows a measurement example during a daytime unstable regime. Large coherent eddies dominated the flow (Fig. 8b), and their signatures on vertical temperature profiles were captured with the DTS measurements. Positive temperature perturbations were correlated with upward vertical air motions and negative with downward motions, indicating upward-directed sensible heat flux due to the unstable stratification. The entire 120 m measurement domain was coupled due to large coherent eddies effectively mixing the air throughout the vertical column. The temporal extent and amplitude of the temperature fluctuations changed with height, corresponding to changes in the dominant turbulence timescale and flux magnitude with height.

Figure 8Example of simultaneous lidar, DTS and 3D sonic anemometer measurements during daytime convective conditions. (a) Lidar vertical wind velocity component (gap is due to non-vertical operation during scanning sequence). (b) Temperature profiles measured with the DTS system (colour) and wind speed fluctuations (arrows) measured with 3D sonic anemometers. For illustration purposes, the wind speed fluctuations were averaged over 10 s, and the vertical wind component was multiplied by 10 prior to plotting. Length of the arrows denote the magnitude of the wind speed fluctuations, and direction is determined by the sign of vertical (w) and horizontal (u) wind speed fluctuations (up – w>0; down – w<0; left – u<0; right – u>0). (c–g) Temperature perturbations (T-T) measured with DTS and sonic anemometers at different heights on the mast. Local time (universal coordinated time – UTC+2), the corresponding stability parameter (ζ) and mean wind speed (U) measured at 27 m height are given in the figure title. Canopy height was approximately 17 m and is highlighted with a white dotted line in panel (b). For illustration purposes, the data gaps from cable holder locations were filled with linear interpolation prior to plotting.


Figure 9Mean potential temperature gradient (a) and standard deviation (b) profiles for the two example periods shown in Figs. 8 and 10. Dots in panel (b) show the σT values estimated with the 3D sonic anemometers at different heights.


By relying on Taylor's frozen turbulence hypothesis, two-dimensional spatial details of the large temperature patterns can be delineated from the vertical DTS measurements. The patterns with low temperatures (related to downward air motions called sweeps) occasionally reached the forest floor but did not always penetrate the whole canopy layer (between 10 and 17 m) since they were likely disrupted by the canopy passage. The patterns with high temperatures, representing upward motions called ejections, typically originated from the forest floor, traversed the forest canopy and continued moving upwards (compare Fig. 8a and b). Ramp cliff patterns connected to the ejection sweep cycle were evident in the time series collected just above the canopy (Fig. 8e and f) but not so clearly at higher levels or below the canopy (Fig. 8c, d and g). These patterns have been previously suggested to be connected to the wind shear and corresponding inflection point instability close to the canopy top (Raupach et al.1996; Finnigan2000; Cava et al.2004; Thomas and Foken2007b; Göckede et al.2007). For this reason, the signatures of ejections would not be expected to reach very high above the canopy. The mean potential temperature profile (Fig. 9a) showed unstable stratification, except below the canopy and in the upper parts of the profile where the stratification was close to neutral.

The patterns observed with the DTS system and the 3D sonic anemometer at the mast top agreed qualitatively with the vertical-pointing wind lidar (Fig. 8), yet quantitative analysis on the agreement was hindered by the fact that the lidar instrument was located approximately 400 m upwind from the measurement mast. While the temporal lag caused by this horizontal displacement was taken into account by shifting the lidar time series in time to maximise the correlation with the mast-top 3D sonic anemometer, the eddies may have already been deformed during advection from the lidar measurement location to the measurement mast, in addition to changes in wind direction. Hence, a direct comparison between the lidar in its current position and the mast measurements is not possible. Nevertheless, these observations demonstrate the capabilities of joint measurements with lidar and DTS in capturing a continuous vertical profile of turbulence from forest floor up to boundary layer top.

Figure 10Same as Fig. 8 but for a strongly stable situation at night.


Figure 11Evolution of potential temperature gradient (dθ/dz(a) for the nighttime period shown in Fig. 10 and concurrent wind speed time series from two heights (b). The black dots in panel (a) highlight the location of the maximum gradient at each time step. The gradient was calculated using spline fit for each time step and estimating the gradient for each height from the first derivative of the fit. White dotted line shows the canopy height.


In contrast to the daytime example, during the nighttime example the air temperature patterns suggested a decoupling caused by the strong thermal stratification of the air (Figs. 9a and 10b), resulting from radiative cooling of the forest (net radiative loss of 51 W m−2) and low mechanical production of turbulence (i.e. weak wind shear). A total of two strong temperature inversion layers were evident, namely one at the canopy top between 12 and 25 m heights and one at 80 m, which descended down to 60 m height during the latter part of the period (Fig. 11a). Vertical movement of these inversion layers resulted in a bimodal vertical profile for σT during this period (Fig. 9b). These inversion layers decoupled the air column into the following three layers: below-canopy airspace (below 10 m height), canopy layer with organised motions (between 10 and 80 m) and a residual layer with slowly varying non-turbulent motions (above 80 m). The maximum potential temperature gradient resided at the canopy height due to radiative cooling of the canopy (Figs. 9a and 11a). Inverted ramp patterns were observed in the temperature time series at the beginning of the example period (Fig. 10), and a canopy wave was observed in the middle of the period (Fig. 11a). By following the evolution of the maximum temperature gradient (Fig. 11a), the amplitude of the canopy wave was estimated to be about 8 m. The switch from inverted ramps to canopy waves was likely due to a decrease in wind shear and, hence, a decrease in turbulence production (Fig. 11b).

These two examples, from contrasting mixing regimes, illustrate the unique observational capabilities of the DTS to analyse the spatial structure of the organised patterns, which can only be guessed at individual heights from classical sonic anemometry (compare, e.g., Fig. 8b and e). In principle, similar results could have been obtained with an array of several in situ instruments (Gao et al.1989; Lee et al.1997; Poulos et al.2002; Horst et al.2004; Mahrt and Vickers2005; Mahrt et al.2009; Patton et al.2010; Bou-Zeid et al.2010; Feigenwinter et al.2010; Thomas2011; Serafimovich et al.2011; Mahrt et al.2014), yet such measurement setups do not provide spatially continuous measurements with the same spatial resolution as the DTS system and are labour-intensive and require careful cross-instrument calibration. Furthermore, continuous profiles are required for following the evolution of, for example, inversion layers (see Fig. 11) and, hence, not easily detectable with arrays of in situ instruments. Note, however, that these two examples relied only on temperature observations; for similar observations of wind vectors, actively heated fibre optics would be needed (Sayde et al.2015).

4 Conclusions and outlook

This study demonstrated the unique observational potential of DTS measurements for capturing the salient features of atmospheric flow within and above-forest canopies across contrasting flow regimes. Despite the fact that these results were obtained with a different DTS machine compared to Thomas et al. (2012) and that different cable suspension techniques were used, these results are in line with their findings on DTS capturing second-order moments of air temperature variability and complement them. In spite of the higher noise floor in the measurements and the limited high-frequency response compared to sonic anemometers, we found the DTS technique to accurately capture the second- and third-order moments, which allows for detection of the spatial structure of coherent air motions dominating the temperature variability at the forest–air interface. In accordance with Thomas et al. (2012), the measurements performed best when either large variability (high heat fluxes) or strong stratification (low mixing and clear skies) was present. Cross-canopy temperature gradients were found to be biased due to radiation-induced errors, especially during morning and evening. These biases could be reduced with radiation shields around the DTS cable; however, such shields would severely disrupt the airflow and, hence, inhibit the estimation of turbulence from DTS measurements. Other approaches have been proposed as well (de Jong et al.2015; Sigmund et al.2017). Despite these shortcomings, DTS measurements can provide the missing spatial details of atmospheric mixing close to the surface which cannot be acquired with conventional in situ or remote sensing methods.

A combination of DTS and conventional in situ instrumentation has been shown to be advantageous when studying low-mixing nocturnal airflows (Thomas et al.2012; Zeeman et al.2015; Pfister et al.2017; Mahrt et al.2019; Pfister et al.2019) and morning transitions (Higgins et al.2018) over smooth surfaces. In the current study, it was shown that the DTS measurements can be used in a wide range of mixing conditions throughout the day, both above and within rough forest canopies. When affixed to a tall mast, DTS measurements can capture coherent motions up to heights of 7hc (Fig. 8b) or more. Hence, observations with different DTS measurement configurations should prove useful in locations where the flow exhibits persistent spatial patterns which are not transported with the mean flow and, thus, cannot be studied with standard point-based in situ instrumentation. These locations include, but are not limited to, flows within and just above forest canopies and street canyons and flows over edges and other spatial discontinuities in the surface characteristics.

Together with aiding the analysis of different flow modes, DTS measurements show potential for investigating the spatial details of atmosphere–ecosystem interactions and gas transfer between these two domains. For instance, the vertical profile of temperature fluctuations can help to separate the canopy and forest floor contributions to above-canopy scalar fluxes (Thomas et al.2008; Klosterhalfen et al.2019), since the measurements allow the tracking of the temperature signal related to coherent ejection and sweep air motions through the forest canopy (Fig. 8b). The measurements can also aid in trying to decipher the persistent problem related to advection at the EC flux measurement sites when trying to close the mass balance of, for example, CO2 for a given target ecosystem (Aubinet et al.2010), since part of this issue is likely to be related to the unresolved horizontal variability of turbulent mixing within the forest (Feigenwinter et al.2010). This issue could be explored with DTS measurements utilising horizontally and vertically distributed fibre-optic cables.

Ultimately, the DTS measurements enable studies to go beyond the typical time series from isolated point measurements, yet, at the same time, provide measurements at the same temporal and spatial scale as conventional in situ instruments. Many topical open questions in the field of ecosystem–atmosphere interactions require spatio-temporal information on the processes involved, and DTS has been shown here to be a suitable tool for gathering such information.

Data availability

Data have been uploaded to the Zenodo open data repository (, Peltola et al.2020).

Author contributions

OP designed the experiment and did the data processing and analysis. KL and CKT supported the DTS data analysis work. IM provided technical support in the field. EO'C maintained the Doppler lidar instrument and provided the data. OP wrote the first version of the paper, and KL, IM, EO'C, CKT and TV provided input.

Competing interests

The authors declare that they have no conflict of interest.


Technical staff at the Hyytiälä research station are acknowledged for their help during the measurement campaign. Olli Peltola has been supported by the postdoctoral researcher project (grant no. 315424) funded by the Academy of Finland. Karl Lapo and Christoph K. Thomas received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant no. 724629; project DarkMix). ICOS-Finland and Atmospheric Mathematics (AtMath) projects by the University of Helsinki are also acknowledged.

Financial support

This research has been supported by the Academy of Finland, Biotieteiden ja Ympäristön Tutkimuksen Toimikunta (grant no. 315424) and the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (DarkMix; grant no. 724629).

Review statement

This paper was edited by Christof Ammann and reviewed by three anonymous referees.


Alekseychik, P., Mammarella, I., Launiainen, S., Rannik, U., and Vesala, T.: Evolution of the nocturnal decoupled layer in a pine forest canopy, Agr. Forest Meteorol., 174–175, 15–27,, 2013. a

Aubinet, M., Feigenwinter, C., Heinesch, B., Bernhofer, C., Canepa, E., Lindroth, A., Montagnani, L., Rebmann, C., Sedlak, P., and Van Gorsel, E.: Direct advection measurements do not help to solve the night-time CO2 closure problem: Evidence from three different forests, Agr. Forest Meteorol., 150, 655–664,, 2010. a

Baldocchi, D.: Measuring fluxes of trace gases and energy between ecosystems and the atmosphere – the state and future of the eddy covariance method, Glob. Change Biol., 20, 3600–3609,, 2014. a

Barlow, J. F.: Progress in observing and modelling the urban boundary layer, Urban Climate, 10, 216–240,, 2014. a

Bohrer, G., Katul, G. G., Walko, R. L., and Avissar, R.: Exploring the Effects of Microscale Structural Heterogeneity of Forest Canopies Using Large-Eddy Simulations, Bound.-Lay. Meteorol., 132, 351–382,, 2009. a

Bou-Zeid, E., Higgins, C., Huwald, H., Meneveau, C., and Parlange, M. B.: Field study of the dynamics and modelling of subgrid-scale turbulence in a stable atmospheric surface layer over a glacier, J. Fluid Mech., 665, 480–515,, 2010. a

Cava, D., Giostra, U., Siqueira, M., and Katul, G.: Organised Motion and Radiative Perturbations in the Nocturnal Canopy Sublayer above an Even-Aged Pine Forest, Bound.-Lay. Meteorol., 112, 129–157,, 2004. a

Cava, D., Katul, G. G., Sempreviva, A. M., Giostra, U., and Scrimieri, A.: On the Anomalous Behaviour of Scalar Flux–Variance Similarity Functions Within the Canopy Sub-layer of a Dense Alpine Forest, Bound.-Lay. Meteorol., 128, 33,, 2008. a

Cheng, Y., Sayde, C., Li, Q., Basara, J., Selker, J., Tanner, E., and Gentine, P.: Failure of Taylor's hypothesis in the atmospheric surface layer and its correction for eddy-covariance measurements, Geophys. Res. Lett., 44, 4287–4295,, 2017. a

Coppin, P. A., Raupach, M. R., and Legg, B. J.: Experiments on scalar dispersion within a model plant canopy part II: An elevated plane source, Bound.-Lay. Meteorol., 35, 167–191,, 1986. a

Couvreux, F., Guichard, F., Masson, V., and Redelsperger, J.-L.: Negative water vapour skewness and dry tongues in the convective boundary layer: observations and large-eddy simulation budget analysis, Bound.-Lay. Meteorol., 123, 269–294,, 2007. a

Dakin, J. P., Pratt, D. J., Bibby, G. W., and Ross, J. N.: Distributed optical fibre Raman temperature sensor using a semiconductor light source and detector, Electron. Lett., 21, 569–570,, 1985. a

de Jong, S. A. P., Slingerland, J. D., and van de Giesen, N. C.: Fiber optic distributed temperature sensing for the determination of air temperature, Atmos. Meas. Tech., 8, 335–339,, 2015. a, b, c, d, e

Egerer, U., Gottschalk, M., Siebert, H., Ehrlich, A., and Wendisch, M.: The new BELUGA setup for collocated turbulence and radiation measurements using a tethered balloon: first applications in the cloudy Arctic boundary layer, Atmos. Meas. Tech., 12, 4019–4038,, 2019. a

Euser, T., Luxemburg, W. M. J., Everson, C. S., Mengistu, M. G., Clulow, A. D., and Bastiaanssen, W. G. M.: A new method to measure Bowen ratios using high-resolution vertical dry and wet bulb temperature profiles, Hydrol. Earth Syst. Sci., 18, 2021–2032,, 2014. a, b

Feigenwinter, C., Mölder, M., Lindroth, A., and Aubinet, M.: Spatiotemporal evolution of CO2 concentration, temperature, and wind field during stable nights at the Norunda forest site, Agr. Forest Meteorol., 150, 692–701,, 2010. a, b

Finnigan, J.: Turbulence in Plant Canopies, Annu. Rev. Fluid Mech., 32, 519–571,, 2000. a, b, c, d

Finnigan, J., Shaw, R. H., and Patton, E. G.: Turbulence structure above a vegetation canopy, J. Fluid Mech., 637, 387–424,, 2009. a, b

Foken, T. and Wichura, B.: Tools for quality assessment of surface-based flux measurements, Agr. Forest Meteorol., 78, 83–105, 1996. a

Foken, T., Leuning, R., Oncley, S. R., Mauder, M., and Aubinet, M.: Corrections and Data Quality Control BT - Eddy Covariance: A Practical Guide to Measurement and Data Analysis, Springer Netherlands, Dordrecht, 85–131,, 2012. a

Franz, D., Acosta, M., Altimir, N., Arriga, N., Arrouays, D., Aubinet, M., Aurela, M., Ayres, E., López-Ballesteros, A., Barbaste, M., Berveiller, D., Biraud, S., Boukir, H., Brown, T., Brömmer, C., Buchmann, N., Burba, G., Carrara, A., Cescatti, A., Ceschia, E., Clement, R., Cremonese, E., Crill, P., Darenova, E., Dengel, S., D'Odorico, P., Filippa, G., Fleck, S., Fratini, G., Fuß, R., Gielen, B., Gogo, S., Grace, J., Graf, A., Grelle, A., Gross, P., Grönwald, T., Haapanala, S., Hehn, M., Heinesch, B., Heiskanen, J., Herbst, M., Herschlein, C., Hörtnagl, L., Hufkens, K., Ibrom, A., Jolivet, C., Joly, L., Jones, M., Kiese, R., Klemedtsson, L., Kljun, N., Klumpp, K., Kolari, P., Kolle, O., Kowalski, A., Kutsch, W., Laurila, T., De Ligne, A., Linder, S., Lindroth, A., Lohila, A., Longdoz, B., Mammarella, I., Manise, T., Jiménez, S., Matteucci, G., Mauder, M., Meier, P., Merbold, L., Mereu, S., Metzger, S., Migliavacca, M., Mölder, M., Montagnani, L., Moureaux, C., Nelson, D., Nemitz, E., Nicolini, G., Nilsson, M., De Beeck, M., Osborne, B., Löfvenius, M., Pavelka, M., Peichl, M., Peltola, O., Pihlatie, M., Pitacco, A., Pokorný, R., Pumpanen, J., Ratié, C., Rebmann, C., Roland, M., Sabbatini, S., Saby, N., Saunders, M., Schmid, H., Schrumpf, M., Sedlák, P., Ortiz, P., Siebicke, L., Šigut, L., Silvennoinen, H., Simioni, G., Skiba, U., Sonnentag, O., Soudani, K., Soulé, P., Steinbrecher, R., Tallec, T., Thimonier, A., Tuittila, E.-S., Tuovinen, J.-P., Vestin, P., Vincent, G., Vincke, C., Vitale, D., Waldner, P., Weslien, P., Wingate, L., Wohlfahrt, G., Zahniser, M., and Vesala, T.: Towards long-Term standardised carbon and greenhouse gas observations for monitoring Europe's terrestrial ecosystems: A review, Int. Agrophys., 32, 439–455, 2018. a

Frehlich, R., Meillier, Y., and Jensen, M. L.: Measurements of Boundary Layer Profiles with In Situ Sensors and Doppler Lidar, J. Atmos. Ocean. Tech., 25, 1328–1340,, 2008. a

Gao, W., Shaw, R. H., and Paw U, K. T.: Observation of organized structure in turbulent flow within and above a forest canopy, Bound.-Lay. Meteorol., 47, 349–377,, 1989. a, b

Garratt, J. R.: Surface influence upon vertical profiles in the atmospheric near-surface layer, Q. J. Roy. Meteorol. Soc., 106, 803–819,, 1980. a

Göckede, M., Thomas, C., Markkanen, T., Mauder, M., Ruppert, J., and Foken, T.: Sensitivity of Lagrangian Stochastic footprints to turbulence statistics, Tellus B, 59, 577–586,, 2007. a

Hari, P. and Kulmala, M.: Station for measuring ecosystem-atmosphere relations (SMEAR II), Boreal Environ. Res., 10, 315–322, 2005. a

Higgins, C. W., Froidevaux, M., Simeonov, V., Vercauteren, N., Barry, C., and Parlange, M. B.: The Effect of Scale on the Applicability of Taylor’s Frozen Turbulence Hypothesis in the Atmospheric Boundary Layer, Bound.-Lay. Meteorol., 143, 379–391,, 2012. a

Higgins, C. W., Katul, G. G., Froidevaux, M., Simeonov, V., and Parlange, M. B.: Are atmospheric surface layer flows ergodic?, Geophys. Res. Lett., 40, 3342–3346,, 2013. a

Higgins, C. W., Wing, M. G., Kelley, J., Sayde, C., Burnett, J., and Holmes, H. A.: A high resolution measurement of the morning ABL transition using distributed temperature sensing and an unmanned aircraft system, Environ. Fluid Mech., 18, 683–693,, 2018. a, b, c, d

Higgins, C. W., Drake, S. A., Kelley, J., Oldroyd, H. J., Jensen, D. D., and Wharton, S.: Ensemble-Averaging Resolves Rapid Atmospheric Response to the 2017 Total Solar Eclipse , Front. Earth. Sci., 7, 198,, 2019. a, b

Hirsikko, A., O'Connor, E. J., Komppula, M., Korhonen, K., Pfüller, A., Giannakaki, E., Wood, C. R., Bauer-Pfundstein, M., Poikonen, A., Karppinen, T., Lonka, H., Kurri, M., Heinonen, J., Moisseev, D., Asmi, E., Aaltonen, V., Nordbo, A., Rodriguez, E., Lihavainen, H., Laaksonen, A., Lehtinen, K. E. J., Laurila, T., Petäjä, T., Kulmala, M., and Viisanen, Y.: Observing wind, aerosol particles, cloud and precipitation: Finland's new ground-based remote-sensing network, Atmos. Meas. Tech., 7, 1351–1375,, 2014. a

Horst, T. W.: A simple formula for attenuation of eddy fluxes measured with first-order-response scalar sensors, Bound.-Lay. Meteorol., 82, 219–233, 1997. a

Horst, T. W. and Lenschow, D. H.: Attenuation of Scalar Fluxes Measured with Spatially-displaced Sensors, Bound.-Lay. Meteorol., 130, 275–300,, 2009. a, b, c

Horst, T. W., Kleissl, J., Lenschow, D. H., Meneveau, C., Moeng, C.-H., Parlange, M. B., Sullivan, P. P., and Weil, J. C.: HATS: Field Observations to Obtain Spatially Filtered Turbulence Fields from Crosswind Arrays of Sonic Anemometers in the Atmospheric Surface Layer, J. Atmos. Sci., 61, 1566–1581,<1566:HFOTOS>2.0.CO;2, 2004. a

Ibrom, A., Dellwik, E., Flyvbjerg, H., Jensen, N. O., and Pilegaard, K.: Strong low-pass filtering effects on water vapour flux measurements with closed-path eddy correlation systems, Agr. Forest Meteorol., 147, 140–156,, 2007. a, b, c

Izett, J. G., Schilperoort, B., Coenders-Gerrits, M., Baas, P., Bosveld, F. C., and van de Wiel, B. J. H.: Missed Fog?, Bound.-Lay. Meteorol.,, 2019. a, b

Kang, Y., Belusic, D., and Smith-Miles, K.: Classes of structures in the stable atmospheric boundary layer, Q. J. Roy. Meteorol. Soc., 141, 2057–2069,, 2015. a

Katul, G., Peltola, O., Grönholm, T., Launiainen, S., Mammarella, I., and Vesala, T.: Ejective and Sweeping Motions Above a Peatland and Their Role in Relaxed-Eddy-Accumulation Measurements and Turbulent Transport Modelling, Bound.-Lay. Meteorol., 169, 163–184,, 2018. a

Keller, C. A., Huwald, H., Vollmer, M. K., Wenger, A., Hill, M., Parlange, M. B., and Reimann, S.: Fiber optic distributed temperature sensing for the determination of the nocturnal atmospheric boundary layer height, Atmos. Meas. Tech., 4, 143–149,, 2011. a, b

Klosterhalfen, A., Graf, A., Brüggemann, N., Drüe, C., Esser, O., González-Dugo, M. P., Heinemann, G., Jacobs, C. M. J., Mauder, M., Moene, A. F., Ney, P., Pütz, T., Rebmann, C., Ramos Rodríguez, M., Scanlon, T. M., Schmidt, M., Steinbrecher, R., Thomas, C. K., Valler, V., Zeeman, M. J., and Vereecken, H.: Source partitioning of H2O and CO2 fluxes based on high-frequency eddy covariance data: a comparison between study sites, Biogeosciences, 16, 1111–1132,, 2019. a

Launiainen, S., Vesala, T., Mölder, M., Mammarella, I., Smolander, S., Rannik, U., Kolari, P., Hari, P., Lindroth, A., and Katul, G.: Vertical variability and effect of stability on turbulence characteristics down to the floor of a pine forest, Tellus B, 59, 919–936,, 2007. a, b, c

Lee, X., Neumann, H. H., Hartog, G., Mickle, R. E., Fuentes, J. D., Black, T. A., Yang, P. C., and Blanken, P. D.: Observation of gravity waves in a boreal forest, Bound.-Lay. Meteorol., 84, 383–398,, 1997. a

Lenschow, D. H., Wulfmeyer, V., and Senff, C.: Measuring second- through fourth-order moments in noisy data, J. Atmos. Ocean. Tech., 17, 1330–1347,<1330:MSTFOM>2.0.CO;2, 2000. a, b, c, d, e, f

Mahrt, L.: Boundary-layer moisture regimes, Q. J. Roy. Meteorol. Soc., 117, 151–176,, 1991. a

Mahrt, L.: Stably Stratified Atmospheric Boundary Layers, Annu. Rev. Fluid Mech., 46, 23–45,, 2014. a

Mahrt, L. and Thomas, C. K.: Surface Stress with Non-stationary Weak Winds and Stable Stratification, Bound.-Lay. Meteorol., 159, 3–21,, 2016. a

Mahrt, L. and Vickers, D.: Boundary-Layer Adjustment Over Small-Scale Changes of Surface Heat Flux, Bound.-Lay. Meteorol., 116, 313–330,, 2005. a

Mahrt, L., Thomas, C. K., and Prueger, J. H.: Space–time structure of mesoscale motions in the stable boundary layer, Q. J. Roy. Meteorol. Soc., 135, 67–75,, 2009. a, b

Mahrt, L., Sun, J., Oncley, S. P., and Horst, T. W.: Transient Cold Air Drainage down a Shallow Valley, J. Atmos. Sci., 71, 2534–2544,, 2014. a

Mahrt, L., Pfister, L., and Thomas, C. K.: Small-Scale Variability in the Nocturnal Boundary Layer, Bound.-Lay. Meteorol., 174, 81–98,, 2019. a, b, c

Mammarella, I., Launiainen, S., Gronholm, T., Keronen, P., Pumpanen, J., Rannik, U., and Vesala, T.: Relative Humidity Effect on the High-Frequency Attenuation of Water Vapor Flux Measured by a Closed-Path Eddy Covariance System, J. Atmos. Ocean. Tech., 26, 1856–1866,, 2009. a, b

Mauder, M., Cuntz, M., Druee, C., Graf, A., Rebmann, C., Schmid, H. P., Schmidt, M., and Steinbrecher, R.: A strategy for quality and uncertainty assessment of long-term eddy-covariance measurements, Agr. Forest Meteorol., 169, 122–135,, 2013. a, b

Mölder, M., Grelle, A., Lindroth, A., and Halldin, S.: Flux-profile relationships over a boreal forest – roughness sublayer corrections, Agr. Forest Meteorol., 98–99, 645–658,, 1999. a

Montagnani, L., Grünwald, T., Kowalski, A., Mammarella, I., Merbold, L., Metzger, S., Sedlák, P., and Siebicke, L.: Estimating the storage term in eddy covariance measurements: the ICOS methodology, Int. Agrophys., 32, 551–567,, 2018. a

Nakai, T., Hiyama, T., Petrov, R. E., Kotani, A., Ohta, T., and Maximov, T. C.: Application of an open-path eddy covariance methane flux measurement system to a larch forest in eastern Siberia, Agr. Forest Meteorol., 282–283, 107860,, 2020. a

Nemitz, E., Mammarella, I., Ibrom, A., Aurela, M., Burba, G., Dengel, S., Gielen, B., Grelle, A., Heinesch, B., Herbst, M., Hörtnagl, L., Klemedtsson, L., Lindroth, A., Lohila, A., McDermitt, D., Meier, P., Merbold, L., Nelson, D., Nicolini, G., Nilsson, M., Peltola, O., Rinne, J., and Zahniser, M.: Standardisation of eddy-covariance flux measurements of methane and nitrous oxide, Int. Agrophys., 32, 517–549,, 2018. a

Newsom, R., Calhoun, R., Ligon, D., and Allwine, J.: Linearly Organized Turbulence Structures Observed Over a Suburban Area by Dual-Doppler Lidar, Bound.-Lay. Meteorol., 127, 111–130,, 2008. a

Pahlow, M., Parlange, M. B., and Porté-Agel, F.: On Monin–Obukhov Similarity In The Stable Atmospheric Boundary Layer, Bound.-Lay. Meteorol., 99, 225–248,, 2001. a

Patton, E. G., Horst, T. W., Sullivan, P. P., Lenschow, D. H., Oncley, S. P., Brown, W. O. J., Burns, S. P., Guenther, A. B., Held, A., Karl, T., Mayor, S. D., Rizzo, L. V., Spuler, S. M., Sun, J., Turnipseed, A. A., Allwine, E. J., Edburg, S. L., Lamb, B. K., Avissar, R., Calhoun, R. J., Kleissl, J., Massman, W. J., Paw U, K. T., and Weil, J. C.: The Canopy Horizontal Array Turbulence Study, B. Am. Meteorol. Soc., 92, 593–611,, 2010. a, b

Pearson, G., Davies, F., and Collier, C.: An Analysis of the Performance of the UFAM Pulsed Doppler Lidar for Observing the Boundary Layer, J. Atmos. Ocean. Tech., 26, 240–250,, 2009. a

Peltola, O., Hensen, A., Helfter, C., Belelli Marchesini, L., Bosveld, F. C., van den Bulk, W. C. M., Elbers, J. A., Haapanala, S., Holst, J., Laurila, T., Lindroth, A., Nemitz, E., Röckmann, T., Vermeulen, A. T., and Mammarella, I.: Evaluating the performance of commonly used gas analysers for methane eddy covariance flux measurements: the InGOS inter-comparison field experiment, Biogeosciences, 11, 3163–3186,, 2014. a

Peltola, O., Lapo, K., Martinkauppi, I., O'Connor, E., Thomas, C. K., and Vesala, T.: Dataset for “Suitability of fiber-optic distributed temperature sensing to reveal mixing processes and higher-order moments at the forest-air interface” [Data set], Zenodo,, 2020. a

Pfister, L., Sigmund, A., Olesch, J., and Thomas, C. K.: Nocturnal Near-Surface Temperature, but not Flow Dynamics, can be Predicted by Microtopography in a Mid-Range Mountain Valley, Bound.-Lay. Meteorol., 165, 333–348,, 2017. a, b, c, d

Pfister, L., Lapo, K., Sayde, C., Selker, J., Mahrt, L., and Thomas, C. K.: Classifying the nocturnal atmospheric boundary layer into temperature and flow regimes, Q. J. Roy. Meteorol. Soc., 145, 1515–1534,, 2019. a, b, c

Poggi, D., Porporato, A., Ridolfi, L., Albertson, J. D., and Katul, G. G.: The Effect of Vegetation Density on Canopy Sub-Layer Turbulence, Bound.-Lay. Meteorol., 111, 565–587,, 2004. a

Poulos, G. S., Blumen, W., Fritts, D. C., Lundquist, J. K., Sun, J., Burns, S. P., Nappo, C., Banta, R., Newsom, R., Cuxart, J., Terradellas, E., Balsley, B., and Jensen, M.: CASES-99: A Comprehensive Investigation of the Stable Nocturnal Boundary Layer, B. Am. Meteorol. Soc., 83, 555–582,<0555:CACIOT>2.3.CO;2, 2002. a, b

Rannik, U.: On the surface layer similarity at a complex forest site, J. Geophys. Res.-Atmos., 103, 8685–8697,, 1998. a, b

Rannik, Ü., Peltola, O., and Mammarella, I.: Random uncertainties of flux measurements by the eddy covariance technique, Atmos. Meas. Tech., 9, 5163–5181,, 2016. a, b

Raupach, M. R., Finnigan, J. J., and Brunei, Y.: Coherent eddies and turbulence in vegetation canopies: The mixing-layer analogy, Bound.-Lay. Meteorol., 78, 351–382,, 1996. a, b

Rebmann, C., Aubinet, M., Schmid, H., Arriga, N., Aurela, M., Burba, G., Clement, R., De Ligne, A., Fratini, G., Gielen, B., Grace, J., Graf, A., Gross, P., Haapanala, S., Herbst, M., Hörtnagl, L., Ibrom, A., Joly, L., Kljun, N., Kolle, O., Kowalski, A., Lindroth, A., Loustau, D., Mammarella, I., Mauder, M., Merbold, L., Metzger, S., Mölder, M., Montagnani, L., Papale, D., Pavelka, M., Peichl, M., Roland, M., Serrano-Ortiz, P., Siebicke, L., Steinbrecher, R., Tuovinen, J.-P., Vesala, T., Wohlfahrt, G., and Franz, D.: ICOS eddy covariance flux-station site setup: a review, Int. Agrophys., 32, 471–494,, 2018. a

Runkle, B. R. K., Wille, C., Gazovic, M., and Kutzbach, L.: Attenuation Correction Procedures for Water Vapour Fluxes from Closed-Path Eddy-Covariance Systems, Bound.-Lay. Meteorol., 142, 401–423,, 2012. a

Sayde, C., Thomas, C. K., Wagner, J., and Selker, J.: High-resolution wind speed measurements using actively heated fiber optics, Geophys. Res. Lett., 42, 10–64,, 2015. a, b, c

Schilperoort, B., Coenders-Gerrits, M., Luxemburg, W., Jiménez Rodríguez, C., Cisneros Vaca, C., and Savenije, H.: Technical note: Using distributed temperature sensing for Bowen ratio evaporation measurements, Hydrol. Earth Syst. Sci., 22, 819–830,, 2018. a, b, c

Schlegel, F., Stiller, J., Bienert, A., Maas, H.-G., Queck, R., and Bernhofer, C.: Large-Eddy Simulation Study of the Effects on Flow of a Heterogeneous Forest at Sub-Tree Resolution, Bound.-Lay. Meteorol., 154, 27–56,, 2015. a

Schotanus, P., Nieuwstadt, F. T. M., and Debruin, H. A. R.: Temperature-Measurement with a Sonic Anemometer and its Application to Heat and Moisture Fluxes, Bound.-Lay. Meteorol., 26, 81–93, 1983. a

Selker, J. S., Thévenaz, L., Huwald, H., Mallet, A., Luxemburg, W., van de Giesen, N., Stejskal, M., Zeman, J., Westhoff, M., and Parlange, M. B.: Distributed fiber-optic temperature sensing for hydrologic systems, Water Resour. Res., 42, W12202,, 2006. a, b

Serafimovich, A., Thomas, C., and Foken, T.: Vertical and Horizontal Transport of Energy and Matter by Coherent Motions in a Tall Spruce Canopy, Bound.-Lay. Meteorol., 140, 429,, 2011. a

Sigmund, A., Pfister, L., Sayde, C., and Thomas, C. K.: Quantitative analysis of the radiation error for aerial coiled-fiber-optic distributed temperature sensing deployments using reinforcing fabric as support structure, Atmos. Meas. Tech., 10, 2149–2162,, 2017. a, b, c, d

Sun, J., Nappo, C. J., Mahrt, L., Belušić, D., Grisogono, B., Stauffer, D. R., Pulido, M., Staquet, C., Jiang, Q., Pouquet, A., Yagüe, C., Galperin, B., Smith, R. B., Finnigan, J. J., Mayor, S. D., Svensson, G., Grachev, A. A., and Neff, W. D.: Review of wave-turbulence interactions in the stable atmospheric boundary layer, Rev. Geophys., 53, 956–993,, 2015. a

Taylor, G. I.: The Spectrum of Turbulence, P. Roy. Soc. A-Math. Phy., 164, 476–490,, 1938. a

Thomas, C. and Foken, T.: Organised Motion in a Tall Spruce Canopy: Temporal Scales, Structure Spacing and Terrain Effects, Bound.-Lay. Meteorol., 122, 123–147,, 2007a. a, b, c

Thomas, C. and Foken, T.: Flux contribution of coherent structures and its implications for the exchange of energy and matter in a tall spruce canopy, Bound.-Lay. Meteorol., 123, 317–337,, 2007b. a, b, c

Thomas, C., Mayer, J.-C., Meixner, F. X., and Foken, T.: Analysis of Low-Frequency Turbulence Above Tall Vegetation Using a Doppler Sodar, Bound.-Lay. Meteorol., 119, 563–587,, 2006. a

Thomas, C., Martin, J. G., Goeckede, M., Siqueira, M. B., Foken, T., Law, B. E., Loescher, H. W., and Katul, G.: Estimating daytime subcanopy respiration from conditional sampling methods applied to multi-scalar high frequency turbulence time series, Agr. Forest Meteorol., 148, 1210–1229,, 2008.  a, b

Thomas, C. K.: Variability of Sub-Canopy Flow, Temperature, and Horizontal Advection in Moderately Complex Terrain, Bound.-Lay. Meteorol., 139, 61–81,, 2011. a, b

Thomas, C. K., Kennedy, A. M., Selker, J. S., Moretti, A., Schroth, M. H., Smoot, A. R., Tufillaro, N. B., and Zeeman, M. J.: High-Resolution Fibre-Optic Temperature Sensing: A New Tool to Study the Two-Dimensional Structure of Atmospheric Surface-Layer Flow, Bound.-Lay. Meteorol., 142, 177–192,, 2012. a, b, c, d, e, f, g, h, i, j

Träumner, K., Wieser, A., Ruck, B., Frank, C., Röhner, L., and Kottmeier, C.: The suitability of Doppler lidar for characterizing the wind field above forest edges, Forestry, 85, 399–412,, 2012. a

Tyler, S. W., Selker, J. S., Hausner, M. B., Hatch, C. E., Torgersen, T., Thodal, C. E., and Schladow, S. G.: Environmental temperature sensing using Raman spectra DTS fiber-optic methods, Water Resour. Res., 45, W00D23,, 2009. a

van de Boer, A., Moene, A. F., Graf, A., Schüttemeyer, D., and Simmer, C.: Detection of Entrainment Influences on Surface-Layer Measurements and Extension of Monin–Obukhov Similarity Theory, Bound.-Lay. Meteorol., 152, 19–44,, 2014. a

van de Giesen, N., Steele-Dunne, S. C., Jansen, J., Hoes, O., Hausner, M. B., Tyler, S., and Selker, J.: Double-ended calibration of fiber-optic Raman spectra distributed temperature sensing data, Sensors-Basel, 12, 5471–5485,, 2012. a, b

Vercauteren, N., Mahrt, L., and Klein, R.: Investigation of interactions between scales of motion in the stable boundary layer, Q. J. Roy. Meteorol. Soc., 142, 2424–2433,, 2016. a

Vickers, D. and Thomas, C. K.: Observations of the scale-dependent turbulence and evaluation of the flux–gradient relationship for sensible heat for a closed Douglas-fir canopy in very weak wind conditions, Atmos. Chem. Phys., 14, 9665–9676,, 2014. a

Wilczak, J. M., Oncley, S. P., and Stage, S. A.: Sonic Anemometer Tilt Correction Algorithms, Bound.-Lay. Meteorol., 99, 127–150,, 2001. a

Zeeman, M. J., Selker, J. S., and Thomas, C. K.: Near-Surface Motion in the Nocturnal, Stable Boundary Layer Observed with Fibre-Optic Distributed Temperature Sensing, Bound.-Lay. Meteorol., 154, 189–205,, 2015. a, b, c, d

Short summary
We evaluated the suitability of fiber-optic distributed temperature sensing (DTS) for observing spatial (>25 cm) and temporal (>1 s) details of airflow within and above forests. The DTS measurements could discern up to third-order moments of the flow and observe spatial details of coherent flow motions. Similar measurements are not possible with more conventional measurement techniques. Hence, the DTS measurements will provide key insights into flows close to roughness elements, e.g. trees.