不同模式扰动方案在风暴尺度集合预报中的对比试验研究

由于适用于中长期集合预报的模式扰动技术在风暴尺度集合预报系统中的影响并不明确,为探究不同模式扰动方案在风暴尺度集合预报中的效果,基于WRF模式设计了3组模式扰动方案:多物理扰动（MP）方案、随机物理倾向扰动（Stochastically Perturbed Parameterization Tendencies,SPPT）方案以及由MP方案与SPPT方案组合构建的一种新混合扰动（SPMP）方案。对2013年7月5—6日发生在江淮流域的一次强对流天气过程进行了数值模拟。结果表明:MP方案在积分前期的降水概率评分较高,对高层大气的扰动效果更为合理;SPPT方案主要作用于积分中后期,对大气低层及近地面的扰动效果最为理想,尤其是对于地面水汽场的模拟;SPMP方案能显著提高大气中高层各预报变量的离散度,降低均方根误差,提升集合成员的可信度,有效弥补降水预报评分在单独使用MP方案和SPPT方案不同积分时段的不足。在扰动水平传播方向上,SPMP方案的扰动形态主要受MP方案主导;垂直方向上,SPMP方案在低层的扰动形态与SPPT方案一致,在高层受MP方案控制。波谱能量分析表明3组方案的扰动能量随积分时间均有向大尺度传播的趋势,SPMP方案能有效补偿两种方案能量在各尺度的耗散。      

大气粒子散射相函数的参数化方案比较及其改进

大气粒子散射相函数的参数化是大气辐射传输参数化的重要组成部分。文中全面比较了大气粒子的HenyeyGreenstein(HG)方案和双Henyey-Greenstein(DHG)方案,并在四流球谐函数展开累加法中,应用这两种相函数参数化方案计算气溶胶、云、霾粒子的反射率、透射率或吸收率。该研究结果表明:HG方案无法表现相函数的后向峰值,因而其计算的大气粒子反射率和透射率精度较差;DHG方案能较好地表征相函数的整体特征,但是该方案计算的相函数易出现后向异常峰值或为负值,并导致计算得到的气溶胶、云、霾粒子的反射率和透射率精度甚至会低于HG方案。对DHG方案进行进一步研究,提出了改进的DHG方案(MDHG)。MDHG方案计算结果稳定,并能很好表征相函数的前向和后向峰值的特征,其计算的大气粒子的反射率和透射率精度也较高。因此,MDHG方案是一种理想的相函数参数化方案。     

复杂地形上气象场对空气质量数值模拟 结果影响的研究

利用WRF模式三种边界层参数化方案(YSU、MYJ、ACM2)产生的气象场分别驱动多尺度空气质量模式CMAQ,对兰州市西固区冬季2005年1月27日至2月2日期间SO₂和NO₂浓度进行了数值模拟,将模拟结果与同期监测的污染物浓度进行对比分析,结果表明:WRF模式不同边界层参数化方案模拟输出的气象场驱动CMAQ模式所模拟的SO₂和NO₂浓度均可以反映出污染物的时空变化特征,CMAQ模式具有模拟复杂下垫面高分辨率污染物输送特征的能力;WRF模式的边界层参数化方案选为局地与非局地闭合方案(ACM2)时,模拟的气象场驱动CMAQ模式得到的空气污染物浓度分布特征最优,这主要是由于ACM2的湍流输送机制较为合理,模拟的边界层低层气象场更接近实际,从而可以较好地模拟污染物的输送特征;当CMAQ模式的垂直混合方案与WRF模式的湍流输送方案一致时(均采用ACM2方案),模式间的兼容性好.     

Parameterization of bedform morphology and defect density with fingerprint analysis techniques

A novel method for parameterizing the morphology of seafloor ripples with fingerprint analysis numerical techniques is presented. This fully automated analysis tool identifies rippled areas in two-dimensional imagery of the seafloor, and returns ripple orientation and wavelength as well as a new morphological parameter, the spatial density of ripple defects. In contrast to widely used manual and spectral parameterization methods, this new technique yields a unique probability distribution for each derived parameter, which describes its spatial variability across the sampled domain. Here we apply this new analysis technique to synthetic and field collected side-scan sonar seafloor images in order to assess the methods capacity to define bed geometry across a wide range of simulated and observed morphological conditions. The resulting orientation and wavelength values compare favorably with those of the existing manual and spectral parameterization methods, and are superior under environmental conditions characterized by low signal to noise ratios as well as high planform ripple sinuosity. Furthermore, the resulting ripple defect density values demonstrate correlation with ripple orientation, wave direction, and the Shields parameter, which is consistent with recent investigations that have theoretically linked this parameter to hydrodynamic forcing conditions. The presented fingerprint analysis method surpasses the capacity of existing methods for ripple parameterization and promises to yield greater insight into theoretical and applied problems associated with the temporal and spatial variability of ripple morphology across a wide spectrum of marine environments.     

湖泊数值模拟研究现状综述

作为陆面过程的重要组成部分,湖泊在天气气候预测中的作用得到了相关研究者的广泛关注,并成为大气科学研究领域中的一个热点.主要综述了当前湖泊模式的基本类型、所考虑的关键物理过程及其参数化方案,并尝试分析了各个方案的优劣以及模拟效果.相关研究表明,现有的湖泊模式对浅湖的模拟比较成熟,而对深湖和冰湖的模拟有待改进.未来的研究亟需构筑适用于大湖、深湖以及冰湖的参数化方案,发展能全面模拟各类型湖泊的数值模式,并且湖泊数值模拟的改进依赖于今后更多高质量的全球湖泊观测结果.     

Parameterized expressions for an improved Rouse equation

The equilibrium profile of suspended sediment concentration （SSC） in steady, uniform turbulent flow has been extensively studied since the proposition of Rouse equation （Rouse, 1937）. Many researchers keep studying and attempting to improve the Rouse equation by modifying the mixing length concept. Bose and Dey （2009） improved Rouse equation based on the modified sediment diffusivity, however, their proposed sediment diffusivity has two parameters, the depth modification factor α and the reciprocal of the sediment Schmidt number β = 1/S, and the corresponding expressions were not given properly in their paper. This paper aims to parameterize the relevant coefficients, the results shows： （1） α is parameterized from a reasonable approximation by the analysis of surface boundary condition of SSC; （2） considering some effect factors on the settling velocity, β can be expressed by an extended and modified Montes-Ippen formula, and its related parameters are calibrated by substantial data about the relationship between the primary Rouse parameter Z. = ω/ku and its counterpart Z＇ = ω/βku.. Through verification with some experimental and field datasets, and comparison with related formula, the results show that the sediment diffusivity and SSC profiles have a good agreement with measured data used the present parameterized expressions of a and,8, so the present improved Rouse equation is reasonable, and can be applied in practical applications. Finally, an approximation and simplified expression of the improved Rouse equation is deduced by using the approach of perturbation, which can be applied easily for the estimation of suspended sediment transport rate.     

Sensitivity of simulated wintertime Arctic atmosphere to vertical resolution in the ARPEGE/IFS model

The current state-of-the-art general circulation models, including several of those used by the IPCC, show considerable biases in the simulated present day high-latitude climate compared to observations and reanalysis data. These biases are most pronounced during the winter season. We here employ ideal vertical profiles of temperature and wind from turbulence-resolving simulations to perform a priori studies of the first-order eddy-viscosity closure scheme employed in the ARPEGE/IFS model. This reveals that the coarse vertical resolution (31 layers) of the model cannot be expected to realistically resolve the Arctic stable boundary layer. The curvature of the Arctic inversion and thus also the vertical turbulent-exchange processes cannot be reproduced by the coarse vertical mesh employed. To investigate how turbulent vertical exchange processes in the Arctic boundary layer are represented by the model parameterization, a simulation with high vertical resolution (90 layers in total) in the lower troposphere is performed. Results from the model simulations are validated against data from the ERA-40 reanalysis. The dependence of the surface air temperature on surface winds, surface energy fluxes, free atmosphere stability and boundary layer height is investigated. The coarse-resolution run reveals considerable biases in these parameters, and in their physical relations to surface air temperature. In the simulation with fine vertical resolution, these biases are clearly reduced. The physical relation between governing parameters for the vertical turbulent-exchange processes improves in comparison with ERA-40 data.     

Numerical Simulation of Effect of Sea Spray on Typhoon ＇Morakot＇ Structure

将海洋飞沫参数化引入到高分辨率、非静力中尺度模式中,并对0908号台风＂Morakot＂进行了数值模拟,研究了海洋飞沫对台风＂Morakot＂结构和强度的影响。结果表明：（1）不论是否考虑海洋飞沫作用,模式均能较好地模拟出台风＂Morakot＂的移动路径,说明海洋飞沫对其移动路径影响不大;（2）引入海洋飞沫参数化后,台风眼墙区域的切向风速、径向风速、垂直速度、涡度、云水混合比、雨水混合比等物理量均增强,表明飞沫对台风结构变化的影响明显;（3）海洋飞沫对台风＂Morakot＂演变的直接影响是在对流层低层,低层风速明显增大,大风速区的影响尤为显著;（4）飞沫的蒸发使台风范围内的潜热和感热通量明显增强,尤其是潜热通量,其大值区对应着台风中心附近的最大风速区。由于水汽和热量输送的增强,使台风眼壁附近的云水量与雨水量增多,因此降水强度明显增加。     

海洋飞沫对台风“Morakot”结构影响的数值模拟研究

将海洋飞沫参数化引入到高分辨率、非静力中尺度模式中,并对0908号台风"Morakot"进行了数值模拟,研究了海洋飞沫对台风"Morakot"结构和强度的影响。结果表明:(1)不论是否考虑海洋飞沫作用,模式均能较好地模拟出台风"Morakot"的移动路径,说明海洋飞沫对其移动路径影响不大;(2)引入海洋飞沫参数化后,台风眼墙区域的切向风速、径向风速、垂直速度、涡度、云水混合比、雨水混合比等物理量均增强,表明飞沫对台风结构变化的影响明显;(3)海洋飞沫对台风"Morakot"演变的直接影响是在对流层低层,低层风速明显增大,大风速区的影响尤为显著;(4)飞沫的蒸发使台风范围内的潜热和感热通量明显增强,尤其是潜热通量,其大值区对应着台风中心附近的最大风速区。由于水汽和热量输送的增强,使台风眼壁附近的云水量与雨水量增多,因此降水强度明显增加。