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

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

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

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

南京夏季城郊湍流统计特性及湍流通量对比分析

利用2010年南京夏季城市热岛三维结构试验中的湍流观测资料,对南京市委党校和南京信息工程大学观测场2个观测点的湍流统计量和湍流通量进行了对比分析。结果表明：城市平均风速大于郊区;城郊2地的湍流强度都呈现出Ix〉Iy〉Iz的趋势,水平方向上的湍强城郊相差不大,而垂直方向上,郊区的Iz要明显低于城市;城、郊的风速归一化标准差σu/u＊、σv/u＊、σw/u＊都符合MoninObukhov相似理论的＂1/3＂定律,城市的垂直风速归一化标准差大于郊区,而水平方向上的风速归一化标准差则明显小于其他下垫面;郊区风速归一化标准差与常熟农田的拟合结果较为相似,而城市风速归一化标准差与长白山森林的拟合结果更为接近;夏季城市以感热通量为主,而郊区湍流能量的输送以潜热通量为主,且城市的湍流热通量受太阳辐射的影响更大。     

鄂东两次暴雨前后近地层物理量场异常特征分析

利用湖北黄石长江南岸50 m高铁塔的风梯度、超声风及气温、湿度观测资料,对2007年5月31日及7月1-2日两场暴雨前后近地层风温湿场和降水前湍流特征进行了计算分析,探索其异常变化特征,为进一步认识黄石地区强降水的近地面层物理过程提供依据。结果表明:(1)降水前温度下降,湿度上升,平均风向发生明显转变,水平风速与垂直气流速度在降水前出现显著增大;降水后,风向平缓变化,平均风速再次增大。(2)垂直气流速度在降水前频繁上下振荡,激发湍流活动。(3)降水前湍流动能增大,湍流动量通量和感热通量增大向上传递,湍流强度逐步增大。湍流动能最大值的出现早于湍流通量,湍流活动在临近降水时进一步加强。可见暴雨发生前后近地层风温湿场有明显的变化,降水前湍流活动增强。     

城市冠层上下大气湍流特征分析

利用兰州市榆中县城市冠层架设的3台涡动相关仪(EC)观测大气湍流资料,分析了城市冠层上下不同下垫面湍流通量和不同风向范围内湍流动能变化特征,之后对城市冠层上下3台EC观测湍流动能最小、最大分布方向上分别对应的最小最大湍流动能的风速谱进行了研究,进一步检验了局地相似理论在城市冠层上下的适用性。结果表明:(1)城市冠层之上水泥、砖石等构成下垫面和城市冠层之内草坪下垫面观测感热通量、摩擦速度较为接近,观测潜热通量、CO2通量在白天差异明显。(2)城市冠层之上的湍流动能总体上大于冠层之内,冠层之上气流来向的上风向较为开阔时湍流动能较大,而冠层之内气流来向的上风向为街道口时湍流动能较大。(3)城市冠层之上的湍涡尺度大于冠层之内,城市冠层小尺度湍涡风速谱在惯性副区基本符合-2/3次方关系,且准各向同性,大尺度湍涡风速谱在惯性副区不符合-2/3次方关系,且各向异性。(4)不稳定层结下,城市冠层上下无量纲速度方差与稳定度基本满足1/3次方局地相似关系,稳定层结下不满足;城市冠层上下无量纲温度、湿度、CO2浓度方差在所有层结下均不满足-1/3次方局地相似关系。(5)近中性层结下,城市冠层上下u、v、w方向无量纲速度方差分别为3.52,3.03,1.49和2.62,2.22,1.50。     

2010年春季北京地区强沙尘暴过程的微气象学特征

利用北京大学校园地区PM10质量浓度观测资料、中国科学院大气物理研究所325m气象塔气象要素梯度和湍流观测资料,分析了北京地区2010年3月20～22日两次强沙尘暴过程微气象学要素和沙尘参量的时空演变以及湍流输送特征,为理解北京地区强沙尘暴天气沙尘输送规律和微气象学特征提供参考。结果表明:3月20～22日强沙尘暴过程前后不同高度温度先升后降,气压和相对湿度则相反。强沙尘暴来临时,高层风速先迅速增大,低层风速增加略有滞后,风切变明显加强,PM10浓度最大值和风速极大值出现时间较吻合。强沙尘暴过境时,不同高度向下的湍流动量输送、向上的湍流热量输送和湍流动能明显加强。与3月21日非沙尘暴日相比,强沙尘暴过程湍流动量通量增加,有利于沙尘粒子的水平和垂直输送过程;由于冷锋过境,水平热通量增大;垂直热通量因白天温度垂直梯度减小而减小,夜间因逆温层被破坏而增加;水平湍流动能对湍流动能占主要贡献,垂直湍流动能仅占水平湍流动能的10%～25%。     

不均匀地表产生的中尺度通量的数值试验

采用北京大学三维的复杂地形中尺度模式,结合陆面过程模式(SiB),模拟了草原和沙漠并存的下垫面的边界层大气运动.利用SiB模式计算了地表辐射、感热、潜热通量,并且预报地表温度.中尺度模式则模拟了沙漠地区受热抬升,形成的辐合运动,垂直速度的分布,不同高度上水平流场的变化以及中尺度动量和热量通量,把中尺度通量跟湍流通量进行了比较,以确定这种中尺度运动在GCM模式的参数化过程中的重要性.试验表明中尺度通量尤其是热量通量要比湍流通量大很多.     

北京地区一次强沙尘天气过程的中尺度通量特征

应用大气所325 m铁塔的湍流资料和梯度资料,对2003年3月20～23日影响北京的一次强沙尘天气过程的中尺度通量和湍流通量进行了分析,结果表明:冷空气的影响是从上层开始的.沙尘爆发前,边界层有强逆温存在,动量以中尺度输送为主,湍流很弱.沙尘爆发时,向下的湍流动量通量明显加大.湍流动量通量是沙尘沉降的主要原因,但沙尘爆发前的中尺度过程不能忽略.沙尘爆发前后,感热通量均以小尺度的湍流输送为主,中尺度输送不明显.     

一次暴雪过程前后近地层物理量场特征分析

利用铁塔风梯度观测资料和超声风温仪观测资料,对2008年1月18—21日暴雪前后,湖北黄石长江岸边近地层风场和湍流作了计算分析,探索其异常变化特征,为认识黄石地区暴雪近地层发生发展的物理过程提供依据。结果表明,暴雪前,风向转变,水平风速和垂直风速明显增大,湍流通量的输送较活跃,湍流动能和湍流强度有显著峰值出现;降雪过程结束后,湍流动能再次增大后缓慢减弱。可见此次暴雪过程前后近地层物理量场有异常变化     

云贵高原西部大理地区近地层湍流特征分析

利用大理国家气候观象台2007年3—12月观测资料,采用涡动相关法等计算方法分析了该地区湍流强度、湍流方差、湍流通量等特征量的日变化和干湿季变化特征。结果表明,湍流强度干季大于湿季;湍流方差与稳定度满足1／3次方定律,风速方差在稳定条件下比不稳定条件下离散,水平方向比垂直方向离散;湍流通量有明显日变化特征,感热、动量通量干季大于湿季,潜热通量湿季大于干季,干湿季热量交换以潜热为主;干季能量闭合率大于湿季;不同风向条件下平均水汽密度存在差异.     

50 Years of the Monin–Obukhov Similarity Theory

In weak wind stable conditions, eddy-correlation fluxes calculated using conventional averaging times of 5 min or longer to define the perturbations are severely contaminated by poorly sampled mesoscale motions. A method is developed to identify the averaging time for each individual data record that captures the turbulence while excluding most of the mesoscale motions. The method is based on multiresolution decomposition of the heat flux, and provides an objective procedure for selecting the averaging time for calculating eddy-correlation fluxes. Eddy-correlation data collected in weak turbulence conditions over grass, snow, a pine forest and the ocean are used to demonstrate the approach.When the small-scale turbulence and mesoscale motions are clearly separated by a gap region in the heat flux cospectra, the variable window width reduces the influence of nonstationarity by more effectively filtering out mesoscale motions compared to traditional methods using constant averaging time. For records where turbulence and mesoscale motions overlap in scale, the method is not well posed, although such records occur infrequently for our datasets. These ambiguous cases correspond to significant nonstationarity at scales that overlap with turbulence scales. The improved turbulence fluxes calculated with the proposed method are the appropriate fluxes for evaluating flux-gradient relationships and Monin–Obukov similarity theory for developing improved model parameterizations of turbulence for weakly turbulent flows     

A Hierarchy of Energy- and Flux-Budget (EFB) Turbulence Closure Models for Stably-Stratified Geophysical Flows

Here we advance the physical background of the energy- and flux-budget turbulence closures based on the budget equations for the turbulent kinetic and potential energies and turbulent fluxes of momentum and buoyancy, and a new relaxation equation for the turbulent dissipation time scale. The closure is designed for stratified geophysical flows from neutral to very stable and accounts for the Earth’s rotation. In accordance with modern experimental evidence, the closure implies the maintaining of turbulence by the velocity shear at any gradient Richardson number Ri, and distinguishes between the two principally different regimes: “strong turbulence” at ${Ri \ll 1}$ typical of boundary-layer flows and characterized by the practically constant turbulent Prandtl number Pr _T; and “weak turbulence” at Ri > 1 typical of the free atmosphere or deep ocean, where Pr _T asymptotically linearly increases with increasing Ri (which implies very strong suppression of the heat transfer compared to the momentum transfer). For use in different applications, the closure is formulated at different levels of complexity, from the local algebraic model relevant to the steady-state regime of turbulence to a hierarchy of non-local closures including simpler down-gradient models, presented in terms of the eddy viscosity and eddy conductivity, and a general non-gradient model based on prognostic equations for all the basic parameters of turbulence including turbulent fluxes.     

On the gap in the spectra of surface-layer atmospheric turbulence

A large set of tower data was used to identify the gap that separates small-scale turbulence and mesoscale structures in the cospectra of surface fluxes. The cospectra were obtained using a multi-resolution decomposition algorithm. The gap time scale τ _g was found by fitting a fifth-order polynomial to the cospectra and identifying special points occurring after the peak at small scales. In unstable conditions (day) τ _g was found to fall as the mean wind speed increased, while no such dependence was observed in stable conditions (night). The gap scale was found to change very weakly with stability both in moderately stable and moderately unstable conditions, with a sharp drop from about 1100 to 250 s occurring in near-neutral conditions. The vertical fluxes computed at different averaging intervals were found to correlate exceptionally well with each other, the scatter being somewhat larger during the night. Although considerable discrepancy may occur for individual records, when averaged over 10 months, the difference in the flux estimated at 7 to 109 min intervals never exceeded 4%, which is comparable or less than the instrumental error.     

Towards Closing the Surface Energy Budget of a Mid-latitude Grassland

Observations for May and August, 2005, from a long-term grassland meteorological station situated in central Netherlands were used to evaluate the closure of the surface energy budget. We compute all possible enthalpy changes, such as the grass cover heat storage, dew water heat storage, air mass heat storage and the photosynthesis energy flux, over an averaging time interval. In addition, the soil heat flux was estimated using a harmonic analysis technique to obtain a more accurate assessment of the surface soil heat flux. By doing so, a closure of 96% was obtained. The harmonic analysis technique appears to improve closure by 9%, the photosynthesis for 3% and the rest of the storage terms for a 3% improvement of the energy budget closure. For calm nights (friction velocity u _* < 0.1 m s⁻¹) when the eddy covariance technique is unreliable for measurement of the vertical turbulent fluxes, the inclusion of a scheme that calculates dew fluxes improves the energy budget closure significantly.     

Time Scales in the Unstable Atmospheric Surface Layer

Calculation of eddy covariances in the atmospheric surface layer (ASL) requires separating the instantaneous signal into mean and fluctuating components. Since the ASL is not statistically stationary, an inherent ambiguity exists in defining the mean quantities. The present study compares four methods of calculating physically relevant time scales in the unstable ASL that may be used to remove the unsteady mean components of instantaneous time signals, in order to yield local turbulent fluxes that appear to be statistically stationary. The four mean-removal time scales are: (t _c ) based on the location of the maximum in the ogive of the heat flux cospectra, () the location of the zero crossing in the multiresolution decomposition of the heat flux, (t ^*) the ratio of the mixed-layer depth over the convective velocity, and () the convergence time of the vertical velocity and temperature variances. The four time scales are evaluated using high quality, three-dimensional sonic anemometry data acquired at the Surface Layer Turbulence and Environmental Science Test (SLTEST) facility located on the salt flats of Utah’s western desert. Results indicate that and , with t _c achieving values about 2–3 times greater than t ^*. The sensitivity of the eddy covariances to the mean-removal time scale (given a fixed 4-h averaging period during midday) is also demonstrated.     

A Simple Model for the Afternoon and Early Evening Decay of Convective Turbulence Over Different Land Surfaces

A simple model to study the decay of turbulent kinetic energy (TKE) in the convective surface layer is presented. In this model, the TKE is dependent upon two terms, the turbulent dissipation rate and the surface buoyancy fluctuations. The time evolution of the surface sensible heat flux is modelled based on fitting functions of actual measurements from the LITFASS-2003 field campaign. These fitting functions carry an amplitude and a time scale. With this approach, the sensible heat flux can be estimated without having to solve the entire surface energy balance. The period of interest covers two characteristic transition sub-periods involved in the decay of convective boundary-layer turbulence. The first sub-period is the afternoon transition, when the sensible heat flux starts to decrease in response to the reduction in solar radiation. It is typically associated with a decay rate of TKE of approximately t ⁻² (t is time following the start of the decay) after several convective eddy turnover times. The early evening transition is the second sub-period, typically just before sunset when the surface sensible heat flux becomes negative. This sub-period is characterized by an abrupt decay in TKE associated with the rapid collapse of turbulence. Overall, the results presented show a significant improvement of the modelled TKE decay when compared to the often applied assumption of a sensible heat flux decreasing instantaneously or with a very short forcing time scale. In addition, for atmospheric modelling studies, it is suggested that the afternoon and early evening decay of sensible heat flux be modelled as a complementary error function.     

Turbulence Characteristics Of The Stable Boundary Layer Over A Mid-Latitude Glacier. Part Ii: Pure Katabatic Forcing Conditions

Observations obtained over a glacier surface in a predominantlykatabatic flow and with a distinctwind maximum below 13-m height are presented. The data werecollected using a 13-m high profilemast and two sonic anemometers (at about 2.5-m and 10-m heights).The spectra at frequencies belowthat of the turbulence range appear to deviate considerably fromthe curves obtained by Kaimal andco-workers during the 1968 Kansas experiment. The characteristicsof these deviations are compared tothe observations of others in surface-layers disturbed by anykind of large-scale outer-layer (orinactive) turbulence. In our case the disturbances arelikely to be induced by the highmountain ridges that surround the glacier. Moreover, the deviationsobserved in the cospectra seemto result from an, as yet, unspecified interaction between theinactive outer-layer turbulenceand the local surface-layer turbulence. Near the distinctwind maximum turbulence production ceasedwhile turbulence itself did not, probably the result ofturbulence transport from other levels. Consequently, we studied thelocal similarity relations using _w instead of u_* as an alternative velocity scale. Wellbelow the wind maximum, and for relatively low stability(0< Ri_g <0.2), the flow behaves accordingto well established local-scaling similarity relationshipsin the stable boundary layer. For higherstability (Ri_g > 0.2), and near or above the wind maximum, the boundary-layer structure conforms tothat of z-less stratification suggesting that the eddy sizeis restricted by the local stability ofthe flow. In line with this we observed that the sensibleheat fluxes relate remarkably well to thelocal flow parameters.     

Using Empirical Mode Decomposition to Filter Out Non-turbulent Contributions to Air–Sea Fluxes

A methodology based on Empirical mode decomposition (EMD) was used to filter out non-turbulent motions from measurements of atmospheric turbulence over the sea, aimed at reducing their contribution to eddy-covariance (EC) estimates of turbulent fluxes. The proposed methodology has two main objectives: (1) to provide more robust estimates of the fluxes of momentum, heat and CO\(_2\); and (2) to reduce the number of flux intervals rejected due to non-stationarity criteria when using traditional EC data processing techniques. The method was applied to measurements from a 28-day cruise (HALOCAST 2010) in the Eastern Pacific region. Empirical mode decomposition was applied to 4-h long time series data and used to determine the cospectral gap time scale, \(T_\mathrm{{gap}}\). Intrinsic modes of oscillation with characteristic periods longer than the gap scale due to non-turbulent motions were assumed and filtered out. Turbulent fluxes were then calculated for sub-intervals of length \(T_\mathrm{{gap}}\) from the filtered 4-h time series. In the HALOCAST data, the gap scale was successfully identified in 89% of the 4-h periods and had a mean of 37 s. The EMD approach resulted in the rejection of 11% of the flux intervals, which was much less than the 68% rejected when using standard filtering methods based on data non-stationarity. For momentum and sensible heat fluxes, the averaged difference in flux magnitude between the traditional and EMD approaches was small (3 and 1%, respectively). For the CO\(_2\) flux, the magnitude of EMD flux estimates was on average 16% less than fluxes estimated from linear detrended 10-min time series. These results provide evidence that the EMD method can be used to reduce the effects of non-turbulent correlations from flux estimates.     

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.     

Velocity and temperature spectra and cospectra in an unstable suburban atmosphere

Boundary-layer flow over very rough surfaces is poorly understood so the applicability of standard micrometeorological theory is uncertain. This study presents observations of the turbulent fluctuations of meteorological parameters over a suburban area. Even though the height of measurement is considered to be close to the junction between the inertial and roughness sub-layers, the wind and temperature spectra and the momentum and sensible heat flux cospectra are in good agreement with reference data from smoother surfaces. Recommendations are made concerning site requirements, height of measurement and averaging times for the study of turbulence and turbulent fluxes over suburban terrain.