数值模式误差对降水四维变分资料同化及预报的影响

利用一个无量纲的水汽发展方程,针对同化时间窗口内出现和不出现降水两种情况,分析了不同模式误差和初始误差对降水四维变分资料同化预报效果的影响。结果表明,应用四维变分资料同化方法进行降水预报前,应该充分考虑数值模式中的误差,才能得到比较满意的同化及预报结果。假定同化窗口内获得的比湿观测场是准确的,当不存在模式误差时,四维变分资料同化方法可以有效地消除初始场误差,找到比湿真实初始场;而存在模式误差时,四维变分资料同化后的模式初始场会偏离真实的比湿初始场,并且模式误差越大,偏离程度越严重。在一些模式误差情况下,由于模式误差在同化窗口及延伸预报时段的作用不同,进行四维变分资料同化处理后,尽管累积降水量的预报结果在同化时间窗口内优于同化前的预报,而在最终预报时刻反而差于同化前。     

变分四维同化方法若干问题的数值试验

资料误差和模式误差都可能影响变分四维同化的结果,针对这一问题利用浅水方程模式进行了变分四维同化的数值模拟试验。试验分三种情况进行：（１）仅修正初始场,（２）仅修正模式;（３）二者同时修正。试验结果表明,当模式有误差时,传统的变分四维同化方法（仅修正初始场）可能将模式的误差混淆到初始场中去,尽管在同化期间可得到较好的拟合,但却不一定能保证同化后有好的预报。如果不修正初始场而修正模式,当模式误差不大时     

四维变分同化方法在暴雨预报中的应用

本文利用PSU/NCAR的MM5数值预报模式及其伴随模式,以我国1999年6月23日~6月24日的一次梅雨锋暴雨过程为个例,作了两组试验:控制试验和同化试验,并对两组试验的降水预报效果以及初始场进行对比分析,结果表明:四维变分资料同化方法可以将各种不同类型、不同时次的观测资料同化到模式中,将这些资料中有用的中尺度信息引入到模式初始场,有效改善初始场,从而提高暴雨预报水平。     

REM模式伴随系统的建立及其四维变分资料同化初步试验

Regional-Eta-Coordinate-Model(REM)中尺度模式对中国区域性降水显示出公认的较高预报能力,建立其四维变分资料同化系统是完善该模式,进一步提高其预报效果的重要工作。本研究编写了REM模式的切线性模式和伴随模式,介绍了建立REM模式伴随系统的过程,并利用实际天气个例资料,分别对REM模式的切线性模式、伴随模式及定义的目标函数梯度进行了正确性检验,检验结果表明对REM模式的切线性模式及伴随模式编写是成功的。利用REM模式的伴随系统,对1998年06月08日00时到09日00时和2000年08月01日00时到02日00时两个实际天气个例进行了四维变分资料同化试验。从数值试验的结果分析可以看到,进行四维变分资料同化后,两个天气个例在预报结束时刻其预报结果对风场和湿度场的预报都有明显改善,对温度场和高度场的预报也有所改善。对于累积降水的预报,两个个例利用四维变分资料同化后得到的初始场进行的预报结果则有较大不同,在个例1中,变分同化后对降水中心的位置和降水强度的预报都有明显改善,预报结果更接近于观测场;个例2中,变分同化后对降水中心位置和强度的预报则没有改善,产生这种现象的原因可能是由于定义的目标函数中没有加进背景场项,也可能是由于采用的观测资料时次比较少,还需要进一步进行研究和试验。     

中尺度WRF数值模式系统本地化业务试验

利用中尺度WRF数值模式及WRF三维变分同化系统,在对比试验的基础上,选取了适合本地的积云过程、微物理过程和辐射过程的方案组合;选择了NCEP/GFS作为模式的背景场;统计计算了以云南为中心的区域背景误差协方差并替换了三维变分同化系统中原有的背景误差协方差;同时,考虑模式底层高度与地面观测站高度的差异,进行了地面资料地形订正.通过上述试验研究,建立了本地化的中尺度WRF数值预报业务系统,该系统能较好地刻画本地下垫面的动力和热力状况,预报能力有明显改善.     

"7·8"特大暴雨过程雷达风场资料的三维变分同化研究

基于中尺度暴雨预报模式AREM和三维变分同化系统GRAPeS-3dvar对2003年7月8日发生在长江中游的一场特大暴雨过程,进行了武汉、宜昌多普勒雷达风场资料的三维变分同化研究,结果表明:(1)增加了武汉、宜昌多普勒雷达反演风场资料后,模式对湘-鄂交界处特大暴雨区的模拟效果改善非常显著.(2)分别单独增加武汉、宜昌多普勒雷达风场资料的三维变分同化后,降水模拟都有明显改进,但要差于对两部雷达风场资料同时进行变分同化的结果,表明同化的中尺度初始风场信息越多,初始场的质量越高,降水模拟效果越好.(3)雷达风场资料的三维变分同化,改善了分析场中尺度结构的描述,从而减轻了spin-up现象,使得模式在积分初期就能模拟出与实况相近的强降水.     

单多普勒雷达资料在伴随模式同化系统中的应用研究

利用变分方法反演单多普勒雷达资料,得到风矢量场。同时,利用MM5伴随模式同化系统,结合一次暴雨过程,设计四种方案,进行数值模拟试验。结果表明,通过变分方法反演的雷达资料的应用对于暴雨的分布预报有明显的改善作用;运用伴随方法同化雷达资料后可改善对暴雨中心的预报;对于各个物理量误差的减少,雷达资料的应用也起到重要作用,尤其是对于风场作用较为明显;雷达资料的应用可加快伴随模式同化系统目标函数的收敛,得到最优初始场。     

多普勒雷达资料同化的稀疏化方式对暴雨过程的影响研究

利用中尺度非静力WRF(Weather Research and Forecasting)模式及其三维变分同化系统,对2007年7月淮河流域的一次强降雨过程进行多普勒雷达径向速度资料的三维变分同化试验,重点考察雷达资料的不同稀疏化方式对同化结果以及对暴雨数值模拟的影响。结果表明:同化多普勒雷达径向速度资料使得模式初始风场包含了更丰富的中尺度特征信息,有效调整了初始场的环流结构,能够改善模式对暴雨过程的模拟效果;以不同的稀疏化处理方式同化多普勒雷达径向速度资料对分析场会产生不同的影响,进而影响模式的降水预报效果,本次试验中当极坐标网格径向分辨率取10 km的时候降水过程的预报效果最好。     

云迹风在热带气旋路径数值预报中的应用研究

通过一系列四维变分同化试验对GMS5卫星资料反演的云迹风资料在西北太平洋热带气旋的初始化及路径数值预报中的作用进行研究,同化资料为中国国家卫星气象中心提供的GMS5水汽和红外云迹风资料,其中70%在400hPa以上,50%集中在200～300hPa。应用美国NCAR/PSU中尺度模式MM5及其四维变分同化系统,同化窗口为6h,对初始时刻和6h后的云迹风进行同化。同化前对云迹风资料进行了简单的类似ECMWF初值检验方法的质量控制。对2002年8个西北太平洋热带气旋共进行了22组试验。结果表明,采用四维变分同化技术同化云迹风对热带气旋路径预报有一定改善,12,24,36和48h预报的平均距离误差分别降低5%,12%,10%和7%,但同化云迹风的作用与初始气旋强度有关。选择初始中心海平面气压960hPa作为强、弱气旋的分类标准,则11个较强气旋平均路径误差12h减小了13%,12h以后的预报误差减小率维持在20%以上。而对于11个较弱气旋,平均路径误差反而略有增加,说明同化云迹风资料对不同初始强度的气旋作用也有所不同。其主要原因是由于强度较强的热带气旋往往具有较为深厚的垂直结构,因此受高层大气流场的影响更明显;同时,较弱热带气旋的云迹风观测相对稀少且凌乱,并且更容易受环境气流的影响,因此对于较弱的热带气旋,当模式变量与模式或变量之间在同化后不够协调的话,就会产生负效应。     

多普勒雷达径向速度同化在淮河暴雨数值模拟中的应用

针对2007年7月淮河流域的一次强降雨过程,利用WRF中尺度数值模式及其三维变分同化系统(WRF-3DVAR),开展了多普勒雷达径向速度的三维变分同化对暴雨过程模拟效果的影响研究。结果表明:WRF-3DVAR能够有效地同化多普勒雷达径向速度资料,同化后使得模式初始场出现了一定的调整,包含更详尽的中尺度特征信息,进而显著改善模式对大暴雨过程前12h降水的模拟效果。在高分辨率中尺度数值模式中有效地利用多普勒天气雷达资料,能较好地提高中尺度降雨预报。     

The definition of mesoscale selective forecast error covariances for a limited area variational analysis

Summary ?The paper deals with an alternative formulation of the so-called NMC (National Meteorological Center, now National Centers for Environmental Prediction) statistics to compute the background error covariance matrix to be used in a mesoscale variational analysis. While the standard method uses differences of forecasts valid for the same time, but starting from different analysis times, the new formulation required the recomputation of the short-term forecast with the initial and lateral boundary data that come from the long-term run. In the frame of a limited-area model, this approach forces the error variances at large scales to decrease drastically, because those scales are controlled by the (constant data) lateral boundary coupling. As a result, the background cost function acts more scale selectively, with an emphasis on medium scales. The analysis increments obtained from the 3D-VAR system show that the analysis increments are sharper and more concentrated with the new formulation, both in single observation and in full observation experiments. This work is part of a wider project for building a variational assimilation system inside the ALADIN model. The complete system should concentrate on mesoscale features and it should not reanalyse those scales that were already treated by the global model (ARPEGE). Some difficulties and perspectives are drawn in the concluding discussion. Received February 12, 2001; revised July 24, 2001     

FY-2C云迹风资料同化应用对台风预报的影响试验研究

针对0505号台风“海棠”, 采用WRF区域中尺度模式进行控制试验和两个同化试验, 利用WRF-3DVAR同化系统同化FY-2C红外和水汽两个通道云迹风反演产品, 同化分云迹风经质量控制和未经质量控制两组同化试验。通过三组试验分析云迹风资料对降水和风场等的预报结果的影响, 并进行24小时降水量分级Ts评分检验以及风场点对点检验。结果表明: 同化经质量控制云迹风资料可以提高降水落区和强度预报的准确度, 不同等级的Ts评分较其它试验都有较明显改进; 风场预报模拟也有所改善。增加两例台风, 使用与“海棠” 相似的处理方法进行模拟试验, 并对模拟结果24小时降水分析与检验, 得到与“海棠”类似结论。因此, 经过合理性选择的云迹风资料的加入, 有利于补充初始场中可能未包含的中尺度信息, 从而提高试验中对于降水、风场等的模拟效果, 提高WRF模式的模拟预报能力。     

多普勒雷达资料同化对江苏一次飑线过程的数值模拟

应用新一代中尺度预报模式WRF模式及其3DVar同化系统, 针对江苏地区2009年6月14日飑线过程进行了多普勒雷达资料的同化试验研究, 在对雷达资料进行严格质量控制的基础上, 设计一系列尺度化因子优化调整及同化频率的敏感性试验。试验结果表明:同化后初始场得到不同程度改善, 适当的尺度化因子设定, 能够有效改进对模式初始场中700 hPa风场和850 hPa温度场以及组合反射率因子等要素的分析, 进而改善短时降水预报和风暴的垂直结构配置;并且同化频率越高, 对初始场的组合反射率因子分布与观测更为接近, 短时降水预报越准确。     

A COMPARISON OF THE BLENDING AND CONSTRAINING METHODS TO INTRODUCE LARGE-SCALE INFORMATION INTO GRAPES MESOSCALE ANALYSIS

To solve the problem of mesoscale analysis error accumulation after a period of continuous cycle data assimilation (CCDA), a blending method and a constraining method are compared to introduce global analysis information into the Global/Regional Assimilation and Prediction Enhanced System mesoscale three-dimensional variational data assimilation system (GRAPES-Meso 3Dvar). Based on a spatial filter used to obtain a blended analysis, the blending method is weighted toward the T639 global analysis for scales larger than the cutoff wavelength of 1,200 km and toward the GRAPES mesoscale analysis for wavelengths below that. The constraining method considers the T639 global analysis data as an extra source of information to be added in the 3DVar cost function. The cloud-resolving GRAPES-Meso system (3 km resolution) with a 3 h analysis cycle update is chosen, and forecast experiments on an extreme precipitation event over the eastern part of China are presented. The comparison shows that the inclusion of large-scale information with both methods has a positive impact on the regional model, in which the 3 h background forecasts are slightly closer to the radiosonde observations. The results also show that both methods are effective in improving large-scale analysis while reserving the well-featured mesoscale information, leading to an enhancement in the balance and accuracy of the analysis. Subjective verification reveals that the introduction of large-scale information has a visible beneficial impact on the forecast of precipitation location and intensity. The methodologies and experiences presented in this paper could serve as a reference for ongoing efforts toward the development of multi-scale analysis in GRAPES-Meso.     

APPLICATION OF BDA SCHEME IN TYPHOON TRACK PREDICTION

The MM5, which is the PSU/NCAR mesoscale nonhydrostatic limited-area model, and its adjoining modeling system are used in this paper. Taking T106 analysis data as background field the authors generate an optimal initial condition of a typhoon by using two bogus data assimilation schemes, and conduct some numerical simulating experiments. The results of No.9608 typhoon (Gloria) show that the optimal initial field have some dramatic improvements, such as inaccurate position of typhoon center, weaker typhoon circulation and incomplete inner structure of the typhoon, which are caused by shortage of data over the sea. Some improvements have been made in the track forecast. Through several comparing experiments, the initialization optimized by BDA scheme is found to be more reasonable than GFDL scheme and its typhoon track forecast is better.     

华北地区一次大雾过程的三维变分同化试验

文中采用WRF非静力数值预报模式及其三维变分同化系统(WRF3D-Var),对2006年1月13—14日发生在华北地区及山东半岛的一次大雾过程进行了包括GTS(Global Telecommunication System)资料、AMDAR(Aircraft Meteorological Data Relay)资料和9210资料的不同资料组合的三维变分同化试验,以及时间间隔分别为6、3和1h不同时间频率的循环同化试验,并以同化分析场为初始场进行了36h的模拟试验。对同化分析场和模拟结果进行了分析,分析结果表明,采用三维变分方法同化AMDAR等多种非常规观测资料后,分析场均有明显的改变,对雾区的模拟结果也有局部不同程度的修正。进一步分析起修正作用的原因得知同化资料后对低层的湿度和层结趋稳性有所改善。同化GTS资料对低层的增湿贡献明显,但对层结趋稳性贡献不大;而同化AMDAR资料主要使层结趋稳性明显,对增湿无贡献;9210资料对低层湿度和层结趋稳性均有贡献。不同时间间隔的循环同化试验表明,多时次的循环同化比单时次的同化分析增量要大,逐时循环同化与6和3h循环同化相比,可明显改善模拟效果。     

Impact of ECMWF,NCEP, and NCMRWF global model analysis on the WRF model forecast over Indian Region

The global model analysis has significant impact on the mesoscale model forecast as global model provides initial condition (IC) and lateral boundary conditions (LBC) for the mesoscale model. With this objective, four operational global model analyses prepared from the European Centre for Medium-Range Weather Forecasts (ECMWF), National Centers for Environmental Prediction (NCEP) Global Data Assimilation System (GDAS), NCEP Global Forecasting System (GFS), and National Centre for Medium Range Weather Forecasting (NCMRWF) are used daily to generate IC and LBC of the mesoscale model during 13th December 2012 to 13th January 2013. The Weather Research and Forecasting (WRF) model version 3.4, broadly used for short-range weather forecast, is adopted in this study as mesoscale model. After initial comparison of global model analyses with Atmospheric Infrared Sounder (AIRS) retrieved temperature and moisture profiles, daily WRF model forecasts initialized from global model analyses are compared with in situ observations and AIRS profiles. Results demonstrated that forecasts initialized from the ECMWF analysis are closer to AIRS-retrieved profiles and in situ observations compared to other global model analyses. No major differences are occurred in the WRF model forecasts when initialized from the NCEP GDAS and GFS analyses, whereas these two analyses have different spatial resolutions and observations used for assimilation. Maximum RMSD is seen in the NCMRWF analysis-based experiments when compared with AIRS-retrieved profiles. The rainfall prediction is also improved when WRF model is initialized from the ECMWF analysis compared to the NCEP and NCMRWF analyses.     

Numerical Simulation of Intense Precipitation Events South of the Alps: Sensitivity to Initial Conditions and Horizontal Resolution

Summary  In this paper we describe the results of several numerical experiments performed with the limited area model LAMBO, based on a 1989 version of the NCEP (National Center for Environmental Prediction) ETA model, operational at ARPA-SMR since 1993. The experiments have been designed to assess the impact of different horizontal resolutions and initial conditions on the quality and detail of the forecast, especially as regards the precipitation field in the case of severe flood events. For initial conditions we developed a mesoscale data assimilation scheme, based on the nudging technique. The scheme makes use of upper air and surface meteorological observations to modify ECMWF (European Centre for Medium Range Weather Forecast) operational analyses, used as first-guess fields, in order to better describe smaller scales features, mainly in the lower troposphere. Three flood cases in the Alpine and Mediterranean regions have been simulated with LAMBO, using a horizontal grid spacing of 15 and 5 km and starting either from ECMWF initialised analysis or from the result of our mesoscale analysis procedure. The results show that increasing the resolution generally improves the forecast, bringing the precipitation peaks in the flooded areas close to the observed values without producing many spurious precipitation patterns. The use of mesoscale analysis produces a more realistic representation of precipitation patterns giving a further improvement to the forecast of precipitation. Furthermore, when simulations are started from mesoscale analysis, some model-simulated thermodynamic indices show greater vertical instability just in the regions where strongest precipitation occurred. Received March 2, 1999/Revised May 30, 1999     

Four-dimensional variational data assimilation for mesoscale and storm-scale applications

Summary ?The status and progress of the four-dimensional variational data assimilation (4DVAR) are briefly reviewed focusing on application to prediction of mesoscale/storm-scale atmospheric phenomena. Theoretical background is provided for each important component of the 4DVAR system – forecast and adjoint models, observations, background, cost function, preconditioning, and minimization. An overview of practical issues specific for mesoscale/storm-scale 4DVAR is then presented in terms of high-resolution observations, nonlinearity and discontinuity problem, model error, errors from lateral boundary condition, and precipitation assimilation. Practical strategies for efficient and simplified 4DVAR are also introduced, e.g., incremental 4DVAR, poor man’s 4DVAR, and inverse 3DVAR. A new concept on hybrid approach is proposed to combine an efficient 4DVAR scheme and the standard 4DVAR scheme aiming at reducing computational demand required by the standard 4DVAR while improving the accuracy of the simplified 4DVAR. Applications to both hydrostatic and nonhydrostatic models are illustrated and our vision on opportunities and directions for future research is provided. Received March 12, 2001; revised July 24, 2001; accepted September 5, 2001