Airborne measurements of turbulent fluxes during LITFASS-98: Comparison with ground measurements and remote sensing in a case study

Summary ?Simultaneous flight measurements with the research aircraft Do 128 and the helicopter-borne turbulence probe Helipod were performed on 18 June 1998 during the LITFASS-98 field experiment. The area-averaged turbulent vertical fluxes of momentum, sensible, and latent heat were determined on a 15 km × 15 km and a 10 km × 10 km flight pattern, respectively. The flights were carried out over heterogeneous terrain at different altitudes within a moderately convective boundary layer with Cumulus clouds. Co-spectra-analysis demonstrated that the small scale turbulent transport was completely sampled, while the comparatively small flight patterns were possibly of critical size regarding the large-scale turbulence. The phygoide of the airplane was identified as a significant peak in some co-spectra. The turbulent fluxes of momentum and sensible heat at 80 m above the ground showed systematic dependence on the location of the flight legs above the heterogeneous terrain. This was not observed for the latent heat flux, probably due to the vertical distribution of humidity in the boundary layer. Statistical error analysis of the fluxes F showed that the systematic statistical error ΔF was one order of magnitude smaller than the standard deviation σ _F . The difference between area-averaged fluxes derived from simultaneous Helipod and Do 128 measurements was much smaller than σ _F , indicating that the systematic statistical error was possibly over-estimated by the usual method. In the upper half of the boundary layer the airborne-measured sensible heat flux agreed well with windprofiler/RASS data. A linear fit was the best approximation for the height dependence of all three fluxes. The linear extrapolations of the latent and sensible heat fluxes to the ground were in good agreement with tower, scintillometer, and averaged ground-station measurements on various surface types. Systematic discrepancies between airborne and ground-based measurements were not found. Received June 18, 2001; revised December 21, 2001; accepted June 3, 2002     

Mean Profiles of Moisture Fluxes in Snow-Filled Boundary Layers

Profiles of moisture fluxes have been examined for convective boundary layers containing clouds and snow, using data derived from aircraft measurements taken on four dates during the 1983/1984 University of Chicago lake-effect snow project. Flux profiles were derived from vertical stacks of aircraft cross-wind flight legs taken at various heights over Lake Michigan near the downwind shore. It was found that, if ice processes are taken into account, profiles of potential temperature and water content were very similar to those presented in past studies of convective boundary layers strongly heated from below. Profiles of total water content and equivalent potential temperature adjusted for ice were nearly invariant with height, except very near the top of the boundary layer, suggesting that internal boundary-layer mixing processes were rapid relative to the rates at which heat and vapour were transported into the boundary layer through entrainment and surface fluxes. Ice was found to play a significant, measurable role in boundary-layer moisture fluxes. It was estimated that 40 to 57% of the upward vapour flux was returned to the surface in the form of snow, converting about 45 to 64% of the surface latent heat flux into sensible heat in the snow-producing process. Assuming advective fluxes are relatively small (thought to be appropriate after the first few tens of km over the lake as suggested by past studies), the boundary layer was found to warm at a rate faster than could be explained by surface heat fluxes and latent heat releases alone, the remainder of the heating presumably coming from radiational processes and entrainment. Discussions of moisture phase change processes throughout the boundary layer and estimates of errors of these flux measurements are presented.     

Estimation Of Atmospheric Water Vapour Flux Profiles In The Nocturnal Unstable Urban Boundary Layer With Doppler Sodar And Raman Lidar

The vertical wind profiles determined by Doppler sodar and the water vapourmixing ratio profiles obtained by Raman lidar are used to estimate the atmosphericwater vapour flux profiles in the nocturnal urban boundary layer under unstableconditions. The experiment was conducted for several nights in the central areaof Rome under a variety of moisture conditions and different urban boundary-layerflow regimes. Despite some scatter in the profiles, the latent heat flux is found tobe positive throughout the depth of the nocturnal urban boundary-layer. Thelayer-averaged flux shows a variation between -4 to +40 W m^-2, whileindividual values of flux in excess of +150 W m^-2 pertain to a case offree convection during cold air advection caused by the sea breeze. The qualityof flux estimates is found to be highly limited by the low sampling rates employedin the experiment resulting in errors to the order of 60%. Therefore, the results mustbe viewed as estimates rather than precise measurements. The skewness profiles ofthe turbulent fluctuations of vertical velocity and water vapour mixing ratio are alsopositive.     

Wind profiles,momentum fluxes and roughness lengths at Cabauw revisited

We describe the results of an experiment focusing on wind speed and momentum fluxes in the atmospheric boundary layer up to 200 m. The measurements were conducted in 1996 at the Cabauw site in the Netherlands. Momentum fluxes are measured using the K-Gill Propeller Vane. Estimates of the roughness length are derived using various techniques from the wind speed and flux measurements, and the observed differences are explained by considering the source area of the meteorological parameters. A clear rough-to-smooth transition is found in the wind speed profiles at Cabauw. The internal boundary layer reaches the lowest k-vane (20 m) only in the south-west direction where the obstacle-free fetch is about 2 km. The internal boundary layer is also reflected in the roughness lengths derived from the wind speed profiles. The lower part of the profile (< 40 m) is not in equilibrium and no reliable roughness analysis can be given. The upper part of the profile can be linked to a large-scale roughness length. Roughness lengths derived from the horizontal wind speed variance and gustiness have large footprints and therefore represent a large-scale average roughness. The drag coefficient is more locally determined but still represents a large-scale roughness length when it is measured above the local internal boundary layer. The roughness length at inhomogeneous sites can therefore be determined best from drag coefficient measurements just above the local internal boundary layers directly, or indirectly from horizontal wind speed variance or gustiness. In addition, the momentum and heat fluxes along the tower are analysed and these show significant variation with height related to stability and possibly surface heterogeneity. It appears that the dimensionless wind speed gradients scale well with local fluxes for the variety of conditions considered, including the unstable cases.     

The structure of the convective atmospheric boundary layer as revealed by lidar and Doppler radars

Results on the structure of the convective atmospheric boundary layer based on the analyses of data from the instrumented NSSL-KTVY tower, airborne Doppler lidar, and ground-based Doppler radars are presented. The vertically averaged wind over the boundary layer was found to be insensitive to baroclinicity, supporting the hypothesis of Arya and Wyngaard (1975). The computed momentum flux profiles were affected by baroclinicity. Horizontal wind spectra from lidar, radar, and tower data compared well with each other both in shape and magnitude. A consistent peak found near 4 km in all the computed spectra might have been caused by horizontally symmetric cells with horizontal wavelength 4 times the boundary-layer height as shown in Kuettner (1971) for the case of weak wind shear.     

Can Water Vapour Raman Lidar Resolve Profiles of Turbulent Variables in the Convective Boundary Layer?

High-resolution water vapour measurements made by the Atmospheric Radiation Measurement (ARM) Raman lidar operated at the Southern Great Plains Climate Research Facility site near Lamont, Oklahoma, U.S.A. are presented. Using a 2-h measurement period for the convective boundary layer (CBL) on 13 September 2005, with temporal and spatial resolutions of 10 s and 75 m, respectively, spectral and autocovariance analyses of water vapour mixing ratio time series are performed. It is demonstrated that the major part of the inertial subrange was detected and that the integral scale was significantly larger than the time resolution. Consequently, the major part of the turbulent fluctuations was resolved. Different methods to retrieve noise error profiles yield consistent results and compare well with noise profiles estimated using Poisson statistics of the Raman lidar signals. Integral scale, mixing-ratio variance, skewness, and kurtosis profiles were determined including error bars with respect to statistical and sampling errors. The integral scale ranges between 70 and 130 s at the top of the CBL. Within the CBL, up to the third order, noise errors are significantly smaller than sampling errors and the absolute values of turbulent variables, respectively. The mixing-ratio variance profile rises monotonically from ≈0.07 to ≈3.7 g² kg⁻² in the entrainment zone. The skewness is nearly zero up to 0.6 z/z _i , becomes −1 around 0.7–0.8 z/z _i , crosses zero at about 0.95 z/z _i , and reaches about 1.7 at 1.1 z/z _i (here, z is the height and z _i is the CBL depth). The noise errors are too large to derive fourth-order moments with sufficient accuracy. Consequently, to the best of our knowledge, the ARM Raman lidar is the first water vapour Raman lidar with demonstrated capability to retrieve profiles of turbulent variables up to the third order during daytime throughout the atmospheric CBL.     

Spectral Structure of Small-Scale Turbulent and Mesoscale Fluxes in the Atmospheric Boundary Layer over a Thermally Inhomogeneous Land Surface

Spectral analysis was performed on aircraft observations of a convective boundary layer (CBL) that developed over a thermally inhomogeneous, well-marked mesoscale land surface. The observations, part of the GAME-Siberia experiment, were recorded between April and June 2000 over the Lena River near Yakutsk City. A special integral parameter termed the ‘reduced depth of the CBL’ was used to scale the height of the mixed layer with variable depth. Analysis of wavelet cospectra and spectra facilitated the separation of fluxes and other variables into small-scale turbulent fluctuations (with scales less than the reduced depth of the CBL, approximately 2 km) and mesoscale fluctuations (up to 20 km). This separation approach allows for independent exploration of the scales. Analyses showed that vertical distributions obeyed different laws for small-scale fluxes and mesoscale fluxes (of sensible heat, water vapour, momentum and carbon dioxide) and for other variables (wind speed and air temperature fluctuations, coherence and degree of anisotropy). Vertical profiles of small-scale turbulent fluxes showed a strong decay that differed from generally accepted similarity models for the CBL. Vertical profiles of mesoscale fluxes and other variables clearly showed sharp inflections at the same relative (with respect to the reduced depth of the CBL) height of approximately 0.55 in the CBL. Conventional similarity models for sensible heat fluxes describe both small-scale turbulent and mesoscale flows. The present results suggest that mesoscale motions that reach up to the relative level of 0.55 could be initiated by thermal surface heterogeneity. Entrainment between the upper part of the CBL and the free atmosphere may cause mesoscale motions in that region of the CBL.     

Entrainment Processes in the Convective Boundary Layer with Varying Wind Shear

Large-eddy simulations (LES) are performed to investigate the entrainment andthe structure of the inversion layer of the convective boundary layer (CBL) withvarying wind shears. Three CBLs are generated with the constant surface kinematicheat flux of 0.05 K m s^-1 and varying geostrophic wind speeds from 5 to 15m s^-1. Heat flux profiles show that the maximum entrainment heat flux as afraction of the surface heat flux increases from 0.13 to 0.30 in magnitude withincreasing wind shear. The thickness of the entrainment layer, relative to the depthof the well-mixed layer, increases substantially from 0.36 to 0.73 with increasingwind shear. The identification of vortices and extensive flow visualizations nearthe entrainment layer show that concentrated vortices perpendicular to the meanboundary-layer wind direction are identified in the capping inversion layer for thecase of strong wind shear. These vortices are found to develop along the mean winddirections over strong updrafts, which are generated by convective rolls and to appearas large-scale wavy motions similar to billows generated by the Kelvin–Helmholtzinstability. Quadrant analysis of the heat flux shows that in the case of strong windshear, large fluctuations of temperature and vertical velocity generated by largeamplitude wavy motions result in greater heat flux at each quadrant than that inthe weak wind shear case.     

BOUNDARY LAYER GROWTH AND LAPSE RATE CHANGES DETERMINED BY LIDAR AND SURFACE HEAT FLUX IN SOFIA*

In this study the results from a boundary layer experiment,conducted in autumn 1991 over a flat,build-up urban area in Southeast Sofia,together with some models for mixed layer growth rates are used to investigate the layered structure of the vertical atmospheric stability distribution in the Sofia Valley.Lidar measurements of aerosol layer heights and morning boundary layer development are combined with surface eddy correlation measurements of kinematic heat and moisture fluxes,profiles of temperature and humidity,wind speed and wind direction.A diagnostic method is presented for determining vertical lapse rates using surface meteorological measurements and lidar returns observed during the transition from nighttime stable stratification to daytime convective boundary layer after the sunrise.     

BOUNDARY LAYER GROWTH AND LAPSE RATE CHANGES DETERMINED BY LIDAR AND SURFACE HEAT FLUX IN SOFIA

In this study the results from a boundary layer experiment,conducted in autumn 1991 over a flat,build-up urbanarea in Southeast Sofia,together with some models for mixed layer growth rates are used to investigate the layered struc-ture of the vertical atmospheric stability distribution in the Sofia Valley.Lidar measurements of aerosol layer heightsand morning boundary layer development are combined with surface eddy correlation measurements of kinematic heatand moisture fluxes,profiles of temperature and humidity,wind speed and wind direction.A diagnostic method is pres-ented for determining vertical lapse rates using surface meteorological measurements and lidar returns observed duringthe transition from nighttime stable stratification to daytime convective boundary layer after the sunrise.     

Validation of the “Lokal-Modell” over heterogeneous land surfaces using aircraft-based measurements of the REEEFA experiment and comparison with micro-scale simulations

Summary An aircraft-based experimental investigation of the atmospheric boundary layer (ABL) structure and of the energy exchange processes over heterogeneous land surfaces is presented. The measurements are used for the validation of the mesoscale atmospheric model “Lokal-Modell” (LM) of the German Weather Service with 2.8 km resolution. In addition, high-resolution simulations using the non-hydrostatic model FOOT3DK with 250 m resolution are performed in order to resolve detailed surface heterogeneities. Two special observation periods in May 1999 show comparable convective boundary layer (CBL) conditions. For one case study vertical profiles and area averages of meteorological quantities and energy fluxes are investigated in detail. The measured net radiation is highly dependent on surface albedo, and the latent heat flux exhibits a strong temporal variability in the investigation area. A reduction of this variability is possible by aggregation of multiple flight patterns. To calculate surface fluxes from aircraft measurements at low altitude, turbulent energy fluxes were extrapolated to the ground by the budget method, which turned out to be well applicable for the sensible heat flux, but not for the latent flux. The development of the ABL is well captured by the LM simulation. The comparison of spatiotemporal averages shows an underestimation of the observed net radiation, which is mainly caused by thin low-level clouds in the LM compared to observed scattered CBL clouds. The sensible heat flux is reproduced very well, while the latent flux is highly overestimated especially above forests. The realistic representation of surface heterogeneities in the investigation area in the FOOT3DK simulations leads to improvements for the energy fluxes, but an overestimation of the latent heat flux still persists. A study of upscaling effects yields more structures than the LM fields when averaged to the same scale, which are partly caused by the non-linear effects of parameter aggregation on the LM scale.     

An Optimal Inverse Method Using Doppler Lidar Measurements to Estimate the Surface Sensible Heat Flux

Inverse methods are widely used in various fields of atmospheric science. However, such methods are not commonly used within the boundary-layer community, where robust observations of surface fluxes are a particular concern. We present a new technique for deriving surface sensible heat fluxes from boundary-layer turbulence observations using an inverse method. Doppler lidar observations of vertical velocity variance are combined with two well-known mixed-layer scaling forward models for a convective boundary layer (CBL). The inverse method is validated using large-eddy simulations of a CBL with increasing wind speed. The majority of the estimated heat fluxes agree within error with the proscribed heat flux, across all wind speeds tested. The method is then applied to Doppler lidar data from the Chilbolton Observatory, UK. Heat fluxes are compared with those from a mast-mounted sonic anemometer. Errors in estimated heat fluxes are on average 18 %, an improvement on previous techniques. However, a significant negative bias is observed (on average $-63\,\%$ ) that is more pronounced in the morning. Results are improved for the fully-developed CBL later in the day, which suggests that the bias is largely related to the choice of forward model, which is kept deliberately simple for this study. Overall, the inverse method provided reasonable flux estimates for the simple case of a CBL. Results shown here demonstrate that this method has promise in utilizing ground-based remote sensing to derive surface fluxes. Extension of the method is relatively straight-forward, and could include more complex forward models, or other measurements.     

Lidar and Sun photometer observations of atmospheric boundary-layer characteristics over an urban area in a mountain valley

We present measurements of the vertical aerosol structure and the aerosol optical depth in the lower troposphere performed above the city of Sofia (an urban area situated in a mountain valley), western Bulgaria by means of a ground-based aerosol lidar operating continuously for a number of years. The lidar measurements were accompanied by measurements of the aerosol optical depth (AOD) in the visible and near infrared regions of the spectrum performed in October 2004 using Microtops II radiometers. The maximum values of the AOD were found to occur 1–2 h before the complete development of the atmospheric boundary layer, i.e. during the residual layer destruction, which confirms our hypothesis concerning the slope circulation effect on the processes taking place in the atmospheric boundary layer. The AOD values obtained by the lidar are lower than those taken by the sun photometer. Further, the AOD exhibits two different types of behaviour. In the case of a ‘clear atmosphere’ (i.e. in the absence of volcanic eruptions and/or dust transport from the Sahara) most of the aerosol accumulated within the atmospheric boundary layer over the urban area considered. The combined use of the two instruments allows the comparison between the optical characteristics of the atmospheric aerosol (e.g. aerosol extinction coefficient, etc.) obtained by the lidar and through an independent method (sun photometer).     

Obtaining Potential Virtual Temperature Profiles,Entrainment Fluxes,and Spectra from Mini Unmanned Aerial Vehicle Data

We present a simple but effective small unmanned aerial vehicle design that is able to make high-resolution temperature and humidity measurements of the atmospheric boundary layer. The air model used is an adapted commercial design, and is able to carry all the instrumentation (barometer, temperature and humidity sensor, and datalogger) required for such measurements. It is fitted with an autopilot that controls the plane’s ascent and descent in a spiral to 1800 m above ground. We describe the results obtained on three different days when the plane, called Aerolemma-3, flew continuously throughout the day. Surface measurements of the sensible virtual heat flux made simultaneously allowed the calculation of all standard convective turbulence scales for the boundary layer, as well as a rigorous test of existing models for the entrainment flux at the top of the boundary layer, and for its growth. A novel approach to calculate the entrainment flux from the top-down, bottom-up model of Wynagaard and Brost is used. We also calculated temperature fluctuations by means of a spectral high-pass filter, and calculated their spectra. Although the time series are small, tapering proved ineffective in this case. The spectra from the untapered series displayed a consistent −5/3 behaviour, and from them it was possible to calculate a dimensionless dissipation function, which exhibited the expected similarity behaviour against boundary-layer bulk stability. The simplicity, ease of use and economy of such small aircraft make us optimistic about their usefulness in boundary-layer research.     

Impact of a sea breeze on the boundary-layer dynamics and the atmospheric stratification in a coastal area of the North Sea

In-situ sodar and lidar measurements were coupled with numerical simulations for studying a sea-breeze event in a flat coastal area of the North Sea. The study’s aims included the recognition of the dynamics of a sea-breeze structure, and its effects on the lower troposphere stratification and the three-dimensional (3D) pollutant distribution. A sea breeze was observed with ground-based remote sensing instruments and analysed by means of numerical simulations using the 3D non-hydrostatic atmospheric model Meso-NH. The vertical structure of the lower troposphere was experimentally determined from the lidar and sodar measurements, while numerical simulations focused on the propagation of the sea breeze inland. The sea-breeze front, the headwind, the thermal internal boundary layer, the gravity current and the sea-breeze circulation were observed and analysed. The development of a late stratification was also observed by the lidar and simulated by the model, suggesting the formation of a stable multilayered structure. The transport of passive tracers inside the sea breeze and their redistribution above the gravity current was simulated too. Numerical modelling showed that local pollutants may travel backward to the sea above the gravity current at relatively low speed due to the shearing between the landward gravity current and the seaward synoptic wind. Such dynamic conditions may enhance an accumulation of pollutants above coastal industrial areas.     

A Study of the Internal Boundary Layer due to a Roughness Change in Neutral Conditions Observed During the LINEX Field Campaigns

Summary As an aspect of the LINEX field studies (1996–1997; Lindenberg near Beeskow, Germany), the characteristics of the internal boundary layer (IBL) that is associated with a step change of the surface roughnesses in neutral constant stress layers was investigated and is reported in this paper. Both smooth to rough (in 1996) and rough to smooth (in 1997) types of flow, have been studied based upon the profiles of mean wind and temperature realised from a 10-m mast and eddy correlation measurements taken at two levels (2 m and 5 m). Depending upon wind direction, the fetch at the site varied between 140 m and 315 m within the wind sector (200° to 340°) used for the field investigations. The height of the IBL, δ, had been determined from the intersect of the logarithmic wind-profiles below (< 2 m) and above (> 6 ) the interface. Values of δ obtained at the experimental site compared fairly well to the existing theoretical/empirical fetch-height relationships of the form: δ=aċx ^b , where a, b, are empirical constants. The ratio for the friction velocities below and above the IBL as measured directly by the eddy correlation techniques showed that for fetches less than 250 m there was an increase (decrease) of about 20% of the momentum flux arising from the smooth to rough (rough to smooth) transitions. Influences of distant obstructions (e.g., bushes, pockets of trees) on the surface flow were markedly important on the examined wind profiles and such can be indicative as multiple IBLs. Received September 1, 1997 Revised August 5, 1998     

étude expérimentale de la couche limite au-dessus d'un relief modéré proche d'une chaîne de montagne

For the heterogeneous site described in the first part, some aspects of the turbulent structure of the planetary boundary layer are studied. Using mixed-layer scaling, the normalized profiles are compared with those obtained over flat terrain during convective conditions. The measurements were made with the same instrumented aircraft at both sites. The dissipative and spectral length scales are smaller over complex terrain within the whole boundary layer. This is due to the shifting of the wavelength peak toward the high frequencies by dynamic turbulence.This last effect can also explain the increase of the dissipation rate over the heterogeneous site during strong wind conditions. The vertical profiles of sensible heat flux and temperature-water vapor correlation show a lack of entrainment process at the top of the boundary layer. This fact suggests that the investigated boundary layer is advected from the neighbouring plain over the complex site (plateau de Lannemezan).     

Determination of boundary-layer parameters using wind profiler/RASS and sodar/RASS in the frame of the LITFASS project

Summary ?Progress in technology as well as signal processing has promoted Wind Profiler Radar (WPR) or sodar with RASS additions to become standard tools in profiling of the atmospheric boundary layer. Apart from these instruments’ basic abilities in profiling mean winds and temperature, this paper will give an emphasis on the profiling of ABL height as well as the turbulent fluxes of sensible heat and momentum both, with respect to methods as well as with respect to realization. The special focus will thereby be laid on the demands for vertical profiling, which were defined within the LITFASS-project of the German Meteorological Service. In the frame of this project, some special measuring campaigns have been performed where remote-sensing systems were used to assess their abilities in profiling ABL parameters. On the base of some case studies from these campaigns comparisons are shown, where results from sodar/RASS and WPR/RASS measurements are compared to measurements from airborne sensor systems and results from numerical models. Regarding turbulent heat fluxes, we found excellent agreement for remotely-sensed flux profiles from WPR/RASS with both, numerical models and airborne in-situ measurements. However, as the inherent errors of the remotely-sensed fluxes are in the order of ± 20 ⋯ 30 W/m² typically, current signal processing does not allow to interpret small-scale vertical structures in the profiles with respect to surface inhomogeneities yet. Received June 16, 2001; revised February 20, 2002; accepted May 30, 2002     

Atmospheric turbulence within and above a douglas-fir stand. Part II: Eddy fluxes of sensible heat and water vapour

This is the second paper describing a study of the turbulence regimes and exchange processes within and above an extensive Douglas-fir stand. The experiment was conducted on Vancouver Island during a two-week rainless period in July and August 1990. Two eddy correlation units were operated in the daytime to measure the fluxes of sensible heat and water vapour and other turbulence statistics at various heights within and above the stand. Net radiation was measured above the overstory using a stationary net radiometer and beneath the overstory using a tram system. Supplementary measurements included soil heat flux, humidity above and beneath the overstory, profiles of wind speed and air temperature, and the spatial variation of sensible heat flux near the forest floor.The sum of sensible and latent heat fluxes above the stand accounted for, on average, 83% of the available energy flux. On some days, energy budget closure was far better than on others. The average value of the Bowen ratio was 2.1 above the stand and 1.4 beneath the overstory. The mid-morning value of the canopy resistance was 150–450 s/m during the experiment and mid-day value of the Omega factor was about 0.20. The daytime mean canopy resistance showed a strong dependence on the mean saturation deficit during the two-week experimental period.The sum of sensible and latent heat fluxes beneath the overstory accounted for 74% of the available energy flux beneath the overstory. One of the reasons for this energy imbalance was that the small number of soil heat flux plates and the short pathway of the radiometer tram system was unable to account for the large horizontal heterogeneity in the available energy flux beneath the overstory. On the other hand, good agreement was obtained among the measurements of sensible heat flux made near the forest floor at four positions 15 m apart.There was a constant flux layer in the trunk space, a large flux divergence in the canopy layer, and a constant flux layer above the stand. Counter-gradient flux of sensible heat constantly occurred at the base of the canopy.The transfer of sensible heat and water vapour was dominated by intermittent cool downdraft and warm updraft events and dry downdraft and moist updraft events, respectively, at all levels. For sensible heat flux, the ratio of the contribution of cool downdrafts to that of warm updrafts was greater than one in the canopy layer and less than one above the stand and near the forest floor.