“00.7”北京特大暴雨模拟中气象资料同化作用的评估

针对2000年7月4～5日北京地区的一次特大暴雨过程(24 h降水量达240 mm),文中利用MM5/WRF三维变分系统和MM5非静力模式,对此次特大暴雨过程中的各种气象监测资料(地基GPS大气柱水汽含量、常规探空、高空测风、地面常规观测和地面自动气象站)的同化作用通过观测系统数值试验进行了评估.结果表明与传统的客观分析方案相比较,MM5/WRF三维变分同化系统可直接引入非常规地基GPS大气柱水汽含量监测资料,提供更好的大气初始分析场.在三维变分同化方案下,各种大气监测资料均对改进此次特大暴雨模拟有不同程度的贡献,其中,常规探空和高空测风监测资料对改进预报结果的影响最大,地面常规观测和地面自动气象站观测资料作用次之,地基GPS大气柱水汽含量资料在与其他大气监测资料相互优势互补后,可很好地改善模式大气的分析质量,通过三维变分同化技术在区域数值天气预报模式初始场中引入地基GPS大气水汽监测网资料,使此次强降水个例的6 h和24 h测站降水预报的TS评分值在1,5,10和20 mm预报检验阈值下分别提高了1%～8%.研究结果对利用三维变分数值系统,评估气象监测网资料在改进高影响天气事件预报中的作用有借鉴意义.     

IMPACT OF DIABATIC PROCESSES IN AGCM ON 4-DIMENSIONAL VARIATIONAL DATA ASSIMILATION*

The impact of diabatic processes on 4-dimensional variational data assimilation (4D-Var) was studied using the 1995 version of NCEP's global spectral model with and without full physics.The adjoint was coded manually.A cost function measuring spectral errors of 6-hour forecasts to "observation" (the NCEP reanalysis data) was minimized using the L-BFGS (the limited memory quasi-Newton algorithm developed by Broyden,Fletcher,Goldfard and Shanno) for optimizing parameters and initial conditions.Minimization of the cost function constrained by an adiabatic version of the NCEP global model converged to a minimum with a significant amount of decrease in the value of the cost function.Minimization of the cost function using the diabatic model, however,failed after a few iterations due to discontinuities introduced by physical parameterizations.Examination of the convergence of the cost function in different spectral domains reveals that the large-scale flow is adjusted during the first 10 iterations,in which discontinuous diabatic parameterizations play very little role.The adjustment produced by the minimization gradually moves to relatively smaller scales between 10-20th iterations.During this transition period,discontinuities in the cost function produced by "on-off" switches in the physical parameterizations caused the cost function to stay in a shallow local minimum instead of continuously decreasing toward a deeper minimum. Next,a mixed 4D-Var scheme is tested in which large-scale flows are first adiabatically adjusted to a sufficient level,followed by a diabatic adjustment introduced after 10 to 20 iterations. The mixed 4D-Var produced a closer fit of analysis to observations,with 38% and 41% more decrease in the values of the cost function and the norm of gradient,respectively,than the standard diabatic 4D-Var,while the CPU time is reduced by 21%.The resulting optimal initial conditions improve the short-range forecast skills of 48-hour statistics.The detrimental effect of parameterization discontinuities on minimization was also reduced.     

On ensemble representation of the observation-error covariance in the Ensemble Kalman Filter

Evensen (2003) presents a modification of the Ensemble Kalman Filter (EnKF), in which the observation-error and background-error covariance matrices are both represented by ensembles, in contrast to the usual practice, where only the background error is so represented. It is shown that this modification can cause the ensemble to collapse to a single member, in the common situation where the number of observations is more than twice the number of ensemble members, and to be rank-deficient when the number of observations is greater than or equal to the ensemble size. It is also shown that some further modifications to the scheme, presented by Evensen as offering numerical efficiencies, can prevent this collapse. However, these latter modifications are shown in some simple numerical examples to require tuning to produce acceptable results, which are nevertheless inferior to those of the standard EnKF.Acknowledgements The author acknowledges useful discussions with Peter Steinle, and other participants at the EnKF workshop held in BMRC in November, 2003.     

雷达径向风速同化对台风麦莎模拟的影响

本文应用WRF-3DVAR系统同化多普勒雷达径向风速资料,并利用WRF模式对台风麦莎（2005年8月6日00时至7日00时）进行数值模拟,以此检验多普勒雷达径向风速资料在改进模式初始场及提高台风路径和降水预报准确度等方面的应用效果及意义,探讨不同同化时间间隔对同化效果的影响。通过对比同化试验和控制试验发现：雷达径向风速资料能通过调整初始风场结构,改进对台风结构的模拟,从而提高模式对台风麦莎的路径、强度和降水的模拟能力;雷达资料信号在进入模式后有一定的时效性,在一定程度上缩小同化时间间隔有助于提高同化效果。     

在NCEPGDAS中同化MSG和GOES资料

首次将MSG-2（Meteosat Second Generation-2）卫星上的旋转增强可见光及红外成像仪（Spin—ning Enhanced Visibleand Infrared Imager,SEVIRI）的观测资料同化到美国国家环境预报中心（Na—tional Centers for Environmental Prediction,NCEP）全球资料同化系统（globaldataassimilationsys—tern,GDAS）中。对当前的地球静止业务环境卫星（Geostationary Operational Environmental Satel—lite,GOES）成像仪资料的同化问题也进行了进一步探讨。利用CRTM（The Community Radiative Transfer Model）模式,对SEVIRI辐射率观测资料进行了模拟。为了对红外辐射率资料进行模拟,CRTM模式中的几个关键部分得到改进,例如：动态更新地面发射率资料以及采用了快速精确的气体吸收模块。为了改进对SEVIRI和GOES成像仪辐射率资料的模拟效果,采用了GSICS（The Global Space—Based Inter—Calibration System）标定订正。初步研究结果表明,包含对SEVIRI辐射率资料的水汽通道（6．25和7．35μm）和二氧化碳通道（13．40μm）的同化对GFS（Global Forecast System）6d预报具有显著的正影响;而对其他5个SEVIRI红外窗口通道资料的同化则减小了这种正影响。通过应用GSICS标定算法,订正了SEVIRI和GOES-12成像仪观测资料的偏差,提高了对GFS预报的影响。此外,还需作进一步研究来提高对SEVIRI红外窗口通道辐射率资料同化的有效性。     

LAPS分析场在一次强对流天气过程尺度分析中的应用

利用LAPS（local analysis and prediction system）同化系统融合多种观测资料,对2010年8月25日发生在上海的一次强对流天气过程进行中尺度分析。结果表明,这次强对流天气主要是由中、低空中尺度辐合系统直接触发形成的。强对流形成阶段,地面有分散辐合形成并逐渐加强,成为触发中尺度垂直环流的主要机制,垂直结构上出现低层辐合高层辐散的有利配置,风暴中心附近出现明显的上升气流区,中高层相对湿度显著增加。成熟阶段,强对流云体中心附近的对流层底层开始出现下沉气流,上升气流在其拖曳作用下明显倾斜。衰减阶段,下沉气流加强使中尺度环流动力结构和水汽供应受到破坏,垂直结构上转为底层辐散高层辐合。因此,与天气尺度分析相比,基于LAPS的中尺度分析能更深刻地揭示中小尺度系统的三维结构和时空演变特征。     

散射计和下投式探空仪资料对台风模拟的影响

为了评估同化QuikSCAT海平面气压场资料和Dropsonde资料对台风"Bilis"和"Kat-rina"模拟的影响,通过WRF-3DVAR变分同化系统将QuikSCAT海平面气压场资料和Dropsonde资料同化到WRF模式中,并以未同化任何资料和同化了QuikSCAT海平面气压场以及同化了Drop-sonde资料的数据为初始场,应用WRF模式进行模拟预报,并对得到的初始分析场和模拟预报得到的结果与台风实况报文进行了对比,得到的结论是:对台风初始时刻海平面气压场分布、初始台风中心定位、台风初始暖心结构的改善、台风路径、海平面最低气压场、海面最大风速和雷达回波的模拟方面,同化QuikSCAT海平面气压场的效果比同化Dropsonde资料明显,其中二者对海平面最低气压场的改进效果比海面最大风速的要好。但是两种资料的同化对初始相对湿度和台风强度细节变化方面都没有产生明显改善作用。同化QuikSCAT海平面气压场试验虽然仅仅同化了单层海平面气压场,但是在三维变分同化系统的动力约束作用下,也可以影响到中高层的各个要素场。     

A comparison of shelf observation platforms for assimilation in an eddy-resolving ocean model

An assessment of the likely benefits of assimilating in situ temperature (T) and salinity (S) observations from repeat glider transects and surface velocity observations from high-frequency radar arrays into an eddy-resolving ocean model is presented. The deployment of new shelf observation platforms around Australia is being undertaken through the Australian Integrated Marine Observing System program. In this study, various options for an observing system along the coast of New South Wales, Australia, are assessed for their benefits to an ocean forecast and reanalysis system. The forecast system considered here uses ensemble optimal interpolation (EnOI) for data assimilation. Using error estimates from the EnOI scheme, estimates of the theoretical analysis errors are calculated for different observing systems that include a range of remotely sensed and in situ observations. The results demonstrate that if HF radar observations are assimilated along with the standard components of the global ocean observing system, the analysis errors are likely to reduce by as much as 80% for velocity and 60% for T, S and sea-level in the vicinity of the observations. Owing to the relatively short along-shore decorrelation length-scales for T and S near the shelf, the glider observations are likely to provide the forecast system with a more modest gain.     

MODIS反照率产品在模拟黄河源区陆面过程和降水中的应用

地表反照率是陆面过程中一个重要的物理量, 其变化直接影响地表能量的收支状况, 进而可以影响气温和降水等其它气象要素。本文利用WRF (Weather Research and Forecasting) 模式, 通过两组数值模拟试验分别探讨了地表反照率改变在黄河源区不同下垫面情况下潜热、 感热的分配关系, 详细分析了地表反照率改变对降水变化的影响机制, 最后应用EOS/MODIS地表反照率产品替代原模式低时空分辨率的地表反照率。研究结果表明: (1)当地表反照率减少(增加)时, 模拟的区域平均地表温度、感热、潜热数值相应增大(减少)。当地表反照率减少0.1时, 地表温度上升约1.0 K, 感热和潜热量增量比约为3∶1。 (2) 地表反照率改变对降水量变化影响最大的区域是黄河源区下游的草场区域, 其次是黄河源头区域, 最小的是黄河源区北部的稀疏植被区域。地表反照率通过对大气动力、 热力以及水汽条件的影响, 使得降水发生的环境改变, 主要体现在: 当地表反照率减少时, 地表气压的减少使得大气低层的辐合气流增强, 有利于上升运动的发生; 2.0 m气温的升高增强了大气近地层的热力不稳定度; 2.0 m比湿的增加表明近地层空气水汽含量增加。 (3) 与实况对比分析发现, 使用卫星遥感产品后在月尺度上能够更准确地模拟降水量的变化过程。