Analysis of Air Pressure Fluctuations and Topsoil Gas Concentrations within a Scots Pine Forest

High-precision differential air pressure measurements were conducted in the below-canopy space of a Scots pine forest and in the forest soil to investigate small air pressure fluctuations and their effect on soil gas flux. In addition to air pressure measurements, tracer gas concentration in the soil and airflow characteristics above and below the canopy were measured. Results suggest that air pressure fluctuations in the frequency range of 0.01 Hz–0.1 Hz are strongly dependent on above-canopy wind speed. While amplitudes of the observed air pressure fluctuations (<10 Pa) increase significantly with increasing above-canopy wind speed, the periods decrease significantly with increasing above-canopy wind speed. These air pressure fluctuations are associated with the pressure-pumping effect in the soil. A pressure-pumping coefficient was defined, which describes the strength of the pressure-pumping effect. During the measurement period, pressure-pumping coefficients up to 0.44 Pa·s −1 were found. The dependence of the pressure-pumping coefficient on mean above-canopy wind speed can be described well with a polynomial fit of second degree. The knowledge of this relation simplifies the quantification of the pressure-pumping effect in a Scots pine forest considerably, since only the mean above-canopy wind speed has to be measured. In addition, empirical modeling revealed that the pressure-pumping coefficient explains the largest fraction of the variance of tracer gas concentration in the topsoil.


Introduction
The occurrence of small air pressure fluctuations (<10 Pa) has been known since the 1960s [1,2].Subsequent research found a relationship between the standard deviation of air pressure fluctuations and airflow variables like the mean wind speed and the friction velocity [3][4][5][6][7][8].In addition, an influence of air pressure fluctuations on soil gas transport was observed [2,[9][10][11].This influence is referred to as turbulence-induced "pressure pumping" and has attracted great interest in recent years [8,[12][13][14][15][16][17][18][19][20].Results from present research suggest that soil gas transport is enhanced up to 100% by air pressure fluctuations induced by airflow, dependent on the investigated soil [18,19,21].This can be an important driver for the soil-atmosphere exchange rates of greenhouse gases.However, there is still need for further in situ quantification of the variables associated with the pressure-pumping effect, such as amplitudes and frequencies of air pressure fluctuations.
The pressure-pumping effect might also influence the soil gas measurements conducted with chamber methods.Studies indicated a modification of soil gas transport under chambers used for soil flux measurements [22][23][24].However, recent results from field studies report both positive and negative correlations between measured soil gas efflux and near-surface turbulence [13,[25][26][27][28][29][30].
The knowledge of pressure-pumping-related pressure fluctuations is still limited, since the measurement of air pressure fluctuations is very challenging.One problem associated with the measurement of air pressure fluctuations is the need for pressure probes that do not add artificial pressure effects, which become apparent at higher wind speeds [31][32][33][34][35]. Another problem is the small amplitudes of air pressure fluctuations.Currently, only differential pressure sensors have the accuracy to measure such small pressure amplitudes.Differential pressure sensors need a reference [5,36,37], which makes simultaneous pressure measurements even more complex.The challenge of simultaneous air pressure fluctuation measurements over larger distances is to develop a system that uses either the same reference or several calibrated references for all involved pressure sensors.Moreover, the measured air pressure signals may be altered due to lengths of tubes [12] and temperature effects [38].Further problems with differential pressure measurements include saturation or damage of the pressure sensors due to diurnal changes of the barometric air pressure [37].
Measurements of atmospheric pressure fluctuations in a forest are rare [5,39,40], and the altering of pressure measurements caused by a forest canopy remains unclear.Some studies focused solely on air pressure fluctuations at the forest floor [2,6], with the advantage that the reference chamber can be buried in the ground to achieve relatively constant conditions in the chamber [41].Other studies tried to retrieve air pressure fluctuations from airflow and air temperature characteristics [7,42], which are easier to measure.
Therefore, this study aims to (1) specify the characteristics of air pressure fluctuations associated with pressure pumping; (2) clarify functional relationships between airflow and air pressure fluctuation characteristics; and (3) verify the effect of air pressure fluctuations on topsoil gas concentrations in a Scots pine forest.Figure 1a shows the normalized plant area density (PAD, in m 2 •m −3 ) profile at the research site.The research site has a lattice tower with a top height of 30 m.It is equipped with a large number of meteorological instruments, including psychrometers to measure dry-bulb and wet-bulb temperature, radiometers to measure the components of the radiation balance, and ultrasonic anemometers [43].Soil temperature and soil moisture were measured at two locations next to the tower at a depth of −0.10 m with Aquaflex probes (Umweltanalytische Produkte GmbH, Ibbenbüren, Germany).Furthermore, a soil temperature ( , in K) was measured with PT100 sensors at the depths −0.01 m ( ), −0.03 m ( ), −0.05 m ( ), −0.10 m ( ), −0.20 m ( ), and −0.40 m ( ).Soil heat flux was measured at two locations ( and , in W•m −2 ) near the tower at a depth of −0.03 m with heat flux plates (HFP01SC, Hukseflux, Delft, The Netherlands).

Measurement Site
According to the World Reference Base (WRB) classification, the soil at the research site is classified as Haplic Regosol [44].Its topsoil is covered by a humus type of mull followed by a silty loam texture over a subsoil consisting partly of sand and gravel [8].

Airflow Measurements
The wind vector components in − ( ), − ( ), and − direction ( ) were measured simultaneously (in m•s −1 ) at five different heights with ultrasonic anemometers (81000VRE, R.M. Young Company, Traverse City, MI, USA) during the period 29 March 2016 to 13 July 2016.The ultrasonic anemometers were mounted on the tower at the heights 2 m ( /ℎ = 0.11), 9 m ( /ℎ = 0.50), 18 m ( /ℎ = 1.00),21 m ( /ℎ = 1.15), and 30 m ( /ℎ = 1.67) above ground level (a.g.l.).The sampling rate of the ultrasonic anemometers was 10 Hz and data were analyzed per 30 min interval.At the measurement site, airflow from northern and southern directions dominated during the measurement period (Figure 1b).To minimize the influence of the tower on the measurements, all ultrasonic anemometers were installed on 1.5 m long supporting booms directed to the west.The maxima of the 30 min mean wind speed during the measurement period at , , , , and were , and , = 6.5 m•s −1 .

Pressure Sensors
Below-canopy and soil air pressures ( , in Pa) were measured with three piezo-resistive differential pressure sensors (GMSD 2.5 MR, Greisinger Electronic GmbH, Regenstauf, Germany) to achieve the highest accuracies currently available (sensitivity 0.1 Pa, accuracy 1%, measurement range −2 hPa to +2.5 hPa).Two sensors were installed at the same heights as the below-canopy airflow measurements ( at , at ).An additional pressure sensor measured the air pressure fluctuations in the soil ( at = −0.03m).The air pressure data were also analyzed per 30 min interval.One inlet of each pressure sensor in the below-canopy space was connected with a tube (inner diameter: 2 mm) to a low-cost pressure head.The pressure head consists of a metal pipe (inner diameter: 2 mm) which was inserted half way into a plastic sphere with 25 holes (outer diameter: 72 mm, hole diameter: 9 mm).The sphere was used to ensure omnidirectionality of the air pressure measurement and to prevent errors due to dynamic pressure effects.The other inlet of the pressure sensor was connected to the reference, which filters all high-frequency air pressure fluctuations (see below).
This design was tested against a state-of-the-art high-performance DigiPort pressure head (DigiPort, Paroscientific Inc., Redmond, WA, USA) used for barometric air pressure measurements at high wind speeds [33].The DigiPort pressure head has been designed to provide accuracy in air pressure measurements of 1 Pa at a wind speed of 5 m•s −1 [45].For all 30 min intervals with mean wind speeds up to 1.5 m•s −1 the measurement values obtained from the DigiPort pressure head correlated well with the measurement values obtained from our pressure head (Pearson correlation coefficient = 0.99 ± 0.01).In the below-canopy space, 30 min mean values of the wind speed never exceeded 1.5 m•s −1 .Therefore, it was believed that the air pressure measurements below the forest canopy were not disturbed by dynamic pressure errors.

Pressure Reference
The pressure reference consists of two nested volumes: a glass jar (volume: 0.225 l inside an insulating box of Styrofoam (volume: 4.2 l) which is in exchange with the atmosphere (Figure 2).The inlet of the pressure sensor for the reference is connected to the glass jar.Air exchange between the glass jar and the Styrofoam box is possible through a capillary, which acts as a pressure equalizer.The capillary has a length = 0.15 m and inner diameter = 2, = 0.124 mm, with being the radius of the capillary.To avoid temperature effects on the air pressure measurement, the temperature inside the reference chamber should only change slowly over time [37].Therefore, thermal packs were put inside the box to increase thermal inertia, and the insulating box was wrapped in a reflecting blanket to minimize heating up due to insulation.In addition, silica gel inside the box was used to maintain a low humidity and to avoid blocking the capillary with water droplets.The high-frequency cutoff of the pressure reference was determined theoretically as follows.The volumetric flux of air (in m 3 •s −1 ) with dynamic viscosity (in kg•m −1 •s −1 ) through the capillary caused by the difference between the air pressure in the reference chamber (in Pa) and the air pressure of the atmosphere (in Pa) is [46]: With the ideal gas law, constant (in m 3 ), and assuming constant air temperature (in K), the resulting change of amount of substance (in mol) in the chamber per change of time ′ can be written as: The universal gas constant is = 8.3144598 J • mol • K [47].With the molar volume of air (in m 3 •mol −1 ) Equations ( 1) and ( 2) can be equated: Solving this differential equation from time = 0 to time with (0) = and ( ) = yields: The time (in s), when the difference between the air pressure in the reference chamber and the air pressure outside has been attenuated to 1/ times the initial value is: The time constant along with the values for and used for the calculation are listed in Table 1.After a time of 3 , the pressure difference between the reference chamber and the outside was smaller than 5% of the initial value.This value was used to calculate the frequency cutoff .Thus, the reference acts as a low-pass filter with = 1/3 ≈ 4 × 10 Hz.The filter behavior was verified by evaluating the quotient of the Fourier spectral energy density of the filtered pressure signal and the Fourier spectral energy density of the unfiltered pressure signal.The resulting filter behavior yielded a cutoff frequency of 5 × 10 Hz, i.e., about 33 min, and was therefore in good agreement with the theoretical considerations.

Soil Gas Measurements
To investigate the effect of air pressure fluctuations on soil gas concentration, a state-of-the-art gas concentration measurement system was used [49].The measurement system consists of a gas sampling pole, a tracer gas injection device, a valve system and a microgas chromatograph (3000 Micro GC Gas 176 Analyzer, Inficon GmbH, Cologne, Germany).The gas sampling pole has a length of 50 cm and features gas sampling membranes at seven depths (from +0.01 m to −0.41 m depth).At a soil depth of −0.21 m, a tracer gas was continuously injected.In this study, helium was used as tracer gas, since it is inert and has a low solubility in water [49].The soil air was sampled through the membranes at several depths and transferred to the microgas chromatograph, which analyzed the soil air and determined the helium concentration.Thus, the helium concentration profile could be monitored.Assuming molecular diffusion as the main soil gas transport process, a steadystate profile of helium concentration is reached.A deviation of the soil helium profile from the steadystate profile indicated an external effect on soil gas transport.
The gas sampling pole was installed approximately 10 m southwest of the tower in a hole drilled in the soil.Helium was continuously injected at a constant rate into the soil during the period 5 July 2016 to 13 July 2016.Helium concentrations (He, in ppm) were measured at the five depths of −0.01 m ( ), −0.06 m ( ), −0.11 m ( ), −0.31 m ( ), and −0.41 m ( ), and the helium concentrations measured at −0.01 m were used to study the effect of air pressure fluctuations on topsoil gas concentrations.

Airflow Data
The half-hourly time series of wind vector data were despiked and a double rotation was applied [50].In addition, the data of heights -were rotated in the mean wind direction at , ensuring that for the time series of all heights the same coordinate system was used.Data of 30 min intervals when wind from eastern directions prevailed were excluded from the analysis (<10% of all 30 min intervals) to rule out influences from the measurement tower.
To distinguish between different atmospheric stabilities, the Obukhov length (in m) was calculated [51,52]: with the virtual potential temperature (in K), the friction velocity * = + / (in m•s −1 ), the von Kármán constant = 0.4 and the gravitational acceleration = 9.81 m•s −1 .The canopytop stability parameter = / was calculated for every 30 min interval [50].All available halfhourly data sets were then assigned to the six stability classes listed in Table 2.For every 30 min interval the mean wind speed at canopy height ( ) and the friction velocity at canopy height ( * ) were calculated.
Five exchange regimes (C1 − C5) were defined based on the correlation of coherent structures derived from momentum flux analysis over all measurement heights [50].Exchange regime C1 is associated with uncoupled momentum flux above and below the forest canopy, while exchange regime C5 is associated with fully coupled momentum flux from to .That means, during times specified as C1, no momentum is transferred from downwards.On the other hand, during times specified as C5, momentum is transferred from down to the forest floor at [50].

Air Pressure Data
All pressure sensors were sampled with a rate of 2.0 Hz.Every 30 min interval influenced by precipitation were excluded from analysis.For each remaining interval the air pressure fluctuations were investigated separately for three different frequency ranges: high frequency range ( , 0.1 Hz -1.0 Hz), medium frequency range ( , 0.01 Hz -0.1 Hz) and low frequency range ( , 0.001 Hz -0.01 Hz).A fourth-order infinite impulse response (IIR) filter was used as a band-pass filter to separate the air pressure signal into separate frequency ranges.The air pressure fluctuations were described statistically by calculating the arithmetic means ̅ and standard deviations .
Mean amplitudes of the air pressure fluctuations were determined through the upper envelope and lower envelope of the air pressure fluctuations: the mean peak-to-peak amplitude was defined as ̂= − .The dominating periods of the air pressure fluctuations (in s) per 30 min interval were determined by investigating the normalized wavelet variance spectrum [53][54][55].The wavelet coefficients were calculated for 100 scales using the Morlet wavelet, corresponding to up to 103 s.Then, the normalized wavelet variance spectrum was calculated.The scale of the first peak in the spectrum is associated with the dominating wavelet scale of the air pressure fluctuations.The Fourier-equivalent frequency (in Hz) of a wavelet scale can be calculated from the center frequency of the used wavelet (in cycles per unit time) and the sampling period (in s) [56]: Thus, the dominating period is = 1/ .The correlation and time lags between two air pressure datasets and was determined by calculating the Pearson correlation coefficient ( , ) and the cross-correlation between and .The Wilcoxon rank sum (WRS) test was used to evaluate if differences between air pressure fluctuation properties for different stability conditions and different exchange regimes were significant (level of significance = 0.05) [57].To describe the strength of the air pressure fluctuations, the pumping coefficient (in Pa•s −1 ) was defined as follows.First, the absolute slope |∆ /∆ | (with ∆ = 1/ = 0.5 s) between two subsequent measurement points is calculated for every measurement point in the 30 min interval.Then, was obtained by calculating the mean of all absolute slopes in the 30 min interval: Therefore, the pumping coefficient describes the mean change in air pressure fluctuation amplitude per second in a 30 min interval.

Empirical Modeling of Topsoil Helium Concentrations
To identify the factors that have an effect on the temporal evolution of , the ensemble learning method random forests (RF) implemented in the Matlab ® Statistics and Machine Learning Toolbox (The MathWorks Inc., Natick, MA, USA, Release 2016a) was applied.Using bootstrap samples, RF combined binary decision trees, which were built on 66% of the available helium concentration data (D1).The remaining 34% of the helium concentration data (D2) were used for model evaluation and prediction of the topsoil helium concentration.
The following predictor variables were considered as potentially informative input to the RF model: dry-bulb and wet-bulb air temperature above and below the canopy, components of the radiation balance at canopy top, soil temperatures, soil moistures, and soil heat fluxes.
Prior to model building, (1) all variables were detrended using a sixth degree polynomial; and (2) the strength of collinearity among the predictor variables was assessed by the variance inflation factor ( < 2) and Belsley collinearity diagnostics [58].If collinearity was detected among two predictor variables, the predictor variable contributing less to the final RF model accuracy was excluded from further model building.The relative contribution of the selected predictor variables to the final RF model output was evaluated by the predictor importance (PI, in %) quantified for D2.
The PI values were used to identify important predictor variables which strongly impact the RF model accuracy after being randomly permuted.After testing for collinearity, various combinations of predictor variables were evaluated for their power to predict .Starting with one predictor variable, further predictor variables were sequentially added to the RF model and retained when the model error decreased.The coefficient of determination ( ), the mean squared error ( ) and the mean absolute error ( ) were used to assess the accuracy with which the RF model simulated the measured concentration values.

Frequency Characteristics of Air Pressure Fluctuations
Figure 3 shows the raw air pressure signals measured below the canopy ( and ) and in the soil ( ) over a six-hour period (22 April 2016, 02:00-08:00) with increasing above-canopy wind speed .The low-frequency fluctuations are visible in the air pressure signals of all measurement heights.However, as increases, the air pressure signals contain fluctuations of higher frequencies that do not occur at low values.Other than results from a previous study [39], the fluctuations did not show a deterministic relationship with , , , , momentum flux ( ′ ′), or air temperature ( ) associated with coherent structures [59].Therefore, subsequent analysis focuses on 30 min statistics.The dependence of the air pressure amplitudes at ( ̂ ) on above-canopy mean wind speed ( ) was used to determine whether the amplitudes of the fluctuations in a frequency range change as a function of wind speed (Figure 4a).In addition, the correlation of the fluctuations between and , and and was determined by analysis of ( , ) and ( , ) with respect to (Figure 4b).Combination of both information yielded the frequency range of interest.While no clear dependence of the mean amplitudes of the low frequency air pressure fluctuations ( ̂ , ) on mean above-canopy wind speed was found, ̂ increased with increasing in the medium and high frequency range (Figure 4a).Therefore, low frequency fluctuations were contained in the air pressure signal independent of wind speed, and were not investigated further.The fluctuations in the medium frequency range (0.01 Hz -0.1 Hz) showed strong correlation below the canopy, while the high frequency fluctuations were mostly uncorrelated (Figure 4b,c).Moreover, the amplitudes of the air pressure fluctuations in the medium frequency range ( ̂ , ) were around five times larger than the amplitudes of the air pressure fluctuations in the high frequency range ( ̂ , ).Based on these results, subsequent analysis of air pressure fluctuations focuses on the medium frequency range.This result is also in agreement with a previous study that suggested air pressure fluctuations with frequencies <0.1 Hz are important [12].

Mean Wind Speed at Canopy Height
Since strongly correlated air pressure fluctuations were measured at to , dependences are only shown for one air pressure signal ( ).While no clear dependence of , on the mean wind speed measured at the same height was found (Figure 5a), , showed a quadratic dependence on (Figure 5b) in accordance with a previous study [7].This emphasizes that these air pressure fluctuations are not produced locally.In contrast to other studies, no clear linear dependence of on the mean wind speed [5,6], or an exponential dependence of on the friction velocity * [8] was found.A similar dependence on was found for the amplitudes ̂ , (Figure 6a).Moreover, a decrease of , with increasing was found (Figure 6b).At the highest values of , the values decreased to 20 s. Results from WRS test showed that ̂ increased significantly with increasing , but that there were no significant differences between ̂ at , and .Moreover, results from WRS test also showed a significant decrease of mean periods with increasing .

Atmospheric Stability
The dependence of the amplitudes ̂ and periods was analyzed under different stability conditions (Table 3).Results from the WRS test showed that ̂ was significantly larger under near-neutral conditions than ̂ under all other stability conditions.This can be attributed to higher wind speeds under near-neutral conditions, which are assumed to be the main generator of the air pressure fluctuations.Since dependences of on * had been found in previous studies [19], and * is related to ′ ′, a dependence of on ′ ′ was expected.In addition, / ( : volumetric mass density) and ′ ′ have the same unit.Therefore, differences between the air pressure fluctuations under the momentum flux-derived exchange regimes C1 (no coupling, 83% of all 30 min intervals) and C5 (completely coupled, 12% of all 30 min intervals) were investigated.Table 4 lists the means for and ̂ under the different exchange regimes.Results from WRS test showed no significant differences between the two exchange regimes.Since mean values for under C1 and C5 are similar, these results emphasize that the occurrence of air pressure fluctuations in the medium frequency range are dependent on .

Pressure-Pumping Coefficient
Since pressure pumping is especially important for soil gas transport, the dependence of the pressure-pumping coefficient in the soil ( , ) on the mean above-canopy wind speed was investigated.Values of , vary between 0 during intervals of low and 0.44 Pa•s −1 during intervals of high (Figure 7).Results from WRS test showed a significant increase of  In contrast, the pressure-pumping coefficient calculated for the low-frequency range The threshold value of = 1.5 m•s −1 for the occurrence of air pressure fluctuations in the medium-frequency range is also a reasonable threshold value for pressure pumping in the soil for this measurement site.

Influence on Soil Gas Transport
The mean helium concentration at the injection depth for the period 5 July 2016 to 13 July 2016 (〈 〉) was 388 ± 69 ppm.Towards the soil surface, the mean helium concentration decreased to 〈 〉 = 115 ± 11 ppm and 〈 〉 = 28 ± 5 ppm.
The evaluation of collinearity among the most important predictor variables revealed that and , , and and (mean of and ) were collinear ( > 3.5) while , did not exhibit any collinearity with other important predictor variables.Therefore, was excluded from building the final RF model.
The predictor variable combination that provided the best discrimination between RF modeling results for D2 data and ( = 0.72) included , , , , and .Results from PI evaluation demonstrate that the relative impact of , on the predictive accuracy of the RF model was greatest ( = 38%).
The second most important predictor variable was ( = 28%), followed by ( = 22%), and ( = 12%) (Figure 8).Considering the results from the collinearity analysis, the most accurate RF model including as a predictor variable yielded = 0.64 .This finding indicates that , is better suited for empirical modelling of than .The RF model results indicate that large proportions of variance can be explained by at least two processes.On the one hand, variance is strongly affected by the soil heat flux.On the other hand, a substantial proportion of variance can be explained by wind-induced air pressure fluctuations.

Conclusions
The results of the present study reveal a strong positive, quadratic relationship between 30 min mean values of wind speed at canopy top and pressure-pumping coefficient in the topsoil.This simplifies the quantification of the pressure-pumping effect considerably, since only wind speed has to be measured.However, it is still unclear whether this relationship is site-dependent or universally applicable.The pressure pumping may be affected by the canopy type, canopy height, and foliage.Therefore, the presented methodology has to be conducted at other sites with different canopy characteristics in the future.
The pressure-pumping coefficient was calculated based on air pressure fluctuations occurring in the frequency range 0.01 Hz -0.1 Hz.It is a measure for the half-hourly intensity of air pressure fluctuations and describes the mean change in pressure per second.Empirical modeling of helium concentration demonstrated that the pressure-pumping coefficient is an important predictor for changes in the topsoil gas concentration, and thus, an important factor for soil gas transport.
Knowledge of the half-hourly amplitudes and frequencies of the air pressure fluctuations can serve as a basis for the investigation of the pressure-pumping effect in the laboratory.By reproducing air pressure fluctuations based on the findings of this study, laboratory studies allow for a clear quantification of air pressure fluctuations on topsoil gas concentrations.
Air pressure fluctuations and airflow characteristics were measured at the forest research site Hartheim, operated by the Chair of Environmental Meteorology of the University of Freiburg.It is located approximately 25 km southwest of Freiburg in the flat southern Upper Rhine Valley (47°56′04″N, 7°36′02″E, 201 m above sea level).The forest at the research site is a single-layered plantation of Scots pines (Pinus sylvestris L.) and was established in the 1960s.The mean tree height (ℎ) in the year 2016 was approximately 18 m, and the mean stand density was 580 trees•ha −1 .The mean plant area index is 1.5.

Figure 1 .
Figure 1.(a) Normalized vertical profile of the plant area density (PAD), and (b) wind rose at height at the measurement site.

Figure 2 .
Figure 2. Schematic picture of the pressure sensor and the pressure reference.

Figure 4 .
Figure 4. Dependence of (a) the amplitudes of the air pressure signal measured at height ( ̂ ); (b) the absolute correlation coefficient between the air pressure signals at heights and (| ( , )|); and (c) the absolute correlation coefficient between the air pressure signals at heights and (| ( , )|) on mean above-canopy wind speed ( ) for different frequency ranges (low, medium, high).

Figure 5 .
Figure 5. Dependences of the standard deviation of the air pressure signal in medium-frequency range ( , ) (a) on the mean below-canopy wind speed ( ); and (b) on the mean above-canopy wind speed ( ).

Figure 6 .
Figure 6.Dependence of (a) amplitudes of the medium frequency part of the air pressure signal measured at height ( ̂ , ); and (b) dominating periods of the medium frequency part of the air pressure signal measured at height ( , ) on the mean above-canopy wind speed ( ).
of , on for > 0.75 m•s −1 can be described by a second-degree polynomial ( = 0.97)

Figure 7 .
Figure 7. Dependence of the pressure-pumping coefficient calculated for the medium-frequency range at height ( , ) on mean above-canopy wind speed ( ).

Table 1 .
Values for and calculated and (time constant) of air for = 298 K (25 °C) and = 1013.25 hPa.

Table 2 .
Classification of atmospheric stability according to the stability parameter .