Observations and parametrization of momentum transfer in heterogeneous terrain consisting of regularly spaced obstacles

Tethered balloon data collected over a heterogeneous, semi-rural area of southern England are presented. The terrain within which the observations were taken is dominated by regularly spaced bluff obstacles that form the boundaries to grass fields. The data are compared with previous surface data to assess how representative these fixed-point measurements are of area-averaged values. It is found that there is a vertical height above which the flow approximates to being homogeneous and observations made above it are representative of the area average. The effective roughness length for momentum is then compared with that derived from a model describing the drag partition around isolated obstacles that are either non-porous or porous to the incident flow. It is shown that the method chosen to partition the drag is of the right order to produce reasonable predictions of the effective roughness length for different roughness densities.     

Measurements Of Thermal Updraft Intensity Over Complex Terrain Using American White Pelicans And A Simple Boundary-Layer Forecast Model

An examination of boundary-layer meteorological and avian aerodynamic theories suggests that soaring birds can be used to measure the magnitude of vertical air motions within the boundary layer. These theories are applied to obtain mixed-layer normalized thermal updraft intensity over both flat and complex terrain from the climb rates of soaring American white pelicans and from diagnostic boundary-layer model-produced estimates of the boundary-layer depth z_i and the convective velocity scale w_*. Comparison of the flatland data with the profiles of normalized updraft velocity obtained from previous studies reveals that the pelican-derived measurements of thermal updraft intensity are in close agreement with those obtained using traditional research aircraft and large eddy simulation (LES) in the height range of 0.2 to 0.8 z_i. Given the success of this method, the profiles of thermal vertical velocity over the flatland and the nearby mountains are compared. This comparison shows that these profiles are statistically indistinguishable over this height range, indicating that the profile for thermal updraft intensity varies little over this sample of complex terrain. These observations support the findings of a recent LES study that explored the turbulent structure of the boundary layer using a range of terrain specifications. For terrain similar in scale to that encountered in this study, results of the LES suggest that the terrain caused less than an 11% variation in the standard deviation of vertical velocity.     

Forecasting weather radar propagation conditions

Summary The increasing use of weather radar quantitative precipitation estimates, particularly in automatic applications such as operational hydrometeorological modelling or assimilation in numerical weather prediction (NWP) models, has promoted the development of quality control procedures on radar data. Anomalous propagation (AP) of the radar beam due to deviation from the standard refractivity vertical profile, is one of the factors that may affect seriously the quality of radar observations because of the increase in quantity and intensity of non-precipitating clutter echoes and consequent contamination of the estimated rainfall field. Another undesired effect of AP is the change in the expected radar echo height, which may be relevant when correcting for beam blockage in radar rainfall estimation in complex terrain. The aim of this paper is to study the use of NWP mesoscale forecasts to predict and monitor AP events. A nested 15-km grid resolution version of the MASS model has been used to retrieve refractivity profiles in the coastal area of Barcelona, near a weather radar and a radiosonde station. Using the refractivity profiles two different magnitudes were computed: the vertical refractivity profile of the lowest 1000 m layer and a ducting index which describes the existence and intensity of the most super-refractive layer contained in the lowest 3-km layer. A comparison between model forecasts and radiosonde diagnostics during a six-month period showed that the model tended to underestimate the degree of super-refraction, with a bias of 4 km⁻¹ and RMSE of 11 km⁻¹ in the 1-km vertical refractivity gradient. Further analysis of the data showed that a combination of previous observations and forecasts allowed to produce modified forecasts improving the original direct model output, decreasing substantially the bias, reducing the RMSE by 20% and improving the skill by 40%, beating also radiosonde observations persistence.     

Scalar-Flux Similarity in the Layer Near the Surface Over Mountainous Terrain

The scaled standard deviations of temperature and humidity are investigated in complex terrain. The study area is a steep Alpine valley, with six measurement sites of different slope, orientation and roughness (i-Box experimental site, Inn Valley, Austria). Examined here are several assumptions forming the basis of Monin–Obukhov similarity theory (MOST), including constant turbulence fluxes with height and the degree of self-correlation between the involved turbulence variables. Since the basic assumptions for the applicability of the MOST approach—horizontally homogeneous and flat conditions—are violated, the analysis is performed based on a local similarity hypothesis. The scaled standard deviations as a function of local stability are compared with previous studies from horizontally homogeneous and flat terrain, horizontally inhomogeneous and flat terrain, weakly inhomogeneous and flat terrain, as well as complex terrain. As a reference, similarity relations for unstable and stable conditions are evaluated using turbulence data from the weakly inhomogeneous and flat terrain of the Cabauw experimental site in the Netherlands, and assessed with the same post-processing method as the i-Box data. Significant differences from the reference curve and also among the i-Box sites are noted, especially for data derived from the i-Box sites with steep slopes. These differences concern the slope and the magnitude of the best-fit curves, illustrating the site dependence of any similarity theory.     

Remote sensing of cloud liquid water

Summary A method is presented to infer cloud liquid water path (LWP in kg/m²) over the ocean from passive microwave measurements of SSM/I. The algorithm to retrieve LWP is based on simulated satellite observations. They are calculated with a radiative transfer model applied to about 3000 radiosonde ascents over the Atlantic Ocean. Since radiosonde observations do not contain direct information about cloud water and ice, these parameters are parameterized based on relative humidity and temperature using modified adiabatic liquid water density profiles. A multiple linear regression is applied to the simulated radiances and the calculated LWP to derive the algorithm. The retrieval accuracy based on the regression analysis including instrumental noise is 0.03 kg/m². Validation of the LWP-algorithm was pursued through a comparison with measurements of a ground-based 33 GHzmicrowave radiometer on board of R.V. Poseidon during the International Cirrus Experiment 1989 at the North Sea (ICE'89). The LWP values agree within the range of uncertainty caused by the different sampling characteristics of the observing systems. The retrieval accuracy for clear-sky cases determined using colocated METEOSAT data over the North Sea is 0.037 kg/m² and confirms the accuracy estimated from regression analysis for the low liquid water cases.The algorithm was used to derive maps of monthly mean LWP over the Atlantic Ocean. As an example the Octobers of the 5 years 1987–1991 were selected to demonstrate the interannual variability of LWP. The results were compared with the cloud water content produced by the climate model ECHAM-T2 from the Max-Planck-Institut Hamburg.Observations during ICE'89 were used to check the accuracy of the applied radiative transfer model. Brightness temperatures were calculated from radiosonde ascents launched during the overpass of DMSP-F8 in cloud-free situations. The channel-dependent differences range from about –2 to 3 K.The possibility to identify different cloud types using microwave and infrared observations was examined. The main conclusion is that simultaneous microwave and infrared measurements enable the separation of dense cirrus and cirrus with underlying water clouds. A classification of clouds with respect to their top heights and LWP was carried out using a combination of SSM/I derived LWP and simultaneously recorded Meteosat IR-data during ICE'89.With 11 Figures     

Stability Dependence of Height Scales and Effective Roughness Lengths of Momentum and Heat Transfer Over Roughness Changes

It is widely accepted that the correct formulation of an effective roughness length, defined as the area average of the roughness length in heterogeneous terrain, relies upon the appropriate de-termination of a height scale. At this height a meteorological quantity is approximately in equilibrium with local surface conditions and independent of horizontal position. This research note determines explicitly the different height scales from the perturbation solutions of flow velocity and temperature, as well as the fluxes of momentum and heat, in a stratified boundary layer. These solutions are derived from the asymptotic approximation theory and shown to capture major characteristics of momentum and heat transfer over heterogeneous terrain with changes of the underlying roughness lengths. The effective roughness lengths can then be computed by use of these height scales. The dependence of height scales and effective roughness lengths upon stratification is also discussed briefly.     

上海徐家汇地区台风“灿鸿”与良态风风速特征

上海徐家汇地区建筑分布密集且高低不一,是典型的具有非均一下垫面特征的城市地貌。本文利用该地区的地理信息研究了粗糙度长度的分布规律,并基于80m高度的风速实测数据,对台风"灿鸿"和良态风作用下的平均风速、湍流强度、阵风因子等参数与粗糙度长度之间的关系进行了分析。结果表明:不同风向角对应的计算扇区内建筑物高度、分布密度的差别导致了粗糙度长度值随风向角发生明显变化,但是变化幅度随着计算扇区的增大而减小;台风"灿鸿"作用下的平均风速最大值大于良态风,两者对应的粗糙度长度变化范围差别甚微,但是台风作用下的粗糙度长度中位数较小且分布相对集中;台风"灿鸿"作用下各向湍流强度均随着平均风速的增加呈明显的减小趋势,但不随粗糙度长度变化;良态风作用下,各向湍流强度不随平均风速变化,而随着粗糙度长度的增加而增加;台风"灿鸿"和良态风作用下,各向阵风因子均随湍流强度的增加而增大,但前者作用下的阵风因子略大于后者。     

Growth Of The Summer Daytime Convective Boundary Layer At Anand

The heights of the daytime convective boundary layer (CBL), computed by a one-dimensional model for a bare soil surface at a semi-arid station,Anand, during the dry and hot summer month of May 1997, are presented. As input, the model requires surface heat flux, friction velocity and air temperature as functions of time. Temperature data at the one-metre level from a tower and sonic anemometer data at 9.5 m collected during the period 13–17 May 1997 in the Land Surface Processes Experiment (LASPEX-97) are used to compute hourly values of surface heat flux, friction velocity and Obukhov length following the operational method suggested by Holtslag and Van Ulden [J. Climate Appl. Meteorol. 22,517–529 (1983)]. The model has been tested with different values for the potential temperature gradient ( ) above the inversion. The model-estimated CBL heights comparefavourably with observed heights obtained from radiosonde ascents.     

Near-Wall Flow Development After A Step Change In Surface Roughness

An experimental study of the initial flow field downstream of a step change in surface roughness is presented. The roughness length of the downstream surface was approximately tenfold that of the upstream roughness and, unlike all previous studies, attention was concentrated on the roughness sublayer region beneath the inertial (log-law) region. The experiments were conducted at a boundary layer Reynolds number of about 6 × 10⁴ (based on layer thickness andfree-stream velocity) and around a longitudinal location where the (downstream) roughness length, z_o2, was about 1% of the boundary-layer thickness atthe roughness change point.The thickness of the roughness sublayer was found for the two roughness. It was observed that the vertical profiles of mean velocity and turbulence characteristics started to show similarity after about 160z₀₂ downstream of the roughness change. The presence of a shear stress overshoot is shown to depend strongly on the precise location (with respect to the roughness elements) at which the measurements are made and the thickness of the equilibrium layer is shown to be very sensitive to the way it is defined. It is demonstrated that the growing equilibrium layer has first to encompass the roughness sublayer before any thickness of inertial sublayer can be developed. It follows that, in somepractical cases, like flows across some urban environments, the latter(log-law) region may never exist at all.     

Airflow from mustard to fallow

Real terrain is complex and its effects on the atmosphere are not well understood. The two‐dimensional change‐of‐roughness problem is investigated in an experimental study of airflow from mustard to fallow. The surface friction velocity is found to overshoot and return slowly to the equilibrium value. The wind profile in the modified region is more nearly logarithmic than the Peterson model predicts. The simple Elliot theory is in fair agreement with the wind profile data.

Values for the ratio of the standard deviation of the vertical wind to the friction velocity reveal no variation with stability and are in close agreement with the value for laboratory flows. Values for the ratio of the standard deviation of the longitudinal wind component to the friction velocity are comparable to other recently determined atmospheric values and show a variation attributable to large‐scale terrain properties.     

Charnock’s Roughness Length Model and Non-dimensional Wind Profiles Over the Sea

An analysis tool for the study of wind speed profiles over the water has been developed. The profiles are analysed using a modified dimensionless wind speed and dimensionless height, assuming that the sea surface roughness can be predicted by Charnock’s roughness length model. In this form, the roughness dependency on wind speed is extracted and the variations on the wind profile are due solely to atmospheric stability. The use of the Charnock’s non-dimensional wind profile is illustrated using data collected from a meteorological mast installed in the Danish North Sea. The best fit with the observed mean non-dimensional wind profile under neutral atmospheric conditions is found using a value of 1.2 × 10⁻² for Charnock’s parameter. The stability correction on the neutral wind profile suggested by the Businger-Dyer relations was found to perform well over the sea.     

Near Wall Flow over Urban-like Roughness

In this study, comprehensive measurements over a number of urban-type surfaces with the same area density of 25% have been performed in a wind tunnel. The experiments were conducted at a free stream velocity of 10 m s^-1 and the main instrumentation was 120 ° x-wire anemometry, but measurement accuracy was checked using laser Doppler anemometry.The results haveconfirmed the strong three-dimensionalityof the turbulent flow inthe roughness sublayer and the depths of the inertial sublayer (log-law region) and roughness sublayer for each surface have been determined. Spatial averaging has been used to remove the variability of the flow in the roughness sublayer due to individual obstacles and it is shown that the spatially averaged mean velocity in the inertial sublayer and roughness sublayer can,together, be described by a single log-law with a mean zero-plane displacement and roughness length for the surface, provided that the proper surface stress is known. The spatially averaged shear stresses in the inertial sublayer and roughness sublayer are compared with the surface stress deduced from form drag measurements on the roughness elements themselves.The dispersive stress arising from the spatial inhomogeneity in the mean flow profiles was deduced from the data and is shown to be negligible compared with the usual Reynolds stresses in the roughness sublayer. Comparisons have been made between a homogeneous (regular element array) surface and one consisting of random height elements of the same total volume. Although the upper limits of the inertial sublayer for both surfaces were almost identical at equivalent fetch, the roughness sublayer was much thicker for the random surface than for the uniform surface, the friction velocity and the roughness length were significantly larger and the `roughness efficiency' was greater. It is argued that the inertial sublayer may not exist at all in some of the more extreme rough urban areas. These results will provide fundamental information for modelling urban air quality and forecasting urban wind climates.     

é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.     

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.     

Roughness lengths for temperature and momentum over heterogeneous terrain

Measurements are presented describing the behaviour of roughness lengths for momentum and temperature for heterogeneous terrain of low fractional rough cover. The momentum roughness lengths were found to be up to several orders of magnitude larger than that for temperature, behaviour which is characteristic of a bluff-rough type of surface and which is consistent with previous experimental and theoretical studies. This contrasts with observations over uniform, vegetated surfaces, which show only an order-of-magnitude increase of momentum roughness length over that for temperature.     

应用印痕函数研究测风塔资料代表性

影响测风塔代表性的因素很多,为了合理地确定测风塔观测结果能够代表的区域范围,需要通过适当的方法综合考虑这些因素。湍流通量代表性研究中常用的印痕函数是风速、稳定度、下垫面粗糙度等的函数,可以定量地刻画单点观测结果的代表性。利用印痕函数方法对中国风能资源观测网中来自不同下垫面的44个测风塔2010年的观测结果进行了测风塔资料的代表性分析。首先,根据两层(10、70 m)的风速和温度计算稳定度参数L_p并确定中性层结,粗糙度的计算选取了中性层结条件下满足对数风廓线的资料。然后,利用中性层结下70 m高度主导风向和风向频率较大方位的资料分析了两种印痕函数方法(Schuepp和Hsieh方法)的适用性及代表性范围与粗糙度的关系。两种方法得到的印痕函数对应距离(即代表性范围)比较接近。但是,在复杂山地Schuepp方法已不适用。当粗糙度小于10 m时,粗糙度与印痕函数对应距离的关系可用对数-线性函数描述。根据经验推荐了利用印痕函数积分90%对应的距离以及下垫面显著变化位置来确定测风塔代表范围。但是,印痕函数积分阈值仍需要通过示踪实验或数值模拟检验。     

Estimation of sector roughness lengths and the effect on prediction of the vertical wind speed profile

An estimate of roughness length is required by some atmospheric models and is also used in the logarithmic profile to determine the increase of wind speed with height under neutral conditions. The choice of technique for determining roughness lengths is generally constrained by the available input data. Here, we compare sets of roughness lengths derived by different methods for the same site and evaluate their impact on the prediction of the vertical wind speed profile.Wind speed and direction data have been collected at four heights over a three-year period at the North Norfolk Wind Monitoring Site. Wind speed profiles were used to generate sector roughness lengths based on the logarithmic profile formula. This is the only direct way of determining roughness lengths. The simplest and cheapest method is to use maps with published tables giving roughness length estimates for different terrain types. Alternatively Wieringa (1976, 1986) and Beljaars (1987) give formulae for determining roughness lengths from wind speed gusts or standard deviations.The four sets of estimated roughness lengths vary considerably. They were used to estimate 34 m wind speeds from 12.7 m observations. The profile-derived roughnesses are used simply as a check on the prediction of the wind speed profiles. The terrain-derived roughness lengths give reasonable results. Gust-derived and standard deviation roughnesses both predict wind speeds which are lower than the observed ones. The error is greater in the case of standard deviation roughnesses. If stability corrections are applied in the prediction of the vertical wind speed profile, the results are considerably improved.     

Airborne measurements of turbulent trace gas fluxes and analysis of eddy structure in the convective boundary layer over complex terrain

Using the new high-frequency measurement equipment of the research aircraft DO 128, which is described in detail, turbulent vertical fluxes of ozone and nitric oxide have been calculated from data sampled during the ESCOMPTE program in the south of France. Based on airborne turbulence measurements, radiosonde data and surface energy balance measurements, the convective boundary layer (CBL) is examined under two different aspects. The analysis covers boundary-layer convection with respect to (i) the control of CBL depth by surface heating and synoptic scale influences, and (ii) the structure of convective plumes and their vertical transport of ozone and nitric oxides. The orographic structure of the terrain causes significant differences between planetary boundary layer (PBL) heights, which are found to exceed those of terrain height variations on average. A comparison of boundary-layer flux profiles as well as mean quantities over flat and complex terrain and also under different pollution situations and weather conditions shows relationships between vertical gradients and corresponding turbulent fluxes. Generally, NO_x transports are directed upward independent of the terrain, since primary emission sources are located near the ground. For ozone, negative fluxes are common in the lower CBL in accordance with the deposition of O₃ at the surface.The detailed structure of thermals, which largely carry out vertical transports in the boundary layer, are examined with a conditional sampling technique. Updrafts mostly contain warm, moist and NO_x loaded air, while the ozone transport by thermals alternates with the background ozone gradient. Evidence for handover processes of trace gases to the free atmosphere can be found in the case of existing gradients across the boundary-layer top. An analysis of the size of eddies suggests the possibility of some influence of the heterogeneous terrain in mountainous area on the length scales of eddies.     

Surface drag and turbulence over an inhomogeneous land surface

Data collected over an inhomogeneous semi-rural area are presented. The data are compared with previous surface-layer data to determine how representative the fixed-point flux measurements are of area averages. Departures from the standard surface-layer results are found to be relatively small (~10–20%), which supports the concept of a blending height above which the flow ceases to respond to variations in the underlying surface and becomes horizontally homogeneous.Effective roughness lengths are derived for different wind directions and the relationship between the effective roughness length and upwind surface is examined in the light of recent ideas on averaging surface roughness lengths. It is found that by averaging drag coefficients, realistic values of the effective roughness length can be calculated which are not very sensitive to the precise choice of the component roughness lengths.     

基于一维变分算法的地基微波辐射计遥感大气温湿廓线研究

为研究地基微波辐射计遥感温、湿度廓线的一维变分算法的反演能力,用北京地区2010—2011年00和12时(世界时)的多通道地基微波辐射计亮温资料进行试验。首先,利用同时次的地面观测资料、红外亮温(由地基微波辐射计自带红外传感器测得)及探空观测数据,给出提取无云样本的方案,得到432个无云样本;再以辐射传输模式计算得到的模拟亮温为参考,对无云条件下的观测亮温进行质量控制;然后利用探空数据进行模拟试验,结果发现,一维变分算法对3 km以下的温度廓线有较大调整。使反演结果更加接近探空,而对湿度廓线在0—10 km都有不同程度的优化;最后利用一维变分算法对地基微波辐射计观测亮温进行大气温湿廓线反演,将结果与探空对比可以看出,温度廓线的均方根误差小于2.9 K,绝对湿度的均方根误差小于0.47 g/m~3;进一步与地基微波辐射计自带神经网络的反演结果比较表明,一维变分的反演结果更接近实际大气。