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.     

A model of turbulence spectra in the atmospheric surface layer

The spectral equations of turbulent kinetic energy and temperature variance have been solved by using Onsager's energy cascade model and by extending Onsager's model to closure of terms that embody the interaction of turbulent and mean flow.The spectral model yields the following results: In a stably stratified shear flow, the peak wave numbers of the spectra of energy and temperature variance shift toward larger wave numbers as stability increases. In an unstably stratified flow, the peak wave numbers of energy spectra move toward smaller wave numbers as instability increases, whereas the opposite trend is observed for the peak wave numbers of temperature variance spectra. Hence, the peak wave numbers of temperature spectra show a discontinuity at the transition from stable to unstable stratification. At near neutral stratification, both spectra reveal a bimodal structure.The universal functions of the Monin-Obukhov similarity theory are predicted to behave as _m ~ _H ~ (- Z/L)^-1/3 in an extremely unstable stratification and as _m ~ _H ~ z/L in an extremely stable stratification. For a stably stratified flow, a constant turbulent Prandtl number is expected.     

The Theoretical Model of Atmospheric Turbulence Spectrum in Surface Layer

It is shown that the slope of energy spectrum obtained from the velocity solution of Kdv-Burgers equation lies between -5/3 and -2 in the dilogarithmic coordinates paper. The spectrum is very close to one of Kolmogorov’s isotropic turbulence and Frisch’s intermittent turbulence in inertial region. In this paper, the Kdv-Burgers equation to describe atmospheric boundary layer turbulence is obtained. In the equation, the 1 / Re, corresponds to dissipative coefficient v, to dispersive coefficient β, then (v/ 2β)2 corresponds to .We prove that the wave number corresponding to maximum energy spectrum decreases with the decrease of stability (i.e., the increase of in eddy-containing region. And the spectrim amplitude decreases with the increase of (i.e., the decrease of stability). These results are consistent with actual turbulence spectrum of atmospheric surface layer from turbulence data.     

大气边界层湍流相干结构的识别

首先利用数字滤波方法对淮河流域试验的大气边界层湍流观测资料进行三项分解，将大气边界层湍流的风速信号分解为近似各项同性的小尺度涡和各向异性的大尺度涡。然后再将大尺度涡信号进行离散正交小波分解，寻求相干结构的主要特征尺度。对于大气边界层湍流垂直脉动风速来说，其相干结构的主要特征尺度为16s；对径向与纬向脉动来说，其相干结构的主要特征尺度为32－64s。在此基础上，利用小波的反变换提取出相干结构的信号与非相干结构的信号，并计算两间的相关系数，最大仅有0．02。此外，对原始大气湍流观测信号不进行数字滤波，直接利用本中子波分析法提取湍流相干结构所得结果作比较研究；并探讨了采用对称或似对称离散正交小波对此研究的影响。     

Resolved-scale turbulence in the atmospheric surface layer from a large eddy simulation

Large eddy simulation has encountered difficulties in handling turbulence in the atmospheric surface layer due to deficiencies in sub-grid scale models. This paper addresses the possibility of resolving the turbulence in the upper part of the surface layer by a low-aspect ratio of grid spacing. Results show that resolved-scale shear stresses dominate over the sub-grid scale components so that effects due to the sub-grid scale model can be ignored in this region. The effects of the lower boundary condition on the resolved-scale turbulence in the upper part of the surface layer are discussed. It is concluded that the normalized mean velocity shear and resolved turbulence in the upper part of the surface layer are not affected by the specification of the lower boundary condition. In addition, the present work proposes a new independent model parameter, the Smagorinsky Reynolds Number (Re_SM), and demonstrates that this number determines the resolved turbulence in the upper part of the surface layer.     

Experimental investigation of atmospheric boundary layer turbulence

This paper describes how to measure turbulence in the atmospheric boundary layer (ABL) in order to address certain problems in modern atmospheric physics. These problems mainly relate to the Earth's energy budget (including the hydrological cycle) and biogeochemical cycles. Starting from the main characteristic numbers and the basic equations of atmospheric turbulent flow, we show what turbulence parameters are important to measure. Special attention is given to the various methods used to compute the turbulent fluxes. We analyse the range of scales which has to be measured to properly capture the eddies contributing to the turbulent transfers. This range of scales determines what sensors can be used in the atmospheric surface layer and in the ABL. We describe the most widely used instruments and their performances. The principal platforms used to deploy these instruments are examined. Aircraft are described in more details, because they allow a thorough exploration of the ABL. In the last section, some examples of ABL turbulence signals measured in various conditions are presented. These examples illustrate horizontally homogeneous turbulence as well as inhomogeneous signals for which standard analysis techniques cannot be used. We show how some recent techniques, like wavelet transforms, can help to investigate this kind of signal. At the end, we present what would be interesting to do in the near future for the study of ABL turbulence.     

An analysis of wind velocity fluctuations in the atmospheric surface layer using an orthonormal wavelet transform

A wavelet analysis can supply information of both the location (time) and the scale of fluctuations. This method is applied to the fluctuations of the natural wind and the turbulent transport of momentum in the atmospheric surface layer. The shapes of both the wavelet spectra and the Fourier spectra of the three components of the wind velocity fluctuations are similar to each other. The quadrant representation of momentum transport shows the scale difference of the transport. The large-scale fluctuations mainly contribute to the downward transport of momentum.     

Thermal plumes and turbulence spectra in the atmospheric boundary layer

Experimental observations on the temperature and wind fields above flat grassy terrain have been obtained with an instrumented 92-m tower during intervals of strong insolation about midday. The turbulence characteristics of the air confirm that free convection prevailed at heights between 16 and 48 m, with some tendency for departure at higher levels. The spectra of temperature and vertical velocity contain gaps at wave numbers in the range 0.01–0.025 m^–1. These are attributed to natural thermal plumes that act as sources of extra energy input to the Kolmogorov-Obukhov-Corrsin scheme of turbulence in or at the low-wave number limit of the inertial subrange. Modified forms of the K-O-C spectral laws for thermally unstable air are derived which agree with the observed spectra over the whole range of wave numbers examined, and which contain the spectral gap at wave numbers corresponding to the thermal plume diameters.     

Analysis of atmospheric flow over a surface protrusion using the turbulence kinetic energy equation

Flow over surface obstructions can produce significantly large wind shears such that adverse flying conditions can occur for aeronautical systems (helicopters, STOL vehicles, etc.) Atmospheric flow fields resulting from a semi-elliptical surface obstruction in an otherwise horizontally homogeneous statistically stationary flow are modelled with the boundary-layer / Boussinesq-approximation of the governing equation of fluid mechanics. The turbulence kinetic energy equation is used to determine the dissipative effects of turbulent shear on the mean flow. Mean-flow results are compared with those given in a previous paper where the same problem was attacked using a Prandtl mixing-length hypothesis. The diffusion and convection of turbulence kinetic energy not accounted for in the Prandtl mixing-length concept cause departures of the mean wind profiles from those previously computed, primarily in the regions of strong pressure gradients. Iso-lines of turbulence kinetic energy and turbulence intensity are plotted in the plane of the flow. They highlight regions of high turbulence intensity in the stagnation zone and sharp gradients in intensity along the transition from adverse to favourable pressure gradient.     

Two-point coherence in the atmospheric surface layer

Two-point, one-dimensional coherence in horizontally homogeneous atmospheric turbulence is studied, both by experiment and analysis. Measurements are carried out using horizontally spaced sensors with the separation perpendicular to the mean velocity. Two-dimensional spectral models and three-dimensional inertial-range spectral tensors are used in the coherence calculations. The one-dimensional coherence for both velocity and scalar fluctuations is found to roll off at a wavenumber much smaller than we would expect from the classical notion of eddy correlation. This is a consequence of the cancellation of Fourier components aliased from the direction of the sensor separation into the streamwise direction. However, the coherence for the three velocity components behaves somewhat differently, reflecting the relative orientations of the velocity component, sensor separation and the mean velocity. These features are well predicted by the calculation. The analysis is also extended to calculate the two-point scalar-vertical velocity cospectrum and the results are in good agreement with our experimental data. The ratio of two- to one-point cospectra decreases at slightly larger wavenumber than the two-point scalar coherence does.     

Temperature structure in the atmospheric surface layer

Production, transport and dissipation terms in the temperature variance equation have been measured in the atmospheric surface layer. The transport term is, within the experimental uncertainty, negligible. The dissipation term, determined by assuming local isotropy, is approximately equal to production under near-neutral conditions. For moderately unstable conditions, the ratio of production to dissipation is 1.4. The resulting imbalance in the budget is attributed to the inequality between the three components of the dissipation term. The Kolmogorov constant for temperature is found to be about 0.8.     

Turbulent dispersion in the atmospheric surface layer

By non-dimensionalizing a trajectory-simulation (TS) model of turbulent dispersion, it is shown that the dimensionless concentration z ₀cu_*/kQ (cu _*/kQ) due to a continuous line (area) source of strength Q in the atmospheric surface layer depends only on z/z ₀, x/z ₀, z ₀/L and z _s/z₀, where z _s is the source height. The TS model is used to tabulate concentration profiles due to ground-level line and area sources. Concentration profiles generated by the TS model for elevated sources are shown to be inconsistent with the Reciprocal Theorems of Smith (1957) and it is suggested that this is because the flux-mean gradient closure scheme inherent in the Reciprocal Theorem is invalid for an elevated source.     

An investigation of the turbulence balance equations in the atmospheric boundary layer

Three Cardington turbulence probes were operated at heights of 61, 91, and 152 m over a 5-h period in an unstable boundary layer. The results have been used to assess (at the middle level) the relative importance of the terms in the turbulence equations for kinetic energy, temperature variance, stress, and the vertical and total horizontal heat fluxes.Now at the Building Research Station, Garston, Watford, England.     

Microscale temperature fluctuations in the atmospheric surface layer

Microscale temperature fluctuations were measured at 2 m above a grassy surface. The temperature-derivative spectrum was in general agreement with earlier results but the bump at nondimensional wavenumbers higher than 0.02 was not as pronounced as has been observed. The Obukhov-Corrsin constant for the one-dimensional temperature spectrum was evaluated to be 0.92 ± 0.05, consistent with recent results. The effects of instability and the vertical variation of temperature variance and kinetic energy dissipation are postulated to explain some of the difference with other spectra.     

On the analysis of wind wave-induced disturbances in the atmospheric turbulent surface layer

Simultaneous measurements of wind velocity, air humidity and sea surface wave-elevation fluctuations obtained on a platform in the open Caspian Sea are analyzed.It is shown that wave noises exist in the spectra and cospectra of the lower part of the atmospheric boundary layer, not only for the frequency of the main energy transporting component of the sea waves, but also in other parts of the spectra, at both lower and higher frequencies. The high frequencies are noncoherent with the sea waves and could be considered as measurement errors due to the existence of the waves. A method of elimination of the coherent wave noises from the spectra and cospectra is suggested and the effectiveness of its application is demonstrated.The essential difference between cases of developing and decaying wind waves is demonstrated.     

The atmospheric boundary layer over a wavy surface

The wavy area of north-west Bohemia (Czechoslovakia) is simulated by a cylindrical surface model. A curvilinear orthogonal system of coordinates along the model surface is introduced. The hydrodynamical equations of motion are transformed into this system of coordinates. By applying boundary-layer assumptions, the equations of motion for the atmospheric boundary layer (ABL) above the model are derived. The equations and boundary conditions show an equivalence of the ABL above the model with that above a flat surface with external pressure gradient.     

Measurements of the surface layer turbulence in the urban area

The turbulence characteristics measured in the surface layer over a real underlying surface are presented. The measurements are carried out at the 3- and 16-m levels using acoustic anemometers; at the lower level, the measurements are doubled. The diurnal cycle is obtained of the temperature flux and friction velocity. To obtain the temperature flux, so-called “acoustic” temperature is used as calculated from the sound speed data with corrections for air humidity and pressure. The normalized characteristics of turbulence are presented as dependent on dimensionless height z/L. The values of temperature flux and friction velocity increase with height: their day-integral values at 16-m level exceed those at 3-m level for about 20 and 35%, respectively.     

An examination of turbulence statistics in the surface boundary layer

Simulated data derived from random numbers are used to show that the process of relating % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiabeo8aZnaaBa% aaleaacaWG3baabeaakiaac+cacaWG1bWaaSbaaSqaaiabgEHiQaqa% baaaaa!3D7C!\[\sigma _w /u_ * \]and similar properties to the stability parameter % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiaadQhacaGGVa% Gaamitaaaa!3A42!\[z/L\]is highly susceptible to error. An alternative method, making use of Ri as a stability index, is not affected in this way and is used to re-examine the data obtained in the 1968 Kansas micrometeorological experiment. The relationship % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiabeo8aZnaaBa% aaleaacaWG3baabeaakiaac+cacaWG1bWaaSbaaSqaaiabgEHiQaqa% baGccqWIdjYocaaIXaGaaiOlaiaaikdacaaI1aaaaa!419F!\[\sigma _w /u_ * \simeq 1.25\] % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaamaabmaabaGaaG% ymaiabgkHiTiaaikdacaWG6bGaai4laiaadYeaaiaawIcacaGLPaaa% daahaaWcbeqaaiaaigdacaGGVaGaaG4maaaaaaa!4087!\[\left( {1 - 2z/L} \right)^{1/3} \]is found to provide a good fit to the unstable data, but it is unclear as to whether a small peak observed in stable conditions is real (perhaps associated with gravity waves) or not (possibly a consequence of measurement errors).The properties % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiabeo8aZnaaBa% aaleaacaWG1baabeaakiaac+cacaWG1bWaaSbaaSqaaiabgEHiQaqa% baaaaa!3D7A!\[\sigma _u /u_ * \]and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiabeo8aZnaaBa% aaleaacaWG1baabeaakiaac+cacaWG1bWaaSbaaSqaaiabgEHiQaqa% baaaaa!3D7A!\[\sigma _u /u_ * \] are found to attain a relatively constant value ( 3) in conditions more unstable than about Ri = -0.4. The shape ratio % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiabeo8aZnaaBa% aaleaacaWG1baabeaakiaac+cacqaHdpWCdaWgaaWcbaGaamODaaqa% baaaaa!3E4F!\[\sigma _u /\sigma _v \] is found to decrease to less than unity in very unstable conditions, possibly as a consequence of some undetected error in measurement of _u . In the case of temperature fluctuations, % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOqaaiabeo8aZnaaBa% aaleaacqaHepaDaeqaaOGaai4laiaadsfadaWgaaWcbaGaey4fIOca% beaakiabg2da9iaaicdacaGGUaGaaGyoaiaaiwdacaGGOaGaeyOeI0% IaamOEaiaac+cacaWGmbGaaiykamaaCaaaleqabaGaeyOeI0IaaGym% aiaac+cacaaIZaaaaaaa!4A30!\[\sigma _\tau /T_ * = 0.95( - z/L)^{ - 1/3} \] is found to provide an excellent fit in unstable conditions (Ri < -0.1), even though this form also agrees well with random behavior.Now With: Atmospheric Turbulence and Diffusion Laboratory, NOAA, P. O. Box E, Oak Ridge Tenn., 37830, U.S.A.