Articles | Volume 19, issue 8
Biogeosciences, 19, 2235–2243, 2022
Biogeosciences, 19, 2235–2243, 2022
Research article
26 Apr 2022
Research article | 26 Apr 2022

Changes of the aerodynamic characteristics of a flux site after an extensive windthrow

Changes of the aerodynamic characteristics of a flux site after an extensive windthrow
Bruna R. F. Oliveira1, Jan J. Keizer1, and Thomas Foken2 Bruna R. F. Oliveira et al.
  • 1Earth surface processes team, Centre for Environmental and Marine Studies (CESAM), Department of Environment and Planning, University of Aveiro, Aveiro, Portugal
  • 2Bayreuth Center of Ecology and Environmental Research (BayCEER), University of Bayreuth, Bayreuth, Germany

Correspondence: Thomas Foken (


A maritime pine plantation in central Portugal that has been continuously monitored using the eddy-covariance technique for carbon fluxes since a wildfire in 2017 was significantly affected by two storms during December 2019 that resulted in a large-scale windthrow. This study analyses the impacts of this windthrow on the aerodynamic characteristics of zero-plane displacement and roughness length and, ultimately, their implications for the turbulent fluxes. The turbulent fluxes were only affected to a minor degree by the windthrow, but the footprint area of the flux tower changed markedly so that the target area of the measurements had to be redetermined.

1 Introduction

Heterogeneous surfaces have an influence on turbulent energy and mass fluxes (Stoy et al., 2013). For example, increased fluxes have been found at forest edges (Klaassen et al., 2002). Large-eddy simulations have suggested that fluxes are at their maximum at a distance of about 10 times the canopy height from the forest edge (Kanani-Sühring and Raasch, 2015; Dupont and Brunet, 2009). An increase in flux could also be associated with an increase in stand heterogeneity (Foken et al., 2021). An intrinsic limitation of the cited studies is that they cannot find a suitable measure of heterogeneity that is correlated with turbulent fluxes.

Recent studies on the role of spatial heterogeneity in flux measurements, as referenced above, either address its influence on the closure of the energy balance at the earth's surface (Mauder et al., 2020) or its relevance for footprints (Göckede et al., 2008; Chu et al., 2021). The typical aerodynamic characteristics such as zero-plane displacement and roughness length – which change following a windthrow – have received little attention. The foundations for the use of aerodynamic characteristics stem from aerodynamic studies in wind tunnels, whose findings were then adopted for meteorology (Prandtl, 1932) and led to the introduction of zero-plane displacement (Paeschke, 1937). These characteristics were a significant research focus from the 1950s to 1980s, when fluxes were determined from wind and temperature profiles measured at more than two heights and sometimes involving elaborate procedures (Nieuwstadt, 1978; Marquardt, 1983). All relevant textbooks of today address the basic principles underpinning aerodynamic characteristics (Stull, 1988; Garratt, 1992; Arya, 2001; Foken, 2017). The fundamental problem, however, is that zero-plane displacement and roughness length can only be determined if both profile and flux measurements are available. In the past, flux measurements were typically missing so that fluxes were estimated using complicated approximation approaches (Kader and Perepelkin, 1984). In recent times, by contrast, profile measurements have been typically missing so that they are replaced by reasoned inferences on profile characteristics. In a strict sense, zero-plane displacement and roughness length can only be assumed uniform across homogeneous surfaces. Recent attempts have tried to incorporate stand structure into empirical relationships for determining zero-plane displacement and roughness length (Nakai et al., 2008; Raupach, 1994), especially by deriving stand structure from remote sensing data, based on the paper by Thom (1971). Maurer et al. (2015) carried out a large-eddy simulation to compare various approaches to estimating aerodynamic characteristics and confirmed a near-linear relationship between canopy height and zero-plane displacement. Spatial heterogeneity and canopy structure can be addressed relatively easily by data processing, except if they change through time. Such temporal changes then lead to erroneous flux calculations. The present study wants to demonstrate this for an abrupt, dramatic change in canopy structure due to a windthrow of trees that had been killed by a wildfire some 2 years earlier.

The post-wildfire flux site in Vila de Rei, central Portugal (Oliveira et al., 2021), offered an opportunity to study the impacts of an abrupt change in aerodynamic characteristics following a windthrow caused by two consecutive storms, Elsa and Fabien, without an apparent concomitant change in stand heterogeneity, virtually like a laboratory experiment. Namely, the storms caused an extensive windthrow of the – dead – burnt maritime pine trees between 19 and 21 December 2019, 28 months after the wildfire, while the tumbled trees remained on the ground afterwards (Fig. 1). Still, the two storms did not throw over the – living – eucalypt trees, i.e., neither the individual eucalypt trees along the pine plantation nor those of the eucalypt plantations adjacent to it. These individual eucalypt trees expectedly have an influence on the aerodynamic characteristics (Jegede and Foken, 1999). This prompted us to examine the changes in the roughness length assuming a linear dependence of zero-plane displacement on stand height before and after windthrow. Furthermore, we decided to analyze if the changes in these two aerodynamic characteristics also affected carbon dioxide fluxes. This question, however, is difficult to answer in this particular case. Roughly 2 years after the wildfire, the pine ecosystem was still in its initial phase of post-fire vegetation recovery so that an increase in CO2 uptake is to be expected between the period before and the period after the windthrow, regardless of the aerodynamic changes. Nevertheless, a change in zero-plane displacement can have a significant impact on the Obukhov–Lettau stability parameter (z-d)/L (z: measurement height, d: zero-plane displacement, L: Obukhov length; Foken and Börngen, 2021), which, in turn, is crucial to the correction and assessment of carbon dioxide fluxes.

2 Materials and methods

2.1 Measurement site and measurements

The measurement area, instrumentation and data processing were comprehensively described by Oliveira et al. (2021), so only the details essential for this study are given below.

The study area is located 8 km to the southeast of the municipality of Vila de Rei, 3937 N, 0806 W, in a Mediterranean climate zone. The wildfire affecting the area burned 1250 ha of woodlands. The measurement site included a plateau of sedimentary sandstone deposits, located at an elevation of 250 m a.s.l. and with slopes of up to 5 over an extension of roughly 10 ha. The crowns of the bulk of the burnt maritime pine (Pinus pinaster Ait.) were fully consumed by the fire, but their trunks of approximately 8 m height remained standing. This canopy height and a zero-plane displacement of 3.8 m were used in all calculations before the windthrow in Oliveira et al. (2021). While the eucalypt plantations were hardly affected by the windthrow, the maritime pine area following the windthrow was a mixture of dead pine trunks fallen on the soil surface or on the recovering vegetation with an estimated canopy height of 2–3 m. The vegetation mainly consisted of shrubs, locally intermixed with 2–3-year-old pine seedlings and a few individual resprouting eucalypt treelets. The localized patches of burnt eucalypt trees and stands had a canopy height of approximately 4 m.

In the relatively open part of the pine area, a 12 m high slim tower was installed at the end of September 2017 and equipped with an eddy-covariance system at 11.8 m, including a CSAT3 sonic anemometer (Campbell Scientific Inc.) and an LI-7500A gas analyzer (LI-COR Biosciences; see Fig. 1). The fluxes of momentum, sensible heat, latent heat and carbon dioxide were analyzed with the eddy-covariance method (Aubinet et al., 2012). The data of the eddy-covariance system were calculated with the Campbell Scientific Inc. EasyFlux DL software for a quick inspection in the field, while all further calculations were done with the software package TK3 (Mauder and Foken, 2015). The processing of the turbulence data is described in full detail in an extensive supplement ( in Oliveira et al. (2021). Special note should be made of the use of the double rotation (Kaimal and Finnigan, 1994) and that no gap filling was applied. The analysis of the aerodynamic characteristics was done for two 9-month periods from 22 December 2018 to 30 September 2019 and from 22 December 2019 to 30 September 2020, only using data from neutral stratification conditions (6846 and 5864 half-hourly data sets for the 2018–2019 and 2019–2020 periods, respectively). By contrast, the analysis of the CO2 flux measurements was limited to the periods from May to August in 2019 and in 2020, involving 5832 and 5759 half-hourly data sets, respectively. Only the 2020 data set was used for assessing the influence of different aerodynamic characteristics on CO2 fluxes because the vegetation cover was markedly different in the 2 years.

Figure 1A 12 m slim tower with an eddy-covariance system at the burnt Pinus pinaster plantation (a: photograph by Jan J. Keizer, 22 September 2017) and the plantation after the windthrow (b: photograph by Jan J. Keizer, 20 June 2020).


2.2 Aerodynamic characteristics

The starting point for the analysis of the aerodynamic conditions was the logarithmic wind profile. In order to exclude influences due to the stability of the stratification, the analysis was limited to neutral cases (-0.05z/L0.1). The profile equation contains the measured quantities wind speed u and friction velocity u and the two unknowns of zero-plane displacement d and roughness length z0 (Arya, 2001; Foken, 2017; Stull, 1988):

(1) u z - d u = 1 κ ln z - d z 0 ,

where the von Kármán constant is κ=0.40. If wind speeds are measured at different heights, both unknowns can be determined iteratively. In the present case, however, measurements were only available at 11.8 m height, so one of the two parameters must be estimated. It is common in such cases to estimate the zero-plane displacement, in particular, as being equal to 2/3 of the stand height (d=0.666 zc), as is implemented in calculation programs for eddy-covariance measurements. In reality, however, this multiplication factor varies in the range from 0.5 to 0.8, depending on stand structure and, hence, on plant development over the course of a year (Maurer et al., 2015) as well as wind speed (Marunitsch, 1971). Often the value attributed to this factor depends strongly on the experience of the observer. Oliveira et al. (2021) used zc=7.6 m as canopy height, d=3.8 m =0.5zc as zero-plane displacement and z0=0.4 m as roughness length for the period before the windthrow. This determination was made because of the sparse canopy with charred trunks without leaves. For the investigations after the windthrow, a canopy height of zc=2.7 m and a zero-plane displacement of d=1.8 m =0.666zc were assumed. This relationship between canopy height and zero-plane displacement corresponds to the classical approach from hydrodynamics. The application of this approach is comprehensively described in Foken (2017, Sect. 3.1.1 and 3.1.2). More recent canopy-structure-dependent approaches (Nakai et al., 2008; Raupach, 1994) lack input parameters for the highly disturbed surface.

The roughness length and the dimensionless wind profile uz-du are typically used as measures of the roughness of surface and the friction on the surface. In addition, it is useful to determine the so-called integral turbulence characteristic from the standard deviation of the vertical wind velocity and friction velocity σw/u, which has a value of about 1.25 in the neutral case (Foken, 2017; Garratt, 1992) and of 1.1 for measurements close above the canopy (Finnigan et al., 2009). It can, however, attain higher values under the influence of high roughness (Foken and Leclerc, 2004).

Roughness length can be determined by two methods that are nearly independent. The first method is through Eq. (1) for a given zd:

(2) z 0 = z - d e κ u z - d / u .

The second method is based on the following relation (Panofsky, 1984), with σw/u=1.25:

(3) σ w u z - d = 1.25 κ ln z - d z 0 .

From Eq. (3), roughness length is then follows as

(4) z 0 = z - d e 1.25 κ σ w / u z - d .

2.3 Influence of the changes in aerodynamic characteristics on carbon dioxide fluxes

All relevant software packages for the calculation of the carbon dioxide fluxes use, besides measurement height, canopy height and zero-plane displacement as input parameters. These parameters are mainly needed to determine the Obukhov–Lettau stability parameter z-dL, with the Obukhov length of

(5) L = - u 3 κ g T Q H ρ c p ,

where gravitational acceleration is g, temperature is T, sensible heat flux is QH, air density is ρ and the specific heat for constant pressure is cp.

The Obukhov–Lettau stability parameter, in turn, is crucial for the (i) spectral correction in the high-frequency spectrum (Garratt et al., 2020; Moore, 1986), (ii) stability-dependent turbulence characteristics in quality control (Foken and Wichura, 1996) and (iii) determination of the footprint (Leclerc and Foken, 2014). Therefore, this study also addresses the implications of the windthrow-induced changes in aerodynamics for these three aspects.

The model spectra for frequency correction (e.g., Kaimal et al., 1972) are stability dependent (i), and the integral characteristics used for data quality analysis (if not limited to the neutral case as in Eq. 3) are also stability dependent (ii); e.g., according to Panofsky et al. (1977),

(6) σ w u = 1.3 1 - 2 z - d L 1 / 3 .

To determine the source area of CO2 flux measurements, footprint models (iii) are used, and the input parameters are mainly wind speed, roughness length and stability. In the present case, the widely used model according to Kormann and Meixner (2001) was applied.

Further explanations of the corresponding equations and models shall be omitted, since they are described comprehensively in the literature (Foken et al., 2012; Mauder et al., 2021). Furthermore, supplementary information was provided in the prior publication (Oliveira et al., 2021).

3 Results

3.1 Aerodynamic characteristics

A climatology of the wind field would require a minimum of 10 years of data, even if a period of 30 years is the standard. In other words, changes in the wind field before and after the windthrow could simply reflect the inter-annual variation. Nevertheless, it was striking that the median wind speed was higher after than before the windthrow in 8 of the 12 wind sectors (Fig. 2), suggesting that the windthrow provoked a generalized decrease in zero-plane displacement.

Figure 2Median wind speed from the 12 wind direction sectors for the periods of May to August 2019 and 2020.


According to Table 1, the ratio uz-du before and after the windthrow changed roughly from 5–6 to 8–9, corresponding to an increase of about 60 %. In the northeastern sector, the change was smaller, but the ratio was already comparatively high before the windthrow. Assuming a value for zd of 8 m before and 10 m after the windthrow, the effective measuring heights were 8 and 10 m, respectively.

The roughness length showed a decrease after the windthrow (Table 2). Especially after the windthrow, the roughness lengths agreed relatively well between the two methods. Generally, the roughness length was greatest in the southeastern sector.

The ratio σw/u showed the expected values (Table 3). The values were slightly lower before than after the windthrow, in line with the smaller distance between measurement height and the top of the canopy. The northeastern and the western–northwestern sectors revealed comparatively large roughnesses both before and after the windthrow. The wind sectors highlighted in Tables 1–3 are illustrated in Fig. 3.

Table 1The median of the ratio uz-du before (z-d=8 m) and after (z-d=10 m) the windthrow for a nearly 9-month period (22 December 2018 to 30 September 2019 and 22 December 2019 to 30 September 2020). For an explanation of the sectors in bold, see the text.

Download Print Version | Download XLSX

Table 2The median roughness lengths z0 in meters before (z-d=8 m) and after (z-d=10 m) the windthrow for the nearly 9-month periods, computed following Eqs. (2) and (4). For an explanation of the sectors in bold and italics, see the text.

Download Print Version | Download XLSX

Table 3The median ratio σw/u before (z-d=8 m) and after (z-d=10 m) the windthrow for the nearly 9-month periods. For an explanation of the sectors in bold and values, see the text.

Download Print Version | Download XLSX

Figure 3UAS-based (unmanned aerial system) ortho-photomap with the most disturbed wind sectors. Image from 23 February 2021 by Bruna R. F. Oliveira.


3.2 Influence on carbon dioxide fluxes

Carbon dioxide fluxes were evaluated for the months of May to August 2019 and 2020 to investigate the impacts of the windthrow-induced changes in the stability parameter. This was done without gap filling. Of course, the comparison between before- and after-windthrow fluxes is not straightforward because of the differences in weather conditions as well as in post-fire ecosystem recovery between the two periods. Therefore, the spectral correction and data quality analysis were only done for the 2020 data set, assuming two different zero-plane displacements.

Examination of the ratio σw/u showed minor and not relevant differences for z-d=8 m before the windthrow and z-d=10 m after the windthrow. The median differed by only 0.4 % and, hence, did not affect the quality flagging. Also, the difference with and without spectral correction was reduced, even if slightly over 1 %.

A similar result was obtained for the CO2 and net ecosystem exchange (NEE) fluxes. As shown in Fig. 4, the differences in NEE fluxes with and without stratification-dependent spectral correction was about 1 % for the two different effective measurement heights of 8 and 10 m. The median values of the spectrally corrected fluxes differed less than 0.5 % between both measurement heights. Furthermore, the scatter around these median values was reduced.

Figure 4Comparison of the net ecosystem exchange (NEE) with and without spectral correction for z-d=8 m (Moore T for used vs. F for not used) and z-d=10 m (Moore T vs. F) and with spectral correction for both effective measuring heights (Moore T, z-d=8 m vs. 10 m) for the data set from May to September 2020.


More pronounced implications were expected with respect to the size of the source area. Therefore, the target area of burnt pine woodland with down wood of charred trunks was calculated with the footprint model (Kormann and Meixner, 2001), analogous to Oliveira et al. (2021). For the 2020 period, this was done for both values of z-d=8 and z-d=10 and the assigned roughness lengths of z0=0.7 m and z0=0.3 m before and after the windthrow, respectively (Fig. 5). For comparison, Fig. 5 also shows the results for the 2019 period using z-d=8 m and z0=0.4 m (Oliveira et al., 2021) as well as z0=0.7 m (this study). The differences between the results obtained using z-d=8 m and z0=0.7 m for both periods were mainly due to the differences in wind regimes during the 2019 and 2020 periods, as shown in Fig. 2. The differences between the results obtained for the 2020 period using the different aerodynamic parameter values of before and after the windthrow were relevant. The footprint was substantially larger for the post- than pre-windthrow parameter values, while the number of cases with the target area in the footprint was clearly smaller. Calculating the 2019 data with different z0 values showed no substantial differences in the footprint.

The enlargement of the footprint can be prevented by reducing the measuring height. To simulate this case, Fig. 5 shows the footprint for a measurement height of 9.8 m with a reduction in the wind speed by 10 % (z-d=8 m, z0=0.3 m). The change in measurement height leads to an almost identical distribution of the footprint as before the windthrow in 2019.

Figure 5Percentage of the target area (burnt maritime pine woodland) classes in the footprint area of the eddy-covariance measurements from May to August 2019 (before windthrow) using z-d=8 m and z0=0.4 m (following Oliveira et al., 2021) as well as using z-d=8 m and z0=0.7 m (this study) and from May to August 2020 using z-d=8 m (z0=0.7 m) and z-d=10 m (z0=0.3 m). Furthermore, the fictive footprint for a measurement height of 9.8 m with a reduction in the wind speed by 10 % (z-d=8 m, z0=0.3 m) is shown.


4 Discussion

4.1 Aerodynamic characteristics

Before the windthrow, the wind profile was lifted up by the displacement height. After the windthrow, the roughness of the pine stands was determined by the vegetation that had recovered after the fire (mainly consisting of shrubs, locally intermixed with 2–3-year-old pine seedlings) and the dead pine trunks that had fallen on top of this vegetation or on the soil surface. The two determination methods produced consistent results. Worth stressing is that the two methods (Eqs. 2 and 4) are not completely independent because they can be transformed into each other.

The 120–150 sector had a comparatively high roughness both before and after the windthrow. This was probably due to the greater slope angle or the influence of the tower. By contrast, the 30–60 and 270–300 sectors were affected by additional mechanical turbulence, as was also found by Oliveira et al. (2021) for the first post-fire year. This is illustrated in Fig. 6, using Eq. (1), and shows the dimensionless wind profile as a function of zd for different roughness lengths. Also the median values for the 30 sectors were included in the plot. Because of the physical correlation between all quantities, the plot cannot be used to determine the zero-plane displacement (as referred earlier), but it does indicate the sectors where the assumed parameter values seem appropriate.

Figure 6Normalized wind profile dependent on the displacement height given as zd (z=11.8 m) and the roughness length z0. Median values for wind sectors are plotted for z-d=8 m (before the wind break) and z-d=10 m (after the wind break).


The ratio σw/u was largely constant at neutral stratification and only revealed a small stability dependence. This parameter can be used to detect obstacles at larger distances (Foken and Leclerc, 2004) or near the anemometer, for example. Even single standing trees can generate noticeable mechanical turbulence (Jegede and Foken, 1999). In the 60–90 sector, the mechanical turbulence was probably due to eucalypt trees, whose crowns were quickly re-established after the wildfire by resprouting. The increased values in the 270–330 sector can only be explained by slope parallel flow (see Fig. 1 of Oliveira et al., 2021).

Evaluation of both dimensionless parameters, together with an assessment of terrain and post-fire vegetation recovery, suggested that the assumptions in Oliveira et al. (2021) of z-d=8 m and z0=0.4 m were adequate for most wind sectors for the post-fire and pre-windthrow period. However, a value for z0=0.7 m would have been more appropriate for the calculation of the footprint. Likewise, the assumptions of z-d=10 m and z0=0.3 m were appropriate for most wind sectors for the post-windthrow period. These assumptions, however, were less appropriate for three sectors highlighted in Table 3, most notably after the windthrow, where no increase in the wind speed (Fig. 2) could be detected. The greater variability in roughness after the windthrow can be explained by the higher wind speeds which, in turn, were due to the increase in effective measurement height resulting from the larger distance to the canopy height.

Before the windthrow, a very low roughness height of z0=0.4 m =0.05zc was assumed because of the very wind-permeable nature of the maritime stands which, in turn, reflected the fact that the fire had consumed the complete crowns of the bulk of the pine trees. This very low roughness height could not be confirmed by either the calculations using Eqs. (2) and (4) or Fig. 6. The obtained value of z0=0.7 m agreed with the simple relationship of z0=0.1zc (Foken, 2017; Monteith and Unsworth, 2013) but not that of z0=0.2zc (Kaimal and Finnigan, 1994). The same applied, mutatis mutandis, after the windthrow. By contrast, the assumed values for the ratio d/zc of 0.5 before the windthrow and 0.666 after the windthrow were confirmed by the observations. These findings further confirmed a linear relationship between d and zc, in line what was found using approaches that explicitly consider parameters describing stand structure (Maurer et al., 2015).

4.2 Influence on carbon dioxide fluxes

The investigations carried out clearly showed that due to the different values of the effective measurement height zd before and after the wind break in the determination of the Obukhov–Lettau stability parameter; no influence on the quality flags and the measurement of the net ecosystem exchange could be detected that would have even come close to the typical error range.

By contrast, however, the footprint area increased markedly due to the change in roughness and possibly also wind speed. The difference would probably not have been as large if a slightly better value of z0=0.7 m had been assumed before the windthrow. The increase in footprint area also implied a decrease in the target areas in the footprint. In case target areas differ markedly from non-target areas in terms of carbon dioxide fluxes, the change in aerodynamic conditions would substantially affect flux measurements. At the present study site, this is probably not the case, as the non-target areas mainly differ from the target areas by their greater slope angles and not the burnt forest. A reduction in the footprint through a reduction in the measurement height is usually not possible with long-term measurement programs with a permanently installed mast because reduction in the measurement height will then typically produce flow distortion problems due to the mast. Hydraulic lifted masts are hardly ever used in such programs (Kolle et al., 2021).

5 Conclusions

The windthrow at the end of 2019 had a significant impact on the aerodynamics of the study area. The present analysis addressed dimensionless turbulence characteristics and focused on the parameters of roughness length and zero-plane displacement. Since both parameters are not independent, either the zero-plane displacement or the roughness length must be specified a priori. For the first post-fire year, Oliveira et al. (2021) selected d=3.8 m =0.5zc (z-d=8 m), justified by the open character of the burnt pine stand due to the complete consumption of most pine crowns by the wildfire. The authors' assumption of z0=0.4 m, however, should be revised to z0=0.7 m. A not very precise choice of roughness length did not have a marked effect on the footprint, as the 2019 calculations showed.

The initial assumption of this study of z-d=10 m and a roughness length of z0=0.3 m for the post-windthrow period continues to seem reasonable as well. This implied that the windthrow drastically changed aerodynamic site conditions. The increase in σw/u provided a clear indication that the distance between measurement height and canopy height had increased significantly after the windthrow.

The present study confirmed the disturbances in specific wind sectors signalled by Oliveira et al. (2021). The disturbances in the northeastern and northwestern sectors could be assigned to terrain characteristics. According to Fig. 6, a change in zero-plane displacement in this sector would not result in an improvement. If coordinate rotation is performed by means of double rotation (Kaimal and Finnigan, 1994), the problem is hardly relevant for the measurements; however, if the authors would have done it by means of planar-fit rotation (Wilczak et al., 2001), the disturbed sector would have to be rotated separately.

The change in aerodynamic conditions due to the windthrow did not have marked impacts on the calculation of the carbon dioxide fluxes, but it did substantially increase the footprint area. In the present case, this increase in footprint area implied the inclusion of sloping terrain but with essentially the same pre- and post-fire vegetation cover as the relatively flat target area so that the implications are expected to be minor. The present windthrow occurred at a relative early stage of post-fire ecosystem recovery so that a direct comparison of pre- and post-windthrow carbon dioxide fluxes was considered unwarranted.

Based on the investigation carried out, we generally recommend determining the effective measurement height and the roughness length as precisely as possible when aerodynamic conditions change in order to be able to determine changes in the footprint area. Since areas outside the target area may have an influence on the fluxes, the quality assessment of the measurement area must be carried out again, taking the footprint into account (Foken et al., 2012; Mauder et al., 2021). Alternatively, the measurement height could be adjusted so that the footprint remains almost identical.

Code and data availability

The program for the calculation of the eddy-covariance data is available at (Mauder and Foken, 2015). The daily CO2 flux data are available in Oliveira et al. (2020,; other data are available on request from the first author (

Author contributions

BRFO was responsible for ModelEco's project funding and management and was the scientist responsible for this post-doctoral study including fieldwork, data analysis (CO2 part), and the preparation of the structure of the paper and of several of its sections. JJK was responsible for FIRE-C-BUDs' project funding and management and data analysis (wind part) and revised the paper. TF was the scientific adviser of the project, mainly responsible for the flux part and special aerodynamic analysis, and supported the writing of the paper. The final version of the paper was prepared by all authors.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Financial support

This work was supported by ModelEco (grant no. PTDC/ASP-SIL/3504/2020), funded by national funds (OE) through FCT MCTES; FIRE-C-BUDs (grant no. PTDC/AGR-FOR/4143/2014 – POCI-01-0145-FEDER-016780), funded by FCT MEC with co-funding by FEDER, within the PT2020 partnership agreement and Compete 2020; and CESAM (grant nos. UIDP/50017/2020 and UIDB/50017/2020), funded by national funds through FCT MCTES. This publication was funded by the German Research Foundation (DFG) and the University of Bayreuth within the funding program Open Access Publishing and by CESAM.

Review statement

This paper was edited by Paul Stoy and reviewed by two anonymous referees.


Arya, S. P.: Introduction to Micrometeorology, Academic Press, San Diego, 415 pp., 2001. 

Aubinet, M., Vesala, T., and Papale, D. (Eds.): Eddy Covariance: A Practical Guide to Measurement and Data Analysis, Springer, Dordrecht, Heidelberg, London, New York, 438 pp.,, 2012. 

Chu, H., Luo, X., Ouyang, Z., Chan, W. S., Dengel, S., Biraud, S. C., Torn, M. S., Metzger, S., Kumar, J., Arain, M. A., Arkebauer, T. J., Baldocchi, D., Bernacchi, C., Billesbach, D., Black, T. A., Blanken, P. D., Bohrer, G., Bracho, R., Brown, S., Brunsell, N. A., Chen, J., Chen, X., Clark, K., Desai, A. R., Duman, T., Durden, D., Fares, S., Forbrich, I., Gamon, J. A., Gough, C. M., Griffis, T., Helbig, M., Hollinger, D., Humphreys, E., Ikawa, H., Iwata, H., Ju, Y., Knowles, J. F., Knox, S. H., Kobayashi, H., Kolb, T., Law, B., Lee, X., Litvak, M., Liu, H., Munger, J. W., Noormets, A., Novick, K., Oberbauer, S. F., Oechel, W., Oikawa, P., Papuga, S. A., Pendall, E., Prajapati, P., Prueger, J., Quinton, W. L., Richardson, A. D., Russell, E. S., Scott, R. L., Starr, G., Staebler, R., Stoy, P. C., Stuart-Haëntjens, E., Sonnentag, O., Sullivan, R. C., Suyker, A., Ueyama, M., Vargas, R., Wood, J. D., and Zona, D.: Representativeness of Eddy-Covariance flux footprints for areas surrounding AmeriFlux sites, Agr. Forest Meteorol., 301–302, 108350,, 2021. 

Dupont, S. and Brunet, Y.: Coherent structures in canopy edge flow: a large-eddy simulation study, J. Fluid Mech., 630, 93–128,, 2009. 

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

Foken, T.: Micrometeorology, 2nd Edn., Springer, Berlin, Heidelberg, 362 pp.,, 2017. 

Foken, T. and Börngen, M.: Lettau's Contribution to the Obukhov Length Scale: A Scientific Historical Study, Bound.-Lay. Meteorol., 179, 369–383,, 2021. 

Foken, T. and Leclerc, M. Y.: Methods and limitations in validation of footprint models, Agr. Forest Meteorol., 127, 223–234,, 2004. 

Foken, T. and Wichura, B.: est. Meteorol., 78, 83–105,, 1996. 

Foken, T., Leuning, R., Oncley, S. P., Mauder, M., and Aubinet, M.: Corrections and data quality in: Eddy Covariance: A Practical Guide to Measurement and Data Analysis, edited by: Aubinet, M., Vesala, T., and Papale, D., Springer, Dordrecht, Heidelberg, London, New York, 85–131,, 2012. 

Foken, T., Babel, W., Munger, J. W., Grönholm, T., Vesala, T., and Knohl, A.: Selected breakpoints of net forest carbon uptake at four eddy-covariance sites, Tellus B, 73, 1–12,, 2021. 

Garratt, J., Wilczak, J., Holtslag, A., Schmid, H. P., Grachev, A., Beljaars, A., Foken, T., Chen, F., Fairall, C., Hicks, B., Kusaka, H., Martilli, A., Masson, V., Mauder, M., Oncley, S., Rotach, M., and Tjernström, M.: Commentaries on Top-Cited Boundary-Layer Meteorology Articles, Bound.-Lay. Meteorol., 177, 169–188,, 2020. 

Garratt, J. R.: The Atmospheric Boundary Layer, Cambridge University Press, Cambridge, 316 pp., 1992. 

Göckede, M., Foken, T., Aubinet, M., Aurela, M., Banza, J., Bernhofer, C., Bonnefond, J. M., Brunet, Y., Carrara, A., Clement, R., Dellwik, E., Elbers, J., Eugster, W., Fuhrer, J., Granier, A., Grünwald, T., Heinesch, B., Janssens, I. A., Knohl, A., Koeble, R., Laurila, T., Longdoz, B., Manca, G., Marek, M., Markkanen, T., Mateus, J., Matteucci, G., Mauder, M., Migliavacca, M., Minerbi, S., Moncrieff, J., Montagnani, L., Moors, E., Ourcival, J.-M., Papale, D., Pereira, J., Pilegaard, K., Pita, G., Rambal, S., Rebmann, C., Rodrigues, A., Rotenberg, E., Sanz, M. J., Sedlak, P., Seufert, G., Siebicke, L., Soussana, J. F., Valentini, R., Vesala, T., Verbeeck, H., and Yakir, D.: Quality control of CarboEurope flux data – Part 1: Coupling footprint analyses with flux data quality assessment to evaluate sites in forest ecosystems, Biogeosciences, 5, 433–450,, 2008. 

Jegede, O. O. and Foken, T.: A study of the internal boundary layer due to a roughness change in neutral conditions observed during the LINEX field campaigns, Theor. Appl. Climatol., 62, 31–41,, 1999. 

Kader, B. A. and Perepelkin, V. G.: Profil skorosti vetra i temperatury v prizemnom sloje atmosfery v uslovijach nejtralnoj i neustojtschivoj stratifikacii (The wind and temperature profile in the near surface layer for neutral and unstable stratification), Izv. AN SSSR, Fiz. Atm. i Okeana, 20, 151–161, 1984. 

Kaimal, J. C. and Finnigan, J. J.: Atmospheric Boundary Layer Flows: Their Structure and Measurement, Oxford University Press, New York, NY, 289 pp.,, 1994. 

Kaimal, J. C., Wyngaard, J. C., Izumi, Y., and Coté, O. R.: Spectral characteristics of surface layer turbulence, Q. J. Roy. Meteor. Soc., 98, 563–589,, 1972. 

Kanani-Sühring, F. and Raasch, S.: Spatial variability of scalar concentrations and fluxes downstream of a clearing-to-forest transition: A Large-Eddy Simulation study, Bound.-Lay. Meteorol., 155, 1–27,, 2015. 

Klaassen, W., van Breugel, P. B., Moors, E. J., and Nieveen, J. P.: Increased heat fluxes near a forest edge, Theor. Appl. Climatol., 72, 231–243,, 2002. 

Kolle, O., Kalthoff, N., Kottmeier, C., and Munger, J. W.: Ground based platforms, in: Springer Handbook of Atmospheric Measurements, edited by: Foken, T., Springer Nature, Switzerland, 155–182,, 2021. 

Kormann, R. and Meixner, F. X.: An analytical footprint model for non-neutral stratification, Bound.-Lay. Meteorol., 99, 207–224,, 2001. 

Leclerc, M. Y. and Foken, T.: Footprints in Micrometeorology and Ecology, Springer, Heidelberg, New York, Dordrecht, London, XIX, 239 pp.,, 2014. 

Marquardt, D.: An algorithm for least-sqares estimation of nonlinear parameters, J. Soc. Ind. Appl. Math., 11, 431–441, 1983. 

Marunitsch, S. V.: Charakteristiki turbulentnosti v yslovijach lesa po gradientnym i strukturnym nabljudenijam (Turbulence characteristics of gradients in a forest from structure observations), Trudy Gos. Gidrometeorol. Inst., 198, 154–165, 1971. 

Mauder, M. and Foken, T.: Eddy-Covariance software TK3, Zenodo [software],, 2015. 

Mauder, M., Foken, T., and Cuxart, J.: Surface Energy Balance Closure over Land: A Review, Bound.-Lay. Meteorol., 177, 395–426,, 2020. 

Mauder, M., Foken, T., Aubinet, M., and Ibrom, A.: Eddy-Covariance Measurements, in: Springer Handbook of Atmospheric Measurements, edited by: Foken, T., Springer Nature, Switzerland, 1473–1504,, 2021. 

Maurer, K. D., Bohrer, G., Kenny, W. T., and Ivanov, V. Y.: Large-eddy simulations of surface roughness parameter sensitivity to canopy-structure characteristics, Biogeosciences, 12, 2533–2548,, 2015. 

Monteith, J. L. and Unsworth, M. H.: Principles of Environmental Physics, 4th Edn., Academic Press, Oxford, 401 pp., 2013. 

Moore, C. J.: Frequency response corrections for eddy correlation systems, Bound.-Lay. Meteorol., 37, 17–35,, 1986. 

Nakai, T., Sumida, A., Daikoku, K. i., Matsumoto, K., van der Molen, M. K., Kodama, Y., Kononov, A. V., Maximov, T. C., Dolman, A. J., Yabuki, H., Hara, T., and Ohta, T.: Parameterisation of aerodynamic roughness over boreal, cool- and warm-temperate forests, Agr. Forest Meteorol., 148, 1916–1925,, 2008. 

Nieuwstadt, F. T. M.: The computation of the friction velocity u and the temperature scale T from temperature and wind velocity profiles by least-square method, Bound.-Lay. Meteorol., 14, 235–246,, 1978.  

Oliveira, B. R. F., Keizer, J. J., and Foken, T.: Daily Carbon Dioxide fluxes measured by an eddy-covariance station in a recently burnt Mediterranean pine stand in Central Portugal, PANGAEA [data set],, 2020. 

Oliveira, B. R. F., Schaller, C., Keizer, J. J., and Foken, T.: Estimating immediate post-fire carbon fluxes using the eddy-covariance technique, Biogeosciences, 18, 285–302,, 2021. 

Paeschke, W.: Experimentelle Untersuchungen zum Rauhigkeitsproblem in der bodennahen Luftschicht, Z. Geophys., 13, 14–21, 1937. 

Panofsky, H. A.: Vertical variation of roughness length at the Boulder Atmospheric Observatory, Bound.-Lay. Meteorol., 28, 305–308,, 1984. 

Panofsky, H. A., Tennekes, H., Lenschow, D. H., and Wyngaard, J. C.: The characteristics of turbulent velocity components in the surface layer under convective conditions, Bound.-Lay. Meteorol., 11, 355–361,, 1977. 

Prandtl, L.: Meteorologische Anwendung der Strömungslehre, Beitr. Phys. freien Atmosphäre, 19, 188–202, 1932. 

Raupach, M. R.: Simplified expressions for vegetation roughness lenght and zero-plane displacement as functions of canopy height and area index, Bound.-Lay. Meteorol., 71, 211–216,, 1994. 

Stoy, P. C., Mauder, M., Foken, T., Marcolla, B., Boegh, E., Ibrom, A., Arain, M. A., Arneth, A., Aurela, M., Bernhofer, C., Cescatti, A., Dellwik, E., Duce, P., Gianelle, D., van Gorsel, E., Kiely, G., Knohl, A., Margolis, H., McCaughey, H., Merbold, L., Montagnani, L., Papale, D., Reichstein, M., Serrano-Ortiz, P., Sottocornola, M., Saunders, M., Spano, D., Vaccari, F., and Varlagin, A.: A data-driven analysis of energy balance closure across FLUXNET research sites: The role of landscape-scale heterogeneity, Agr. Forest Meteorol., 171–172, 137–152,, 2013. 

Stull, R. B.: An Introduction to Boundary Layer Meteorology, Kluwer, Dordrecht, 666 pp.,, 1988. 

Thom, A. S.: Momentum absorption by vegetation, Q. J. Roy. Meteorol. Soc., 97, 414–428,, 1971. 

Wilczak, J. M., Oncley, S. P., and Stage, S. A.: Sonic anemometer tilt correction algorithms, Bound.-Lay. Meteorol., 99, 127–150,, 2001. 

Short summary
This study analyzes the impacts of this windthrow on the aerodynamic characteristics of zero-plane displacement and roughness length and, ultimately, their implications for the turbulent fluxes. The turbulent fluxes were only affected to a minor degree by the windthrow, but the footprint area of the flux tower changed markedly so that the target area of the measurements had to be redetermined.
Final-revised paper