A spectral and lag-correlation analysis of turbulence in a deciduous forest canopy

The processes influencing turbulence in a deciduous forest and the relevant length and time scales are investigated with spectral and cross-correlation analysis. Wind velocity power spectra were computed from three-dimensional wind velocity measurements made at six levels inside the plant canopy and at one level above the canopy. Velocity spectra measured within the plant canopy differ from those measured in the surface boundary layer. Noted features associated with the within-canopy turbulence spectra are: (a) power spectra measured in the canopy crown peak at higher wavenumbers than do those measured in the subcanopy trunkspace and above the canopy; (b) peak spectral values collapse to a relatively universal value when scaled according to a non-dimensional frequency comprised of the product of the natural frequency and the Eulerian time scale for vertical velocity; (c) at wavenumbers exceeding the spectral peak, the slopes of the power spectra are more negative than those observed in the surface boundary layer; (d) Eulerian length scales decrease with depth into the canopy crown, then increase with further depth into the canopy; (e) turbulent events below crown closure are more correlated with turbulent events above the canopy than are those occurring in the canopy crown; and (f) Taylor's frozen eddy hypothesis is not valid in a plant canopy. Interactions between plant elements and the mean wind and turbulence alter the processes that produce, transport and remove turbulent kinetic energy and account for the noted observations.     

Estimating the components of the sensible heat budget of a tall forest canopy in complex terrain

Ultrasonic wind measurements, sonic temperature and air temperature data at two heights in the advection experiment MORE II were used to establish a complete budget of sensible heat including vertical advection, horizontal advection and horizontal turbulent flux divergence. MORE II took place at the long-term Carbo-Europe IP site in Tharandt, Germany. During the growing period of 2003 three additional towers were established to measure all relevant parameters for an estimation of advective fluxes, primarily of CO₂. Additionally, in relation to other advection experiments, a calculation of the horizontal turbulent flux divergence is proposed and the relation of this flux to atmospheric stability and friction velocity is discussed. In order to obtain a complete budget, different scaling heights for horizontal advection and horizontal turbulent flux divergence are tested. It is shown that neglecting advective fluxes may lead to incorrect results. If advective fluxes are taken into account, the sensible heat budget based upon vertical turbulent flux and storage change only, is reduced by approximately 30%. Additional consideration of horizontal turbulent flux divergence would in turn add 5–10% to this sum (i.e., the sum of vertical turbulent flux plus storage change plus horizontal and vertical advection). In comparison with available energy horizontal advection is important at night whilst horizontal turbulent flux divergence is rather insignificant. Obviously, advective fluxes typically improve poor nighttime energy budget closure and might change ecosystem respiration fluxes considerably.     

Influence of non-flat terrain and wind direction shear on canopy turbulence

We have analyzed eddy covariance data collected within open canopy to investigate the influence of non-flat terrain and wind direction shear on the canopy turbulence. The study site is located on non-flat terrain with slopes in both south-north and east-west directions. The surface elevation change is smaller than the height of roughness element such as building and tree at this site. A variety of turbulent statistics were examined as a function of wind direction in near-neutral conditions. Heterogeneous surface characteristics results in significant differences in measured turbulent statistics. Upwind trees on the flat and up-sloping terrains yield typical features of canopy turbulence while upwind elevated surface with trees yields significant wind direction shear, reduced u and w skewness, and negligible correlation between u and w. The directional dependence of turbulence statistics is due that strong wind blows more horizontally rather than following terrain, and hence combination of slope related momentum flux and canopy eddy motion decreases the magnitude of Sk _w and r _uw for the downslope flow while it enhances them for the upslope flow. Significant v skewness to the west indicates intermittent downdraft of northerly wind, possibly due to lateral shear of wind in the presence of significant wind direction shear. The effects of wind direction shear on turbulent statistics were also examined. The results showed that correlation coefficient between lateral velocities and vertical velocity show significant dependence on wind direction shear through change of lateral wind shear. Quadrant analysis shows increased outward interaction and reduced role of sweep motion for longitudinal momentum flux for the downslope flow. Multi-resolution analysis indicates that uw correlation shows peak at larger averaging time for the upslope flow than for the downslope flow, indicating that large eddy plays an active role in momentum transfer for the upslope flow. On the other hand, downslope flow shows larger velocity variances than other flows despite similar wind speed. These results suggest that non-flatness of terrain significantly influences on canopy-atmosphere exchange.     

On the effect of clearcuts on turbulence structure above a forest canopy

Summary The modification of the flow structure arising from the removal of large patches of trees in a managed forest plantation near Gainesville, Florida is described. Using wavelet analysis of turbulence measurements taken above a forest canopy hundreds of meters downwind from the forest gap and well outside the footprint, the present paper examines changes in flow characteristics and demonstrates that the presence of the nearby clearcut introduces extraneous coherent events passing by the eddy-covariance flux measurement system.     

Statistics of atmospheric turbulence within a natural black spruce forest canopy

Turbulence statistics were measured in a natural black-spruce forest canopy in southeastern Manitoba, Canada. Sonic anemometers were used to measure time series of vertical wind velocity (w), and cup anemometers to measure horizontal wind speed (s), above the canopy and at seven different heights within the canopy. Vertical profiles were measured during 25 runs on eight different days when conditions above the canopy were near-neutral.Profiles of s and of the standard deviation ( _w ) of w show relatively little scatter and suggest that, for this canopy and these stability conditions, profiles can be predicted from simple measurements made above the canopy. Within the canopy, a negative skewness and a high kurtosis of the w-frequency distributions indicate asymmetry and the persistence of large, high-velocity eddies. The Eulerian time scale is only a weak function of height within the canopy.Although w-power spectra above the canopy are similar to those in the free atmosphere, we did not observe an extensive inertial subrange in the spectra within the canopy. Also, a second peak is present that is especially prominent near the ground. The lack of the inertial subrange is likely caused by the presence of sources and sinks for turbulent kinetic energy within our canopy. The secondary spectral peak is probably generated by wake turbulence caused by form drag on the wide, horizontal spruce branches.     

Turbulence measurements over complex terrain

Horizontal turbulence measurements obtained from 22 wind sensors located on 9 towers in a mountainous coastal area are described and categorized by stability and terrain. Vector wind time series are high-pass filtered, and lateral and longitudinal wind speed variance is calculated for averaging times ranging from 15 s to 2 h. Parameterizations of the functional dependence of variance on averaging time are discussed, and a modification of Panofsky's (1988) uniform terrain technique applicable to complex terrain is presented. The parameterization is applied to the data and shown to be more realistic than a less complicated power law technique. The parameter values are shown to be different than the flat terrain cases of Kaimal et al. (1972), and are primarily a function of sensor location within the complex terrain. The parameters are also examined in terms of their dependence upon season, stability, marine boundary-layer height, and measurement height.     

Wind profiles over complex terrain

A large set of routine wind profiles has been analyzed at three towers, located in various types of complex terrain in New England. After allowing for effects of roughness change, uphill flow and stability, roughness lengths have been estimated. In general, the profiles and roughness lengths could be explained by differences in terrain features and stability.     

复杂地形低风速气象特征分析

利用13个地面站和1个100 m气象塔的全年实测资料分析湖南中北部丘陵地区中尺度范围的低风速气象特征。结果表明：（1）该区域低风速出现频率高,小于15 m/s的风速出现频率年均达416%;（2）低风速与地形情况密切相关,区域内西部山地的低风速出现频率明显高于东部平原地区;（3）低风速的出现具有日变化特点,典型出现时间段是夜间至次日上午;（4）低风速条件下10 min平均风向多出现大角度变化的情况,60°以上的风向变化频率可达20%～30%。（5）水平风向标准差在低风速情况下明显增大,但观测仪器性能会对小风结果产生影响。     

Boundary-layer pollutant concentrations over complex terrain

Given the incident profiles of wind velocity and pollutant concentration, we seek to determine the 3-dimensional concentration field of a pollutant upon a region with complex terrain. The analytic solution of the wind velocity in a 3-dimensional boundary-layer model by Walmsley et al. (1980) is utilized as a forcing function in the simplified concentration perturbation equation for a pollutant. The resulting solution applies to an isolated cosine-squared hill in a neutrally stratified boundary-layer flow with a surface type which absorbs the pollutant totally. The solution shows that the concentration perturbation field is organized in accordance with the wind field. In particular, the east-west cross-section is 180° out of phase with the velocity perturbation field. The vertical profiles of the concentration perturbations for selected grid points approach the value of the upper boundary condition very rapidly.     

The intermittent vertical heat flux over a spruce forest canopy

Vertical heat and momentum fluxes were measured by the eddy correlation method under near-neutral conditions during both day and night above a spruce forest canopy. Results show that 50% of the heat transported to the spruce canopy during night and away from the canopy during the day occurs in extreme magnitude events, the majority of less than one second duration. Extreme-magnitude events were more frequent and lasted longer during the day than during the night. The distributions of the duration of extreme events in the same direction as the net heat flux and the turbulence intensity for both day and night were similar.During the night, the mean horizontal windspeed was about 1 m s^-1 and the measured coincident transport of momentum and heat accounted for 36% of the total heat flux. The predominant mechanism of forced convection during the extreme nightime heat transport events was excess heat sweeps in which the duration of the event is usually less than 1 s. During the day, the mean horizontal windspeed was about 2 m s ^-1 and the measured coincident transport of momentum and heat accounted for only 16% of the total heat flux. Local free convection was suggested to account for 27% of the total heat flux. The predominant mechanism of mixed convection during the extreme daytime heat transport events is deficit heat inward interactions. During both night and day, about 10% of the total heat transport occurred in extreme events working against the thermal gradient.Now atPurdue University, West Lafayette, IN, U.S.A.     

Seasonal and diurnal variations of coherent structures over a deciduous forest

Coherent structures in turbulent flow above a midlatitude deciduous forest are identified using a wavelet analysis technique. Coupling between motions above the canopy (z/h=1.5, whereh is canopy height) and within the canopy (z/h=0.6) are studied using composite velocity and temperature fields constructed from 85 hours of data. Data are classified into winter and summer cases, for both convective and stable conditions. Vertical velocity fluctuations are in phase at both observation levels. Horizontal motions associated with the structures within the canopy lead those above the canopy, and linear analysis indicates that the horizontal motions deep in the canopy should lead the vertical motions by 90°. On average, coherent structures are responsible for only about 40% of overall turbulent heat and momentum fluxes, much less than previously reported. However, our large data set reveals that this flux fraction comes from a wide distribution that includes much higher fractions in its upper extremes. The separation distanceL _s between adjacent coherent structures, 6–10h, is comparable to that obtained in previous observations over short canopies and in the laboratory. Changes in separation between the summer and winter (leafless) conditions are consistent withL _s being determined by a local horizontal wind shear scale.     

复杂地形下冰雹预报的数值试验研究

根据威宁地区的地形、气候特征，建立三维冰雹云数值模式，对冰雹研究有一定参考价值。     

Turbulent dispersion of pollutant over complex terrain

A two-dimensional turbulent diffusion equation is derived in a streamline coordinate system, defined for rotational flow over complex terrain and limited aloft by an elevated, impenetrable inversion. In the first instance, the steady-state equation is solved for an inner region of the boundary layer, in which the effect of curvature is negligible and, for simplicity, it is assumed that vorticity has a power-law dependence upon stream function. A variational method of solution is also discussed, in which vorticity may have a more general representation. A numerical calculation is performed for a special case of symmetrical flow over an isolated hill. The dependence of pollutant concentration upon the flow field, downwind distance and source is examined and the effect of wind acceleration in the neighbourhood of the top of the hill is discussed. It is pointed out that the diffusion model can be applied to any realistic flow field, provided that the streamlines are specified.     

On mountain wave drag over complex terrain

Summary Mountain wave drag is calculated for rotating, stratified, nonhydrostatic Boussinesq flow over a mountain ridge using linear theory for a variety of mountain profiles representing complex/irregular terrain. The inclusion of a sinusoidal corrugation to the familiar witch-of-Agnesi profile creates a stegosaurus profile. The associated drag is greatly enhanced for mesoscale mountains when the corrugation wave-number matches that for the dominant inertia-gravity wave contribution to the cross-mountain surface pressure gradient. Similarly, increasing the jaggedness (by decreasing the exponentb) increases the drag for mesoscale mountains whose topographic spectral intensity,M(k), has the form of a power law:M(k)=mk ^–b wherek is the zonal wavenumber.Spectral analysis of one-kilometer resolution topographic data for the Appalachian Mountains suggests that a power law profile withb=1.7 accurately represents the topographic spectral intensity and that it yields good estimates of the drag.The application of these results to the parameterization of mountain wave drag in general circulation models is discussed.With 7 Figures     

Estimates of effective surface roughness over complex terrain

Turbulence data collected in an area of three-dimensional complex terrain using instruments atteched to the tether cable of a captive balloon together with radiosonde ascents are presented. In addition, data collected using only radiosonde ascents in an area of two-dimensional complex terrain of large slope are also shown. Eddy correlation measurements of the turbulent momentum flux and wind velocity profiles are used to deduce the magnitude of the effective roughness from the drag coefficient and normalised velocity profiles. A relationship connecting the terrain characteristics and the roughness length is compared with the experimental data for both types of terrain plus previous experimental estimates of the roughness length over complex terrain. The formula taken from previous work by Grant and Mason (1990) is found to agree with the data when representing an area of order 100 km².     

Prediction of surface potential temperature over complex terrain

Based on a statistical approach, the surface potential temperature at seven observing stations in complex terrain has been examined. It is shown that the surface potential temperature depends primarily on the rate of change of slope wind and on the geostrophic-level potential temperature.     

Some turbulence characteristics of convectively mixed layers over rugged and homogeneous terrain

Turbulence characteristics of convectively mixed layers over rugged and homogenous terrain are presented, and their differences are sought. The data used in the study were obtained from aircraftborne instruments. Apart from statistics related to the mixing ratio q, there does not seem to be a marked difference between the two sets of data. This lack of difference is attributed to the degree of ruggedness of the terrain and the dominant wind direction. The behaviour of some statistics related to q is attributed to the entrainment effect of the unstable layer of the moist, cool air above the mixed layer.This work was initiated and a large part of it was completed while all three authors were with CSIRO (Commonwealth Scientific Industrial Research Organisation), Australia.Formerly at: Division of Atmospheric Research, CSIRO, Epping, NSW, Australia.Formerly at: Division of Mathematics and Statistics, CSIRO, Lindfield, NSW, Australia.     

Influence of foliar density and thermal stability on profiles of Reynolds stress and turbulence intensity in a deciduous forest

Observations were made of turbulence in an extensive deciduous forest on level terrain using a vertical array of seven three-dimensional sonic anemometer/thermometers within and above the canopy. Data were collected through the period of leaf fall and over a range of thermal stabilities. A bulk canopy drag coefficient was nearly independent of the density of the forest but decreased greatly with the onset of nocturnal stability. The depth of penetration of momentum into the forest increased with leaf fall but, again, was greatly curtailed by stable conditions. Turbulent velocities decreased with increasing depth in the forest but relative turbulence intensities increased to mid-canopy levels. Leaf density influenced turbulence levels but not as strongly as did thermal stability. Thermal effects were adequately described by the single parameter h/L, where h is the canopy height and L is the Monin-Obukhov length. The longitudinal and vertical velocity correlation coefficient was larger in magnitude than expected in the upper layers of the forest but decreased to a small value in the lowest layers where the Reynolds stress was small. The ratio _w /u ^*, where u ^* is the local friction velocity, reflected changes in the uw correlation, becoming smaller than usual in the upper canopy layers. It is believed that these effects result from the intermittent, spatially coherent structures that are responsible for a large fraction of the momentum flux to the forest.     

A case study of the nocturnal boundary layer over a complex terrain

A case study of the structure of the nocturnal boundary layer (NBL) over complex terrain is presented. Observations were made during the third night of Project STABLE (Weber and Kurzeja, 1991), whose main goal was to study turbulence and diffusion over the complex terrain of the Savannah River Site (SRS) near Augusta, Georgia.The passage of a mesoscale phenomenon, defined as a turbulent meso-flow (TMF) with an explanation of the nomenclature used, and a composite structure of the lowest few hundred meters over complex terrain are presented. The spatial extent of the TMF was at least 30–50 km, but the forcing is not well understood. The TMF occurred without the presence of a synoptic-scale cold front, under clear conditions, and with no discernible discontinuity in a microbarograph pressure trace. The structure of the NBL over the complex terrain at SRS differed from the expected homogeneous terrain NBL. The vertical structure exhibited dual low level wind maxima, dual inversions, and a persistent elevated turbulent layer.The persistent elevated turbulent layer, with a spatial extent of at least 30 km, was observed for the entire night. The persistent adiabatic layer may have resulted from turbulence induced by shear instability.