湍流能量闭合法用于大气模式的初探

湍流能量（TKE）闭合方案是近年来大气边界层（PBL）模式中技推崇的方法,它具有11－2阶精度。大气环流模式的计算容量和垂直网络限制,不允许PBL方案过份精细。本文对TKE方案进行简化及网络稀疏化试验。并用Wangara资料验证。验证结果表明：此方案引入PBL内仅设置5层的大气环流模式是可行的。地面以上400米内模式结果与实测吻合。但中、上部位温计算值偏高。关于本方案在大尺度模式中的应用尚须进一步研究。     

Application of the E-ε closure model to simulations of mesoscale topographic effects

A mesoscale Planetary Boundary Layer (PBL) model with a simple turbulence closure scheme based on the turbulence kinetic energy (TKE) equation and the dissipation () equation is used to simulate atmospheric flow over mesoscale topography. Comparative studies with different parameterizations suggest that with a proper closure assumption for turbulence dissipation, the E-model can simulate the circulation induced by the mesoscale topography with results similar to those obtained using the E- model. On the other hand, the first-order closure using O'Brien's cubic interpolation for eddy diffusivities (K) generally produces much larger K profiles in the stable and the unstable regions, which is believed to be due to the overprediction of the height of the PBL. All models with the TKE equation yield quite similar ensemble mean fields, which are found to be little sensitive to the closure assumption for turbulence dissipation, though their predicted magnitudes of TKE and K may differ appreciably. A discussion on the diurnal evolution of the mesoscale topography-induced circulation and the spatial variations of the turbulence fluxes in the surface layer is also given based on the E- model results.     

Frontal substructures within the planetary boundary layer

A two-dimensional mesoscale model, extended by a TKE closure for the subgrid-scale terms and coupled with a soil model, is used to investigate the role of the Planetary Boundary Layer (PBL) for the development and the substructures of two different types of cold fronts. The effects of turbulent friction, large-scale (geostrophic) forcing and the diurnal variation of the terms of the surface energy balance (SEB) equation on the frontal development are studied by 10 different model runs. The ageostrophic cross-frontal circulation in the lowest two kilometres of a cold front results from friction as well as from large-scale forcing. The first one dominates the PBL processes and causes a special boundary-layer structure, which becomes apparent through the existence of seven characteristic zones defined for the x-z cross sections of potential temperature. The arrangement of these characteristic zones depends on the sense of rotation of the frictionally induced part of the ageostrophic circulation and hence on the direction of the along-front jet within the boundary layer. The daytime increase of the terms of the SEB equation for a midlatitude midsummer case leads to a strong enhancement of the frictionally induced cross-frontal circulation. The arrangement of the seven characteristic zones, however, is approximately conserved.     

不同边界层参数化方案对一次梅雨锋暴雨过程湍流交换特征模拟的影响

利用WRF模式结合不同的边界层参数化方案,对2007年7月3—5日发生在江淮流域的一次梅雨锋暴雨过程进行多组数值模拟试验。结果发现,边界层方案的选取对于降水的落区和强度模拟会产生较显著的影响;在降水率及地面要素的模拟上,各方案在降水中后期的模拟差异明显大于降水发生阶段;不同边界层方案的选取对于降水时段内的水平风场、垂直运动和假相当位温的垂直分布都产生影响,直接影响降水时空分布的模拟;不同方案都模拟出了在降水发生之后不同于晴空日变化的湍流动能垂直分布,经分析发现与局地较强的垂直风切变和近地面强湍流气团被抬升有关,而浮力项起着耗散作用;各方案的湍流交换特征与湍流动能特征基本吻合,相比于其他方案,MYJ方案在降水区域的湍流动能及湍流交换强度明显偏弱,对热通量的输送也偏弱;GBM方案在边界层内的湍流混合偏弱而在边界层以上湍流混合显著偏强,热通量输送在边界层以上的高度上误差明显,影响了对降水区域气象要素的模拟能力,仍需要进一步改进。      

A Study On Atmospheric Boundary-Layer Characteristics At Anand, India Using Lsp Experimental Data Sets

An attempt is made to study the planetary boundary layer (PBL) characteristics during the winter period at Anand (22.4°N, 72.6°E), a semi-arid region, which is located in the western part of India. A one-dimensional turbulent kinetic energy (TKE) closure model is used for the study. The structure of the PBL,which consists of profiles of zonal and meridional components of wind,potential temperature and specific humidity, is simulated. A one-dimensional soil heat and moisture transport parameterization scheme is incorporated for the accurate representation of the energy exchange processes at the soil-atmosphere interface. The diurnal variation of fluxes of sensible heat, latent heat, shortwave radiation, net radiation and soil flux, soil temperature at different depths, Richardson number and TKE at the height of the constant flux layer is studied. The model predictions are compared with the available observations obtained from a special Land Surface Processes (LSP) experiment.     

Sensitivity of PBL model predictions to model design and uncertainties in environmental inputs

A simple time-dependent one-dimensional model of the planetary boundary layer (PBL) is described and used to examine the degree to which model design decisions affect model output variables. The model's sensitivity to changes in the environmental conditions is also explored. Averages of the surface fluxes, near-ground wind speeds and other PBL properties from 48 h simulations are compared to control runs. The model-calculated surface fluxes are most sensitive, in decreasing order of importance, to the vertical grid spacing, the form of closure between the surface temperature and the atmosphere, the use of vertical diffusivity smoothing, the choice of maximum time step and choice of turbulence closure scheme. These fluxes are relatively insensitive to mixing-length scaling or choice of implicit time step weighting factor. Sensitivity to changes in soil type exceeds any of the design criteria tested. The modeled fluxes are moderately sensitive to small variations in the horizontal pressure gradient, to unsteadiness in the geostrophic wind and to variations in surface roughness. They are relatively insensitive to uncertainties in local vertical velocities and small (25%) variations applied separately to soil thermal diffusivity or heat capacity. The sensitivity of the average PBL depth (Z _i ) to model and environmental changes are similar to those of surface fluxes except thatZ _i is more sensitive to changes in mixing length, albedo and imposed vertical velocity then are the surface fluxes.     

A note on surface layer parameterizations in the weather research and forecast model

Two surface layer parameterization schemes along with five planetary boundary layer (PBL) schemes in the Weather Research and Forecasting model (WRF) are analyzed in order to evaluate the performance of the WRF model in simulating the surface variables and turbulent fluxes over an Indian sub-continent region. These surface layer schemes are based on the fifth-generation Pennsylvania State University—National Center for Atmospheric Research Mesoscale Model (MM5) parameterization; (a) Old MM5 scheme having Businger-Dyer similarity functions and (b) revised MM5 scheme utilizing the functions that are valid for full ranges of atmospheric stabilities. The study suggests that each PBL scheme can reproduce the diurnal variation of 2 m temperature, momentum flux and sensible heat flux irrespective of the surface layer scheme used for the simulations. However, a comparison of model-simulated values of surface variables and turbulent fluxes with observed values suggests that each PBL scheme is found to systematically over-estimate the nocturnal 2 m temperature and 10 m wind speed with both the revised and old schemes during stable conditions.     

Feedbacks in the Land-Surface and Mixed-Layer Energy Budgets

A mixed-layer model of the surface energy budget and the planetary boundary layer (PBL) is developed, based on the prognostic equations for soil temperature, mixed layer potential temperature and specific humidity and the growth and abrupt collapse of the PBL. Detailed parameterizations of the longwave radiative fluxes are included. The feedbacks in the uncoupled (i.e. surface energy budget with non-responding PBL) and coupled land surface and atmospheric mixed-layer energy budgets are examined. A simplified, time continuous, version of the model, in which the specific humidity budget is the balance of evapotranspiration and dry-air entrainment, and the PBL height is given by the lifted condensation level, is shown to be in good agreement with the complete model. By forcing the simplified model with daily mean rather than periodic solar radiation, an equilibrium model state is achieved where the fluxes are in close agreement with the daily mean fluxes corresponding to the periodic forcing. The model also agrees favorably with measurements from the FIFE field experiment. Feedbacks are examined using the equilibrium model state. The uncoupled and coupled model sensitivities with respect to the minimal stomatal resistance and the atmospheric specific humidity not only differ in magnitude, but in sign as well. This results puts into question the extent to which uncoupled land-surface models that are forced with atmospheric variables may be used in sensitivity studies.     

Nesting Turbulence in an Offshore Convective Boundary Layer Using Large-Eddy Simulations

The applicability of the one-way nesting technique for numerical simulations of the heterogeneous atmospheric boundary layer using the large-eddy simulation (LES) framework of the Weather Research and Forecasting model is investigated. The focus of this study is on LES of offshore convective boundary layers. Simulations were carried out using two subgrid-scale models (linear and non-linear) with two different closures [diagnostic and prognostic subgrid-scale turbulent kinetic energy (TKE) equations]. We found that the non-linear backscatter and anisotropy model with a prognostic subgrid-scale TKE equation is capable of providing similar results when performing one-way nested LES to a stand-alone domain having the same grid resolution but using periodic lateral boundary conditions. A good agreement is obtained in terms of velocity shear and turbulent fluxes, while velocity variances are overestimated. A streamwise fetch of 14 km is needed following each domain transition in order for the solution to reach quasi-stationary results and for the velocity spectra to generate proper energy content at high wavelengths, however, a pile-up of energy is observed at the low-wavelength portion of the spectrum on the first nested domain. The inclusion of a second nest with higher resolution allows the solution to reach effective grid spacing well within the Kolmogorov inertial subrange of turbulence and develop an appropriate energy cascade that eliminates most of the pile-up of energy at low wavelengths. Consequently, the overestimation of velocity variances is substantially reduced and a considerably better agreement with respect to the stand-alone domain results is achieved.     

Aircraft-Based Measurements Of Turbulence Structures In The Katabatic Flow Over Greenland

Turbulence structures in the katabatic flow in the stable boundary layer (SBL) over the ice sheet are studied for two case studies with high wind speeds during the aircraft-based experiment KABEG (Katabatic wind and boundary layer front experiment around Greenland) in the area of southern Greenland. The aircraft data allow the direct determination of turbulence structures in the katabatic flow. For the first time, this allows the study of the turbulence structure in the katabatic wind system over the whole boundary layer and over a horizontal scale of 80 km.The katabatic flow is associated with a low-level jet (LLJ), with maximum wind speeds up to 25 m s^-1. Turbulent kinetic energy (TKE) and the magnitude of the turbulent fluxes show a strong decrease below the LLJ. Sensible heat fluxes at the lowest level have values down to -25 W m^-2. Latent heat fluxes are small in general, but evaporation values of up to +13 W m^-2 are also measured. Turbulence spectra show a well-defined inertial subrange and a clear spectral gap around 250-m wavelength. While turbulence intensity decreases monotonously with height above the LLJ for the upper part of the slope, high spectral intensities are also present at upper levels close to the ice edge. Normalized fluxes and variances generally follow power-law profiles in the SBL.Terms of the TKE budget are computed from the aircraft data. The TKE destruction by the negative buoyancy is found to be very small, and the dissipation rate exceeds the dynamical production.     

一次暴雨过程大气边界层结构的数值分析

1979年7月28日河北唐山地区的强暴雨,10小时总降水量达430mm,降水强度大,从时间和空间上都非常集中。这次太平洋副热带高压北侧暖区发生的强暴雨引起国内气象界的普遍重视。游景炎,陆一强等对这次暴雨过程的大尺度环境、中尺度结构进行了详细分析并分别讨论了边界层急流、强对流性云团以及地形等的作用。游景炎根据地面天气图的分析,发现明显的雷暴高压和中尺度低压。这种中尺度系统在大气边界层内如何表现呢?这次暴雨系统延伸多高呢?各种物理量输送和分布的特性     

Sensitivity of Tropical Cyclone Intensity and Structure to Planetary Boundary Layer Parameterization

The advanced weather research and forecasting model is used to investigate the influence of planetary boundary layer (PBL) processes on intensity and structure of the storm Phailin (2013). Five simulations are conducted with five PBL schemes; Yonsei University (YSU), Mellor?Yamada?Nakanishi?Niino order2.5 (MYNN2), Assymetric Convective Model2 (ACM2), Medium Range Forecast (MRF), and Bougeault and Lacarrere (BouLac). The simulation duration includes the pre???intensification and rapid intensification phase of Phailin before landfall. Results indicate that during the pre???intensification phase, storm’s track and intensity are not much sensitivity to PBL but structural changes are noted. A significant sensitivity of track and intensity to PBL parameterizations are found during rapid intensification phase. BouLac and MRF produced two extremes with 39 hPa intense and 16 km compact storm for BouLac than MRF. Further analysis reveals an outward movement of air parcel just above the boundary layer which causes spin-down for YSU and MYNN2. BouLac is associated with stronger eddy diffusivity and moisture fluxes within the boundary layer and stronger cyclonic vorticity just above the boundary layer than other experiments. Stronger cyclonic vorticity above the boundary layer in BouLac favors intense updraft, facilitating more moisture transport from the boundary layer to upper layers aiding stronger secondary circulation and robustly intensifying the storm. A relatively deeper and drier inflow layer associated with weaker cyclonic vorticity just above the boundary layer reduces the moisture transport and weaken the secondary circulation for MRF than others.     

Turbulence Characteristics of the Shear-Free Convective Boundary Layer Driven by Heterogeneous Surface Heating

Large-eddy simulations (LESs) are employed to investigate the turbulence characteristics in the shear-free convective boundary layer (CBL) driven by heterogeneous surface heating. The patterns of surface heating are arranged as a chessboard with two different surface heat fluxes in the neighbouring patches, and the heterogeneity scale Λ in four different cases is taken as 1.2, 2.5, 5.0 and 10.0 km, respectively. The results are compared with those for the homogeneous case. The impact of the heterogeneity scale on the domain-averaged CBL characteristics, such as the profiles of the potential temperature and the heat flux, is not significant. However, different turbulence characteristics are induced by different heterogeneous surface heating. The greatest turbulent kinetic energy (TKE) is produced in the case with the largest heterogeneity scale, whilst the TKE in the other heterogeneous cases is close to that for the homogeneous case. This result indicates that the TKE is not enhanced unless the scale of the heterogeneous surface heating is large enough. The potential temperature variance is enhanced more significantly by a larger surface heterogeneity scale. But this effect diminishes with increasing CBL height, which implies that the turbulent eddy structures are changed during the CBL development. Analyses show that there are two types of organized turbulent eddies: one relates to the thermal circulations induced by the heterogeneous surface heating, whilst the other identifies with the inherent turbulent eddies (large eddies) induced by the free convection. At the early stage of the CBL development, the dominant scale of the organized turbulent eddies is controlled by the scale of the surface heterogeneity. With time increasing, the original pattern breaks up, and the vertical velocity eventually displays horizontal structures similar to those for the homogeneous heating case. It is found that after this transition, the values of λ/z _i (λ is the dominant horizontal scale of the turbulent eddies, z _i is the boundary-layer height) ≈1.6, which is just the aspect ratio of large eddies in the CBL.     

Intercomparison of Planetary Boundary-Layer Parametrizations in the WRF Model for a Single Day from CASES-99

This study compares five planetary boundary-layer (PBL) parametrizations in the Weather Research and Forecasting (WRF) numerical model for a single day from the Cooperative Atmosphere-Surface Exchange Study (CASES-99) field program. The five schemes include two first-order closure schemes—the Yonsei University (YSU) PBL and Asymmetric Convective Model version 2 (ACM2), and three turbulent kinetic energy (TKE) closure schemes—the Mellor–Yamada–Janjić (MYJ), quasi-normal scale elimination (QNSE), and Bougeault–Lacarrére (BouLac) PBL. The comparison results reveal that discrepancies among thermodynamic surface variables from different schemes are large at daytime, while the variables converge at nighttime with large deviations from those observed. On the other hand, wind components are more divergent at nighttime with significant biases. Regarding PBL structures, a non-local scheme with the entrainment flux proportional to the surface flux is favourable in unstable conditions. In stable conditions, the local TKE closure schemes show better performance. The sensitivity of simulated variables to surface-layer parametrizations is also investigated to assess relative contributions of the surface-layer parametrizations to typical features of each PBL scheme. In the surface layer, temperature and moisture are more strongly influenced by surface-layer formulations than by PBL mixing algorithms in both convective and stable regimes, while wind speed depends on vertical diffusion formulations in the convective regime. Regarding PBL structures, surface-layer formulations only contribute to near-surface variability and then PBL mean properties, whereas shapes of the profiles are determined by PBL mixing algorithms.     

An efficient PBL model for global circulation models-design and validation

An efficient, pianetary boundary layer (PBL) model is developed and validated with empirical data for applications in general circulation models (GCMs). The purpose of this PBL model is to establish the turbulent surface fluxes as a function of the principal external PBL parameters in a numerically efficient way. It consists of a surface layer and a mixed layer matched together with the conditions of constant momentum and heat flux at the interface. An algebraic solution to the mean momentum equations describes the mixed-layer velocity profile and thus determines the surface wind vector. The velocity profile is globally valid by incorporating the effect of variable Coriolis force without becoming singular at the equator. Turbulent diffusion depends on atmospheric stability and is modeled in the surface layer by a drag law and with first-order closure in the mixed layer. Radiative cooling in the stably stratified PBL is considered in a simple manner. The coupled system is solved by an iterative method. In order to preserve the computational efficiency of the large-scale model, the PBL model is implemented into the GISS GCM by means of look-up tables with the bulk PBL Richardson number, PBL depth, neutral drag coefficient, and latitude as independent variables.A validation of the PBL model with observed data in the form of Rossby number similarity theory shows that the internal feedback mechanisms are represented correctly. The model, however, underpredicted the sensible heat-flux. A subsequent correction in the turbulence parameterization yields better agreement with the empirical data. The behavior of the principal internal PBL quantities is presented for a range of thermal stabilities and latitudes.     

荒漠绿洲边界层结构的数值模拟

陆-气相互作用和中小尺度天气系统的研究中,水平不均匀边界层和水平不均匀地表的强迫作用都是重要的物理过程.本文用已建立的陆面过程与大气边界层耦合模式(BLCM),较详细地研究了草地周围为荒漠(半沙漠)的地表植被不均匀而造成的边界层结构特征和局地环流及其昼夜变化.通过边界层顶影响自由大气不同尺度的运动.模式结果揭示出最强的上升和下沉运动是发生在荒漠-草地间动力和热力不连续的界面附近,且呈现出不对称性.     

A subgrid-scale model for large-eddy simulation of planetary boundary-layer flows

A long-standing problem in large-eddy simulations (LES) of the planetary boundary layer (PBL) is that the mean wind and temperature profiles differ from the Monin-Obukhov similarity forms in the surface layer. This shortcoming of LES has been attributed to poor grid resolution and inadequate sub-grid-scale (SGS) modeling. We study this deficiency in PBL LES solutions calculated over a range of shear and buoyancy forcing conditions. The discrepancy from similarity forms becomes larger with increasing shear and smaller buoyancy forcing, and persists even with substantial horizontal grid refinement. With strong buoyancy forcing, however, the error is negligible.In order to achieve better agreement between LES and similarity forms in the surface layer, a two-part SGS eddy-viscosity model is proposed. The model preserves the usual SGS turbulent kinetic energy formulation for the SGS eddy viscosity, but it explicitly includes a contribution from the mean flow and a reduction of the contributions from the turbulent fluctuations near the surface. Solutions with the new model yield increased fluctuation amplitudes near the surface and better correspondence with similarity forms out to a distance of 0.1–0.2 times the PBL depth, i.e., a typical surface-layer depth. These results are also found to be independent of grid anisotropy. The new model is simple to implement and computationally inexpensive.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Velocity profiles,the resistance law and the dissipation rate of mean flow kinetic energy in a neutrally and stably stratified planetary boundary layer

The logarithmic + polynomial approximation is suggested for vertical profiles of velocity components in a planetary boundary layer (PBL) at neutral and stable stratification. The resistance law functions A and B are determined on the basis of this approximation, using integral relations derived from the momentum equations, the Monin-Obukhov asymptotic formula for the wind profile in a stably stratified near-surface layer and the known expressions for the PBL depth. This result gives a realistic and convenient method for calculating the surface friction velocity and direction and the total dissipation rate of mean flow kinetic energy in terms of geostrophic velocity, buoyancy flux at the surface, the roughness parameter and the Coriolis parameter. In the course of these derivations a review is given of current views on the main problems of the neutral and stable PBL.