On the Potential of Using Air Pressure Fluctuations to Estimate Wind-Induced Tree Motion in a Planted Scots Pine Forest

: This paper reports statistical relationships between measured airﬂow, air pressure ﬂuctuations, and the wind-induced motion of planted Scots pine trees ( Pinus sylvestris L.). The results presented illustrate the potential of low-cost, ground-based air pressure measurements for monitoring wind-induced tree response. It is suggested that air pressure ﬂuctuations can be used as surrogate information for above-canopy airﬂow, often used to estimate wind loads on forest trees. We demonstrate that air pressure ﬂuctuations can be measured representatively at the forest ﬂoor and correlate very well with wind speed and direction at mean canopy-top (18 m a.g.l.) and above the 18 m high, 56-year-old forest. Their strong correlation (coefﬁcient of determination R 2 > 0.77) allows a good approximation of airﬂow conditions above the canopy, and, with some limitations, in the below-canopy space. Air pressure ﬂuctuations also correlate very well with wind-induced tree motion with a similar correlation to that between wind speed and tree motion. Furthermore, the main directions of wind-induced tree motion agree very well with the propagation direction of air pressure waves. Above-canopy airﬂow measurements in forests with a large vertical extent are rare, and often require tall wind measurement towers. Therefore, we consider the estimation of airﬂow conditions over forests using ground-based air pressure measurements a promising option for monitoring the airﬂow conditions of relevance for predicting wind-induced tree response over large areas using a minimum of measurement infrastructure.


Introduction
The characteristics of near-ground wind fields are important physical site factors that contribute to the development and survival of trees and forests [1,2]. Local wind fields affect trees' physiological, morphological, and biomechanical adaptations to their atmospheric environment [2][3][4][5][6]. However, information on above-and below-canopy wind characteristics is limited because the standard weather stations operated in the networks of the national weather services are located outside forests. Within stand wind characteristics are mainly known from specialized sites of forest research institutes [7,8], research networks such as the European Integrated Carbon Observation System [9], and the International Cooperative Programme on Assessment and Monitoring of Air pollution Effects on Forests [10]. Otherwise, little systematic information and data are available on the wind fields in and over forests.
The reasons for the lack of wind measurement sites in forests are the great effort required to make such measurements, the lack of infrastructure, and the low spatial representativeness of below-canopy measurements. Due to large vertical extent of forests, it is a link between the strength of pressure pumping and the mean above-canopy wind speed was determined. A subsequent study at the same forest site showed that the propagation direction and speed of air pressure fluctuations are also strongly correlated with the abovecanopy wind speed and direction [35].
Since strong functional relationships can be established between wind conditions over forests and small, wind-induced air pressure fluctuations measured near the ground, an alternative to tower-based wind speed and direction measurements in forests is presented in this paper. The proposed method allows the assessment and monitoring of tree reaction to local mean and turbulent airflow. For this, we first assess the functional relationship between the above-canopy wind speed and air pressure fluctuations measured near the ground in a planted Scots pine forest. Then, the air pressure fluctuations were directly related to the wind-induced tree motion. It is shown that the method is suitable for monitoring the means and extremes of the response of forest trees to the wind in situations where no above-canopy wind speed and direction measurements are available.

Study Site
Airflow, air pressure, and tree motion measurements (sampling frequency 10 Hz) were conducted from 3 to 28 July 2020 at the Hartheim forest research site operated by Environmental Meteorology of the University of Freiburg, Germany [44]. The research site is located in a planted Scots pine forest established in 1964, in the flat, north-south aligned southern Upper Rhine valley (47 • 56 04 N, 7 • 36 02 E, 201 m a.s.l.) on the border of Germany with France. During the measurement period, the stand density was 550 trees per hectare, the mean stand height was 18.0 m, and the mean plant area index of the sparse canopy was 1.9, indicating that the forest is thinly stocked ( Table 1). The Scots pine forest extends 3.2 km in N-S and 2.0 km in W-E directions. The nearest forest edge that turns into a larger canopy gap is 200 m west of the measurement site.

Airflow Measurements
Airflow was measured at five scaffold towers (To1 to To5) arranged in a cross approximately aligned along the longitudinal axis of the Rhine Valley ( Figure 1). The distance between the towers erected at the endpoints of the cross was 48 m. At heights of 21 m (z 1 , 3 m above the canopy), 18 m (z 2 , at canopy top), 9 m (z 3 , 9 m below the crown space), and 2 m (z 4 , 2 m above the forest floor), 20 ultrasonic anemometers (R.M. Young Company, USA, type 81000VRE) were mounted to To1 to To5 on 1.5 m long booms to measure the wind vector components in east-west (u), north-south (v), and vertical (w) directions.  (To1 to To5). The seven sample trees (B1 to B7) are marked with green dots. The five red measurement points (M1, M3, M9, M15, M17) were used to assess the propagation direction of air pressure waves. The orange lines indicate the distance between individual measurement points. (b) Side view of one of the five measuring towers where wind speed (WS), wind direction (WD), and P raw were measured at multiple heights.
Since northern and southern wind directions dominate at the research site [45], the booms were positioned pointing west to minimize the influence of the towers on the airflow measurements. Wind speed (WS) and wind direction (WD) were calculated from the horizontal wind vector components. For the analysis presented and discussed below, WS was averaged to 10-min mean values, WD was calculated from 10-min mean values of u and v. Unless otherwise specified, the WS and WD values always refer to z 1 , i.e., to the above-canopy measurements.

Air Pressure Measurements
Small, wind-induced air pressure fluctuations (P raw ) were measured using amplified differential piezoresistive pressure sensors (Honeywell, USA, type SSCSNBN001NDAA5) with a range of ±248.84 Pa. The two-port sensors are part of a low-cost, highly precise measuring system developed and evaluated in previous studies [25,46]. One port of the sensors was connected to a pressure reference to attenuate high-frequency P raw fluctuations. The pressure head used in this study was of the quad-disc type [47,48], made of white 3 mm PVC sheets and in-house 3D-printed parts.
At the forest floor (z 5 = 0 m), 49 air pressure sensors were installed in a 7 × 7 regular grid with a 12 m distance between the measurement points (M1 to M49). The dimensioning of the measuring grid was based on the findings on the correlation of air pressure fluctuations obtained in an earlier investigation at the research site [35]. They reported a strong correlation of air pressure fluctuations up to a horizontal distance of 15 m above the forest floor. In addition, To1 to To5 carried sensors to measure air pressure fluctuations at z 1 -z 4 . The air pressure measuring systems were mounted on the same booms carrying the sonic anemometers.
The radiometric analog outputs of all air pressure sensors were read using an Arduino Nano based microcontroller. An ethernet shield (type Wiznet W5100 Ethernet-Chip) was connected to the Arduino Nano to send P raw data to a central microcomputer (type Raspberry Pi 3 Model B+, Raspberry Pi Foundation, UK) using the message queuing telemetry transport (MQTT) communication protocol.
The P raw time series were first band-pass filtered for the frequency range 0.01-0.10 Hz (P) which controls pressure pumping [25,35]. To assess the aggregated strength of the P raw fluctuations, PPC was calculated over 10-min intervals: with sf being the sampling frequency, N is the number of P raw measurements in a 10-min interval, and ∆P is the difference between two successive P values.
Since wind-tree interactions are commonly analyzed using wind speed at or above canopy top [22,23,49,50], the functional dependence in 10-min intervals of above-canopy WS measured at z 1 and PPC estimated at z 5 was established with a power law. The fit of the power law to the data was always better than the fit of other simple linear regression models: where a 1 , a 2 , and a 3 are coefficients. The power law was fitted using WS measured at To1 to To5 and the PPC values computed at M1 to M49 (PPC 1 to PPC 49 ), yielding 20 tower-related and 49 measurement grid-related empirical models. The direction of P propagation (PD) was calculated using the approach reported in a previous study [35]. Instead of using maxima of cross-correlation functions calculated between P time series, the simultaneous occurrence of P maxima was utilized to determine propagation characteristics of air pressure waves moving across M1, M3, M9, M15, and M17 aligned in the shape of a cross in the southwestern corner of the ground-level air pressure measurement grid. The location of the sensor cross was chosen because the southwestern corner of the measurement grid is the first to be impacted. Due to the similarity of the air pressure measurements, which will be illustrated later, the evaluated sensor cross could also have been positioned at any other location of the measurement grid. The air pressure waves' motion was analyzed only when all five P time series at M1, M3, M9, M15, and M17 showed at least one peak in the investigated 10-min intervals. The peak selection was based on Matlab's (The MathWorks Inc., Natick, Massachusetts) findpeaks function [51] using a data-specific combination of peak prominence of 0.6 and a peak width of 120 values as selection criteria.
Starting from the central measurement point M9, the time differences to the four measurement points at the corners of the measurement cross were quantified. Based on the time differences, the direction of propagation of all examined air pressure waves was then determined and averaged.

Tree Motion Measurements
Wind-induced stem tilt of seven Scots pines (B1 to B7) was measured in x (t x ) and y (t y ) directions over the height range 2.4 to 2.8 m using the Tree Response Sensor (TRS) which is part of the Tree Motion Monitoring System (TreeMMoSys) developed in an earlier study [52]. The height range corresponds to 1/7 of the tree height (HT), where the first antinodal point of vibration of a clamped-free beam was assumed [53]. Tree response data were collected wirelessly and stored on a ground receiver using TreeMMoSys.
To investigate the integral response of the aerial tree parts to the wind, as in a previous study [21], the stem tilt vector (T) was calculated from t x and t y for B1 to B7: The T values were used to calculate 10-min mean interval values (T mean ) and 98th percentiles (T 98 ). While T mean was used to assess the mean response of B1 to B7 to the wind, T 98 was used to quantify the impact of strong and short-term wind loading from turbulent airflow components.

Bi-Orthogonal Decomposition
To reduce the dimensions of the spatiotemporal dynamics of WS and P, the biorthogonal decomposition (BOD) was applied [54]. BOD can capture important common time-space characteristics in high-dimensional data. The five 10-min above-canopy WS time series and the 49 below-canopy P time series were compiled into space-time signals (C) each and decomposed [46]: where m is the number of BOD components, α m = √ λ m are weighting factors, λ m are eigenvalues, µ m (t) are temporal modes, ν m (q) are spatial modes, and (µ m (t), ν m (q)) form a set of normalized orthogonal functions. The Kaiser criterion [55], defined as the mean of all eigenvalues, was used to determine the number of informative BOD components in the decomposed WS (WS BOD ) and PPC (PPC BOD ) time series useful for further analysis.
As a surrogate for the actual mean wind load on the sampled trees, which is still unknown, squared WS BOD WS 2 BOD and PPC BOD were used. Both variables potentially provide information about the temporal dynamics of tree motion but are not an expression of the effective wind load defined in a previous study [19]. The 98th percentiles of P (P 98 ) were used as surrogate values to simulate gust speed in the 10-min intervals. Gust speed was previously used to represent the impact of strong and short-term wind loads on trees in forests [56][57][58][59].

Wavelet Analysis
As demonstrated below, PPC BOD is used to approximate wind speed at various heights above ground, including below-canopy wind speed. However, this article focuses on the above-canopy wind speed. For a general illustration of the approach's potential, wavelet coherence (WCoh) is used as a local correlation measure. It was calculated to measure the 10-min localized correlation between WS BOD and PPC BOD in the study period over 106 scales [60]: where F WS (b,c) and F PPC (b,c) are continuous wavelet transforms of WS BOD at z 1 to z 4 and PPC BOD with the analytic Morlet wavelet at scales b and positions c along the timeline. The asterisk denotes the complex conjugate, and S indicates smoothing in space and time over 12 scales.

Modeling of Tree Reactions to Short-Term Wind Loading
The ensemble learning method bootstrap aggregating (bagging) was used to model T 98 (T 98,mod ). The bagging algorithm averages predictions from multiple regression models to reduce the variance in the prediction by minimizing the mean squared error [61]. The bootstrap procedure was applied to H = 10 randomly sampled unweighted subsets of the predictors T mean and P 98 per Scots pine tree to train ten tree-specific regression models (G h , with h = 1, . . . , 10): Examination of the predictive accuracy based on the coefficient of determination (R 2 ) of different predictor combinations revealed that T mean and P 98 are best suited for predicting T 98 . It is plausible to assume that high T 98 is coupled with high T mean . Since the relationship between T mean and T 98 is not strictly linear, P 98 was used as a surrogate for considering the short-term impact of gust speed on stem tilt. To measure the prediction accuracy and prevent overfitting, out-of-bag subsets that were not used to train the regression models were evaluated.

Wind Speed above and below the Canopy
During the study period, WS averaged across To1 to To5 varied between 0.2 and 5.7 m/s. Above-canopy WS mean = 1.63 m/s ± 0.11 (mean ± standard deviation) indicate very similar wind speed conditions measured over the whole study period at the research site ( Figure 2a). At the forest floor, averaged WS values were much lower, mostly ranging between 0.1 and 0.7 m/s with WS mean = 0.23 ± 0.07 m/s (Figure 2b). Although the standard deviation of WS at the ground obviously exhibits larger relative variability, it was in a similar range in absolute terms as the wind speed variation at z 1 .
At both heights, WS is subject to pronounced daily variations on many days. The variations occur more regularly at z 4 than at z 1 . The highest WS values mainly happen during the day. The daily component in the WS time series can be attributed to the location of the research site in the Upper Rhine valley, where tertiary circulation systems develop recurrently [44]. A secondary, below-canopy wind speed maximum that could result from thermal stability conditions [62] and mean horizontal air pressure gradients [63,64] was observed in 3.8% of all 10-min intervals. A systematic wind direction shear resulting from air pressure conditions in the below-canopy space as reported in a previous study [63] was not found. It is likely that the below-canopy wind direction distributions determined at To1 to To5 reflect the local wind field influenced by the surrounding trees and understorey.

Air Pressure Fluctuations below the Canopy
Below-canopy P (z 3 to z 5 ) values under different above-canopy wind conditions are shown over two different 10-min intervals in Figure 3. With an above-canopy WS = 0.6 m/s measured at To1, corresponding to the 10th percentile of the To1 wind speed distribution, below-canopy P is small, mostly varying between −0.30 and 0.30 Pa. With an above-canopy WS = 2.7 m/s, representing the 90th percentile of the To1 wind speed distribution, below-canopy P measurably increases and fluctuates between −6.0 and 4.0 Pa. In addition to the magnitude of the P fluctuations, their similarity across heights z 3 to z 5 also increased with increasing above-canopy wind speed.
If all 10-min intervals of the study period are included, the mean values of the Pearson correlation coefficient (R) calculated between below-canopy P time series increased with increasing above-canopy WS from R = 0.76 at WS < 1.0 m/s to R = 0.96 at WS > 4.0 m/s. The same tendency was found in the correlation of the 10-min mean P values in dependence of WS. With increasing above-canopy WS measured at To1, the correlation strength rose from R = 0.78 to R = 0.98.
Previous studies also reported a correlation that increases with higher WS between below-canopy P and above-canopy WS [25,26]. In a further study, it was demonstrated that the propagation of P could also be represented as a function of above-canopy WD [35].
However, although an earlier study narrowed the frequency range of wind-induced pressure fluctuations over forest [25], their origin is still unclear. This is not at least because only a few studies reported results from direct measurements of wind-induced air pressure fluctuations [25,26,35,37,38,65].
Air pressure fluctuations at the forest floor can result from turbulence induced by local airflow-forest-interactions when shear instability at the canopy top generates coherent structures [30]. In turbulent canopy airflow, these coherent structures are responsible for efficient transport of heat and momentum between forests and the atmosphere [31][32][33]. However, the results of previous studies are ambiguous. One study was able to demonstrate that air pressure fluctuations measured below a deciduous forest canopy are strongly correlated with the passage of coherent turbulent structures [27,28]. Another study investigating coherent structures in a Scots pine forest did not find such a correlation [45].
The temporal and spatial similarity of 10-min PPC 1 to PPC 49 was strong. The 10-min PPC values were averaged to hourly values for clarity in Figure 4. Based on these hourly values, the daily PPC 1 to PPC 49 cycles in the measurement grid can be better visualized. All points in the measurement grid showed a pronounced mean daily PPC cycle with the highest hourly values of 0.08 Pa/s and more from 13:00 to 15:00 CET. During these hours, above-and below-canopy WS recurrently reached the highest mean values during the day. The lowest PPC values of 0.04 Pa/s and less regularly occurred from 22:00 to 7:00 CET. The intra-hourly spatial PPC differences in the measurement grid are negligible, varying mainly in the narrow range −0.007 Pa/s to 0.008 Pa/s throughout the day.
The low spatiotemporal variability of the above-canopy WS and below-canopy PPC can be represented by only one BOD component, as evaluated by the Kaiser criterion ( Figure 5). Due to the similarity of the above-canopy WS measured at To1 to To5, the first BOD component (WS BOD ) explains 98% of the variance (EV) in the above-canopy WS time series. The high similarity between PPC 1 to PPC 49 values is also demonstrated by the value of EV = 98% associated with the first BOD component (PPC BOD ). For the analyzed WS and PPC data, three separate BOD components each are sufficient to explain almost all of each variable's spatiotemporal variability at the research site. The mean daily cycles of WCoh show the highest values from 12:00 to 20:00 CET ( Figure 6). This reflects the daily variation of WS BOD at the research site, which, as indicated earlier, has a maximum in the afternoon. With higher WS BOD , the amplitudes of P and the PPC values increase, as has been demonstrated in previous studies [25,26]. Over the day, the WCoh levels at z 1 and z 2 are similar, with strong above-canopy WCoh = 0.79 representing the highest 10-min overall value at 15:10 CET. There is a height-dependent lag between the maximum WCoh values, which cannot yet be explained. The highest WCoh values are reached latest in the day at z 4 . At night, when WS BOD is low, WCoh is also low. The lowest WS BOD values are found at all measurement heights between 2:00 and 6:00 CET in the morning. When WCoh is averaged over the entire study period, its decrease with reduced height above the ground becomes apparent. From z 1 to z 4 , the study period mean values of WCoh are 0.69, 0.69, 0.62, and 0.55. Although the level of correlation decreases from above the forest toward the forest floor, the association between WS BOD and PPC BOD is sufficient to model wind speed from z 3 to z 1 with reasonable certainty using PPC. The consistent lower level of wavelet coherence toward the forest floor undermines the findings of previous studies [27,28,34] of a dominant modulation of the near-surface air pressure pattern by coherent structures.

Approximation of Wind-Induced Tree Motion
The dependence of PPC BOD on WS 2 BOD is strong. As known from previous studies at the research site [25,26,35], below-canopy PPC is high when above-canopy WS is high and vice versa. This relationship is very well reproduced (R 2 = 0.95) by the power law (Figure 7a). The high EV values associated with the first BOD components, together with the strong correlation between WS BOD and PPC BOD , show that a single below-canopy P measurement would have been sufficient to calculate above-canopy WS. In future investigations, such an elaborate measurement system as presented in this paper will no longer be necessary. Since the functional dependence of PPC BOD on WS 2 BOD is strong, PPC BOD was used as a predictor for T mean of trees B1 to B7 and compared to the dependence of T mean on WS 2 BOD . Figure 7b shows the response of the sample trees' stems to WS 2 BOD together with the power law curves. Although B1 to B7 responded slightly different to WS 2 BOD , the functional dependence of T mean on WS 2 BOD is always strong, which is expressed by the high R 2 values, all exceeding R 2 = 0.75. Stem tilt was lowest for the thickest tree (B6), and largest for the tree with the thinnest stem (B7).
The dependence of T mean on PPC BOD exhibits different tree-specific power law curves compared to the curves calculated between T mean and WS 2 BOD (Figure 7c). However, the strength of the relationship is no less than between T mean and WS 2 BOD , as indicated by R 2 > 0.75 (p = 0.00). Thus, with the comparable strength of the functional relationships, it is possible to use ground-based P to approximate the wind conditions above the Scots pine forest that drive stem tilt of B1 to B7 and to therefore model T mean .
In addition to T mean determined by mean wind loads, large stem tilt far from the central tilt tendencies at the upper tail of the T distribution represented by T 98 are also of great interest for the analysis of wind-induced tree response. They are indicators of the trees' maximum resistance to wind loads caused by the turbulent airflow. Therefore, Figure 8 illustrates the potential of using the combination of T mean and P 98 as predictors for modeling T 98 .
The displayed point clouds cluster around the 1:1-lines. The scatter of the points is very low. Measured by the coefficient of determination, the results of the regression-based modeling are convincing for all sample trees (R 2 ≥ 0.88). The simple regression models can reproduce tree-individual differences in T 98 as a function of T mean and P 98 . While the small diameter tree B7 responded most strongly to the extreme wind loads, the large diameter tree B6 shows the smallest T 98 values. This is a response pattern that was already described when discussing T mean . For modeling of T 98 , the relative importance of T mean was 71 ± 6%. The relative importance of P 98 equaled 29 ± 6%. The magnitude of the predictor importance ratio results from the low mean wind and gust speed levels during the study period. At low wind loading, T 98 is closely coupled to T mean . Under these wind conditions, the impact of sporadic, intense short-term wind loads due to high gust speed as approximated by P 98 is not expected.
The distribution of the local wind direction observed above the canopy at To1 and To2 is shown in Figure 9a in the form of wind roses. The wind roses displaying the wind direction distributions were created using the 10-min mean horizontal wind vector components. The wind direction distributions are bimodal but different. Southern and north-northwestern directions along the longitudinal axis of the Upper Rhine valley dominate the distributions. Few values were observed for easterly southeasterly directions, i.e., more across the lateral axis of the Upper Rhine valley.
A similar distribution, but more sharply divided into two sectors, is found from the P propagation direction values (Figure 9b). The wind roses displaying the P propagation direction distributions were created using the 10-min mean values of the time differences of P maxima during the passage of P waves across M1, M3, M15, M17, and M9. The values are more concentrated in the two ranges 135 to 180 • and 315 to 360 • . Almost no values were determined from westerly and easterly directions.
In Figure 9c,e,g, the T 98 shares of B1, B2, and B7 are shown per wind direction sector by radial bar plots. The radial bar plots were chosen instead of the wind rose plots because the directional distribution of T 98 and not T mean is shown as a function of the WD distribution. The radial bar plots illustrate the T 98 distribution as a function of WD.
Depending on WD at To1 and To2, T 98 shares occur across all sectors, with tree-specific differences in magnitude. There is hardly any dependence on wind direction. In addition, the T 98 shares differ per sector. While B2 shows the weakest response to wind loading, the wind-induced response of B7 is most vigorous. The wind-induced motion patterns of B1, B2, and B7 as a function of mean P propagation direction are shown in Figure 9d,f,h. The radial bar plots show the T 98 distribution as a function of PD. The sectoral distribution of the highest T 98 values differs from the distributions shown as a function of WD. The values are more closely distributed within the two directional sectors where the mean P propagation directions occurred most frequently.
At this point, it can only be speculated about the reasons for the differences in the representations of the direction-dependent tree motion. A plausible assumption is that the wind speed measurement is more influenced by the local measurement environment than the air pressure measurement. Taking this reasoning further, it is necessary to plan the position of the airflow measurement even more precisely in future studies of wind-tree interactions than has been the case in the past.
The results obtained for PPC BOD , P 98 , and PD are closely related to WS BOD and WD. This similarity can be used to estimate the temporal dynamics of T mean and T 98 through ground-based air pressure measurements, the assessment of which would otherwise only be possible through more complex, tower-based wind speed measurements above the forest canopy.
Since the functional relationship between air pressure fluctuations measured on the ground and above-canopy wind speed is not due to a known local dependence, it is possible to calibrate the pressure fluctuations and wind-induced tree motion with wind speed measurements made outside of forests, e.g., in the networks of the national meteorological services. For example, the air pressure fluctuations measured in this study can be correlated (R = 0.70) with the nearest wind speed measurement made by the German Meteorological Service at 10 m a.g.l. in the city of Freiburg, 20 km northeast of the measurement site. Thus, we see potential that in the future, the available measuring stations located outside of forests can be used for an improved estimation of wind speed in forested areas.

Conclusions
The presented results show that ground-based measurements of air pressure fluctuations are suitable for estimating the temporal dynamics of mean and extreme tree response in a planted Scots pine forest. The significant advantages of the low-cost air pressure measurements are their ease of installation and the potential for deployment in large numbers over extensive areas. The results demonstrate that one air pressure measurement is sufficient to obtain spatially representative estimates of wind conditions over the forest in the presented measurement grid. This feature of the air pressure measurement can complement and potentially replace more complex tower-based wind measurements in remote forest areas with poorly developed infrastructure and research facilities. If a vertical profile of wind speed measurements is available, at least over a short period of time, then even a parameterization of the mean vertical wind speed profile could be achieved.
To evaluate the air pressure sensors' general suitability, a simplified set up is required to be tested in more complex terrain in future studies, e.g., at existing forest research sites with tall meteorological towers. The air pressure measurements could be used in combination with tree response sensors to capture better the small-scale variability of airflow, which to date remains very difficult to describe and quantify. The proposed approach to use low-cost, ground-based air pressure measurements for airflow approximation opens up opportunities to contribute to improved small-scale storm hazard assessment through the large-scale deployment of air pressure measurement networks located on the ground in forests. It may even be considered that a large-scale measurement network is used to monitor the adaptation of trees to the local wind climate. This could provide new insights into the interactions of wind-induced sway, tree acclimation, air pressure fluctuations, and canopy airflow.