风切变对边界层对流影响的大涡模拟研究

利用"西北干旱区陆-气相互作用野外观测实验"加密观测期间在敦煌站的观测资料以及大涡模式,模拟了对流边界层的发展,以及示踪物从混合层向残留层传输的时空变化。模拟的对流边界层的结构及演变特征与实测结果基本一致。进一步通过有风切变和无风切变的敏感性数值试验,研究了风切变对垂直速度、位温和示踪物浓度的水平分布以及示踪物传输高度的影响。研究结果表明,在有风切变的试验中(甚至风切变仅存在于近地层中),对流边界层的增长加强,而且示踪物被传输的高度也较高。与浮力驱动的对流边界层相比,由浮力和风切变共同驱动的边界层中上升气流较弱而下沉气流较强,但前者的上升气流与下沉气流的分布在垂直方向上更为倾斜。由于夹卷作用的增强,浮力和风切变共同驱动的对流边界层较浮力驱动的对流边界层暖。在夹卷层,浮力和风切变共同驱动的边界层对流的上升气流和下沉气流都比浮力驱动的边界层对流中的强,而且垂直速度的概率密度函数分布也较对称,其位温和示踪物浓度的概率密度函数分布也比浮力驱动的边界层中的平直。对湍流动能收支的分析也表明风切变对湍流动能有重要影响,尤其对夹卷层中的湍流动能切变产生项影响较大。示踪物浓度的概率密度函数垂直分布显示,浮力驱动的边界层中示踪物浓度随高度变化较小,而浮力和风切变共同驱动的边界层中示踪物浓度随高度递减,但是示踪物传输的高度比较高。     

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

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

A TKE-dissipation model for the atmospheric boundary layer

The dissipation, , of turbulent kinetic energy (TKE) is a key parameter in atmospheric boundary-layer (ABL) models. Besides being a sink for momentum, it is often used together with the TKE to define an internal turbulence time scale for closure relations. A prognostic formulation for the dissipation of TKE is formulated, based on isotropic tensor modeling methods. The formulation is coupled to a level 2.5 second-order closure model and evaluated against measurements taken in horizontally homogeneous conditions, as well as against a tailored length-scale formulation. A formulation suitable for convective as well as neutral and stable ABLs is suggested.On leave from Department of Meteorology, Uppsala University, P.O. Box 516, S-751 20 Uppsala, Sweden.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Diurnal Variations of the Flux Imbalance Over Homogeneous and Heterogeneous Landscapes

It is well known that the sum of the turbulent sensible and latent heat fluxes as measured by the eddy-covariance method is systematically lower than the available energy (i.e., the net radiation minus the ground heat flux). We examine the separate and joint effects of diurnal and spatial variations of surface temperature on this flux imbalance in a dry convective boundary layer using the Weather Research and Forecasting model. Results show that, over homogeneous surfaces, the flux due to turbulent-organized structures is responsible for the imbalance, whereas over heterogeneous surfaces, the flux due to mesoscale or secondary circulations is the main contributor to the imbalance. Over homogeneous surfaces, the flux imbalance in free convective conditions exhibits a clear diurnal cycle, showing that the flux-imbalance magnitude slowly decreases during the morning period and rapidly increases during the afternoon period. However, in shear convective conditions, the flux-imbalance magnitude is much smaller, but slightly increases with time. The flux imbalance over heterogeneous surfaces exhibits a diurnal cycle under both free and shear convective conditions, which is similar to that over homogeneous surfaces in free convective conditions, and is also consistent with the general trend in the global observations. The rapid increase in the flux-imbalance magnitude during the afternoon period is mainly caused by the afternoon decay of the turbulent kinetic energy (TKE). Interestingly, over heterogeneous surfaces, the flux imbalance is linearly related to the TKE and the difference between the potential temperature and surface temperature, ΔT; the larger the TKE and ΔT values, the smaller the flux-imbalance magnitude.     

Analytical Model of Roll Vortices in the Convective Boundary Layer

Analytical solutions of convective waves in the convective boundary layer (CBL) were obtained with two-layer linearized atmospheric equations including Rayleigh friction, which represents the turbulent viscosity in the lower CBL. The analytical model shows that the interaction between the convection in the lower layer and gravity waves in the upper layer is one of the causes for the formation of convective bands. The flow and temperature fields obtained by the analytical model present the main characteristics of convective bands found in field observations. We have also investigated the influences of atmospheric conditions on the characteristics of the bands. Results accord with previous knowledge about these phenomena.     

一次梅雨锋暴雨的边界层热通量输送与湍流动能收支分析

利用中尺度WRF模式对于2007年7月一次典型的梅雨锋暴雨过程进行了高分辨率数值模拟,对于边界层内的热通量输送和湍流动能的时空变化特征,以及湍流动能各收支项的分布及变化特征进行了分析。结果表明,降水发生时段内边界层热通量和湍流动能的时空分布特征与晴空日变化特征表现出显著不同,潜热通量随高度自下而上呈现"正—负—正"的分布,感热通量以负值为主,负值中心高度与潜热通量由负转正的高度相对应,湍流动能的发展高度与持续时间都有所增加,降水区近地面湍流动能弱于其他区域,但是在468 m以上高度则显著强于其他区域。降水区湍流动能的来源主要是平均风切变所产生的机械湍流,浮力作用与粘性耗散在降水期间消耗湍流动能,湍流输送作用将低层的湍流动能输送至较高的高度,使低层减小而高层增大,临界高度与湍流动能的大值中心高度对应。     

Evaluation of the Diurnal Cycle in the Atmospheric Boundary Layer Over Land as Represented by a Variety of Single-Column Models: The Second GABLS Experiment

We present the main results from the second model intercomparison within the GEWEX (Global Energy and Water cycle EXperiment) Atmospheric Boundary Layer Study (GABLS). The target is to examine the diurnal cycle over land in today??s numerical weather prediction and climate models for operational and research purposes. The set-up of the case is based on observations taken during the Cooperative Atmosphere-Surface Exchange Study-1999 (CASES-99), which was held in Kansas, USA in the early autumn with a strong diurnal cycle with no clouds present. The models are forced with a constant geostrophic wind, prescribed surface temperature and large-scale divergence. Results from 30 different model simulations and one large-eddy simulation (LES) are analyzed and compared with observations. Even though the surface temperature is prescribed, the models give variable near-surface air temperatures. This, in turn, gives rise to differences in low-level stability affecting the turbulence and the turbulent heat fluxes. The increase in modelled upward sensible heat flux during the morning transition is typically too weak and the growth of the convective boundary layer before noon is too slow. This is related to weak modelled near-surface winds during the morning hours. The agreement between the models, the LES and observations is the best during the late afternoon. From this intercomparison study, we find that modelling the diurnal cycle is still a big challenge. For the convective part of the diurnal cycle, some of the first-order schemes perform somewhat better while the turbulent kinetic energy (TKE) schemes tend to be slightly better during nighttime conditions. Finer vertical resolution tends to improve results to some extent, but is certainly not the solution to all the deficiencies identified.     

珠江三角洲附近地区强风暴的中尺度环境场若干特征分析

用滤波方法对两次强风暴天气过程进行尺度分离,得到与强风暴天气区相对应的中尺度扰动流场,分析结果表明:(1)强风暴天气—般出现在前期的层结不稳定区域内非热成风的大值中心附近。(2)在强对流的活跃阶段,在行星边界层的较下部和对流层的中部,分别有一无辐散层;从地面往上,辐散辐合依次相间出现;在强对流区的上空,存在明显的中尺度辐合和辐散中心。 此外,文中对中尺度扰动和强对流天气的成因作了初步讨论。      

A comparison of parameterized vs. measured transilient mixing coefficients for a convective mixed layer

The parameterization for transilient turbulence coefficients suggested by Stull and Driedonks (SD. 1987) is tested against the large-eddy simulations (LES) of Ebert et al. (ESS, 1989) for the special case of an idealized convective boundary layer. The SD parameterization is based on a nonlocal approximation to the turbulence kinetic energy (TKE) equation, and requires turbulent exchange (i.e., the matrix of transilient mixing coefficients is assumed to be symmetric) and dominance of the smaller eddies (i.e., elements closer to the main diagonal of the matrix are greater). Measurements from the LES model, however, show that the transilient matrix is asymmetric in convective situations, with larger eddies dominating.Mean-state conditions such as the deep convective mixing and mixed-layer growth are satisfactorily described by the parameterization, but the surface layer is too deep and the entrainment zone thickness is poorly defined. Turbulence properties such as skewed vertical velocity distributions are not possible within the constraint of a symmetric matrix, and partial convective overturning is also not possible because of the restriction that small eddies dominate. Future improved parametcrizations might continue to be based on the TKE equation, but should allow transilient matrix asymmetry.Research was performed while the author worked at the Institute of Atmospheric Physics, German Aerospace Research Establishment (DLR), D-8031 Weling. Germany.     

A one-dimensional model for the parameterization of deep convection in the ocean

A one-dimensional penetrative plume model has been constructed to parameterize the process of deep convection in ocean general circulation models (OGCMs). This research is motivated by the need for OGCMs to better model the production of deep and intermediate water masses. The parameterization scheme takes the temperature and salinity profiles of OGCM grid boxes and simulates the subgrid-scale effects of convection using a one-dimensional parcel model. The model moves water parcels from the surface layer down to their level of neutral buoyancy, simulating the effect of convective plumes. While in transit, the plumes exchange water with the surrounding environment; however, the bulk of the plume water mass is deposited at e level of neutral buoyancy. Weak upwelling around the plumes is included to maintain an overall mass balance. The process continues until the negative buoyant energy of the one-dimensional vertical column is minimized. The parameterized plume entrainment rate, which plays a central role in the parameterization, is calculated using modified equations based on the physics of entraining buoyant plumes. This scheme differs from the convective adjustment techniques currently used in OGCMs, because the parcels penetrate downward with the appropriate degree of mixing until they reach their level of neutral stability.     

一次东北冷涡衰退阶段暴雨成因分析

东北冷涡暴雨一般发生在冷涡发展阶段, 冷涡衰退阶段发生暴雨的几率不高, 在日常业务预报中往往容易忽视冷涡衰退阶段暴雨的预报, 容易出现漏报的现象。作者结合观测资料和数值模拟结果, 对2005年6月21日发生在吉林省中东部地区的一次暴雨过程进行分析。结果表明: 这次暴雨发生在东北冷涡衰退阶段, 强降水区处于非常弱的对流不稳定区, 低空不同层次上急流的耦合与脉动, 在时间和空间上形成“接力” 之势, 给雨区带来特定的水汽配置并在雨区上空形成合适的不稳定环境, 导致雨区上空湿位涡的强烈增长, 在弱对流不稳定和条件对称不稳定的共同作用下, 层状云中形成对流并发展, 产生了暴雨天气。重力惯性波是暴雨产生的可能触发机制。     

The Neutral, Barotropic Planetary Boundary Layer, Capped by a Low-Level Inversion

The presence of a low-level, capping inversion layer will affect the height and structure of the planetary boundary layer (PBL). Results from models of varying levels of sophistication, including analytical, turbulent kinetic energy (TKE), second-order closure (SOC), large-eddy simulation (LES) and direct numerical simulation (DNS) models, are used to investigate this influence for the neutral, barotropic PBL. Predicted and observed profiles of stress and geostrophic departure components, and integral measures, such as the parameters of Rossby-number similarity theory, are compared for the KONTUR, Marine Stratocumulus, JASIN, Leipzig, Pre-Wangara and Upavon field experiments.Analytical models of the equilibrium value of inversion height z_i, which depend on the surface friction velocity u_*, and both the Coriolis parameter f and the free-flow Brunt-Väisälä frequency N, are found to give reasonable estimates of the PBL height. They also indicate that only the KONTUR and Marine Stratocumulus experiments were strongly influenced by N. More quantitative comparisons would require larger, more comprehensive datasets. The effects of the presence of a capping inversion on the profile structure were found to be insignificant for h_* = |f|z_i/u_* > 0.15.The simple analytical model performed quite well over all values of h_*; it predicted the profiles of the longitudinal stress component (in the direction of the surface stress) better than the lateral component. The more advanced models performed well for small values of h_* (for flow over the sea), but systematically underestimated the cross-isobaric angle for flow over land. These models predicted the profiles of the lateral stress component better than the longitudinal component. The profiles of the analytical model agreed with those of the advanced models when the constant eddy viscosity of the outer layer was increased.Agreement with DNS was achieved by increasing the eddyviscosity of the analytical model by a factor of 5.Zilitinkevich and Esau(2002, Boundary-Layer Meteorology, 104, 371–379)suggest that the neutral, barotropic values of A and B of Rossby-numbersimilarity theory are not universal constants, but depend on the ratio N/|f|. The dependence for A and B is calculated using the analytical model and TKE models. Over the sea (h_* 0.1; N/|f| 100, where we have used the Zilitinkevich-Esau relation to convert between h_* and N/|f|) there is agreement between the model predictions and observations; however over land where the equilibrium boundary-layer height is greater (h_* 0.35; N/|f| 10) the inconsistency between the advanced model predictions (TKE, SOC, LES, and DNS) and observations, as noted previously by Hess and Garratt, still exists. We attribute this disagreement to violations of the strict assumptions of steady, horizontally homogeneous, neutral, barotropic conditions implicit in the observations. At small values of z_i and a strongly stable background stratification (h_* 0.04; N/|f| 1000) both the TKE and analytical models predict that A and B depend significantly on h_*, however observations are unavailable to confirm these predictions. Zilitinkevich and Esau call this case the `long-lived near-neutral PBL', and state that it is found in cold weather at high latitudes.     

Observational analysis of moisture evolution and variability in the boundary layer during the dryline on 22 May 2002

An observational analysis of boundary layer moisture evolution during the dryline on 22 May 2002 is presented. This dryline occurred during the International H₂O Project (IHOP) and was well observed by a variety of instruments at the intensive observing region (IOR), Homestead site. Although the observed strong upward air motion and the well-mixed boundary layer favored convection, the dryline did not trigger a convective storm. Several operational and research forecast models predicted deep convection at the IOR. High spatial and temporal resolution observational data from National Aeronautics and Space Administration (NASA) lidar instruments, Atmospheric Emitted Radiance Interferometer (AERI), sounding profiles and simulation results are used to investigate the role of moisture during this dryline. It is hypothesized that in addition to convection and lifting, abundant, deep and persistent moisture is required for a dryline to trigger convection. The possible reason why the dryline failed to trigger convection over the IOR is discussed.     

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.     

复杂山区上空垂直速度场和热力对流活动的理想数值模拟

利用英国气象局高分辨率的边界层数值模式BLASIUS,针对中国西北一个复杂山区进行了一系列的理想数值模拟,分析了在不同天气条件下山区上空的垂直速度场分布和对流特征以及地形对热力对流活动的影响,同时讨论了与地形有关的对流触发机制。模式结果表明,复杂山区的垂直运动在稳定层结和风速较大的情况下较易预测,而在中性层结下,山区上空的垂直运动分布随机性强。在Froude数小于0.5的条件下,气流往往被山峰阻塞而在迎风坡造成地形强迫和辐合性抬升,从而易在迎风坡触发深对流活动;在背风坡则由于迎风坡的绕流重新辐合而造成垂直运动。绕流的辐合是触发深对流活动的另一重要因子。在大风或Froude数较大的条件下,地形重力波容易在山地下游被激发。地形重力波与对流活动的相互影响在模式中清楚可见。在适当的条件下,重力波除了可以与对流活动相耦合从而使气团上升到更高的高度外,重力波的走向很可能会影响到深对流系统的传播路径。研究还发现稳定度对相邻两条对流线之间的距离长短也有影响。稳定度较小时,相邻两条对流线之间的平均宽度趋向变大而单个对流线的强度也相应变大。定量化的结论和理论升华值得进一步的数值模拟研究。     

An ensemble of dynamically identical thermals and vertical profiles of turbulent moments in the convective surface layer of atmosphere

An ensemble of convective thermals is distinguished from the surface layer of penetrative turbulent convection over a heated horizontally uniform surface. For an isolated convective element, an integral model of a quasi-stationary spontaneous jet is proposed which admits an exact analytical solution. A simple statistical model is constructed for an ensemble of dynamically identical thermals. In this model, convective thermals ascend in a static environment, their dynamic parameters are described by the equations of an isolated quasi-stationary jet, and their diameters are stochastic. It is shown that the ensemble of thermals rising in a horizontally homogeneous environment forms surface-layer turbulent moments. The analytical relationships for higher turbulent moments of vertical velocity and temperature are compared with experimental data from the second to the forth order inclusive.     

A special MJO event with a double Kelvin wave structure

The second Madden–Julian Oscillation (MJO) event during the field campaign of the Dynamics of the MJO/Cooperative Indian Ocean Experiment on Intraseasonal Variability in the Year 2011 (DYNAMO/CINDY2011) exhibi ted an unusual double rainband structure. Using a wavenumber-frequency spectral filtering method, we unveil that this double rainband structure arises primarily from the Kelvin wave component. The zonal phase speed of the double rainbands is about 7.9 degree per day in the equatorial Indian Ocean, being in the range of convectively coupled Kelvin wave phase speeds. The convection and circulation anomalies associated with the Kelvin wave component are characterized by two anomalous convective cells, with low-level westerly (easterly) and high (low) pressure anomalies to the west (east) of the convective centers, and opposite wind and pressure anomalies in the upper troposphere. Such a zonal wind–pressure phase relationship is consistent with the equatorial free-wave dynamics. While the free-atmospheric circulation was dominated by the first baroclinic mode vertical structure, moisture and vertical motion in the boundary layer led the convection.The convection and circulation structures derived based on the conventional MJO filter show a different characteristic. For example, the phase speed is slower (about 5.9 degree per day), and there were no double convective branches. This suggests that MJO generally involves multi-scales and it is incomplete to extract its signals by using the conventional filtering technique.     

Spatial representativeness of single tower measurements and the imbalance problem with eddy-covariance fluxes: results of a large-eddy simulation study

A large-eddy simulation (LES) study is presented that investigates the spatial variability of temporal eddy covariance fluxes and the systematic underestimation of representative fluxes linked to them. It extends a prior numerical study by performing high resolution simulations that allow for virtual measurements down to 20 m in a convective boundary layer, so that conditions for small tower measurement sites can be analysed. It accounts for different convective regimes as the wind speed and the near-surface heat flux are varied. Moreover, it is the first LES imbalance study that extends to the stable boundary layer. It reveals shortcomings of single site measurements and the necessity of using horizontally-distributed observation networks. The imbalances in the convective case are attributed to a locally non-vanishing mean vertical advection due to turbulent organised structures (TOS). The strength of the TOS and thus the imbalance magnitude depends on height, the horizontal mean wind and the convection type. Contrary to the results of a prior study, TOS cannot generally be responsible for large energy imbalances: at low observation heights (corresponding to small towers and near-surface energy balance stations) the TOS related imbalances are generally about one order of magnitude smaller than those in field experiments. However, TOS may cause large imbalances at large towers not only in the case of cellular convection and low wind speeds, as found in the previous study, but also in the case of roll convection at large wind speeds. In the stably stratified boundary layer for all observation heights neither TOS nor significant imbalances are observed. Attempting to reduce imbalances in convective situations by applying the conventional linear detrending method increases the systematic flux underestimation. Thus, a new filter method is proposed.     

青奥开幕式期间南京地区浅层对流成因及逆向移动特征分析

利用常规观测资料、加密自动站资料、多普勒雷达资料及NCEP 1°×1°再分析资料,针对第二届青年奥运会开幕式期间形成发展于对流层低层偏东气流中的对流引发南京局地短时强降水天气进行了分析,探讨了低层偏东气流中浅层弱对流的形成机制、降水期间伴随在降水云系中逆向移动的两类对流单体的雷达特征和移动原因。结果表明,降水期间,南京上空对流层低层维持偏东风,中层以上则为西南风,两者之间的过渡区对应一干层,该干层由对流层中层的一支偏北气流携至的干空气形成,叠置于对流层低层较浅薄偏东暖湿气流之上,促进了对流层中低层对流不稳定层结的发展。地面风场扰动形成的局地辐合及地面非锋性斜压带激发了对流层低层偏东气流中对流的形成。镶嵌于降水回波中逆向移动的两类对流单体结构差异明显,两类对流单体质心高度、垂直伸展厚度与所在高度层中所盛行的背景风场决定了对流单体的传播与移动。     

低压倒槽影响下积层混合云形成过程模拟研究

利用多普勒雷达资料和中尺度天气预报模式wRF（Weather Research and Forecasting）模拟结果,对2009年5月9—10日发生在太原及其周边地区的一次积层混合云降水形成过程进行分析。结果表明,在积层混合云的形成初期,局地对流云得到发展,随着其强度不断增强,与周围云发生并合过程（包括局地单体对流的并合、积云团的并合和积层混合云内强中心的并合）,形成范围较大的积层混合云云系。局地单体对流和积云团的并合可带来云体的爆发性增长,霰含量、雨水含量大幅增加。积层混合云内强中心的并合对降水强度影响不大,但有利于降水面积扩大。低压倒槽和弱冷锋是此次积层混合云形成和维持的主要影响因素。低压倒槽有利于低层大范围不稳定能量的积累,风向切变有利于近距离云团的发展和并合,山地动力和热力作用有利于局地对流单体、积层混合云内强中心的形成和加强。