The impact of logging on the surrounding flow in a managed forest

The influence of a freshly logged area in a managed pine forest on the flow field is investigated by comparing sodar wind profile data over the forest canopy with the synoptic wind field extracted from North American Regional Reanalysis, National Centers for Environmental Prediction. As a consequence of the pressure gradient arising from the sharp temperature difference between the clearcut and the surrounding uncut forests, the local wind direction over the forest measured with the sodar departs dramatically from the prevailing synoptic wind direction when the latter is transverse to the clearcut-sodar direction. Sodar measurements also indicate systematic strong updrafts during daytime followed by nighttime downdrafts with wind coming from the logged area. This suggests the presence of horizontal advection carrying daytime warm air (or nighttime cool air) from the clearcut to the forested area. This paper also examines the influence of wind velocity, clearcut fetch, and solar radiation on locally generated circulations and advection. The presence of local circulations arising from contrasting neighboring surface characteristics well outside the footprint is of particular relevance for atmospheric flux sites where robust surface?Catmosphere exchange values are sought. This study highlights the high level of circumspection required at the time of identifying locations for flux sites. It also suggests vigilant monitoring of the surrounding landscape during eddy?Cflux measurements particularly in actively managed landscapes.     

Characterizing the dispersive state of convective boundary layers for applied dispersion modeling

Estimates from semiempirical models that characterize surface heat flux, mixing depth, and profiles of temperature, wind, and turbulence are compared with observations from atmospheric field studies conducted in Colorado, Illinois, Indiana, and Minnesota. Sodar observations are compared with tower measurements at the Colorado site, for wind and turbulence profiles. The median surface heat flux, as calculated using surface-layer flux-profile relationships and an energy budget model, was consistently overestimated by 20 to 80%. Several mixing-depth models were evaluated: (1) integration of the hourly surface heat flux and friction velocity, (2) solving for the time rate of change of profiles of virtual potential temperature, and (3) an interpolation scheme used by the U.S. Environmental Protection Agency in regulatory dispersion models. For the late afternoon, 80 to 90% of the estimates from the first and third models were within 40% of the observed values. For the morning hours after sunrise, all were less accurate. Temperature estimates from surface-layer flux-profile relationships compared well with observations within the mixed layer, but were too low for the inversion layer aloft. Wind profiles were derived using surface-layer flux-profile relationships, a windprofile power-law based on Pasquill stability category, and sodar measurements. The sodar measurements were superior to both types of model estimates. Turbulence profiles were derived from sodar measurements and from semiempirical similarity relationships based on mixing depth and Obukhov length. The scatter in the comparisons with the sodar observations is twice that seen in the comparisons with empirical profile relationships. Overall, it appears that uncertainty of as low as 20 to 30% in the characterization of the diffusion meteorology is the exception rather than the rule.On assignment from the National Oceanic and Atmospheric Administration, U. S. Department of Commerce.Disclaimer: Although the research described in this article has been supported by the United States Environmental Protection Agency, it has not been subjected to Agency review and therefore does not necessarily reflect the views of the Agency and no official endorsement should be inferred. Mention of trade names or commercial products does not constitute endorsement or recommendation for use.     

Detection of Inversions and Mixing Height by REMTECH PA2 Sodar in Comparison with Collocated Radiosonde Measurements

Summary A REMTECH PA2 Doppler Sodar is operated regularly at the Czech Hydrometeorological Institute (CHMI) observatory in Prague, collocated with a routine rawinsonde sounding system. The Air Pollution Control Division of CHMI utilises the sodar data in air pollution studies and as an information support for the smog warning system operated in Prague. Besides of the basic software for echo strength and wind profile evaluation, optional routines for deriving parameters such as inversion and mixing height, stability class etc. were delivered by the sodar manufacturer. Based on a sufficiently large data set (more than one year) of synchronous sodar and rawinsonde measurements, an analysis and comparison of inversion and mixing heights provided by both sounding systems have been accomplished in order to evaluate the correctness and accuracy of sodar estimates of these parameters. In contrast to the wind speed and wind direction data, for which a satisfactory agreement with other kind of measurements has been reported by many studies, the results for inversion and mixing height detection were totally disappointing. A direct applicability of inversion height and mixing height data provided by the REMTECH’s automatic routines in air pollution studies or smog warning systems is quite problematical with the present “state of the art”. Received November 3, 1998 Revised April 20, 1999     

The diurnal wind variation in a variable eddy viscosity semi-geostrophic Ekman boundary-layer model: Analytical study

Summary ?A time-dependent semi-geostrophic Ekman boundary-layer model (SG), including slowly varying eddy diffusivity with height and inertial term effects, is developed to investigate the diurnal wind variation in the planetary boundary layer (PBL). An approximate analytical solution of this model is derived by using the WKB method, which extends the Tan and Farahani (1998)’s solution by including the vertical variable eddy viscosity. The features of the diurnal wind variation in the PBL mainly depend on three factors: the latitude, horizontal momentum advection and eddy viscosity. The vertical variable eddy viscosity has little influence on diurnal wind variation in the PBL at the low latitude, however its effect may be exacerbated in the mid- and high latitudes. In comparing with the constant eddy viscosity case, the decreasing (increasing) with height eddy viscosity produces a large (small) maximum wind speed (MWS) in the PBL, however, the eddy viscosity that has a mid-layer peak in the vertical gives rise to a higher height of occurrence of MWS. For the boundary-layer wind structure, there is a singular point when the modified SG inertial oscillation frequency η equals the forcing frequency ω. The isotachs of boundary-layer wind speed have almost no tilt to left or right relative to time evolution and the occurrence time of the MWS is the earliest at the singular point. The feature will be enhanced in the decreasing with height eddy viscosity case and weakened in the eddy viscosity initially increasing with height case. Received April 6, 2001; accepted December 27, 2001     

Evaluation of the profile and the resistance method for estimation of surface fluxes of momentum,sensible and latent heat

Summary Hourly lysimetric and micrometeorological data taken over a grass surface at the Meteorological Research Unit, Cardington U.K. have been analysed. A temperature difference and measurements of wind speed at only one height, combined with an independently estimated effective roughness length allowed sensible heat and momentum fluxes determination by the profile method on an hourly basis. The estimates are compared with direct measurements of sensible heat and friction velocity obtained by the eddy correlation method. The sensible and latent heat fluxes are also modelled by the resistance method. Equations based on the Monin—Obukhov similarity theory are used to account for stability effects through various forms of parameterization Aerodynamic and surface resistances, necessary for the Penman—Monteith equation are calculated from routinely measured meteorological data. The profile method for estimation of sensible heat flux and friction velocity is found to work excellently on the discussed daytime experimental data which correspond mainly to near neutral or slightly unstable conditions.Surface latent and sensible heat fluxes can also be described very well by the resistance method. A slightly better estimate of the sensible heat flux is achieved when stability corrections are taken into account. On the contrary Penman-Monteith equation for estimating latent heat flux is insensitive to adjustments for atmospheric stability.The comparison of the various methods leads to the establishment of empirical relationships which correlate various quantities such as soil heat flux, resistances, evapotranspiration etc. to routinely measured meteorological data.With 8 Figures     

Use of a High-Resolution Sodar to Study Surface-layer Turbulence at Night

Measurements in the atmospheric surface layer are generally made with point sensors located in the first few tens of metres. In most cases, however, these measurements are not representative of the whole surface layer. Standard Doppler sodars allow a continuous display of the turbulent thermal structure and wind profiles in the boundary layer up to 1000 m, with a few points, if any, in the surface layer. To overcome these limitations a new sodar configuration is proposed that allows for a higher resolution in the surface layer. Because of its capabilities (echo recording starting at 2 m, echo intensity vertical resolution of approximately 2 m, temporal resolution of 1 s) this sodar is called the surface-layer mini-sodar (SLM-sodar). Features and capabilities of the SLM-sodar are described and compared with the sodar. The comparison of the thermal vertical structure given by the SLM-sodar and the sodar provides evidence that, in most cases, the surface layer presents a level of complexity comparable to that of the entire boundary layer. Considering its high vertical resolution, the SLM-sodar is a promising system for the study of the nocturnal surface layer. The nocturnal SLM-sodar measurements have shown that, depending on wind speed, the structure of the surface layer may change substantially within a short time period. At night, when the wind speed is greater than 3 m s⁻¹, mechanical mixing destroys the wavy structure present in the nocturnal layer. Sonic anemometer measurements have shown that, in such cases, also the sensible heat flux varies with height, reaching a peak in correspondence with the wind speed peak. Under these conditions the assumption of horizontal homogeneity of the surface layer and the choice of the averaging time need to be carefully treated.     

Foehn signals detected by sodar wind and turbulence measurements in the Rhine Valley, Austria, during the MAP field phase

Summary Within the Mesoscale Alpine Programme MAP conducted in autumn 1999, the vertical structure and the temporal evolution of the planetary boundary layer (PBL) in the Rhine Valley 2km south of Lake Constance were observed with a Remtech PA2 sodar (sound-detection-and-ranging instrument) rendering half-hour averages of the three-dimensional wind profile within the lowest kilometre above ground. During Foehn events, tethered balloon soundings and wind profiler measurements were conducted in addition to the rawinsonde network which was built up for the MAP field campaign.The remote sensing instrument renders a surprisingly high number of valid data during south Foehn. Due to the frequent formation of a cold air pool with stable conditions below the Foehn flow with near-neutral static stability, even more sodar data is valid during Foehn periods than during no Foehn periods. A significant reduction of the sodar data quality is only observed during Foehn events with grounding of the Foehn at the sodar site due to high background noise. At higher levels, a Foehn signal can be detected from the sodar wind and turbulence intensitiy information. With Foehn, higher wind speeds and larger turbulence intensities occur than without Foehn. Comparisons to rawinsonde and tethersonde soundings and wind profiler measurements at sites nearby reveal the spatial inhomogeneity of the Foehn flow within this part of the valley as well as instrumental short-comings. Different methods to determine the mixing height using the vertical sounding devices lead to some uncertainty of mixing height estimates which however can reasonably be explained.     

Study of near-surface models for large-eddy simulations of a neutrally stratified atmospheric boundary layer

In large-eddy simulations (LES) of the atmospheric boundary layer (ABL), near-surface models are often used to supplement subgrid-scale (SGS) turbulent stresses when a major fraction of the energetic scales within the surface layer cannot be resolved with the temporal and spatial resolution at hand. In this study, we investigate the performance of both dynamic and non-dynamic eddy viscosity models coupled with near-surface models in simulations of a neutrally stratified ABL. Two near-surface models that are commonly used in LES of the atmospheric boundary layer are considered. Additionally, a hybrid Reynolds- averaged/LES eddy viscosity model is presented, which uses Prandtl’s mixing length model in the vicinity of the surface, and blends in with the dynamic Smagorinsky model away from the surface. Present simulations show that significant portions of the modelled turbulent stresses are generated by the near-surface models, and they play a dominant role in capturing the expected logarithmic wind profile. Visualizations of the instantaneous vorticity field reveal that flow structures in the vicinity of the surface depend on the choice of the near-surface model. Among the three near-surface models studied, the hybrid eddy viscosity model gives the closest agreement with the logarithmic wind profile in the surface layer. It is also observed that high levels of resolved turbulence stresses can be maintained with the so-called canopy stress model while producing good agreement with the logarithmic wind profile.     

Corrections for Wind-Speed Errors from Sodar and Lidar in Complex Terrain

The quality of lidar and sodar wind estimates is generally judged through comparisons with mast-mounted instruments, and the resulting regressions. Evaluation of the relative merits of lidars versus sodars is complicated by the fact that lidars are generally placed close to a mast whereas sodars are generally placed some distance from a mast so that acoustic reflections off the mast are reduced. This leads to the two technologies, lidar and sodar, not being compared in similar situations. Differences arising from the two geometries can be expected to be larger in complex terrain, where the wind regime can vary significantly spatially. The current work explores these differences in moderately complex terrain. Lidar–mast comparisons are performed with the lidar close to an 80 m mast, and sodar–mast comparisons performed with the sodar 300 m from the mast. Systematic variations in estimated wind speed are found to occur with height, consistent with predictions from a simple flow model. When the lidar was moved to the sodar location, further from the mast, there were significant changes in the estimated wind speeds and a reduction in correlation with the mast-based wind speeds, as expected. However, the correlation between collocated lidar and sodar winds was high. This finding emphasizes that any comparison of two remote sensing instruments needs to be through similar experiments, and that differences in accuracy often reported for the lidar and sodar technologies are likely to be contaminated due to poor comparison configurations. A method was devised to simulate the sodar being collocated with the mast, by using the lidar–sodar measurements and the lidar–mast measurements. It was found that there was then no statistically detectable difference between lidar–mast regressions and sodar–mast regressions for the particular lidar and sodar tested. Both remote sensing instruments were also found to be good estimators of Weibull parameters, as compared with those derived from mast data. The conclusion is that the sodar measured the winds above the sodar with a similar accuracy to the lidar measuring winds above the lidar.     

Comparison of three methods for determining strong wind stress over Lake Flevo

The drag coefficient C _d (10 m) at the center of shallow Lake Flevo (20-km diam) is evaluated for wind speeds u between 5 and 15 m s^?1 independently by three methods. Trivane measurements of eddy-correlation fluxes agree with eddy flux data available for moderate wind speeds from other sites, and can together be represented by C _d(10 m) = 0.0007 μ^0.3. Additional evaluations of water-surface slope give C _d(10 m) ≈ 0.0024, indicating that the stress at the water-surface level may not be entirely accounted for by eddy-correlation measurements well above the waves. Neither the eddy-correlation stress, nor the water-surface stress appears to be accurately estimable from profile measurements of wind, temperature and humidity analyzed without regard to sea state, if u > 10 m s^?1.     

An Approximate Analytical Solution For The Baroclinic And Variable Eddy Diffusivity Semi-Geostrophic Ekman Boundary Layer

The WKB method has been used to develop an approximate solutionof the semi-geostrophic Ekman boundary layer with height-dependenteddy viscosity and a baroclinic pressure field. The approximate solutionretains the same simple form as the classical Ekman solution. Behavioursof the approximate solution are discussed for different eddy viscosityand the pressure systems. These features show that wind structure inthe semi-geostrophic Ekman boundary layer depends on the interactionbetween the inertial acceleration, variable eddy viscosity and baroclinicpressure gradient. Anticyclonic shear has an acceleration effect on theair motion in the boundary layer, while cyclonic shear has a decelerationeffect. Decreasing pressure gradient with height results in a super-geostrophicpeak in the wind speed profile, however the increasing pressure gradient withheight may remove the peak. Anticyclonic shear and decreasing the variableeddy viscosity with height has an enhanced effect on the peak.Variable eddy viscosity and inertial acceleration has an important role in thedivergence and vorticity in the boundary layer and the vertical motion at the top of the boundary layer that is called Ekman pumping. Compared to the constanteddy viscosity case, the variable eddy diffusivity reduces the absolute value ofEkman pumping, especially in the case of eddy viscosity initially increasing with height. The difference in the Ekman pumping produced by different eddy diffusivity assumptions is intensified in anticyclonic flow and reduced in cyclonic flow.     

Sodar Detected Top-Down Convection in a Nocturnal Cloud-Topped Boundary Layer: A Case Study

Nocturnal convection, originating in a well-mixed marine cloud-topped boundary layer, advected onshore, was observed using a Doppler sodar on the Tyrrhenian coast in Italy. The horizontal and vertical dimensions of the downdrafts were evaluated. The oscillation frequency triggered by the downdrafts at the inversion layer, derived from the harmonic analysis of the sodar measured vertical velocity (w), is compared with the Brunt-Vaisala frequency, obtained from the rawinsonde temperature profile. A similarity function for the ²_w vertical profile was used to fit the sodar experimental data and to retrieve the depth of the mixing layer and the sensible heat flux at the top of the cloud layer. The results are in agreement with the convection layer depth observed in the sodar echoes facsimile record, and with the energy budget evaluated at the top of the cloud layer using the rawinsonde profiles.     

The CRPE minisodar: Applications in micrometeorology and in physics of precipitations

The technical characteristics of the CRPE minisodar are described and it is compared with the classical Doppler sodar.The results of surface layer analysis are given and simple methods are demonstrated for estimating the sensible heat flux and friction velocity in the surface layer.The system is shown to be feasible for working on the birth of gravity waves and the analysis of wind profile and turbulence evolution.Hydrometeorological studies with the minisodar are discussed with special emphasis on the problem of the relationship between wind shear in the surface layer and hydrometeor fall speed.In the conclusion, the prospective uses of the minisodar are presented.     

Regional shear stress of broken forest from radiosonde wind profiles in the unstable surface layer

Mean wind speed profiles were measured by tracking radiosondes in the unstable atmospheric boundary layer (ABL) over the forested Landes region in southwestern France. New Monin-Obukhov stability correction functions, recently proposed following an, analysis by Kader and Yaglom, as well as the Businger-Dyer stability formulation were tested, with wind speeds in the surface sublayer to calculate the regional shear stress. These profile-derived shear stresses were compared with eddy correlation measurements gathered above a mature forest stand, at a location roughly, 4.5 km from the radiosonde launch site. The shear stress values obtained by means of the newly proposed stability function were in slightly better agreement with the eddy correlation values than those obtained by means of a Businger-Dyer type stability function. The general robustness of the profile method can be attributed in part to prior knowledge of the regional surface roughness (z ₀=1.2 m) and the momentum displacement height (d ₀=6.0 m), which were determined from neutral wind profile analysis. The 100 m drag coefficient for the unstable conditions above this broken forest surface was found to beu _* ² /V ₁₀₀ ² =0.0173.     

Case Studies of the Wintertime Convective Boundary-Layer Structure in the Urban Area of Milan, Italy

The paper describes some aspects of the convective boundary-layer structure based on simultaneous sodar and tethersonde measurements during a field experiment in the urban area of Milan in the period 8 to 20 February, 1993. During this period, fog episodes and strong low-level elevated inversions (with lower boundaries < 400 m) were observed most of the time. A close agreement in the mixing height values, derived from the sodar and tethersonde profiles, has been achieved under these conditions. The validity of the similarity relationships, which have been originally derived to describe the vertical velocity variance and heat flux profiles over horizontally homogeneous terrain under quasi-stationary conditions, was evaluated when applied to the urban boundary layer.     

Composing the diurnal cycle of mixing height from simultaneous sodar and wind profiler measurements

Mixing heights have been determined independently from sodar and boundary-layer wind profiler measurements carried out at the Meteorological Observatory Lindenberg of the German Weather Service between 1 June 1994 and 6 July 1994. Good agreement between both systems was found for mixing height values between about 250 m and 700 m, i.e., in an evolving convective boundary layer. Considering the typical sounding ranges of the sodar (50 m up to 800 m) and of the wind profiler (200 m up to higher than 3000 m), simultaneous operation of the two systems is demonstrated to be a promising tool for continously monitoring the mixing height throughout its complete diurnal cycle.     

Measurements Of Turbulence Transfer In The Near-Surface Layer Over The Southeastern Tibetan Plateau

Turbulent flux measurements at Qamdo site over the Tibetan Plateau during TIPEX from May 18 to June 30, 1998 are presented. Sensible heat dominated,accounting for about 66% of the available energy (the sum of net radiation and soil heat flux) prior to the monsoon(dry period), reducing to about 31%, with latent heat increased to about 56% of available energy,in the monsoon season (wet period). Surface energy budget closure on average was about 0.80 (0.85)prior to the monsoon and 0.89 (0.76) during the monsoon using eddy correlation (profile) methods. The sum of latent and sensible heat fluxes calculated from the flux-profilemethod was smaller by about 15% than that from eddy correlation. Martano's method is used toestimate the surface aerodynamic roughness length z₀ and zero plane displacement d from singlelevel sonic anemometer data, giving d = 0.12 m and z₀ = 0.08 m. The overall neutral dragcoefficient (C_DN) and scalar coefficient (C_HN) were found to be C_DN = 0.0055and C_HN = 0.0059 in the southeastern area of Tibet. Their variations with the mean wind speed at 10 m are discussed.     

The influence of nonstationarity on the turbulent flux–gradient relationship for stable stratification

Extensive eddy-correlation datasets are analyzed to examine the influence of nonstationarity of the mean flow on the flux–gradient relationship near the surface. This nonstationarity is due to wavelike motions, meandering of the wind vector, and numerous unidentified small-scale mesoscale motions. While the data do not reveal an obvious critical gradient Richardson number, the maximum downward heat flux increases approximately linearly with increasing friction velocity for significant stability. The largest of our datasets is chosen to more closely examine the influence of stability, nonstationarity, distortion of the mean wind profile and self-correlation on the flux-gradient relationship. Stability is expressed in terms of z/L, the gradient Richardson number or the bulk Richardson number over the tower layer. The efficiency of the momentum transport systematically increases with increasing nonstationarity and attendant distortion of the mean wind profile. Enhancement of the turbulent momentum flux associated with nonstationarity is examined in terms of the nondimensional shear, Prandtl number and the eddy diffusivity.     

Comparison of wind measurements by means of industrially produced sodar and high-altitude meteorological mast in Obninsk

The real characteristics of the MFAS sodar industrially produced by Scintec AG (Germany) are given. The paper presents a wider range of altitudes as compared with the papers published before as a result of using the meteorological mast with the height of 310 m and longer duration of practically continuous testing embracing different weather conditions in the period from July to November 2008. During the comparison, the 10-min averaging time for the wind speed and direction is used. The comparison of the data of acoustic sounding and contact measurement methods indicates that the sodar provides reliable information on the vertical profiles of average wind speed within the considered altitude range up to 300 m except the cases of heavy rainfalls and low wind speed.     

How do Stability Corrections Perform in the Stable Boundary Layer Over Snow?

We assess sensible heat-flux parametrizations in stable conditions over snow surfaces by testing and developing stability correction functions for two alpine and two polar test sites. Five turbulence datasets are analyzed with respect to, (a) the validity of the Monin–Obukhov similarity theory, (b) the model performance of well-established stability corrections, and (c) the development of new univariate and multivariate stability corrections. Using a wide range of stability corrections reveals an overestimation of the turbulent sensible heat flux for high wind speeds and a generally poor performance of all investigated functions for large temperature differences between snow and the atmosphere above (>10 K). Applying the Monin–Obukhov bulk formulation introduces a mean absolute error in the sensible heat flux of \(6\,\hbox {W m}^{-2}\) (compared with heat fluxes calculated directly from eddy covariance). The stability corrections produce an additional error between 1 and \(5\,\hbox {W m}^{-2}\), with the smallest error for published stability corrections found for the Holtslag scheme. We confirm from previous studies that stability corrections need improvements for large temperature differences and wind speeds, where sensible heat fluxes are distinctly overestimated. Under these atmospheric conditions our newly developed stability corrections slightly improve the model performance. However, the differences between stability corrections are typically small when compared to the residual error, which stems from the Monin–Obukhov bulk formulation.