An Improvement of the Louis Scheme for the Surface Layer in an Atmospheric Modelling System

The Louis scheme for calculating the vertical eddy fluxes within the atmospheric surface layer is improved by broadening the original assumptions. In our approach, the momentum and heat transfer roughness lengths (z₀ and z_T respectively) can be different, and z₀ need not be negligibly small compared with the lowest height (z) in modelling. For these conditions, we choose more consistent wind and potential temperature profile forms, then derive new algorithms for calculating fluxes. Improvement is demonstrated for a wide range of z/L (L is the Obukhov length), z/z₀ and z₀ z_T, by comparing these fluxes with those derived from a theoretical surface-layer model. The improved algorithms can be used in atmospheric modelling systems for more varied surfaces and a wide range of atmospheric stability.     

A Study of Correlations of Scalar Quantities in the Atmospheric Surface Layer

Standard deviations for vertical velocity and scalar quantities, such as temperature, T, and specific humidity, q, were analyzed on the basis of Monin-Obukhov (M-O) similarity theory in the atmospheric surface layer. The correlation coefficient between scalar quantities T and q, R_Tq, was derived from the similarity functions and can be expressed as the ratio of B_T/B_q (B_T≤ B_q), where parameter B is the value of the normalized standard deviation of any scalar quantity at neutral conditions.     

Atmospheric Turbulence Characteristics Over a Temperature-Inhomogeneous Land Surface. Part I: Statistical Characteristics of Small-Scale Spatial Inhomogeneities of Land Surface Temperature

The results of our 1993 and 1994 experiments at the Tsimlyansk experimental site are given. Spatial variations of small-scale (0.1 m < l < 100 m) temperature inhomogeneities of a steppe surface have been measured. Considerable variations in surface radiation temperature T_r, with variances up to 5 K², have been found. The spectra of T_r variations are of a one-mode character. The position of the spectral maximum on the wavenumber scale depends on the surface characteristics and may change from year-to-year. Probability density distributions for steppe surface radiation temperature variations are close to a normal distribution.     

Aspects Of The Atmospheric Surface Layers On Mars And Earth

The structures of mean flow and turbulence in the atmospheric surface boundary layer have been extensively studied on Earth, and to a far less extent on Mars, where only the Viking missions and the Pathfinder mission have delivered in-situ data. Largely the behaviour of surface-layer turbulence and mean flow on Mars is found to obey the same scaling laws as on Earth. The largest micrometeorological differences between the two atmospheres are associated with the low air density of the Martian atmosphere. Together with the virtual absence of water vapour, it reduces the importance of the atmospheric heat flux in the surface energy budget. This increases the temperature variation of the surface forcing the near-surface temperature gradient and thereby the diabatic heat flux to higher values than are typical on the Earth, resulting in turn in a deeper daytime boundary layer. As wind speed is much like that of the Earth, this larger diabatic heat flux is carried mostly by larger maximal values of T_*, the surface scale temperature. The higher kinematic viscosity yields a Kolmogorov scale of the order of ten times larger than on Earth, influencing the transition between rough and smooth flow for the same surface features.The scaling laws have been validated analysing the Martian surface-layer data for the relations between the power spectra of wind and temperature turbulence and the corresponding mean values of wind speed and temperature. Usual spectral formulations were used based on the scaling laws ruling the Earth atmospheric surface layer, whereby the Earth's atmosphere is used as a standard for the Martian atmosphere.     

Effects Of Changing Surface Heat Flux On Atmospheric Boundary-Layer Flow Over Flat Terrain

We examine the unsteady response of a neutral atmospheric boundary layer (ABL) of depth h and friction velocity u _* when a uniform surface heat flux is applied abruptly or decreased rapidly over a time scale t<inf>θ</inf> less than about h /(10u _*). Standard Monin–Obukhov (MO) relationships are used for the perturbed eddy viscosity profile in terms of the changes to the heat flux and mean shear. Analytical solutions for changes in temperature, mean wind and shear stress profile are obtained for the surface layer, when there are small changes in h /|L_MO| over the time scale t_MO~|L _MO|/(10u_*) (where L _MO and t _MO are the length and time scales, respectively). They show that a maximum in the wind speed profile occurs at the top of the thermal boundary layer for weak surface cooling, i.e. a wind jet, whereas there is a flattening of the profile and no marked maximum for weak surface heating. The modelled profiles are approximately the same as those obtained from the U.K. Met Office Unified Model when operating as a mesoscale model at 12-km horizontal resolution. The theoretical model is modified when strong surface heating is suddenly applied, resulting in a large change in h /|L _MO| (>>1), over the time scale t _MO. The eddy structure is predicted to change significantly and the addition of convective turbulence increases the shear turbulence at the ground. A low-level wind jet can form, with convective turbulence adding to the mean momentum of the flow. This was verified by our laboratory experiment and direct numerical simulations. Additionally, it is shown that the effects of Coriolis acceleration diminish (rather than as suggested in the literature, amplify) the formation of the wind jets in the situations considered here. Hence, only when the surface heat flux changes over time scales greater than 1/f (where f is the Coriolis parameter) does the ABL adjust monotonically between its equilibrium states. These results are also applicable to the ABL passing over spatially varying surface heat fluxes.     

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.     

南京地区大气颗粒物影响近地面臭氧的个例研究

通过对2008年4月2～7日南京地区地面气象观测数据以及两个站点空气质量(O₃、NO_x、PM10)监测资料的分析, 发现O₃和PM10之间存在一定程度的反相关。利用一个光化学箱模式对该个例中大气颗粒物影响近地面臭氧的过程进行模拟, 结果发现大气颗粒物浓度的升高使得气溶胶光学厚度增加20％～40%, 导致NO₂和O₃近地面光解率下降20%～30%, OH和HO₂自由基浓度分别减少20%～50%, 造成O₃净生成率下降30%～40%。研究表明, 颗粒物对光化学过程的抑制造成了大气氧化能力的降低, 是近地面臭氧浓度减少的可能原因。     

Numerical Studies with a Regional Atmospheric Climate Model Based on Changes in the Roughness Length for Momentum and Heat Over Antarctica

A regional atmospheric climate model is used toexamine the effect of changes in the roughnesslengths of momentum (_z0m) and heat (_z0h)on the structure of the lower atmosphere and on thesurface energy fluxes over Antarctica. Fourexperiments were carried out in which _z0mand/or _z0h were altered with respect to acontrol experiment. The changes consisted of (1) alowering of _z0m from a field aggregated froma vegetation map with an orographic correction basedon the European Centre for Medium-Range WeatherForecasts _z0m field, to a constant value of10^-3 m; and (2) a lowering of _z0h from a valueequal to _z0m to a constant value of 10^-3 mor a value dependent on the wind speed via a surfacerenewal model. A reduction of _z0m results in theexpected increase in near-surface wind speed. It alsoresults in an increase in the depth of the layer in whichsouth-easterly near-surface winds prevail, and in adecrease in the strength of the large-scale flow overthe continent, in particular in summer. In theescarpment region a decrease of _z0m is foundto result in too high wind speeds. Surface temperatureson average decrease while atmospheric temperaturesincrease, resulting in an increase of near-surfacestatic stability. Changes in roughness lengths donot significantly change the temperature profiles.The surface fluxes, on average found reduced, aremodelled best by using the _z0h based on thesurface renewal method.     

Variance Method to Determine Turbulent Fluxes of Momentum And Sensible Heat in The Stable Atmospheric Surface Layer

Evidence is presented that in the stable atmospheric surface layer turbulent fluxes of heat and momentum can be determined from the standard deviations of longitudinal wind velocity and temperature, σ _u and σ _T respectively, measured at a single level. An attractive aspect of this method is that it yields fluxes from measurements that can be obtained with two-dimensional sonic anemometers. These instruments are increasingly being used at official weather stations, where they replace the standard cup anemometer–wind vane system. With methods such as the one described in this note, a widespread, good quality, flux network can be established, which would greatly benefit the modelling community. It is shown that a ‘variance’ dimensionless height (ζ _σ) defined from σ _u and σ _T is highly related to the ‘conventional’ dimensionless stability parameter ζ=z/L, where z is height and L is the Obukhov length. Empirical functions for ζ _σ are proposed that allow direct calculation of heat and momentum fluxes from σ _u and σ _T. The method performs fairly well also during a night of intermittent turbulence.     

Model of convective heat exchange due to isolated thermals in the atmospheric boundary layer

A numerical model of convective heat transfer due to isolated thermals in the atmospheric boundary layer is used to describe the temperature profile transformation in undisturbed conditions as a result of intensive dry free convection. Based on some assumptions, the heat transfer Equation (2) is transformed to the form (14) in which the coefficients and the function F are expressed by (d/dz)(ln ) and by parameters of thermals. Equation (14) has been solved numerically with the help of Equation (15) obtained from the statics equation because of Equation (8). The size distribution function f(z, r, t) of the thermals is discrete (Table I), according to Vulf'son (1961). On Figures 1 and 2 are plotted successive temperature profiles for a ground inversion, transformed due to free convection and turbulence (Figures 1a and 2a), and due to turbulence only (Figures 1b and 2b). The profiles are computed from Equation 14 (Figures 1a and 2a) and Equation 16 (Figures 1b and 2b) for k _z= 1 m² s^–1 (Figure 1) and k _z= 10 m² s^–1 (Figure 2). On Figure 3 the real temperature profiles in Sofia for June 22nd 1976 are compared with the profiles computed using the real initial profile for 4.30 h local time. Good qualitative agreement can be seen.     

山谷城市的近地层大气湍流谱特征

山谷城市地区的大气湍流结构研究至今报道很少.本文根据兰州市一次观测实验,对各风速分量谱、温度谱以及动量和感热通量的协谱进行了分析.结果表明,尽管仪器的安装高度较低(6.45m),湍谱特征仍与平坦均匀下垫面得到的一些典型结果基本一致,从而认为,在类似地区进行湍流通量等的直接观测时,对观测高度和地点的选择,可以不必象一些文献要求的那样高.     

Turbulence Structure of the Unstable Atmospheric Surface Layer and Transition to the Outer Layer

We present a new model of the structure of turbulence in the unstable atmospheric surface layer, and of the structural transition between this and the outer layer. The archetypal element of wall-bounded shear turbulence is the Theodorsen ejection amplifier (TEA) structure, in which an initial ejection of air from near the ground into an ideal laminar and logarithmic flow induces vortical motion about a hairpin-shaped core, which then creates a second ejection that is similar to, but larger than, the first. A series of TEA structures form a TEA cascade. In real turbulent flows TEA structures occur in distorted forms as TEA-like (TEAL) structures. Distortion terminates many TEAL cascades and only the best-formed TEAL structures initiate new cycles. In an extended log layer the resulting shear turbulence is a complex, self-organizing, dissipative system exhibiting self-similar behaviour under inner scaling. Spectral results show that this structure is insensitive to instability. This is contrary to the fundamental hypothesis of Monin--Obukhov similarity theory. All TEAL cascades terminate at the top of the surface layer where they encounter, and are severely distorted by, powerful eddies of similar size from the outer layer. These eddies are products of the breakdown of the large eddies produced by buoyancy in the outer layer. When the outer layer is much deeper than the surface layer the interacting eddies are from the inertial subrange of the outer Richardson cascade. The scale height of the surface layer, z _s, is then found by matching the powers delivered to the creation of emerging TEAL structures to the power passing down the Richardson cascade in the outer layer. It is z _s = u _* ³ /k_s, where u _* is friction velocity, k is the von Kármán constant and _s is the rate of dissipation of turbulence kinetic energy in the outer layer immediately above the surface layer. This height is comparable to the Obukhov length in the fully convective boundary layer. Aircraft and tower observations confirm a strong qualitative change in the structure of the turbulence at about that height. The tallest eddies within the surface layer have height z _s, so z _s is a new basis parameter for similarity models of the surface layer.     

卫星遥感结合气象资料计算的青藏高原地面感热特征分析

本文选取1981年7月至2012年12月美国国家航空和航天局（NASA）制作的归一化的动态植被指数（NDVI）资料、根据NDVI值计算地表热力输送系数（C_H）的参数化关系式（C_H-I_NDV）和青藏高原70个常规气象观测资料,计算了青藏高原全区的逐月地表热力输送系数（C_H）,讨论了其时空分布特征,并在此基础上计算了高原70个常规台站的感热通量（SSHF）序列,并与已有感热资料进行了对比。随后,探讨了地面感热通量的气候特征及其年际变化与气候因子的关系。结果表明:高原地区的C_H值具有明显的空间差异和季节差异,表现为东高西低、夏季大、冬季小的特点。感热的年际变化在冬季主要响应于地气温差的变化,夏季则受地面风速影响较大;由于风速减小趋缓,地气温差增大,变化趋势在2003年前后由减弱趋势转变为增强趋势,这种趋势的转变最早发生在2001年秋季,且在高原全区具有较好的一致性。     

我国东西部地区地气温差的年代际变化特征

利用1960~2006年我国地温、气温逐日4个时次[02:00(北京时间,下同)、08:00、14:00和20:00]的台站观测资料,计算并分析了我国东南、西北地区各季地气温差的年代际变化特征。分析结果表明:我国东南部地区各季地气温差在20世纪70年代末以前,大部分年份偏高,高于平均值,而在20世纪70年代末以后,我国东南部地区各季地气温差偏低,在夏季和冬季表现尤为明显。我国西北地区春季和夏季地气温差在20世纪70年代末以前大部分年份偏低,低于平均值;而在20世纪70年代末以后,地气温差则大部分年份明显偏高。我国西北地区秋季地气温差的年代际变化特征不明显,而冬季地气温差的年代际变化趋势与春夏季相反,在20世纪70年代末以前大部分年份偏高,高于平均值,而在20世纪70年代末以后偏低。另外,发现地温和气温对我国东南、西北地区各季地气温差的年代际变化在各季所起的贡献作用不同。     

An Approximate Analytical Solution Of Flux-Profile Relationships For The Atmospheric Surface Layer With Different Momentum And Heat Roughness Lengths

An approximate method for calculating the relationship between z/L(z = reference height, L = Obukhov length) and the bulk Richardsonnumber is presented. If this relationship is known, the momentum andheat fluxes can be computed easily without any iteration. The avoidance of iteration can speed up computationsin large-scale models considerably (up to 10 times) and cases which do not converge or converge very slowly cannot occur. The proposed formulae take into account the difference between momentum (z_0M) and heat roughnesslengths (z_0H). Because the roughness lengths are not neglected at any step of the derivation, the resulting analytical formulae can be used not only between the surface and the reference height but also between two finite levels z₁ andz₂ (by replacing z_0M and z_0H by z₁ and z by z₂). Theequations remain correct even in the limit z₁ z₂.The formulae are based upon the (partially modified) Businger–Dyer flux–profile relationships and,consequently, they are restricted to predominantly homogeneous terrain.These new approximations are an improvement over the existing solutions because they are simpler than most of the formulae in the literature and are able to match the numerical exact solution for different parameter sets (Businger, Dyer, Högström) with an maximum error of about 2% for a wide range of z/L, z/z_0M and z_0M/z_0H.Furthermore, in stable conditions, schemes with and without a finitecritical bulk Richardson number can be approximated. The possibleambiguity of the exact solution =f(RI_B) in (moderately) stable conditions is discussed briefly. The performance of the new formulae is compared to the exact numerical solution and to different formulae proposed in the literature.     

Turbulence Characteristics of the Stable Boundary Layer Over a Mid-Latitude Glacier. Part I: A Combination of Katabatic and Large-Scale Forcing

Profile and eddy-correlation (heights of 4 and 10 m) measurements performed on the Pasterze glacier (Austria) are used to study the characteristics of the stable boundary layer under conditions of katabatic and large-scale forcing. We consider cases where large-scale forcing results in a downslope (or following) ambient wind. The analysis of averaged spectra and cospectra reveals low frequency perturbations that have a large influence on the variances of temperature and horizontal wind components and also alter the cospectra of momentum and sensible heat flux. Only the spectrum of the vertical wind speed is comparable to universal spectra. The low frequency perturbations occur as brief intermittent events and result in downward entrainment of ambient air thereby producing enhanced downward sensible heat fluxes and downward as well as upward momentum fluxes with various magnitudes and timescales. After the variances were high pass filtered, the normalised standard deviations of wind speed and temperature compare favourably to findings in the literature within the range 0>z/L>0.5. For larger z/L they deviate as a result of an increased influence from low frequency perturbations and thus non-stationarity. In line with this, the turbulent kinetic energy budget (at 4 m height) indicates that production (shear) is in balance with destruction (buoyancy and dissipation) within the range 0>z/L>0.3. Non-dimensional gradients of wind speed within the range 0>z/L>0.3 have a slope of about 3.5. The scatter for the dimensionless temperature gradient is quite large, and the slope is comparable to that for wind speed gradients. For z/L>0.3 the imbalance in the turbulent kinetic energy budget grows and non-dimensional gradients for wind speed and temperature deviate considerably from accepted values as a result of increased non-stationarity. Average roughness lengths for momentum and sensible heat flux derived from wind speed and temperature profiles are respectively 1 × 10^-3 m and 6 × 10^-5 m, consistent with the literature. The ratio (z_0h/z_0m) compares to those predicted by surface renewal models. A variation of this ratio with the roughness Reynolds number is not indicated by our data.     

Detecting Near-Surface Coherent Structure Characteristics Using Wavelet Transform with Different Meteorological Elements

In the present study, three wavelet basis functions (Mexican-hat, Morlet, and Wave) were used to analyze atmospheric turbulence data obtained from an eddy covariance system in order to determine effect of six meteorological elements (three-dimensional wind speed, temperature, and CO2 and H2O concentrations) on the time scale of coherent structures. First, we used the degree of correlation between original and reconstructed waveforms to test the three wavelets’ performance when determining the time scale of coherent structures. The Wave wavelet’s reconstructed coherent structure signal best matched the original signal; thus, it was used for further analysis of the time scale, number, and time cover of the meteorological elements. We found similar results for all elements, though there was some internal variation, suggesting that coherent structures are not inherently dependent on these elements. Our results provide a basis for proper coherent structure detection in atmospheric turbulence and improve the understanding of similarities and differences between coherent structure characteristics of different meteorological elements, which is helpful for further research into atmospheric turbulence and boundary layers.     

Effects of Topographical Slope Angle and Atmospheric Stratification on Surface-Layer Turbulence

The effect of topographical slope angle and atmospheric stratification on turbulence intensities in the unstably stratified surface layer have been parameterized using observations obtained from a three-dimensional sonic anemometer installed at 8 m height above the ground at the Seoul National University (SNU) campus site in Korea for the years 1999–2001. Winds obtained from the sonic anemometer are analyzed according to the mean wind direction, since the topographical slope angle changes significantly along the azimuthal direction. The effects of the topographical slope angle and atmospheric stratification on surface-layer turbulence intensity are examined with these data. It is found that both the friction velocity and the variance for each component of wind normalized by the mean wind speed decrease with increase of the topographical slope angle, having a maximum decreasing rate at very unstable stratification. The decreasing rate of the normalized friction velocity (u _* /U) is found to be much larger than that of the turbulence intensity of each wind component due to the reduction of wind shear with increase in slope angle under unstable stratification. The decreasing rate of the w component of turbulence intensity (σ _w /U) is the smallest over the downslope surface whereas that of the u component (σ _u /U) has a minimum over the upslope surface. Consequently, σ _w /u _* has a maximum increasing rate with increase in slope angle for the downslope wind, whereas σ _u /u _* has its maximum for the upslope wind. The sloping terrain is found to reduce both the friction velocity and turbulence intensity compared with those on a flat surface. However, the reduction of the friction velocity over the sloping terrain is larger than that of the turbulence intensity, thereby enhancing the turbulence intensity normalized by the friction velocity over sloping terrain compared with that over a flat surface.     

Variability Of The Stable And Unstable Atmospheric Boundary-Layer Height And Its Scales Over A Boreal Forest

Radiosondes releases during the NOPEX-WINTEX experiment carried out in late winter in Northern Finland were analysed for the determination of the height h of the atmospheric boundary layer. We investigate various possible scaling approaches, based on length scales using micrometeorological turbulence surface measurements and the background atmospheric stratification above h. Under stable conditions, the three previously observed turbulence regimes delineated by values of z/L (L is the Obukhov length) appears as a blueprint for understanding the departures found for the suitability of the Ekman scaling based on L_E = u_/f (u is the friction velocity and f the Coriolis parameter). The length scale L_N = u_/N (where N is the Brunt–Väisälä frequency) appears to be a useful scale under most stable conditions, especially in association with L. Under unstable conditions, shear production of turbulence is still significant, so that the three scales L, L_N and L_E are again relevant and the dimensionless ratios _N = L_N/L and L_N/L_E = N/f describe well the WINTEX data. Furthermore, in the classical scaling framework, the unstable domain may also be divided into three regimes as reflected by the dependence ofu_/f on instability (z/L).     

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.