由地基遥感资料确定大气边界层特征

中国科学院大气物理研究所大气边界层物理和大气化学国家重点实验室拥有用于边界层垂直结构探测的相控阵声雷达（PA2）、低层大气廓线仪（LAP-3000）以及无线电声雷达（RASS）。本文研究如何利用这些设备获取的遥感观测资料,确定边界层湍流热通量的垂直分布、地面热通量以及边界层高度等问题。     

北京地区一次典型大雾天气过程的边界层特征分析

利用中国科学院大气物理研究所325 m铁塔15层风、温、湿梯度观测资料和3层超声资料,对2002年12月1~4日发生在北京地区的持续大雾天气过程进行近地面层大气边界层特征分析。结果表明,近地面大气边界层较大的相对湿度(70%)、较小的风速(3.0 m.s-1)和风速垂直切变(0.02 s-1)、稳定的层结结构以及较低的气温是北京持续大雾天气形成的主要原因。冷空气的侵入使得边界层相对湿度迅速减小,风速和风速垂直切变增大,破坏近地面大气边界层的结构,导致大雾的消散。分析还发现,大雾的维持与消散主要受风场等动力因素的影响,热力层结是大雾维持和消散的必要条件。冷空气的侵入自上而下影响平均风场,而对湍流风场的影响则是自下而上的。尺度分析结果表明,大雾期间,近地面边界层内中尺度动量通量和感热通量都大于湍流尺度的,中尺度动量通量与平均风速基本呈反相关;冷空气的侵入使得湍流通量显著加强,是导致大雾天气消失的主要原因。     

2009年“莫拉克”台风登陆过程阵风特征分析

利用上海台风研究所移动观测车获取的“莫拉克”台风登陆过程中超声风、温等观测资料对地面阵风特性进行了诊断分析.结果表明,在风速时间序列中叠加有周期为3-7 min的阵风扰动,显现出明显的相干结构,即沿顺风方向阵风风速峰期有下沉运动,谷期有上升运动;阵风扰动的各向异性特征明显,沿顺风方向的阵风扰动能量最大,其次是沿侧向和垂直方向的扰动能量;沿顺风方向的阵风垂直动量通量向下传播,而沿侧风方向阵风扰动动量垂直通量总体贡献接近于0.阵风扰动沿顺风方向的积分空间尺度和时间尺度最大,沿侧风方向和垂直方向其次,均明显大于湍流的积分空间和时间尺度.此外,阵风扰动的其他特征还包括:感热垂直通量极小;当平均风速较大时阵风风向变化幅度较小,而风速较小时阵风风向变化幅度则较大;动力学分析表明,阵风扰动主要表现出重力内波的一些特性.     

非均匀对流边界层的地转强迫流动和动量输送

用大涡模拟方法研究地转强迫下的对流边界层流动和地表热力非均匀性影响.模拟重现了典型对流边界层的平均风廓线和动量通量垂直分布.地表热力非均匀性对区域平均风速和动量通量分布无明显影响,但边界层内的局地流动性状和湍流动量输送情况有系统性的改变.下风较热区近地面风速增强而高空流动受到阻塞,上风较冷区之上情况则正好相反.对应于平均流动场的畸变,地表较热区之上边界层大部可以出现动量向上输送的情况,较冷区成为大气动量下传的主要通道.地面应力在较热区增强、较冷区减弱的趋势明显.     

鄱阳湖地区复杂地表条件下一次强降水过程的近地面边界层特征

使用鄱阳湖北部70 m气象塔湍流和梯度观测数据,分析了2011年6月6日夜间一次暴雨过程中近地面边界层特征.结果表明,此次过程是在高空低槽和西南急流的天气背景下,受鄱阳湖复杂地表影响产生的局地性强降水.强降水发生前受东南暖平流影响,近地面边界层中水汽累积,不稳定性增加;强降水过程中,近地层感热、潜热通量迅速增加,同时,近地面层湍流动量通量下传和水平输送增加,鄱阳湖的水汽输送加强降水强度.另外,强降水过程中,近地面湍流动能迅速增大并达到最大值,而平均动能的增大发生在强降水结束后,表明地表作用明显,近地面边界层的湍流场为暴雨提供动力条件.尺度分析表明,强降水前,中尺度动量通量占主要地位,降水过程中湍流动量通量显著加强.     

夏季北京城市大气边界层湍流通量的塔层观测

利用2003年8月13～26日期间中国科学院大气物理研究所大气边界层气象塔47、120和280 m 3层高度上三维风速和温度脉动的观测资料,结合涡动相关方法,分析了夏季北京市城市下垫面的动量和感热等湍流通量。结果表明：3层高度上动量通量和感热通量的日变化彼此非常接近, 且3层高度上最大感热通量达到200～250 W·m-2。观测期间盛行南风,120 m高度的水平风速小于10 m·s-1。同时分析了由慢响应传感器采集的15个高度上的气温和8个高度上的相对湿度。结果表明：受地表加热影响,320 m高度上气温和相对湿度的日变化非常显著。不同高度上动量通量和感热通量的比较表明：城市下垫面的空气动力不均匀性比热力不均匀性显著。用廓线法计算了32和65 m高度之间的湍流通量,结果表明：1）夏季晴好天气情况下,北京城市下垫面能量输送以感热为主;2）由廓线法计算得到的动量通量、感热通量和由涡动相关法直接测量结果相近。此外还初步讨论了大气边界层温度、湿度垂直廓线随时间的演变。     

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

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

日间对流边界层中的非局地动量混合

日间对流边界层最显著的结构特征是在热力作用下所形成的组织化对流。与小尺度湍涡不同的是，组织化对流具有边界层尺度的垂直相干性，可实现垂直贯穿边界层的非局地物质和能量传输。本文针对对流边界层中的动量混合，探究组织化对流对动量输送的贡献。以高精度大涡模拟数据为研究资料，通过傅里叶变换、本征正交分解和经验模态分解3种滤波方法，分离组织化对流和背景湍涡，计算与两者相关的非局地和局地动量通量，发现与组织化对流相关的非局地动量通量是总通量的重要组成部分，并主导混合层中的垂直动量输送。而后，基于协谱和相位谱分析，探究组织化对流的空间结构对动量传输的影响，发现在热力主导的不稳定环境中，单体型环流结构对动量的传输效率较低。而在风切较强的近中性环境中，滚涡型组织化结构可使垂直和水平流向扰动速度的相位差减小，从而提升动量传输效率。研究结果表明，边界层方案需要包含非局地动量通量项，其参数化应考虑整体稳定度对传输效率的影响。     

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

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

北京地区一次PM2.5重污染过程的边界层特征分析

利用北京市环境保护监测中心和美国大使馆的细颗粒物（PM2.5）逐时监测数据,中国科学院大气物理研究所325 m气象梯度塔资料以及实况天气图和探空资料,对2015年11月27日至12月1日北京的PM2.5重污染过程的边界层特征进行了分析。研究发现:这次重污染过程持续时间长、强度大,其中PM2.5浓度超过75 μg/m3的时次共计126 h,超过150 μg/m3共计116 h,小时最高PM2.5浓度为522 μg/m3。在高低空环流场配置的影响下,近地面静风和多层逆温结构抑制了污染物在水平和垂直方向上的输送,加上边界层内的深厚湿层,使得其中气溶胶不断吸湿增长,高PM2.5浓度得以维持。在重污染期间,湍流动能较低,不利于污染物的水平和垂直扩散。垂直方向的湍流动能一直占水平方向的15%~20%左右,水平湍流动能占主要贡献。摩擦速度与湍流动能呈现出相似的变化趋势,不同高度之间的摩擦速度差别不大。超出前后时次一个数量级的湍流强度尖峰的出现是湍流场发生调整的一个信号,是PM2.5浓度发生剧烈转变的前兆,预示着污染状况更加糟糕。重污染过程中感热通量的输送方向为从地面向大气输送,感热通量和潜热通量都大幅减少,并且表现出明显的日变化特征。对湍流功率谱计算和分析表明,在重污染过程期间,时间尺度为5 min至6 h的中尺度过程对从地面到大气方向的动量和热量通量输送做出了重要贡献。     

Observations of Coherent Turbulence Structures in the Near-Neutral Atmospheric Boundary Layer

Turbulence structures of high Reynolds number flow in the near-neutral atmospheric boundary layer (ABL) are investigated based on observations at Shionomisaki and Shigaraki, Japan. A Doppler sodar measured the vertical profiles of winds in the ABL. Using the integral wavelet transform for the time series of surface wind data, the pattern of a descending high-speed structure with large vertical extent (from the surface to more than 200-m level) is depicted from the Doppler sodar data. Essentially this structure is a specific type of coherent structure that has been previously shown in experiments on turbulent boundary-layer flows. Large-scale high-speed structures in the ABL are extracted using a long time scale (240 s) for the wavelet transform. The non-dimensional interval of time between structures is evaluated as 3.0–6.2 in most cases. These structures make a large contribution to downward momentum transfer in the surface layer. Quadrant analyses of the turbulent motion measured by the sonic anemometer (20-m height) suggest that the sweep motion (high-speed downward motion) plays a substantial role in the downward momentum transfer. In general, the contribution of sweep motions to the momentum flux is nearly equal to that of ejection motions (low-speed upward motions). This contribution of sweep motions is related to the large-scale high-speed structures.     

Coherent Structures Detected in Atmospheric Boundary-layer Turbulence using Wavelet transforms at Huaihe River Basin, China

The presence of coherent structures in turbulent shear flows suggests order in apparently random flows. These coherent structures play an important dynamical role in momentum and scalar transport. To develop dynamical models describing the evolution of such motion, it is necessary to detect and isolate the coherent structures from the background fluctuations. In this paper, we decomposed atmospheric turbulence time series into large-scale eddies, which include coherent structures and small eddies, which are stochastic by using Fourier digital filtering. The wavelet energy computed for the three components of the velocity fluctuations in the large-scale eddies appears to have local maximum values at certain time scales, which correspond to the scales or frequencies of coherent structures. We extract coherent signals from large-scale vortices at this scale by inverse wavelet transform formulae. This method provides an objective technique for examining the turbulence signal associated with coherent structures in the atmospheric boundary layer. The average duration of coherent structures in three directions based on Mexican hat wavelets are 33 s, 34 s and 25 s respectively. Symmetric andanti-symmetric wavelet basis functions give almost the same results. The main features of the structures during the day and night have little difference. The dimensionless durations for u, v and w have linear correlations with each other. These relationships are insensitive to the wavelet basis.     

1998年南海季风爆发期间近海面层大气湍流结构和通量输送的观测研究

1998年5月14日至6月22日,在西沙永兴岛近海铁塔上进行了一次海-气通量观测试验,观测期包括了西南季风爆发前、爆发、爆发后风速加强等几个阶段。这次试验获得该年西南季风爆发的天气特征以及由涡度相关法、廓线法计算的动量、感热、潜热通量及湍流强度等一些统计量分布。分析结果表明,观测期湍流强度σu、σv、σw与平均风速之比为0.096、0.066、0.045;在近中性条件下（z/L≈0）,各相似函数基本为常数,σu/u＊≈3,σv/u＊≈2,σw/u＊≈1.25;在稳定条件下u、v、w三方向σ/u＊近似相等;在稳定与不稳定条件下σt/t＊随稳定度参数z/L的变化趋势相似,但符号相反。摩擦速度u＊随风速的变化接近于关系式u＊=0.029U10+0.006（U10为10m高度30min平均风速）,空气动力粗糙度长度z0变化在0.01～0.35mm之间。观测期动量通量变化在0.05～0.30N/m2之间,季风爆发期明显增大;曳力系数CD平均为1.12×10-3,它随平均风速变化可以表示为关系式103CD=0.003U210+0.020U10+0.836或103CD=0.056U10+0.732。感热通量由海洋输向大气,平均值为7.8W/m2,数值在0～15W/m2之间变化,季风爆发后明显增大;潜热通量数值一般变化于50～200W/m2,白天中午数值较高,夜间较低。季风爆发后也明显增大。观测期间鲍恩比（BowenRatio）日平均为0.05左右。     

Flux decomposition into coherent structures

This study examines the intermittency of the momentum flux near the surface and the relation of such intermittency to coherent structures. Toward this goal, variances and covariances are decomposed into coherent structures and less coherent activity. The sampled structures are identified using the Haar transform and then decomposed into eigenvectors of the lagged covariance matrix.The methodology is applied to the momentum flux for a relatively stationary 50-h period of strong winds measured from a 45 m tower in the Lammefjord Experiment. Events of sinking motion with strong horizontal momentum account for the majority of the flux. Such sweeping motions arrive as gust microfronts. The large momentum flux is associated with strong coherent fluctuations of the longitudinal wind component and high correlation with relatively modest fluctuations of vertical motion. In the heated case (HAPEX), a phase lag between the vertical and horizontal velocity fluctuations leads to less efficient momentum transport by the main coherent structures.The event nature of the flux is used to formulate an expression for the flux error due to sampling problems. Estimation of the momentum flux requires a significantly longer record than for the heat flux. Modulation of the flux by mesoscale variations also affects the sampling strategy.     

用连续子波变换提取城市冠层大气湍流的相干结构

切变湍流的相干结构是湍流研究中的重大发现,它表明湍流在表面上看来不规则运动中具有可检测的有序运动,这种相干结构在切变湍流的脉动生成和发展中起着主宰作用.因此识别和提取相干结构对于认识和研究湍流是非常重要的.用数字滤波法将包含相干结构的大尺度信号提取出来以后,再用子波分析,根据子波能量极大值的判别方法,分别确定出大气湍流三个方向上的速度脉动信号相干结构的频率或时间尺度,然后由确定尺度上的连续子波反演公式,提取出大气湍流三个方向上的速度脉动信号相干结构所对应的波形.     

Large-Eddy Simulation of Coherent Turbulence Structures Associated with Scalar Ramps Over Plant Canopies

Large-eddy simulations were performed of a neutrally-stratified turbulent flow within and above an ideal, horizontally- and vertically-homogeneous plant canopy. Three simulations were performed for shear-driven flows in small and large computational domains, and a pressure-driven flow in a small domain, to enable the nature of canopy turbulence unaffected by external conditions to be captured. The simulations reproduced quite realistic canopy turbulence characteristics, including typical ramp structures appearing in time traces of the scalar concentration near the canopy top. Then, the spatial structure of the organised turbulence that caused the scalar ramps was examined using conditional sampling of three-dimensional instantaneous fields, triggered by the occurrence of ramp structures. A wavelet transform was used for the detection of ramp structures in the time traces. The ensemble-averaged results illustrate that the scalar ramps are associated with the microfrontal structure in the scalar, the ejection-sweep structure in the streamwise and vertical velocities, a laterally divergent flow just around the ramp-detection point, and a positive, vertically-coherent pressure perturbation. These vertical structures were consistent with previous measurements made in fields or wind tunnels. However, the most striking feature is that the horizontal slice of the same structure revealed a streamwise-elongated region of high-speed streamwise velocity impacting on another elongated region of low-speed velocity. These elongated structures resemble the so-called streak structures that are commonly observed in near-wall shear layers. Since elongated structures of essentially similar spatial scales were observed in all of the runs, these streak structures appear to be inherent in near-canopy turbulence. Presumably, strong wind shear formed just above the canopy is involved in their formation. By synthesis of the ensemble-averaged and instantaneous results, the following processes were inferred for the development of scalar microfronts and their associated flow structures: (1) a distinct scalar microfront develops where a coherent downdraft associated with a high-speed streak penetrates into the region of a low-speed streak; (2) a stagnation in flow between two streaks of different velocities builds up a vertically-coherent high-pressure region there; (3) the pressure gradients around the high-pressure region work to reduce the longitudinal variations in streamwise velocity and to enhance the laterally-divergent flow and lifted updrafts downstream of the microfront; (4) as the coherent mother downdraft impinges on the canopy, canopy-scale eddies are formed near the canopy top in a similar manner as observed in conventional mixing-layer turbulence.     

Analysis of aircraft measurements of momentum flux in the subcloud layer over the tropical atlantic ocean during gate

The momentum flux data obtained by the gust probe aboard the NOAA DC-6 aircraft during GATE are analyzed. Vertical profiles are obtained for Phases I and III and correlated with vertical wind velocity profiles using the geostrophic departure method. Reasonable agreement is obtained using the horizontal equations of motion with negligible advective acceleration. The vertical profiles of momentum flux and wind speed variance compare well with the numerical model results of Deardorff (1972) and Wyngaard et al. (1974). Vertical distributions of power spectra for vertical eddy motion and cospectra corresponding to the momentum flux components are obtained along with the height variation of the dominant length scales of vertical eddy motion and the dissipation rate of turbulence kinetic energy. When normalized by mixed-layer similarity, these results agree well with previous determinations in the boundary layer over tropical oceans and over land.