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

Doppler sodar derived values of the temperature structure parameter C _infT ^sup2 , the vertical velocity variance ¯′w ², and the rate of dissipation of turbulent energy ?, were measured during unstable conditions above the Lannemezan heterogeneous site. The vertical profiles of these turbulent parameters, normalized by the classical convective scales are compared with those obtained using the same acoustic sounder above an homogeneous site during convective conditions. The typical decrease of C _infT ^sup2 as Z ^-4/3 is partially verified on the heterogeneous site: for the lower levels, C _infT ^sup2 exhibits an increase with Z whereas for the intermediate levels C _infT ^sup2 . decreases as Z ^-4/3. For the upper levels, C _infT ^sup2 increases with Z due to a signal-to-noise ratio lower than 1. The vertical profiles of ¯′w ² above the two sites are rather similar. However, near the base of the convective inversion Z _i , the values measured on the heterogeneous site are more scattered. The same scattering is also observed with the ? values; moreover, for the lower levels (Z<0.17Z _i ) the increase of ? as Z decreases is more important at the homogeneous site than at the heterogeneous one. It is suggested that these particular features observed at lower levels above the heterogeneous site are mainly related to a complex local boundary layer induced by the near environment of the sodar (vegetation and relief).     

Similarity scaling applied to sodar observations of the convective boundary layer above an irregular hill

Nine profiles of the temperature structure parameter C _T ² and the standard deviation of vertical velocity fluctuations ( _w) in the convective boundary layer (CBL) were obtained with a monostatic Doppler sodar during the second intensive field campaign of the First ISLSCP Field Experiment in 1987. The results were analyzed by using local similarity theory. Local similarity curves depend on four parameters: the height of the mixed layer (z _i ), the depth of the interfacial layer (), and the temperature fluxes at the top of the mixed layer (Q _i ) and the surface (Q _o). Values of these parameters were inferred from sodar data by using the similarity curve for C _T ² and observations at three points in its profile. The effects of entrainment processes on the profiles of C _T ² and _wnear the top of the CBL appeared to be described well by local similarity theory. Inferred estimates of surface temperature flux, however, were underestimated in comparison to fluxes measured by eddy correlation. The measured values of _wappeared to be slightly smaller than estimates based on available parmeterizations. These discrepancies might have been caused by experimental error or, more likely, by the distortion of turbulence structure above the site by flow over the nonuniform terrain at the observation site.     

SODAR-based estimation of TKE and momentum flux profiles in the atmospheric boundary layer: Test of a parameterization model

Summary A simple parameterization for the estimation of turbulent kinetic energy (TKE) and momentum flux profiles under near-neutral stratification based on sodar measurements of the vertical velocity variance has been tested using data from the LINEX-2000 experiment. Measurements included operation of a phased-array Doppler sodar DSDPA.90 and of a sonic anemometer USA-1 mounted at a meteorological tower at a height of 90m. Good agreement has been found between the TKE and momentum flux values derived from the sonic and sodar data (with correlation coefficients r>0.90 and a slope of the regression lines of about 1.01.1) suggesting the possible use of sodar measurements of _w ² to derive turbulence parameter profiles above the tower range.     

Similarity of structure-function parameters in the stably stratified boundary layer

The structure-function parametersC _T ² andC _v ² of temperature and velocity, respectively, from the 1973 Minnesota experiments and from large-eddy and direct numerical simulations show a smooth transition from M–O similarity to the local scaling hypothesized by Nieuwstadt for the outer regions of the stable boundary layer. Under that hypothesis, turbulence statistics aloft depend on the local vertical fluxes of momentum and temperature, so these results suggest that remote-sensing measurements ofC _T ² andC _v ² could be used to infer vertical profiles of those fluxes. We argue that the sensitivity of the fluxes to unsteadiness, baroclinity, terrain slope, and breaking gravity waves precludes the universality of the vertical profiles of structure-function parameters in the stable PBL. We find that theC _T ² profile is particularly sensitive to these effects, which is consistent with observations that it varies considerably from case to case.     

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.     

Vertical velocity field in the convective boundary layer as observed with an acoustic Doppler sodar

The morning development of the daytime convective boundary layer (CBL) during fine weather has been observed with an acoustic Doppler sodar of the C.R.P.E. In particular, the vertical profile of the vertical velocity third-order statistic W* ³ has been obtained. This quantity is a maximum near 0.3z _I where z _I, is the height of the CBL. The histogram of vertical velocity in the CBL shows a relationship between W ³ and the convective velocity W _* and is useful for convective plume determination.     

An aircraft study of turbulence dissipation rate and temperature structure function in the unstable marine atmospheric boundary layer

The dissipation rate of turbulent kinetic energy, , and the temperature structure function parameter, C _T ², have been measured over water from the near surface (Z = 3 m) to the top of the boundary layer. The near surface values of and C _T ² were used to calculate the velocity and temperature Monin-Obukhov scaling parameters u _* and T _*. The data collected during unstable lapse rates were used to evaluate the feasibility of extrapolating the values of and C _T ² as a function of height with empirical scaling formulae. The dissipation rate scaling formula of Wyngaard et al. (l971 a) gave a good fit to an average of the data for Z < 0.8 Z _i. In the surface layer the scaling formula of Wyngaard et al. (1971b) disagreed with the C _T ² values by as much as 50%. This disagreement is due to an unexpected reduction in the measured values of C _T ² forZ < 30 m. At this point it is not clear if the discrepancy is a unique property of the marine boundary layer or if it is simply some unknown instrumental or analytical problem. The mixed layer scaling results were similar to the overland results of Kaimal et al. (1976).     

Rapid vertical sampling in stably stratified turbulent shear flows

A new method for obtaining instantaneous vertical profiles of two components of velocity and temperature in thermally stratified turbulent shear flows is presented. In this report, the design and construction of the traversing system will be discussed and results to date will be presented. The method is based on rapid vertical sampling whereby probe sensors are moved vertically at a high speed such that the measurement is approximately instantaneous. The system is designed to collect many measurements for the calculation of statistics such as vertical wave number spectra, mean square vertical gradients, and Thorpe scales. Results are presented for vertical profiles of temperature and compared to vertical profiles measured by single-point Eulerian time averages. The quality of the vertical profiles is found to be good over many profiles. Some comparisons are made between vertical measurements and standard single-point Eulerian measurements for three cases of stably stratified turbulent shear flow in which the initial microscale Reynolds number, Re_λ≈30. In case 1, the mean conditions are characterized by a gradient Richardson number, Ri_g=0.015, for which the flow is “unstable”, meaning the spatially evolving turbulent kinetic energy (E_k) grows. In case 2, Ri_g=0.095, for which the evolving turbulent kinetic energy is almost constant. In case 3, the flow is highly stable, where Ri_g=0.25 and E_k decays with spatial evolution. The measurements indicate anisotropy in the small scales for all cases. In particular, it is found that the ratio grows initially to a maximum and then decays with further evolution. Maximum Thorpe displacements are measured and compared to single-point measures of the vertical scales. It is found that vertical length scales derived from single-point measurements, such as the Ozmidov scale, L_O=(ε/N³)^1/2 and the overturn scale, L_t=θ′/(dT/dz), do not represent well the wide range of overturning scales which are actually present in the turbulence.     

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     

Markov-chain simulation of particle dispersion in inhomogeneous flows: The mean drift velocity induced by a gradient in Eulerian velocity variance

The Langevin equation is used to derive the Markov equation for the vertical velocity of a fluid particle moving in turbulent flow. It is shown that if the Eulerian velocity variance _wE is not constant with height, there is an associated vertical pressure gradient which appears as a force-like term in the Markov equation. The correct form of the Markov equation is: w(t + t) = aw(t) + b _wE + (1 – a)T _L ( _wE ²)/z, where w(t) is the vertical velocity at time t, a random number from a Gaussian distribution with zero mean and unit variance, T _L the Lagrangian integral time scale for vertical velocity, a = exp(–t/T _L), and b = (1 – a ²)^1/2. This equation can be used for inhomogeneous turbulence in which the mean wind speed, _wE and T _L vary with height. A two-dimensional numerical simulation shows that when this equation is used, an initially uniform distribution of tracer remains uniform.     

A sodar study of the temperature structure parameter in the convective boundary layer

From sodar measurements gathered during the Voves experiment (France, summer 1977), the variations of the temperature structure parameter C _T ² were studied in the morning planetary boundary layer. Dimensionless profiles of C _T ² are consistent with the mixed-layer scaling of Kaimal et al. (1976); however, for z < 0,5 z _i, the decrease of C _T ² as z ^4/3 should be weighted according to Frisch and Ochs (1975).When the final breakup of the nocturnal inversion is achieved, the variations of the maximum of the C _T ² profile are in good agreement with those predicted by Wyngaard and Le Mone (1980). Discrepancies are observed mainly when the mixed layer is shallow and mechanical turbulence is important compared with buoyancy-driven turbulence.     

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.     

Displaced-Beam Small Aperture Scintillometer Test. Part Ii: Cases-99 Stable Boundary-Layer Experiment

The performance of the Scintec displaced-beam small aperture scintillometer (DBSAS) in the stable boundary layer (SBL) is investigated using data gathered during the CASES-99 experiment in Kansas, U.S.A. The DBSAS is superior to the eddy-covariance method in determining vertical fluxes of sensible heat and momentumclose to the ground and/or over short (< 1 min) averaging intervals. Both aspects are of importance in the shallow and non-stationary SBL.The friction velocity, u_*, the temperature scale, _*, and from these the sensible heat flux, H, were calculated from the indirectly determined dissipation rate, , and the structure parameter of temperature, C_T ², by the DBSAS, which was operated over a path length of 112 m. All these variables are compared with eddy-covariance data for 10-minute time averages. Previously reported systematic errors in the DBSAS, overestimation of u_* for low u_* values and underestimation of u_* for high u_* values, have in part been dealt with by adjusting the beam displacement distance from 2.7 mm to 2.6 mm in the calculations. The latter adjustment is presented as a working hypothesis, not a general solution.     

étude expérimentale de la couche limite de montagne

The characteristics of the boundary layer over complex terrain (Lannemezan - lat.: 43.7° N and, long.: 0.7 ° E) are analyzed for various scales, using measurements obtained during the COCAGNE Experiment. In this first part, the dynamic characteristics of the flow are studied with respect to atmospheric stability and the relief at small (~20 km) and medium scales (~100 km). These relief scales depend on the topographical profile of the Lannemezan Plateau along the dominant axis of the wind (E-W) and the Pyrénées Mountains located at the south of the experimental site. The terrain heterogeneities have a standard deviation of ~48 m and a wavelength of ~2 km.The averaged vertical profiles of wind speed and direction over the heterogeneous terrain are analyzed. The decrease of wind speed within the boundary layer is greater than over flat terrain (WANGARA Experiment). However, a comparison between ETTEX (complex terrain) and COCAGNE vertical wind speed profiles shows good agreement during unstable conditions. In contrast, during neutral conditions a more rapid increase with normalized height is found with COCAGNE than with ETTEX and WANGARA data. The vertical profiles of wind direction reveal an influence of the Pyrénées Mountains on the wind flow. The wind rotation in the BL is determined by the geostrophic wind direction-Pyrénées axis angle (negative deviation) as the geostrophic wind is connected with the Mountain axis.When the geostrophic wind does not interact with the Pyrénées axis, the mean and turbulent wind flow characteristics (drag coefficient C _D, friction velocity u _*) depend on the topography of the plateau. When the wind speed is strong (>6 m s ^-1), an internal boundary layer is generated from the leading edge of the Plateau.     

Deriving characteristic parameters of the convective boundary layer from sodar measurements of the vertical velocity variance

This paper deals with the derivation of the convective mixing height and the characteristic convective velocity w _* from profiles of _w measured by sodar. The parameters were obtained by fitting an analytical profile to the observed data. Results were compared with values obtained by the meteorological preprocessor of a dispersion model and from noon radiosoundings. In addition, a Monte Carlo method was applied to study the influence of measurement errors. It turned out that it is inherently difficult to determine the depth of deep mixed layers from sodar measurements with a limited range, although the determination of w _* should be possible. However, a significant underestimation of _w , and thus w _*, was found, which is probably due to disproportional sampling of updrafts and downdrafts.     

Characteristics of atmospheric turbulence in the surface layer over Antarctica

This paper presents meteorological measurements made during the antarctic summer period, on two 9 m and 3 m towers, on the rocky and ice shelf terrains of the Indian antarctic stations Maitri and Dakshin Gangotri, respectively. The measurements of fluctuations in temperature and wind speed made with relatively lesser precision instrumentation pertain to smaller wave numbers ~10^-2 m^-1 appropriate to outer scale L ₀ of the atmospheric turbulence spectrum. Autocorrelation analysis of the fluctuations in temperature and wind speed has been performed. A new autoregressive scheme has been developed to represent the computed autocorrelation functions by a Yule statistical model, and to estimate the correlation period T ₀ of the turbulent medium. Height profiles of outer scale L ₀ of turbulence may be given in terms of T ₀ and mean wind speed u. Further, the similarity theory of Monin-Obukhov has been used to compute height profiles of temperature structure parameter C _T ². At Maitri, values of L ₀ and C _T ² are higher between 03–22 h local time than between 22–03 h. Values of L ₀ and C _T ² are smaller over the ice shelf terrain of the Dakshin Gangotri station, compared to those over the rocky terrain of the Maitri station.     

The response of a propeller anemometer to turbulent flow with the mean wind vector perpendicular to the axis of rotation

The system transfer function ¦H(v)¦² at frequencyv (units of Hz) for a vertical velocity propeller anemometer in a statistically stationary and horizontally homogeneous turbulent flow is determined from: (1) experimental estimates of propeller velocity spectra; and (2) estimates of Eulerian vertical velocity spectra based on the hypothesis that degradation of the input vertical velocity Fourier components occurs in the inertial subrange. The experimental estimates of ¦H(v)¦² were adequately summarized with the mathematical expression for the system transfer function of a first-order system with parameterT which has units of time and is analogous to the time constant of a horizontal velocity propeller anemometer. Dimensional analysis techniques and the Monin-Obukhov similarity hypothesis were used to construct a model for the system parameterT which yielded the result that _w /D ₁ ( _w /)^1/3, where _w , andD ₁ denote the standard deviation of the input vertical velocity fluctuations, the horizontal mean wind speed, and the diameter of the propeller, respectively. The system parameterT is interpreted in terms of the time required for the propeller velocity statistics to become asymptotically independent of time upon being released from rest in a statistically stationary turbulent flow.Currently on leave of absence from the Indian Institute of Technology, New Delhi, India.     

Multiscaling properties of concentration fluctuations in dispersing plumes revealed using an orthonormal wavelet decomposition

The dynamical characteristics of concentration fluctuations in a dispersing plume over the energetic and inertial-convective range of scales of turbulent motion are studied using a multiscale analysis technique that is based on an orthonormal wavelet representation. It is shown that the Haar wavelet concentration spectrum is similar to the Fourier concentration spectrum in that both spectra exhibit an extensive inertial-convective subrange spanning about two decades in frequency, with a scaling exponent of -5/3. Analysis of the statistical properties (e.g., fluctuation intensity, skewness, and kurtosis) of the concentration wavelet coefficients (i.e., the concentration discrete detailed signal) suggests that the small scales are always more intermittent than the large scales. The degree of intermittency increases monotonically with decreasing scale within the inertial-convective subrange, reaching a plateau at the very small scales associated with the beginning of the near-dissipation subrange. The probability density function (pdf) of the concentration discrete detailed signal displays stretched exponential tails with an intermittency exponent (tail slope) q that increases as ^a , where is the scale or dilation and a is a power-law exponent that is dependent on downwind distance, plume height, and stratification strength with typical values in the range from about 0.25 to 0.35. It is shown that the concentration variance cascade process requires a phase coherency of eddies between different scales at the small-scale end of the inertial-convective subrange.The variation of the concentration wavelet statistics with height above the ground is investigated. The increased mean shear near the ground smooths the fine-scale plume structure for scales within the inertial-convective subrange, producing a weaker spatiotemporal intermittency in the concentration field compared to that measured higher up in the plume. The pdf of the concentration detailed signal at a fixed scale possesses less elongated tails with decreasing height z. The intermittency exponent q is found to decrease roughly linearly with increasing z.Finally, the results of the wavelet decomposition are combined to provide a conceptual model of the turbulent transport, stirring, and mixing regimes in a dispersing plume. The implications of the results for contaminant texture in a plume are discussed.     

Estimation of abl Parameters Using the Vertical Velocity Measurements of an Acoustic Sounder

The friction velocity, the surface heat flux and the height of the Atmospheric Boundary Layer (ABL) are important parameters. In this work, vertical velocity variance ( _w ² ) and wind velocity structure parameter (C _v ² ) profiles estimated by acoustic sounder measurements are used, along with similarity relations, to estimate these parameters in the unstable Atmospheric Boundary Layer and the friction velocity in the stable one. The data were collected by two acoustic sounders with different height range and resolution under various atmospheric conditions (stability) and at two experimental sites in different terrain. The C _v ² profiles are estimated using gate difference of the vertical velocity measurements and the assumption of local isotropy. The vertical velocity data are corrected for the significant effects of noisy measurements and sampling volume averaging on the _w ² and C _v ² estimations using original techniques that are presented in this work. The results of the similarity method using acoustic sounder data are compared against estimates of the corresponding atmospheric parameters obtained from direct measurements. The comparison confirms the ability of the method to provide reasonably accurate estimates of these parameters especially in the middle of the day.     

Eddy flux measurements over the ocean and related transfer coefficients

Eddy correlation measurements of vertical turbulent fluxes made during AMTEX 1975 are used to assess the reliability of flux prediction from established bulk transfer relations, using both surface-layer and planetary boundary-layer formulations. The surface-layer formulae predict momentum and latent heat fluxes to an accuracy comparable to the direct eddy correlation method, using transfer coefficients of C _DN (at 10m and in neutral conditions) increasing with wind speed, and a constant C _EN - 1.5 × 10 ^–3 . The data suggest C _CHN , for sensible heat, increases significantly with wind speed and is on average 30% lower than C _CEN The boundary-layer drag coefficient, C _GD , agrees within about 40% of recently published values using a vertically averaged geostrophic wind to the height of the lowest temperature inversion, corrected for trajectory curvature. Values of _* / from which C _CGH is derived, are in excellent agreement if the published values are modified to account for inappropriate surface temperatures used in their derivation. Preliminary values of C _GE are also presented.