台风“珍珠”登陆期间动量通量的多尺度分析

利用台风"珍珠"登陆前后的近地层湍流观测资料,分析了该台风经过观测场地前后地面气象要素的变化及其动量通量特征。结果表明,台风"珍珠"经过观测场地前后的近地层气象要素发生了急剧的变化,并且在台风前部存在强的中尺度对流系统,反映在风速能谱密度结构上,频率f在3×10-4~2×10-3Hz之间的中尺度信号对能谱的贡献比平稳天气形势下的能谱贡献大很多,尤其是顺风方向风速的能谱密度的峰值与湍流信号的峰值相当;动量通量分析结果表明,台风中心经过观测场地的前后三小时,近地层通量以向下输送的中尺度通量为主,湍流通量的贡献相对于中尺度通量较小,也是向下输送的;而在其他时段,近地层通量主要以向上输送的湍流通量为主,中尺度通量量值很小,可以忽略。     

Variation of turbulence in the marine boundary layer over the Arabian Sea during indian southwest monsoon (MONEX 79)

Analysis of high frequency (20 Hz) turbulence data collected from low level flights by the National Center for Atmospheric Research (NCAR) Electra aircraft over the central Arabian Sea on 20 and 24 June, 1979 as part of MONEX 79 indicates the influence of the Somali Jet on boundary-layer turbulence. Different stages of monsoon development were evident on the two observation days from mean boundary-layer profiles. However, turbulence statistics of wind speed components and temperature in the monsoon boundary layer for both days are generally greater than those observed in laboratory experiments or tropical and trade wind boundary layers in which a strong jet was not present. Analysis of high frequency wind, temperature and humidity data was made to obtain fluxes of momentum and heat. Magnitudes of the sensible and latent heat fluxes are three to five times larger than the values observed over the monsoon boundary layer over the Bay of Bengal. The turbulent kinetic energy budget over the Arabian Sea for 24 June indicates the importance of buoyancy, and to a lesser extent shear as the dominant term. Dissipation serves as the primary sink term.     

南京雾-霾天气个例湍流运动特征的对比研究

利用超声风速仪、能见度仪等探测资料,采用Fortran、Matlab、Origin等软件处理数据并绘图分析,对比了南京郊区2013年12月4日一次霾过程、2013年12月7—9日一次雾过程和2013年12月3日一个晴天的湍流运动特征。结果表明,霾天和晴天平均水平风速、平均动能、湍流动能、湍流强度、摩擦速度、动量通量和热量通量都有明显的日变化,而这些参量在雾天没有明显的日变化,不同天气湍流参量大小有差异;霾、雾、晴天近中性层结下,u、v、w三个方向风速归一化标准差近似为常数,霾天分别为3.15、2.72、1.17,雾天分别为3.11、2.45、1.25,晴天分别为3.40、3.45、1.50;不稳定条件下霾、雾、晴天风速归一化标准差和湍流动能归一化标准差符合1/3幂次律,稳定条件下霾、雾、晴天无因次湍流动能均满足1/3幂次律;不稳定条件下霾、雾、晴天温度和湿度归一化标准差满足-1/3幂次律,稳定条件下霾、雾、晴天温度归一化标准差符合-2/3幂次律;雾天归一化u、v、w谱与Kansas谱吻合度比霾天和晴天高,霾天u、v、w谱峰值频率约为0.01 Hz,雾天u谱峰值频率约为0.004 Hz,v、w谱峰值频率约为0.1 Hz,晴天u、v、w谱峰值频率约为0.01 Hz。得出结论:雾天的边界层结构与霾天和晴天有所不同;不稳定条件下霾、雾、晴天风速、湍流动能、温度、湿度归一化标准差均符合莫宁-奥布霍夫相似理论;雾天归一化u、v、w谱比霾天和晴天更加符合莫宁-奥布霍夫局地各向同性理论,且霾天和晴天以机械湍流为主,雾天既有机械湍流也有热力湍流。     

Multifractal Characteristics of Intermittent Turbulence in the Urban Canopy Layer

The multifractality of energy and thermal dissipation of fully developed intermittent turbulence is investigated in the urban canopy layer under unstable conditions by the singularity spectrum for the fractal dimensions of sets of singularities characterizing multifractals. In order to obtain high-order moment properties of smallscale turbulent dissipation in the inertial range, an ultrasonic anemometer with a high sampling frequency of 100 Hz was used. The authors found that the turbulent signal could be singular everywhere. Moreover, the singular exponents of energy and thermal dissipation rates are most frequently encountered at around 0.2, which is significantly smaller than the singular exponents for a wind tunnel at a moderate Reynolds number. The evidence indicates a higher intermittency of turbulence in the urban canopy layer at a high Reynolds number, which is demonstrated by the data with high temporal resolution. Furthermore, the temperature field is more intermittent than the velocity field. In addition, a large amount of samples could be used for verification of the results.     

Vertical Structure Of Turbulence In Offshore Flow During Rasex

The adjustment of the boundary layer immediately downstream froma coastline is examined based on two levels of eddy correlation data collected on a mast at the shore and six levels of eddy correlation data and profiles of mean variables collected from a mast 2 km offshore during the Risø Air-Sea Experiment. The characteristics of offshore flow are studied in terms of case studies and inter-variable relationships for the entire one-month data set. A turbulent kinetic energy budget is constructed for each case study.The buoyancy generation of turbulence is small compared to shear generation and dissipation. However, weakly stable and weakly unstable cases exhibit completely different vertical structure. With flow of warm air from land over cooler water, modest buoyancy destruction of turbulence and reduced shear generation of turbulence over the less rough sea surface cause the turbulence to rapidly weaken downstream from the coast. The reduction of downward mixing of momentum by the stratification leads to smaller roughness lengths compared to the unstable case. Shear generation at higher levels and advection of stronger turbulence from land often lead to an increase of stress and turbulence energy with height and downward transport of turbulence energy toward the surface.With flow of cool air over a warmer sea surface, a convective internal boundary layer develops downstream from the coast. An overlying relatively thick layer of downward buoyancy flux (virtual temperature flux) is sometimes maintained by shear generation in the accelerating offshore flow.     

Scale Properties of Anisotropic and Isotropic Turbulence in the Urban Surface Layer

The scale properties of anisotropic and isotropic turbulence in the urban surface layer are investigated. A dimensionless anisotropic tensor is introduced and the turbulent tensor anisotropic coefficient, defined as C, where \(C = 3d_{3}\,+\,1 (d_{3}\) is the minimum eigenvalue of the tensor) is used to characterize the turbulence anisotropy or isotropy. Turbulence is isotropic when \(C \approx 1\), and anisotropic when \(C \ll 1\). Three-dimensional velocity data collected using a sonic anemometer are analyzed to obtain the anisotropic characteristics of atmospheric turbulence in the urban surface layer, and the tensor anisotropic coefficient of turbulent eddies at different spatial scales calculated. The analysis shows that C is strongly dependent on atmospheric stability \(\xi = (z-z_{\mathrm{d}})/L_{{\textit{MO}}}\), where z is the measurement height, \(z_{\mathrm{d}}\) is the displacement height, and \(L_{{\textit{MO}}}\) is the Obukhov length. The turbulence at a specific scale in unstable conditions (i.e., \(\xi < 0\)) is closer to isotropic than that at the same scale under stable conditions. The maximum isotropic scale of turbulence is determined based on the characteristics of the power spectrum in three directions. Turbulence does not behave isotropically when the eddy scale is greater than the maximum isotropic scale, whereas it is horizontally isotropic at relatively large scales. The maximum isotropic scale of turbulence is compared to the outer scale of temperature, which is obtained by fitting the temperature fluctuation spectrum using the von Karman turbulent model. The results show that the outer scale of temperature is greater than the maximum isotropic scale of turbulence.     

Subgrid-Scale Dynamics of Water Vapour, Heat, and Momentum over a Lake

We examine the dynamics of turbulence subgrid (or sub-filter) scales over a lake surface and the implications for large-eddy simulations (LES) of the atmospheric boundary layer. The analysis is based on measurements obtained during the Lake-Atmosphere Turbulent EXchange (LATEX) field campaign (August–October, 2006) over Lake Geneva, Switzerland. Wind velocity, temperature and humidity profiles were measured at 20 Hz using a vertical array of four sonic anemometers and open-path gas analyzers. The results indicate that the observed subgrid-scale statistics are very similar to those observed over land surfaces, suggesting that the effect of the lake waves on surface-layer turbulence during LATEX is small. The measurements allowed, for the first time, the study of subgrid-scale turbulent transport of water vapour, which is found to be well correlated with the transport of heat, suggesting that the subgrid-scale modelling of the two scalars may be coupled to save computational resources during LES.     

HUBEX试验区近地面层的湍流输送

1998年淮河流域能量和水循环试验(HUBEX)期间曾进行了1个月的近地面层湍流观测.分析不稳定条件下湍流的统计特性和谱特征,并与Monin-Obukhov相似理论进行了比较.结果表明,不稳定的时候各湍流量的统计特征与相似理论的预期相符.虽然不稳定条件下温度和湿度涨落的相关系数很高,谱的式样也相近,但湿度谱的峰值频率高于温度谱.协谱曲线的形状显示感热通量的谱峰较宽,表现出w和T在较宽范围的强相关性,而水汽通量谱在高频段下降很快,说明水汽的输送更多地出现在低频部分.从谱相关系数可见,在近中性的时候,各尺度湍流涡的热量输送效率普遍较低,随着不稳定性增强而显著提高.分析还发现,不论不稳定性的程度如何,小尺度湍流的水汽输送效率都较低.水汽通量谱的相关系数随稳定度的变化不如热通量的谱相关系数大,表明近中性时除小尺度湍流外其他湍流涡的水汽输送效率高于热量输送.     

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.     

The Effect of Wind-Turbine Wakes on Summertime US Midwest Atmospheric Wind Profiles as Observed with Ground-Based Doppler Lidar

We examine the influence of a modern multi-megawatt wind turbine on wind and turbulence profiles three rotor diameters ( $D$ D ) downwind of the turbine. Light detection and ranging (lidar) wind-profile observations were collected during summer 2011 in an operating wind farm in central Iowa at 20-m vertical intervals from 40 to 220 m above the surface. After a calibration period during which two lidars were operated next to each other, one lidar was located approximately $2D$ 2 D directly south of a wind turbine; the other lidar was moved approximately $3D$ 3 D north of the same wind turbine. Data from the two lidars during southerly flow conditions enabled the simultaneous capture of inflow and wake conditions. The inflow wind and turbulence profiles exhibit strong variability with atmospheric stability: daytime profiles are well-mixed with little shear and strong turbulence, while nighttime profiles exhibit minimal turbulence and considerable shear across the rotor disk region and above. Consistent with the observations available from other studies and with wind-tunnel and large-eddy simulation studies, measurable reductions in wake wind-speeds occur at heights spanning the wind turbine rotor (43–117 m), and turbulent quantities increase in the wake. In generalizing these results as a function of inflow wind speed, we find the wind-speed deficit in the wake is largest at hub height or just above, and the maximum deficit occurs when wind speeds are below the rated speed for the turbine. Similarly, the maximum enhancement of turbulence kinetic energy and turbulence intensity occurs at hub height, although observations at the top of the rotor disk do not allow assessment of turbulence in that region. The wind shear below turbine hub height (quantified here with the power-law coefficient) is found to be a useful parameter to identify whether a downwind lidar observes turbine wake or free-flow conditions. These field observations provide data for validating turbine-wake models and wind-tunnel observations, and for guiding assessments of the impacts of wakes on surface turbulent fluxes or surface temperatures downwind of turbines.     

北京地区一次强沙尘暴过程的大气边界层结构和湍流通量输送特征

利用中国科学院大气物理研究所北京325 m气象塔的风速和温度平均场观测资料和湍流资料,以及北京市气象台地面常规气象资料和逐日08:00和20:00(北京时间)的探空资料,分析了2002年3月18～22日沙尘暴过境前后北京城市边界层结构特征和湍流输送特征,结果表明:1)在沙尘暴爆发前,边界层中水平风速一直较小;气温较高,大气层结稳定,在边界层上部有强大的逆温层.随着冷锋过境,沙尘暴爆发,边界层中水平风速和平均湍流速度急剧增强;温度也突然变化,先迅速增强后又持续下降,逆温层迅速被破坏.2)沙尘暴初期,280 m上为系统性上升气流,而47和120 m则为系统性的下沉气流.随着沙尘暴爆发,湍流动能、向下传输的动量以及向上传输的感热也迅速增大,并且120 m高度的湍能、动量通量以及感热通量明显高于47和280 m,这与北京的局地环流有关.3)本次沙尘暴过程中,120和47 m层的摩擦速度都明显超过了北京的临界摩擦速度,表明局地起沙也是本次沙尘暴过程中北京沙尘的一个重要沙源.     

A study of intermittent turbulence with cases-99 tower measurments

Characteristics of intermittent turbulence events in the stably stratified nocturnal boundary layer are investigated with data collected in the CASES-99 tower array of 300-m radius. The array consists of a central 60-m tower with eddy covariance measurements at eight levels and six satellite towers with eddy covariance measurements at 5 m. A significant increase in the magnitude of vertical wind velocity () and spectral information are used to define the onset of an intermittent turbulence event. Normally, only a subset of 5 m-levels in the tower network experience an intermittent turbulence event concurrent with one at the 5 m-level on the main tower. This behaviour reveals the small horizontal extent of most events. Intermittent turbulence events at the main tower 5-m level are normally confined to a layer much thinner than the 60-m tower height. The turbulent kinetic energy budget is evaluated for intermittent turbulence events observed at the 5-m level on the main tower. Generally, the onset of an intermittent turbulence event is not closely related to the reduction of the gradient Richardson number below 0.25, the critical Richardson number of turbulence generation for linear instability. Possible explanations including the influence of advected turbulence patches are discussed.     

应用E-ε湍流动能闭合湖泊热力学过程模型对东太湖湖-气交换的模拟

采用考虑沉水植物影响的E-ε湍流动能闭合湖泊热力学过程模型,模拟2013年8月东太湖湖-气交换过程,并利用太湖的站点观测数据对模型进行了验证。太湖水温的模拟值与观测值吻合较好,模型计算的各层水温与观测值相比,均方根误差均未超过1℃。同时模型也较好地模拟出太湖表面感热通量和潜热通量,潜热通量的模拟值与观测值的标准差为54.7 W/m2。由于湖水较浅,太湖的水温层结会明显受到天气状况的影响。晴朗小风条件下的湖水呈现显著的热分层现象,当风速为0.8 m/s,高层和底层的温差达到7.9℃。大风天气条件驱动较强的水体湍流混合,水温的热分层消失,风速为12 m/s,湖泊上层与底层的水温差仅0.12℃。此外,模拟结果较好地呈现出了东太湖沉水植物的存在通过增大湖体消光系数,减小到达湖体内部的热量,并增加对湖水的阻力,影响湖体中湍流动能的分布,并进而影响湖水温度的分布。综上所述,该模型能够较好地模拟出浅水大湖湖-气交换的过程。      

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.     

北京市郊区及城区边缘的大气湍流结构特征

为了比较城市边缘地区和郊区大气边界层中湍流结构的差异,本文分析了1986年春、夏两季在北京市城市边缘地区和郊区两处大气湍流观测结果.结果表明,由于城市边缘地区的下垫表面具有较高数值的粗糙度,使其近地面层大气受到较强的动力影响,在其湍流速度分量谱曲线;无量纲化湍流速度分量标准偏差,σ_u/u*,σ_v/u*和σ_w/u*;以及无量纲化湍流耗散率φ_ε等湍流特征量在近中性层结条件下显示出与大气稳定度的关系比较微弱.     

The Near-Calm Stable Boundary Layer

For the near-calm stable boundary layer, nominally 2-m mean wind speed <0.5 ms⁻¹, the time-average turbulent flux is dominated by infrequent mixing events. These events are related to accelerations associated with wave-like motions and other more complex small-scale motions. In this regime, the relationship between the fluxes and the weak mean flow breaks down. Such near-calm conditions are common at some sites. For very weak winds and strong stratification, the characteristics of the fluctuating quantities change slowly with increasing scale and the separation between the turbulence and non-turbulent motions can become ambiguous. Therefore, a new analysis strategy is developed based on the scale dependence of selected flow characteristics, such as the ratio of the fluctuating potential energy to the kinetic energy. In contrast to more developed turbulence, correlations between fluctuating quantities are small, and a significant heat flux is sometimes carried by very weak vertical motions with large temperature fluctuations. The relation of the flux events to small-scale increases of wind speed is examined. Large remaining uncertainties are noted.     

Impressions of the turbulence variability in a weakly stratified,flat-bottom deep-sea ‘boundary layer’

The character of turbulent overturns in a weakly stratified deep-sea is investigated in some detail using 144 high-resolution temperature sensors at 0.7 m intervals, starting 5 m above the bottom. A 9-day, 1 Hz sampled record from the 912 m depth flat-bottom (<0.5% bottom-slope) mooring site in the central-north Alboran Sea (W-Mediterranean) demonstrates an overall conservative temperature range of only 0.05 °C, a typical mean buoyancy period as large as 3 h and a 1 Hz-profile-vertically-averaged turbulence dissipation rate maximum of only 10⁻⁸ m² s⁻³. Nonetheless, this ‘boundary layer’ varies in height between <6 and >104 m above the bottom and is thus not homogeneous throughout; the temperature variations are seldom quiescent and are generally turbulent in appearance, well exceeding noise levels. The turbulence character is associated with small-scale internal waves; examples are found of both shear- and convection-driven turbulence; particular association, although not phase-locked, is found between turbulence variations and tidal rather than with inertial motions; the mean buoyancy frequency of a few times the inertial frequency implies the importance of ‘slantwise convection’ in the direction of the earth rotational vector rather than in the direction of gravity. Such convection is observed both in near-homogeneous and weakly stratified form.     

The Askervein hill project: A finite control volume prediction of three-dimensional flows over the hill

A three-dimensional nonlinear numerical model, that has been extensively used previously to predict environmental water flows, was applied to predict the flow over an isolated hill, Askervein. Predictions are reported for winds approaching the hill from 210 ° and 180 ° clockwise from north, both under almost neutral atmospheric conditions.The model predictions were compared with data collected during a major field study in 1983. From the comparisons it was concluded that the model predicts the mean flow variables with good accuracy. Larger discrepancies were found for quantities related to the turbulence, pointing to deficiencies in the turbulence model, and perhaps in some of the measurements.     

The partitioning of attached and detached eddy motion in the atmospheric surface layer using Lorentz wavelet filtering

Townsend's attached eddy hypothesis states that the turbulent structure in the constant stress layer can be decomposed into attached and detached eddy motion. This paper proposes and tests a methodology for separating the attached and detached eddy motion from time series measurements of velocity and temperature. The proposed methodology is based on the time-frequency localization and filtering capabilities of the orthonormal wavelet transforms. Using a relative entropy statistical measure, the optimal wavelet basis is identified first. The turbulence time series measurements are then transformed into the wavelet domain where the contribution of specific events in the time-frequency domain is identified. The filtering scheme utilizes a recently constructed Lorentz thresholding methodology that successfully eliminates all wavelet coefficients associated with the detached eddy motion. While this filtering scheme lacks the compression efficiency of the classical Donoho and Johnstone's universal thresholding model, it conserves the higher-order statistics and important turbulence interactions related to the Reynolds stresses. Following the filtering scheme, the attached eddy motion time series is re-constructed by an inverse wavelet transform of the non-zero wavelet coefficients. The proposed partitioning methodology for attached and detached eddy motion is tested using 56 Hz triaxial sonic anemometer velocity and temperature measurements above a uniform dry lake bed in Owens valley, California, for a wide range of atmospheric stability conditions. Validation that the wavelet filtered time series represents the attached eddy motion is also discussed in the context of conservation of turbulence energy and surface fluxes.     

大气边界层阵风相干结构的产生条件

壁湍流相干结构的发现是近代湍流研究的重大进展之一,从20世纪50年代开始,在大气边界层湍流中也发现了相干结构——对流云街,并进行了系统的研究。近些年来,人们发现在近地层湍流中也存在相干结构。利用北京325 m气象塔对城市下垫面中大风和小风天气的风速分析,发现较有规律的周期3~6 min的阵风,且有明显的相干结构,而对不同下垫面的阵风研究,均发现存在这种相干结构,这种阵风相干结构对通量输送有不可忽视的作用。本文利用2012年4月甘肃省民勤县巴丹吉林沙漠观测塔的超声风速和平均场风速、温度观测资料,对阵风相干结构的产生条件进行了分析。采用傅立叶变换,将三维超声风速按频率分成基流(周期10分钟以上)、阵风扰动(周期1到10分钟)、湍流脉动(周期小于1分钟)三部分,结合平均场的资料分析发现:阵风相干结构出现在静力中性、不稳定甚至略微稳定的条件下,或者说机械作用主导的大气边界层,阵风区就会出现相干结构,热力作用对其有抑制和干扰的作用。从而,阵风的相干结构和壁面相干结构都出现在中性条件下,是机械湍流的现象,都主导着动量能量的输运。阵风区的相干结构并不等同于对流云街,他们出现在不同的大气稳定度条件下且尺度不同。