中尺度涡旋混合参数空间变化对物理场以及CFC-11模拟的影响

本文选用中国科学院大气物理研究所全球海洋模式(LICOM),对中尺度涡旋参数化方案(GM90方案)中等密度扩散系数和等密度面厚度扩散系数(统称为涡旋扩散系数A_ρ)对物理场及CFC-11(一氟三氯甲烷)分布的影响进行了研究。本文做了两个试验,即涡旋扩散系数采用常系数方式(对照试验)和采用在非绝热层以下A_ρ随海洋浮力频率垂直变化的参数化方案(浮力试验)。模拟结果表明,依浮力频率垂直变化的方案对模式物理场的模拟能力有一定程度的提升,如南极绕极流的输送强度比常系数方案增大了约20%~30%,与观测事实更接近;浮力试验对对照试验中过强的南极中层水有一定的削弱作用,使得模式对南大洋高纬次表层位密度的模拟有一定的改善。稍有不足的是,浮力试验对南极底层水也有一定的削弱,使得2000~3000 m深度位密度模拟较常系数方案偏低。通过对CFC-11分布、存储以及输送的研究发现,次网格参数取值的不同对南大洋CFC-11模拟情况有较大影响。浮力试验加大了南北高纬CFC-11海表的吸收通量,对南极大陆周边海域向南大洋主储藏区(34°S~60°S)的CFC-11输送能力有一定的增强,使得南大洋对CFC-11储藏量增大,大部分海区与观测资料更接近。通过CFC-11断面分析,浮力试验对南大洋上层海洋位密度模拟的改善使得CFC-11垂直结构与观测更接近。     

The Role of Shear in the Morning Transition Boundary Layer

We use large-eddy simulation (LES) to better define the early stages of the morning transition boundary layer. Previous LES studies relating to the morning transition boundary layer focus on the role of the entraining convective boundary layer (CBL). By using a combination of different domain sizes and grid lengths, the full evolution from the stable boundary layer (SBL) to the CBL is modelled here. In the early stages of the morning transition the boundary layer is shown to be a combination of a shallow mixed layer capped by a significant shear driven stable boundary layer (the so-called mixed CBL–SBL state). The mixed CBL–SBL state is the key to understanding the sensitivity to shear. Turbulent kinetic energy budgets also indicate that it is shear driven. The negative flux from the mixed CBL–SBL state extends much further above the minimum than is typically found for the CBL later in the day, and the depth of penetration scales as w _m /N _i , where w _m is the combined friction and convective velocity scale and N _i the static stability at the inversion top.     

当前对流参数化研究存在问题的讨论

基于对流参数化的发展过程,讨论了当前对流参数化方案中存在的问题,着重分析了大尺度水汽辐合型和质量通量型方案的缺陷,同时对东亚季风区和我国的研究现状进行了分析,提出了通用对流参数化方案的设计思想和应该考虑的原则,具有重要的参考借鉴意义。     

利用多普勒天气雷达估算对流云火箭增雨防雹用弹量的方案

利用新一代多普勒天气雷达产品资料,结合人工影响天气原理,引用人工影响天气数值实验的人工冰核浓度等有关参数,将垂直累积液态含水量(VIL)及回波顶高度(ET)等导出产品合理地应用于火箭增雨防雹作业用弹量计算公式中,解决了以前人为或间接估计对流云作业区体积及其含水量的难题,提出客观定量计算人工增雨防雹作业用弹量的初步方案,并以其计算的用弹量为主要参数设计火箭作业参数指挥界面.经2005年一次增雨作业应用检验,证明结果切实可行,对节约人工影响天气成本,有效指导对流云火箭增雨防雹作业具有一定的意义.     

Doppler Lidar Measurements of Vertical Velocity Spectra in the Convective Planetary Boundary Layer

We utilized a Doppler lidar to measure spectra of vertical velocity w from 390m above the surface to the top of the daytime convective boundary layer (CBL). The high resolution 2μm wavelength Doppler lidar developed by the NOAA Environmental Technology Laboratory was used to detect the mean radial velocity of aerosol particles. It operated continuously during the daytime in the zenith-pointing mode for several days in summer 1996 during the Lidars-in-Flat-Terrain experiment over level farmland in central Illinois, U.S.A. The temporal resolution of the lidar was about 1 s, and the range-gate resolution was about 30m. The vertical cross-sections were used to calculate spectra as a function of height with unprecedented vertical resolution throughout much of the CBL, and, in general, we find continuity of the spectral peaks throughout the depth of the CBL. We compare the observed spectra with previous formulations based on both measurements and numerical simulations, and discuss the considerable differences, both on an averaged and a case-by-case basis. We fit the observed spectra to a model that takes into account the wavelength of the spectral peak and the curvature of the spectra across the transition from low wavenumbers to the inertial subrange. The curvature generally is as large or larger than the von Kármán spectra. There is large case-to-case variability, some of which can be linked to the mean structure of the CBL, especially the mean wind and the convective instability. We also find a large case-to-case variability in our estimates of normalized turbulent kinetic energy dissipation deduced from the spectra, likely due for the most part to a varying ratio of entrainment flux to surface flux. Finally, we find a relatively larger contribution to the low wavenumber region of the spectra in cases with smaller shear across the capping inversion, and suggest that this may be due partly to gravity waves in the inversion and overlying free atmosphere.     

台风“云娜”在近海强度变化及结构特征的数值研究Ⅰ:云微物理参数化对云结构及降水特征的影响

首先对AREM模式模拟的台风基本结构和云结构进行验证,检验了模拟结果的可靠性.在此基础上,设计了5组试验来研究云微物理参数化方案对台风"云娜"云结构及降水特征的影响.试验设计主要突出冰相云微物理过程、云微物理特征引发的冷却效应以及霰下落速度的重要性.结果表明:云微物理参数化过程对云的发展和降水特征的影响更为显著.各试验的水凝物分布和强度不同,降水类型和强度存在较大差异,由此引起的云中热力结构也有较大区别;在所有试验方案中,24 h降水率最大差异为52.5 mm/h.云微物理过程对云和降水特征的具体影响表现在:(1)如果不考虑雨水蒸发冷却效应,此时台风内核上升运动强度最强(达到-19 Pa/s),雨水和霰粒子增长最明显,相对于对照试验增量分别为1.8和2.5 g/kg.(2)霰和雪的融化对于螺旋雨带中雨滴的增长十分重要,但他们可能不是云墙中雨水形成的主导因子.(3)不同方案的降水模拟特征也存在较大差别,采用暖云参数化后,降水区域最小,但其中对流降水比例最大(63.19%);霰落速减半后,降水区域最大,其中非对流降水比例也最大(51.15%).     

The Vertical Structures of Atmospheric Temperature Anomalies Associated with El Nio Simulated by the LASG/IAP AGCM： Sensitivity to Convection Schemes

The vertical structures of atmospheric temperature anomalies associated with El Nio are simulated with a spectrum atmospheric general circulation model developed by LASG/IAP (SAMIL). Sensitivity of the model’s response to convection scheme is discussed. Two convection schemes, i.e., the revised Zhang and Macfarlane (RZM) and Tiedtke (TDK) convection schemes, are employed in two sets of AMIP-type (Atmospheric Model Intercomparison Project) SAMIL simulations, respectively. Despite some deficiencies in the up...     

对流有效位能计算的新方案

对流有效位能是强对流天气分析预报的重要参数。本文通过理论推导,提出了载水气块和非载水气块两种情况下,对流有效位能的两个新的计算方案,分别记作CAPEw和CAPE。该方案与以往的方案有很大的不同和根本的区别,更加合理和符合实际。编制了相应的计算程序,可以计算两种情况下的对流抑制能量CINw、CIN,对流有效位能CAPEw、CAPE,抬升凝结高度LCL(二者相同),自由对流高度LFCw、LFC,对流凝结高度CCL(二者相同)以及平衡高度ELw、EL;可以预报对流温度Tg(二者相同);可以分析地面和高空温度、湿度、等压面高度发生变化时,CINw、CIN和CAPEw、CAPE等的数值变化情况,便于业务应用和理论研究。还讨论了影响对流有效位能局地变化的因子和预报思路。     

Performance assessment of a GCM land surface scheme using a fine‐scale calibrated hydrological model: an evaluation of MOSES for the Nile Basin

Land surface schemes (LSSs) represent the interface between land surface and the atmosphere in general circulation models (GCMs). Errors in LSS‐simulated heat and moisture fluxes can result from inadequate representation of hydrological features and the derivation of effective surface parameters for large heterogeneous GCM gridboxes from small‐scale observations. Previous assessments of LSS performance have generally compared simulated heat and moisture fluxes to observations over a defined experimental domain for a limited period. A different approach has been evaluated in this study, which uses a fine‐resolution calibrated hydrological model of the study basin to provide a quasi‐observed runoff series for direct comparison with simulated runoff from a selected LSS at GCM scale. The approach is tested on two GCM gridboxes covering two contrasting regions within the Nile Basin. Performance is mixed; output from the LSS is generally compatible with that of the fine‐resolution model for one gridbox while it cannot reproduce the runoff dynamics for the other. The results also demonstrate the high sensitivity of runoff and evapotranspiration to radiation and precipitation inputs and show the importance of subtle issues such as temporal disaggregation of climatic inputs. We conclude that the use of a fine‐resolution calibrated model to evaluate a LSS has several advantages, can be generalized to other areas to improve the performance of global models and provides useful data that can be used to constrain LSS parameterizations. Copyright © 2009 John Wiley & Sons, Ltd.     

A new computational strategy for solving two‐phase flow in strongly heterogeneous poroelastic media of evolving scales

We develop a new computational methodology for solving two‐phase flow in highly heterogeneous porous media incorporating geomechanical coupling subject to uncertainty in the poromechanical parameters. Within the framework of a staggered‐in‐time coupling algorithm, the numerical method proposed herein relies on a Petrov–Galerkin postprocessing approach projected on the Raviart–Thomas space to compute the Darcy velocity of the mixture in conjunction with a locally conservative higher order finite volume discretization of the nonlinear transport equation for the saturation and an operator splitting procedure based on the difference in the time‐scales of transport and geomechanics to compute the effects of transient porosity upon saturation. Notable features of the numerical modeling proposed herein are the local conservation properties inherited by the discrete fluxes that are crucial to correctly capture the fingering patterns arising from the interaction between heterogeneity and nonlinear viscous coupling. Water flooding in a poroelastic formation subject to an overburden is simulated with the geology characterized by multiscale self‐similar permeability and Young modulus random fields with power‐law covariance structure. Statistical moments of the poromechanical unknowns are computed within the framework of a high‐resolution Monte Carlo method. Numerical results illustrate the necessity of adopting locally conservative schemes to obtain reliable predictions of secondary recovery and finger growth in strongly heterogeneous deformable reservoirs. Copyright © 2011 John Wiley & Sons, Ltd.