Coherent Momentum Exchange above and within a Scots Pine Forest

Biorthogonal decomposition (BOD) is used to detect and study synchronous coherent structures occurring at multiple levels in the vertical momentum flux (u1w1) within and above a planted Scots pine forest during a 12-week continuous measurement period. In this study, the presented method allowed for the simultaneous detection and quantification of the number of coherent structures (N), their duration (D) and separation (S) at five measurement heights (z1–z5) covering the range z1/h = 0.11 to z5/h = 1.67, with h being the mean stand height at the measurement site. Results presented for five different exchange regimes (C1–C5) and for four different atmospheric stability conditions (stable, transition to stable, near-neutral, forced convection) demonstrate that during the measurement period, above-canopy momentum flux was only to a limited extent involved in the evolution of spatiotemporal momentum flux patterns found within the below-canopy space. Fully-coupled turbulent momentum exchange over the investigated height range occurred during 19% of all analyzed half-hourly datasets. Across the analyzed exchange regimes, the median contribution of strong sweeps and ejections to total momentum transfer above the canopy varied between 30% and 39% while covering 28%–32% of the time. In the below-canopy space, the contribution of coherent structures varied between 19% and 21% while covering the same amount of time. This suggests that momentum transfer through synchronous coherent structures is very efficient above the forest canopy, but attenuated in the below-canopy space. Since the majority of the presented results agrees well with the results from previous studies that analyzed coherent structures at single levels, the BOD is a promising tool for the consistent investigation of synchronous coherent structures at multiple measurement heights.


Introduction
Various methods have been applied to detect and investigate strong, organized turbulent structures that intermittently exchange substantial amounts of heat, mass, momentum and other scalar quantities between plant canopies and the roughness sublayer located directly above them [1]. These structures are often referred to as coherent structures, because their characteristics appear coherently at several levels within and above plant canopies [2,3] and can be separated from weaker, small-scale fluctuations in turbulent flow [4]. Two quantitative methods, which have been most commonly applied to investigate coherent structure characteristics in data measured in the field, are quadrant analysis [5,6] and wavelet analysis [7][8][9]. As can be expected, there are considerable differences between these methods in their coherent structure detection capabilities [7,10]. Furthermore, the definition of coherent structures varies in previous studies, which introduces differences in coherent structure properties [11].
Although the same coherent structures are detected at several levels within and above forest canopies, results from previous studies that applied quadrant and/or wavelet analysis were often reported from single-level evaluations [12][13][14]. Wavelet analysis requires the determination of peak scales, which are derived from the measurements at each height separately [15]. This can pose a problem when investigating coherent structures that occur simultaneously at several heights, since different peak scales introduce differences in the appearance of coherent structures. In a recent study, the wavelet analysis was used to evaluate synchronous coherent structures at two different measurement heights by introducing time offsets between the heights [16]. Thus, for the investigation of coherent structure characteristics, it is more desirable to have a method that extracts structures from time series that are occurring at multiple heights simultaneously. Then, the spatiotemporal evolution and propagation of coherent structures could consistently be analyzed over the height range covered by the airflow measurements.
Since tower-based micrometeorological measurements in plant canopies are mainly carried out in the canopy layer and in the roughness sublayer, which extends up to two or three canopy heights [17], there is a great interest in improving the interpretability of the dynamics of coherent structures when moving through plant canopies and in quantifying the contribution of coherent structures to total turbulent canopy-atmosphere exchange. This is important when chemical concentration gradients [18,19] and turbulent carbon and energy fluxes are measured within and above plant canopies [20] to study and model matter and energy cycles. A better understanding of coherent structure propagation into forest canopies is also of importance when tree response to wind loading is investigated [21,22].
To gain a deeper understanding of the spatiotemporal evolution and propagation of coherent structures above and within tall plant canopies, it is necessary to detect and separate coherent structures from small-scale turbulence [14]. A method that allows for this separation, as well as for the simultaneous analysis of several time series is the biorthogonal decomposition (BOD) [23]. The BOD is able to capture important common time-space characteristics in high-dimensional datasets. In the field of plant canopy turbulence, the BOD has already successfully been applied in previous studies focusing on the detection of coherent patterns in the interactions of airflow and plant canopies [22,24,25].
The objectives of this study are to investigate the spatiotemporal evolution of synchronous coherent structures occurring in the vertical turbulent momentum flux within and above a planted Scots pine forest canopy and to quantify their typical characteristics under different atmospheric stability conditions and exchange regimes. For that purpose, the BOD is applied on wind vector component time series, which were measured at five levels within and above a Scots pine plantation canopy.

Measurement Site
Airflow measurements were carried out at the forest research site Hartheim, which is located approximately 25 km southwest of Freiburg (southwest Germany) in the flat southern Upper Rhine Valley (47˝56 1 04"N, 7˝36 1 02"E, 201 m above sea level). The research site was operated for more than 40 years by the Meteorological Institute of the University of Freiburg. The forest at the measurement site is a single-layered Scots pine (Pinus sylvestris L.) plantation, which was established in the 1960s. In the year 2014, the forest's mean height (h) was approximately 18 m, and its mean stand density around the measurement site was 580 trees¨ha´1. The mean plant area index (PAI) is 1.5. Figure 1a shows the normalized plant area density (PAD) profile at the measurement site.

Airflow Measurements
During the period from 3 June 2014-25 August 2014, the wind vector components in the x-(u), y-(v) and z-direction (w) were measured (sampling rate: 10 Hz) above and within the forest canopy using five ultrasonic anemometers (R.M. Young Company, Type 81000VRE). The ultrasonic anemometers were mounted on a 30 m high lattice tower at 2 m (z1/h = 0.11), 9 m (z2/h = 0.50), 18 m (z3/h = 1.00), 21 m (z4/h = 1.15) and 30 m (z5/h = 1.67) above ground level (a.g.l.). Since airflow from northern and southern directions dominates at the measurement site (>90% during the measurement period; Figure 1b), the ultrasonic anemometers were oriented to the west on 1.5 m-long supporting booms to minimize the influence of the tower on the wind measurements. Furthermore, all data from eastern directions were excluded from the analysis. During the measurement period, half-hourly wind speed was lower than 5.0 m·s −1 over 98% of the time. The maximum 30-min mean wind speed value at z5 was 9.1 m·s −1 .

Data Processing
First of all, the 1750 available half-hourly time series of wind vector data were despiked with a one-step procedure. The despiking method analyzed the mean and standard deviation for each 300-s window. Spikes were defined as values with an amplitude of at least 3.5 standard deviations away from the 300-s median values. If spikes were detected, then they were replaced by the corresponding values of a median filtered version of the time series. For a better comparison with previous studies, a double rotation was applied on the wind vector data to align the x-axis of the wind vector coordinate system into mean wind direction and to set the mean lateral and vertical wind vector components to zero [26,27]. The data of all measurement heights were rotated according to the rotation angles derived from the data at z5 to ensure that the same coordinate system was used for all time series. After that, Reynolds decomposition was used to split the rotated wind vector components into mean (denoted by an overbar) and turbulent parts (denoted by a prime).
To highlight coherent, organized structures in the wind vector component time series, a band pass infinite impulse response (IIR) filter was used to smooth out short-term fluctuations and eliminate low-frequency parts. Based on previous results on the typical duration of coherent structures in and above the Scots pine forest [22], the cut-off frequencies chosen to separate coherent structures from long-term variations and short-term fluctuations were 0.003 Hz and 0.05 Hz. Thus, the filter removes all variations longer than 300 s and also all fluctuations smaller than event durations of 20 s from the u and w time series. The cut-off frequencies were found to be suitable to include all relevant temporal scales associated with coherent w u   flux structures.

Airflow Measurements
During the period from 3 June 2014-25 August 2014, the wind vector components in the x-(u), y-(v) and z-direction (w) were measured (sampling rate: 10 Hz) above and within the forest canopy using five ultrasonic anemometers (R.M. Young Company, Type 81000VRE). The ultrasonic anemometers were mounted on a 30 m high lattice tower at 2 m (z 1 /h = 0.11), 9 m (z 2 /h = 0.50), 18 m (z 3 /h = 1.00), 21 m (z 4 /h = 1.15) and 30 m (z 5 /h = 1.67) above ground level (a.g.l.). Since airflow from northern and southern directions dominates at the measurement site (>90% during the measurement period; Figure 1b), the ultrasonic anemometers were oriented to the west on 1.5 m-long supporting booms to minimize the influence of the tower on the wind measurements. Furthermore, all data from eastern directions were excluded from the analysis. During the measurement period, half-hourly wind speed was lower than 5.0 m¨s´1 over 98% of the time. The maximum 30-min mean wind speed value at z 5 was 9.1 m¨s´1.

Data Processing
First of all, the 1750 available half-hourly time series of wind vector data were despiked with a one-step procedure. The despiking method analyzed the mean and standard deviation for each 300-s window. Spikes were defined as values with an amplitude of at least 3.5 standard deviations away from the 300-s median values. If spikes were detected, then they were replaced by the corresponding values of a median filtered version of the time series. For a better comparison with previous studies, a double rotation was applied on the wind vector data to align the x-axis of the wind vector coordinate system into mean wind direction and to set the mean lateral and vertical wind vector components to zero [26,27]. The data of all measurement heights were rotated according to the rotation angles derived from the data at z 5 to ensure that the same coordinate system was used for all time series. After that, Reynolds decomposition was used to split the rotated wind vector components into mean (denoted by an overbar) and turbulent parts (denoted by a prime).
To highlight coherent, organized structures in the wind vector component time series, a band pass infinite impulse response (IIR) filter was used to smooth out short-term fluctuations and eliminate low-frequency parts. Based on previous results on the typical duration of coherent structures in and above the Scots pine forest [22], the cut-off frequencies chosen to separate coherent structures from long-term variations and short-term fluctuations were 0.003 Hz and 0.05 Hz. Thus, the filter removes all variations longer than 300 s and also all fluctuations smaller than event durations of 20 s from the Atmosphere 2016, 7, 61 4 of 20 u 1 and w 1 time series. The cut-off frequencies were found to be suitable to include all relevant temporal scales associated with coherent u 1 w 1 flux structures.

Classification of Atmospheric Stability
The classification of atmospheric stability was based on half-hourly values of the Obukhov length Λ [28], which was calculated according to [29]: with the virtual potential temperature θ v , the friction velocity u˚" , the von Kármán constant κ = 0.4 and the gravitational acceleration g = 9.81 m¨s´1. Following [30], the variations of heat and momentum flux were analyzed as a function of the canopy-top stability parameter ζ 3 " z 3 {Λ z3 to determine the boundary values for ζ 3 for the different atmospheric stabilities. Then, all available half-hourly datasets were assigned to the six ζ 3 derived stability classes listed in Table 1. It has to be noted that theoretically, Λ is not applicable over heterogeneous surfaces. However, Λ calculated above forest canopies still approximates the prevalent atmospheric stability conditions reasonably well and is therefore widely used for non-ideal surfaces, as well.
An interval length of 30 min might have an effect on within-canopy turbulent exchange under stable conditions [6,31]. Therefore, frequency characteristics of u 1 and w 1 time series were compared under different stability conditions using Fourier analysis. Since no stability-dependent, deterministic differences in Fourier spectra were found, 30 min was used as the interval length under all stability conditions. Following [32], all datasets classified as very stable and free convection were excluded from further analysis.

Wavelet Analysis
In a number of previous studies that focused on the detection and analysis of coherent structures as part of turbulent exchange processes through plant canopies, wavelet analysis was applied [7,8,33,34]. In comparison to Fourier transform, the wavelet transform shows good resolution in both frequency and time domains, which is achieved by decomposing a time series with a family of wavelet functions. The wavelet functions are produced by dilatation and translation of a mother wavelet ψ n,s (t) [35]: ψ n,s " ψˆt´n s˙ ( 2) with n being the translation parameter, which corresponds to the position of the wavelet on the time axis, and s is the scale dilatation parameter, which corresponds to the width of the wavelet function. Wavelet functions that are commonly used in plant canopy turbulence analysis include the Mexican hat wavelet [7,8,11,12,34,36], the complex Morlet wavelet [32,33,37] or the Haar wavelet [13,36]. The continuous, one-dimensional wavelet transform of a time series f (t) is defined as the convolution of f (t) with ψ n,s (t) [35]: where C n (s) are the wavelet coefficients. The wavelet variance W(s), i.e., the amount of signal energy contained at s, can be calculated as follows [11]: Wpsq "`8 Calculating W(s) over all scales yields the wavelet variance spectrum. From analysis of its maxima (W(s peak )), the peak scale s peak corresponding to the typical duration of coherent structures can be detected [38]. The wavelet coefficients C peak associated with W(s peak ) are used for further analysis. The zero-crossings of these coefficients allow the detection of abrupt changes in time series, which are commonly used to detect individual coherent structures [7,8,11,12,39].
To each s, a corresponding Fourier-equivalent frequency (F a ) and, thus, a corresponding time scale, can be calculated from the center frequency F c of the used wavelet and the sampling frequency ∆ [40]: The detection method for coherent structures presented in this paper was compared to wavelet analysis. For this comparison, the Mexican hat (MH) wavelet was used to detect coherent structures in half-hourly u 1 and w 1 time series, because it is well-localized in time [8]. MH-coefficients were calculated for the scales 50-750, which correspond to time periods from 20-300 s. If W(s) showed several prominent peaks, then the peak scale corresponding to the shortest period was assumed to be associated with coherent structures [12].

Biorthogonal Decomposition
To investigate coherent structures in the momentum flux that are common to more than one measurement height, a method is desired that detects coherent structures in airflow time series over several heights simultaneously. The comparison of results from wavelet analysis over several heights is not always straightforward because height-specific s peak may not always be the same, leading to a different coherent structure appearance at every height.
The method that is proposed for simultaneous detection of coherent structures over several heights is the biorthogonal decomposition (BOD) [23]. For the application of the BOD, the time series are compiled into space-time signals Upx, tq on (XˆT) with X being the set of measurement points in space and T being the corresponding measurement times to each measurement point. It can be expressed as: with k being the number of BOD components, α k " a λ k are weighting factors, λ k are roots of the eigenvalues, µ k (t) are temporal modes, ν k (x) are spatial modes and (ν k (x), µ k (t)) form a set of normalized orthogonal functions. It was demonstrated that the eigenmodes of the spatial operator: are spatial modes with the corresponding eigenvalues λ k [23]. The asterisk indicates the complex conjugate of Upx, tq. The eigenmodes of the temporal correlation operator: are temporal modes with the same set of eigenvalues λ k . By multiplication with ? α k , the spatial information of the signal can be separated from its temporal information [23,41]. Since the norm of the series Upx, tq converges, it is possible to truncate the series to the first k terms [23].
The total energy of the analyzed signals is defined as [41]: The explained variance of each k-th term in the series Upx, tq, called the relative eigenvalue λ rk , can be calculated according to [41]: Information about phase lags in time and space between the time series of Upx, tq are contained in the different behavior of the chronos and topos [24,25,41].
In this study, the BOD is applied to the turbulent wind vector components γ " u 1 , w 1 ( : Since wind vector components were measured at five levels, the BOD yields a maximum of five components after the decomposition of U γ .
The main purpose of the application of the BOD is to reduce the dimensionality of U γ . Normally, not all five BOD components contain information that is useful for the detection of coherent structures in the u 1 w 1 time series. Therefore, the number of useful BOD components derived from the u 1 and w 1 time series was determined for every 30-min interval by analyzing scree plots [42]. Here, scree plots show BOD component-specific eigenvalues as a function of the number of BOD components. The Kaiser criterion [43], which is defined as the mean value of all eigenvalues, was used to determine the number of useful BOD components m γ of U γ . The example shown in Figure 2 leads to one relevant BOD component that is used for the reconstruction. Every BOD component with an eigenvalue larger than the Kaiser criterion was considered to contain useful information and was used to reconstruct decomposed U γ : The reconstructed half-hourly U γ,r time series were then used to calculate u 1 r w 1 r time series: Based on the u 1 r w 1 r time series, characteristics of coherent structures in turbulent momentum exchange through the forest canopy were analyzed.

Classification of Exchange Regimes
Following [33], five exchange regimes (C1-C5) were defined that enable the case-by-case detection of u 1 r w 1 r patterns that form when coherent structures propagate through the height range covered by the airflow measurements. Overall, 16 combinations of momentum flux information are possible ( Table 2). These combinations were assigned to C1-C5 based on the absolute value of the Pearson correlation coefficient (R), which was calculated between u 1 r w 1 r time series of adjacent heights for all half-hourly intervals. After evaluating different R-thresholds, reconstructed momentum flux was defined to be coupled between two adjacent measurement heights when R ě 0.8. The measurement height z 5 is considered as the origin of the penetration depth of coherent structures.

Classification of Exchange Regimes
Following [33], five exchange regimes (C1-C5) were defined that enable the case-by-case detection of       If only one BOD component is used to build u 1 r and w 1 r , then it follows that |R| = 1 between all heights. This is due to the fact that only the absolute values of u 1 r and w 1 r change. If the calculation of either u 1 r or w 1 r involves more than one BOD component, then |R| < 1.
Applying this correlation criterion, C1 corresponds to all half-hourly u 1 r w 1 r time series at z 1 -z 4 that are uncoupled with u 1 r w 1 r time series at z 5 . If C5 is detected, then u 1 r w 1 r is fully coupled from z 5 -z 1 .

Detection and Classification of Coherent Structures
Similar to quadrant analysis, coherent structures were detected when local u 1 r w 1 r variations fell below -Hˇˇu 1 r w 1 rˇ, with H = 3 to consider only strong events [6,44]. Detected coherent structures were classified according to the conventions of quadrant analysis based on the signs of u 1 and w 1 [5,45,46]: sweep (+u 1 ,´w 1 ) and ejection (´u 1 , +w 1 ).
After coherent structure detection, the signs of u 1 r and w 1 r associated with individual coherent structures were determined after zero-crossings of u 1 r w 1 r and classified as sweep or ejection. The number (N meas ) of coherent structures, their duration (D meas ) and separation (S meas ) were determined for all half-hourly intervals. To account for bias in N meas , D meas and S meas values obtained under atmospheric stability conditions with generally lower wind speed (stable, transition to stable, forced convection), all values are normalized. The normalization factor β is calculated as the ratio between the mean wind speed M and the mean wind speed under near-neutral conditions M NN : β " M{ M NN . Since N meas values are smaller due to lower wind speed, N meas is normalized as follows: D meas and S meas values are multiplied by β since D meas and S meas are larger due to lower wind speed: S " β S meas (16) Typical coherent structure characteristics are presented by median values ( r N, r D, r S) and the interquartile range (IQR) in order to provide robust estimation of the central tendencies in the total number of coherent structures. The Wilcoxon rank sum (WRS) test [47] was used to evaluate whether there is a significant difference (level of significance α = 0.05) in coherent structure properties determined for different heights, coherent structure types, exchange regimes and stability conditions. Using Taylor's hypothesis of frozen turbulence [48], the length scale L of a coherent structure can be obtained by: L " M C D with M C " 1.8 M z 3 being the convective velocity [5] and M z 3 being the mean wind speed at canopy height; D is the mean normalized duration of a detected coherent structure. The contribution of sweeps to the overall turbulent momentum exchange (F coh,sw ) and the contribution of ejections to the overall turbulent momentum exchange (F coh,ej ) were calculated as follows [13]: where n sw is the number of detected sweeps, n ej is the number of detected ejections, m i is the number of continuous 10-Hz time series of u 1 r w 1 r during the i-th detected coherent structure and N tot is the total number of measurement points in a 30-min interval. Then, the total contribution of both, sweeps and ejections, to the overall turbulent momentum exchange was calculated: The exchange efficiency E was calculated from half-hourly F coh [11]: with C being the net time cover of all coherent structures, which was calculated by: where T = N tot /10 Hz and D i is the duration of the individual coherent structure i. Furthermore, the attenuation A of the amplitude at z i was calculated for every coherent structure n detected at height z 5 : A n pz i q " maxˇˇu 1 r,n w 1 r,n pz i qˇˇ{maxˇˇu 1 r,n w 1 r,n pz 5 qˇˇ (23) A flowchart (Figure 3) summarizes all steps involved in the data processing. Furthermore, the attenuation A of the amplitude at zi was calculated for every coherent structure n detected at height z5: A flowchart (Figure 3) summarizes all steps involved in the data processing.

Wavelet Analysis vs. BOD
To highlight the differences in the temporal behavior of C peak and reconstructed momentum flux u 1 r w 1 r , Figure 4 shows a comparison of u 1 w 1 with u 1 r w 1 r and C peak over a half-hourly interval (12:30-13:00 CET, 2 August 2014). One useful BOD component explaining 79% of total U u 1 variance was used to build u 1 r , and one BOD component explaining 82% of total U w 1 variance was used to build w 1 r . speak vary from height to height, which poses a problem when comparing similarities between w u   of adjacent heights.  The calculation of the absolute differences between the correlation coefficients obtained from BOD and wavelet analysis over all half-hourly intervals leads to the result that the correlation between r u , r w and r r w u   calculated for any two adjacent heights is 21%, 32% and 46% stronger for BOD than for wavelet analysis (Table 3). The colors indicate momentum flux at z 1 -z 5 . Variance explained by u r ' using one BOD component is 79%; variance explained by w r ' using one BOD component is 82%. To improve visibility, u r 'w r ' calculated at z 2 and z 1 was enhanced five-and 20-fold. For comparison; (c) shows the wavelet coefficients C peak related to z 1 -z 5 .
The wavelet analysis is carried out separately at every height. As a consequence, the peak scales s peak vary from height to height, which poses a problem when comparing similarities between u 1 w 1 of adjacent heights. The calculated peak scales for z 1 -z 5  On the other hand, the BOD includes z 1 -z 5 simultaneously and produces similar temporal u 1 r w 1 r behavior at all measurement heights. Thus, u 1 r w 1 r extracts the features in u 1 w 1 that are common to all heights.
The calculation of the absolute differences between the correlation coefficients obtained from BOD and wavelet analysis over all half-hourly intervals leads to the result that the correlation between u 1 r , w 1 r and u 1 r w 1 r calculated for any two adjacent heights is 21%, 32% and 46% stronger for BOD than for wavelet analysis (Table 3).

Exchange Regimes
In the measurement period, C1 covers 42.4% of the analyzed half-hourly intervals. The exchange regime C2 was not detected in the analyzed dataset. The exchange regimes C3-C5, which indicate increasing penetration depth of coherent structures into the below-canopy space, were detected in 23.6%, 15.2% and 18.8% of the half-hourly intervals.
In Table 4, fractions of C1-C5 are summarized as a function of atmospheric stability. From that, it can be inferred that C1 is the dominant exchange regime under stable conditions and transition to stable, covering 92.7% and 64.0% of all cases. Furthermore, it is obvious that with increasing instability, the occurrence of coherent structures penetrating into the canopy clearly increases, as the fractions of C4 and C5 increases. Under forced convection, C4 and C5 dominate, i.e., momentum is transferred from above the canopy into the below-canopy space. Table 3. First quartile (Q1), median and third quartile (Q3) of the difference between the correlation coefficients (R) calculated for adjacent heights (z 1 -z 5 ) obtained from the reconstructed time series using BOD (u r ', w r ', u r 'w r ') and peak wavelet coefficients (C peak,u' , C peak,w' , C peak,u'w' ).
Quantity z 1 -z 2 z 2 -z 3 z 3 -z 4 z 4 -z 5 all Rpu 1 r pz i q, u 1 r pz j qq´RpC peak, u1 pz i q, C peak, u1 pz j qq Q1 0.  A comparison of the results obtained by BOD with the results reported in previous studies shows general agreement on the proportions of exchange regimes. A study investigated exchange regimes in a spruce forest based on wavelet analysis during a period of five days with |ζ| ă 1 [9]. The definition of exchange regimes used in that study was based on the vertical profiles of buoyancy exchange of coherent structures [33]: wave motion (Wa), decoupled canopy (Dc), decoupled sub-canopy (Ds), coupled sub-canopy by sweeps (Cs) and fully-coupled canopy (C).
In terms of coupling between above-and below-canopy airflow, C1 and C2 can be compared to Wa and Dc, C3 to Ds and C4 and C5 to Cs and C. The proportions of exchange regimes reported by [9]  Although a comparison between both studies is not straightforward, because [9] used a classification of exchange regimes based on buoyancy exchange, some common results are found. First, most of the time, there is no coupling between above-canopy air and below-canopy air (Wa + Dc = 51.6%; C1 + C2 = 42.4%). Second, the proportions of exchange regimes associated with a coupling between above-and below-canopy air are in a similar range (Cs + C = 30.4%; C4 + C5 = 34.0%). Furthermore, specific exchange regimes seem to occur at certain times throughout the day ( Figure 5). While C1 was observed more often at night, C3-C5 occurred mostly during daytime with C4 occurrence peaking in the morning and C5 occurrence peaking at noon, which is in good agreement with previous findings [33]. With increasing solar radiation, however, stability decreases, which increases the vertical transport. Thus, coupling between air above and below the canopy becomes more likely.

Detection of Coherent Structures
After applying the BOD, sweeps and ejections were detected separately, and their properties were quantified. In Figure 6, the detection is exemplified for z 5 . It shows u 1 r w 1 r being calculated over an arbitrarily-chosen half-hourly interval (15:00-15:30 CET, 3 June 2014). The reconstruction of u 1 w 1 is based on two BOD components that explain 85.3% of U u 1 variance and 81.6% of U w 1 variance. decreases, which increases the vertical transport. Thus, coupling between air above and below the canopy becomes more likely.

Detection of Coherent Structures
After applying the BOD, sweeps and ejections were detected separately, and their properties were quantified. In Figure 6, the detection is exemplified for z5. It shows  (14)- (16). Table 5 provides comprehensive information on the central tendencies of N for z1-z5 under different exchange regimes. In the entire measurement period Ñ varied between seven and 21, and IQR represented 6-27 detected coherent structures, which is comparable to the range of N found in previous studies using wavelet analysis [11,16].
Under fully-coupled conditions, the results of the WRS test demonstrated that at z5, Ñ related to sweeps and ejections were significantly larger than Ñ below the canopy. Thus, it can be inferred that not all coherent structures detected at z5 contributed to the momentum exchange into the forest. Moreover, for C3-C5, the number of sweeps was significantly larger than the number of ejections at all measurement heights. Therefore, sweeps were the dominating coherent structure type.    Counting sweeps and ejections gives N meas = 14 and N meas = 16. For each sweep and ejection, its duration, as well as the separation of a subsequent event of the same type is determined for z 1 -z 5 by monitoring the crossings of -Hˇˇu 1 r w 1 rˇ. For example, the first ejection shown in Figure 6 has a duration of D meas = 24 s and is followed by a second ejection event after S meas = 41 s. After determination of N meas , D meas and S meas for all 30-min intervals, the values were converted to normalized N, D and S according to Equations (14)- (16). Table 5 provides comprehensive information on the central tendencies of N for z 1 -z 5 under different exchange regimes. In the entire measurement period r N varied between seven and 21, and IQR represented 6-27 detected coherent structures, which is comparable to the range of N found in previous studies using wavelet analysis [11,16]. Under fully-coupled conditions, the results of the WRS test demonstrated that at z 5 , r N related to sweeps and ejections were significantly larger than r N below the canopy. Thus, it can be inferred that not all coherent structures detected at z 5 contributed to the momentum exchange into the forest. Moreover, for C3-C5, the number of sweeps was significantly larger than the number of ejections at all measurement heights. Therefore, sweeps were the dominating coherent structure type.

Typical Event Duration and Separation
The median values of D varied between 19.2 s and 26.0 s with IQR spanning from 14.4 s-36.6 s ( Table 6). Results from the WRS test showed only significant height-dependent differences under fully-coupled conditions, with D of sweeps and ejections being significantly larger at z 1 than at z 5 . Furthermore, D values of sweeps were significantly shorter than D values of ejections at all heights under C3-C5.
Under all exchange regimes, r S values varied between 75.7 s and 106.6 s at z 1 -z 5 with IQR spanning from 37.0-201.2 s (Table 7). Results from the WRS test showed no height-specific differences. However, results indicated that for C1-C4, S values for sweeps were shorter than S values for ejections at all measurement heights.
Although there is general agreement between the tendencies obtained for N, D and S in this work and the results obtained in previous studies, a comparison of absolute values is not straightforward because in previous studies: (i) different methods have been applied to detect and study coherent structures for only one or several levels separately; (ii) investigation periods were often much shorter; (iii) the reported results do not refer to the same central tendencies in the data; (iv) a diversity of coherent structure definitions have been used, leading to different definitions of D and S; and (v) a possible bias due to different wind speed at different stability conditions was often not considered.  A study investigated coherent structures during a period of five days with |ζ| ă 1 conditions in a spruce forest and found D values in the range 10-30 s [9]. Another study determined a median coherent structure duration of 90-115 s with wavelet analysis and 1.3-1.8 s with quadrant analysis during stable and unstable conditions [7]. Based on visual inspection, D was reported to vary between 24 and 39 s below a deciduous forest canopy and between 20 and 23 s above the canopy during three 30-min intervals with stable, neutral and unstable conditions [2]. Furthermore, based on visual analysis, a mean duration of coherent structures above a pine forest of 33-40 s was determined, with separations of 97-124 s during two 100-min intervals under unstable conditions [49]. Two studies used wavelet analysis to detect coherent structures at a pine forest site during slightly unstable conditions and found characteristic durations and separations around 4 s and 29 s [8,15].

Further Coherent Structure Characteristics
Median values obtained for L ( r L), F coh ( r F coh ), E ( r E), C ( r C) and A ( r A) for different exchange regimes are summarized in Table 8. The r L values spanned between r L = 43.6 m determined for C1 and r L = 96.4 m determined for C4. The median values of F coh varied between 18.7% and 21.3% below the canopy and between 29.8% and 39.4% at and above the canopy. Median transport efficiency included values between r E = 0.6 and r E = 0.7 below the canopy and between r E = 1.1 and r E = 1.3 at the canopy height and above the canopy. Median total time cover of coherent structures varied between r C = 27.5% and r C = 29.2%. Results from the WRS test showed that the values of r L for C4 and C5 were significantly larger than the r L values calculated for C1 and C3, suggesting that only large coherent structures reached the subcanopy space. For C3-C5, the values for r F coh and r E at z 5 were significantly smaller than the values at canopy height, but significantly larger than the values below the canopy. This suggests that momentum flux was enhanced through shear-generated momentum flux at canopy height, but weakened through dissipation in the subcanopy space. Height-specific differences in r E indicate that coherent structures transferred momentum efficiently above the canopy. However, below the canopy, the transport efficiency of momentum flux through coherent structures was weakened. This is also visible in the calculated values for r A. The amplitude of coherent structures involved in momentum transport from z 5 -z 1 was strongly attenuated below the canopy.
Median values obtained for L, F coh , E and C related to different stability conditions are summarized in Table 9.
Results from the WRS test indicated a dependence of r L on atmospheric stability. For all heights, r L was significantly larger during near-neutral conditions and forced convection than under stable conditions and transition to stable. Thus, in agreement with findings from a previous study [11], an increase of r L with increasing instability was found. Under stable stratification, the formation of large coherent eddies was inhibited, which led to the dominance of eddies with smaller sizes. Table 9. Stability-related r L (m), r F coh (%), r E (%) and r L (%) calculated at z 1 -z 5 . Although the magnitude of F coh largely depends on the method of analysis [33], the results above the canopy obtained in this study agree well with the results from previous studies: One study [7] reported median values of F coh between 40% and 48% by analyzing coherent structures in a deciduous forest at canopy height. Another study [11] found mean values of F coh between 38% and 51% at 10 m and 30 m above a mixed surface. Over a pine forest, F coh of 92.7%˘3.2% was reported [49]. In an urban area, a study [39] found that F coh ranged between 41% and 107% at z/h = 1.5 under unstable conditions. Values of F coh between 55% and 95% were derived in a walnut orchard during different stability conditions from quadrant analysis using a hole parameter of H = 3 [6].

Variable
Results from the analysis of the contribution of sweeps and ejections to F coh under different atmospheric stabilities are shown in Figure 7. Values for F coh,sw peaked at canopy height and decreased with height above the canopy, whereas F coh,ej still increased above the canopy. This behavior was also observed in a previous study [33]. Increasing F coh,ej above the canopy indicates the transition from a flow dominated by shear-driven coherent structures to a flow associated with rough-wall boundary layers [6].
The authors of a previous study [3] concluded from their results that the magnitude of sweeps are larger than the magnitude of ejections under unstable conditions and that the strength of both decreases with increasing stability. The larger magnitude of sweeps was also observed in this study. However, a decrease in the strength of sweeps and ejections with increasing stability was not found.
An investigation [49] of coherent structures above a pine forest under unstable conditions (-1 < ζ < 0) came to the conclusion that the proportions of sweeps are higher compared to ejections at canopy height and that the importance of ejections increases with increasing height from the canopy top. Numerous studies [1,2,6,33,45] found similar results.
An increase of ejections with height was also observed in this study. Furthermore, the authors of some studies [2,50] concluded from their results that sweeps and ejections contribute most to the overall transport in the canopy and that the sweeps-related contribution dominates. This statement can be confirmed from the findings obtained in this study.
Results from the analysis of the contribution of sweeps and ejections to Fcoh under different atmospheric stabilities are shown in Figure 7. Values for Fcoh,sw peaked at canopy height and decreased with height above the canopy, whereas Fcoh,ej still increased above the canopy. This behavior was also observed in a previous study [33]. Increasing Fcoh,ej above the canopy indicates the transition from a flow dominated by shear-driven coherent structures to a flow associated with rough-wall boundary layers [6]. The authors of a previous study [3] concluded from their results that the magnitude of sweeps are larger than the magnitude of ejections under unstable conditions and that the strength of both decreases with increasing stability. The larger magnitude of sweeps was also observed in this study. However, a decrease in the strength of sweeps and ejections with increasing stability was not found.
An investigation [49] of coherent structures above a pine forest under unstable conditions (-1 < ζ < 0) came to the conclusion that the proportions of sweeps are higher compared to ejections at canopy height and that the importance of ejections increases with increasing height from the canopy top. Numerous studies [1,2,6,33,45] found similar results.
An increase of ejections with height was also observed in this study. Furthermore, the authors of some studies [2,50] concluded from their results that sweeps and ejections contribute most to the overall transport in the canopy and that the sweeps-related contribution dominates. This statement can be confirmed from the findings obtained in this study.
It was also reported that the contribution of ejections increased in the below-canopy space and that the contributions of sweeps and ejections become equal further down to the ground [33]. An increase of ejections with decreasing height was not found in this work, whereas the contributions of sweeps and ejections did approach each other near the ground. During stable conditions, ejections dominated in the sub-canopy space. It was also reported that the contribution of ejections increased in the below-canopy space and that the contributions of sweeps and ejections become equal further down to the ground [33]. An increase of ejections with decreasing height was not found in this work, whereas the contributions of sweeps and ejections did approach each other near the ground. During stable conditions, ejections dominated in the sub-canopy space.
The r E values determined above the canopy were similar to the r E values (1.05-1.66) found by previous studies [11,16]. Under different atmospheric stabilities, r C varied between 24.4% and 30.8%. Results from the WRS test showed significant height-specific differences, as well as stability-specific differences. During transition to stable, r C at z 5 was significantly larger than r C at z 4 and z 3 , but also significantly smaller than r C below the canopy. Under stable conditions, r C at z 5 was larger than directly at z 4 and z 3 . Since F coh always peaked at the canopy top, this result indicates the occurrence of strong, but short sweeps and ejections during transition to stable and stable conditions. Overall, the r C values are in good agreement with previous studies ( [11]: r C = 25%-45%; [49]: r C = 28%-38%).

Conclusions
The contribution of synchronous coherent structures to the exchange of momentum was investigated in a planted Scots pine forest. Coherent structures are known to substantially contribute to the transfer of momentum in vegetation canopies. To quantify their contribution to the total momentum transfer, single-level detection methods have been applied in previous studies. The disadvantage of these methods, however, is that detection results differ between the investigated levels. Therefore, the BOD was presented as a tool for the synchronous detection of coherent structures in the vertical momentum exchange between the Scots pine forest and the atmosphere. The BOD simultaneously treats the momentum flux measured at multiple levels, which is a major advantage over commonly-applied single-level methods, such as wavelet analysis. Furthermore, the BOD allowed for the analysis of coupling between above-and below-canopy momentum fluxes.
It was found that fully-coupled turbulent momentum exchange over all measurement heights occurred during 19% of all analyzed half-hourly datasets, mostly during daytime. At night, under stable conditions, the subcanopy layer was most frequently decoupled from the above-canopy layer. From the synchronous coherent structures extracted from the momentum flux, sweep and ejection characteristics, such as the number of occurrences, duration and separation, were quantified. The median values for the number of occurrence per 30-min interval, duration and separation were 7-21, 14.4-36.6 s and 75.7-106.6 s. The contribution of the detected sweeps and ejections to total momentum flux, their transport efficiency, as well as their time cover were 16.1%-40.8%, 0.6-1.5 and 24.4%-30.8%.
Results suggest that momentum transfer through synchronous coherent structures is very efficient above the forest canopy, but attenuated in the below-canopy space, which might result from the generation of small-scale turbulent structures by the Scots pine trees, which induces a short-circuiting of momentum transfer. This implies that canopy characteristics, such as density and vertical distribution of biomass, have considerable effects on the characteristics of coherent structures.
Although there is general agreement between the majority of the presented results and results obtained in previous studies that used single-level analysis methods, this work is considered as a first case study for coherent structure detection using the BOD. Therefore, there is potential for improvement of the methodology and the reduction of uncertainty associated with the choice of the required thresholds. Further comparisons with wavelet-based detection of coherent structures will show whether the proposed methodology is suitable not only for momentum exchange, but also for the exchange of scalar quantities, such as temperature or CO 2 .