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中,变分同化后对降水中心位置和强度的预报则没有改善,产生这种现象的原因可能是由于定义的目标函数中没有加进背景场项,也可能是由于采用的观测资料时次比较少,还需要进一步进行研究和试验。     

GRAPES全球切线性和伴随模式的调优

伴随技术是四维变分同化（4DVar）系统中计算代价函数梯度的最佳办法,切线性和伴随模式的效果和效率直接影响着4DVar系统的发展。基于GRAPES（Global and Regional Assimilation PrEdiction System）全球切线性和伴随模式1.0版本,利用GRAPES全球模式2.0版本在并行框架和性能等方面的改善,重新优化和设计了GRAPES全球切线性伴随模式2.0版本,提高了GRAPES全球切线性和伴随模式的效果和效率,优化了切线性模式程序结构,使其计算时间最优可控制在非线性模式的1.2倍以内;采用在切线性模式中保存基态的方法,重构了伴随模式的程序结构,使其计算时间最优控制在非线性模式的1.5倍以内;在GRAPES全球切线性物理过程的设计中,将线性物理过程的轨迹基态计算和切线性扰动计算解耦,提高了GRAPES全球切线性和伴随模式的计算效果和效率。     

改进的切线性模式对一个边界层模式变分资料同化的改善

大气、海洋中边界层物理过程的强非线性会对切线性模式近似（尤其是当积分时间长时）有较大的影响，从而给相应的4维变分资料同化问题的求解造成困难。本构造了Mellor-Yamada湍流封闭模式（level 2.5)的一个改进的切线性模式，相比通常的切线性模式和简化的切线性模式可以提高对非线性扰动的逼近。利用这个改进的切线性模式的伴随模式分别进行了1到7天的变分资料同化试验，得到了满意的结果。而用通常的伴随模式和简化的伴随模式都无法得到令人满意的结果。     

Adjoint sensitivity study on idealized explosive cyclogenesis

The adjoint sensitivity related to explosive cyclogenesis in a conditionally unstable atmosphere is investigated in this study.The PSU/NCAR limited-area,nonhydrostatic primitive equation numerical model MM5 and its adjoint system are employed for numerical simulation and adjoint computation,respectively.To ensure the explosive development of a baroclinic wave,the forecast model is initialized with an idealized condition including an idealized two-dimensional baroclinic jet with a balanced three-dimensional moderateamplitude disturbance,derived from a potential vorticity inversion technique.Firstly,the validity period of the tangent linear model for this idealized baroclinic wave case is discussed,considering different initial moisture distributions and a dry condition.Secondly,the 48-h forecast surface pressure center and the vertical component of the relative vorticity of the cyclone are selected as the response functions for adjoint computation in a dry and moist environment,respectively.The preliminary results show that the validity of the tangent linear assumption for this idealized baroclinic wave case can extend to 48 h with intense moist convection,and the validity period can last even longer in the dry adjoint integration.Adjoint sensitivity analysis indicates that the rapid development of the idealized baroclinic wave is sensitive to the initial wind and temperature perturbations around the steering level in the upstream.Moreover,the moist adjoint sensitivity can capture a secondary high sensitivity center in the upper troposphere,which cannot be depicted in the dry adjoint run.     

伴随系统及非线性优化方法在REM模式可预报性研究中的实际个例应用

利用REM模式的伴随系统和非线性优化方法，通过三个实际天气个例，对REM模式的可预报性问题进行了研究。结果表明，REM模式在给定的实际应用中可接受的预报误差范围内，对三个天气个例都具有预报能力。对于个例一，利用现有的常规报文初始观测场，进行简单的插值处理（最优插值等），REM数值模式就可以得到比较满意的预报结果； 对于个例二和个例三，对现有的报文初始观测场进行处理（如四维变分资料同化）后，REM模式在给定的误差允许范围内，对这两个天气个例仍得到满意的预报。研究结果不仅对改进数值模式具有一定的指导意义，而且对如何改进数值模式的初值问题，特别是在中尺度天气预报中如何改进具有一定的参考价值。     

The linearization and adjoint of radiation transfer processes in the NMC spectral model part I: Solar radiative transfer

Summary Linearization and adjoint-model derivation for the solar radiation transfer codes in the NMC spectral model have been carried out. Verification of the validity of resulting tangent linear model and the correctness of the corresponding adjoint have been performed. Applications of derived adjoint model are considered, including parameter estimation for inputs to solar radiation codes with aid of the physics (i.e., the solar radiation codes) and a sensitivity study of the downward solar radiation flux at the earth surface with respect to water vapor amount at various heights.With 5 Figures     

包含“开关”变量的MM5四维变分资料同化系统

数值模式中的湿物理过程对于模拟各种尺度的天气现象非常重要,但它会使得模式方程的状态变量不连续,这种不连续在模式程序中表现为“开关”变量。本文对包括Kuo降水参数化方案的非静力中尺度数值模式MM5的四维变分资料同化系统进行研究,利用对实际降水的模拟,比较详细地讨论了有关的“开关”变量以及切向线性化和伴随问题。结果表明：对于一个离散的数值模式,保持“开关”变量与基态一致,用伴随码方法构造的伴随模式计算的梯度值能够为最小化过程提供较好的下降方向;Kuo方案中对流每隔一个积分步的交替发生并不影响目标函数最小化的收敛速度;“开关”变量的存在也不影响将风、温度、气压和比湿结合起来同化对MM5降水预报准确性的提高。     

CMA-GFS 4DVar边界层过程线性化的改进北大核心CSCD

持续发展和优化切线性模式的线性化物理过程,保持与非线性模式一致是改善四维变分同化(4DVar)分析和预报效果的有效方法之一。目前业务系统的CMA-GFS模式采用基于Charney-Phillips(C-P)跳点的边界层参数化方案,而CMA-GFS 4DVar系统中采用基于Lorenz跳点的边界层线性化方案。为改善CMA-GFS 4DVar系统的边界层分析和预报效果,基于C-P跳点的边界层参数化方案研发了新边界层线性化方案,并通过对方案中地表热量通量和水汽通量扰动、自由大气的理查逊系数扰动、边界层的热量和动量交换系数扰动等进行更加精细地规约化约束,在确保CMA-GFS切线性和伴随模式稳定运行的情况下,减少线性化过程对切线性模式预报精度的影响。切线性近似试验检验表明:相较于原方案,新边界层线性化方案可以减少边界层位温和比湿的相对误差,最大可减少10%。批量4DVar循环同化试验表明:新边界层线性化方案可以有效改善切线性模式对低层位温、风场和比湿扰动的预报精度,减少4DVar内外循环目标泛函的相对差异,并提高700 hPa位势高度的可预报时效。     

基于总能量模的GRAPES全球模式奇异向量扰动研究

以发展基于奇异向量技术为初值扰动的GRAPES全球集合预报系统为目的,在GRAPES模式及其干动力框架下的切线性、伴随模式基础上开展了以总能量模为权重算子的奇异向量计算技术研究,建立奇异向量的计算求解模块,并通过奇异向量检验方法和切线性近似方法验证了奇异向量求解的正确性.通过对中高纬度的GRAPES奇异向量水平结构的线性演变分析,证实了在最优时间间隔内GRAPES奇异向量能够快速增长,并能描述中高纬度大气的斜压不稳定特征.分析在初始时刻和最优化时间间隔时刻的GRAPES奇异向量总能量及其分量(动能和势能)的垂直分布特征,发现在中高纬度区域,GRAPES奇异向量能够描述对流层不同层次的斜压不稳定增长特征.     

Adjoint-based Sensitivity Analysis of a Mesoscale Low on the Mei-yu Front and Its Implicationsfor Adaptive Observation

An adjoint sensitivity analysis of one mesoscale low on the mei-yu Front is presented in this paper.The sensitivity gradient of simulation error dry energy with respect to initial analysis is calculated. And after verifying the ability of a tangent linear and adjoint model to describe small perturbations in the nonlinear model, the sensitivity gradient analysis is implemented in detail. The sensitivity gradient with respect to different physical fields are not uniform in intensity, simulation error is most sensitive to the vapor mixed ratio. The localization and consistency are obvious characters of horizontal distribution of the sensitivity gradient, which is useful for the practical implementation of adaptive observation. The sensitivity region tilts to the northwest with height increasing; the singular vector calculation proves that this tilting characterizes a quick-growing structure, which denotes that using the leading singular vectors to decide the adaptive observation region is proper. When connected with simulation of a mesoscale low on the mei-yu Front, the sensitivity gradient has the following physical characters: the obvious sensitive region is mesoscale, concentrated in the middle-upper troposphere, and locates around the key system; and the sensitivity gradient of different physical fields correlates dynamically.     

Adjoint-based Sensitivity Analysis of a Mesoscale Low on the Mei-yu Front and Its Implications for Adaptive Observation

An adjoint sensitivity analysis of one mesoscale low on the mei-yu Front is presented in this paper. The sensitivity gradient of simulation error dry energy with respect to initial analysis is calculated. And after verifying the ability of a tangent linear and adjoint model to describe small perturbations in the nonlinear model, the sensitivity gradient analysis is implemented in detail. The sensitivity gradient with respect to different physical fields are not uniform in intensity, simulation error is most sensitive to the vapor mixed ratio. The localization and consistency are obvious characters of horizontal distribution of the sensitivity gradient, which is useful for the practical implementation of adaptive observation. The sensitivity region tilts to the northwest with height increasing; the singular vector calculation proves that this tilting characterizes a quick-growing structure, which denotes that using the leading singular vectors to decide the adaptive observation region is proper. When connected with simulation of a mesoscale low on the mei-yu Front, the sensitivity gradient has the following physical characters： the obvious sensitive region is mesoscale, concentrated in the middle-upper troposphere, and locates around the key system; and the sensitivity gradient of different physical fields correlates dynamically.     

Application and Characteristic Analysis of the Moist Singular Vector in GRAPES-GEPS

The singular vector (SV) initial perturbation method can capture the fastest-growing initial perturbation in a tangent linear model (TLM). Based on the global tangent linear and adjoint model of GRAPES-GEPS (Global/Regional Assimilation and Prediction System—Global Ensemble Prediction System), some experiments were carried out to analyze the structure of the moist SVs from the perspectives of the energy norm, energy spectrum, and vertical structure. The conclusions are as follows: The evolution of the SVs is synchronous with that of the atmospheric circulation, which is flow-dependent. The moist and dry SVs are located in unstable regions at mid-to-high latitudes, but the moist SVs are wider, can contain more small- and medium-scale information, and have more energy than the dry SVs. From the energy spectrum analysis, the energy growth caused by the moist SVs is reflected in the relatively small-scale weather system. In addition, moist SVs can generate perturbations associated with large-scale condensation and precipitation, which is not true for dry SVs. For the ensemble forecasts, the average anomaly correlation coefficient of large-scale circulation is better for the forecast based on moist SVs in the Northern Hemisphere, and the low-level variables forecasted by the moist SVs are also improved, especially in the first 72 h. In addition, the moist SVs respond better to short-term precipitation according to statistical precipitation scores based on 10 cases. The inclusion of the large-scale condensation process in the calculation of SVs can improve the short-term weather prediction effectively.     

An Economical Approach to Four-dimensional Variational Data Assimilation

Four-dimensional variational data assimilation (4DVar) is one of the most promising methods to provide optimal analysis for numerical weather prediction (NWP). Five national NWP centers in the world have successfully applied 4DVar methods in their global NWPs, thanks to the increment method and adjoint technique. However, the application of 4DVar is still limited by the computer resources available at many NWP centers and research institutes. It is essential, therefore, to further reduce the computational cost of 4DVar. Here, an economical approach to implement 4DVar is proposed, using the technique of dimension-reduced projection (DRP), which is called ``DRP-4DVar." The proposed approach is based on dimension reduction using an ensemble of historical samples to define a subspace. It directly obtains an optimal solution in the reduced space by fitting observations with historical time series generated by the model to form consistent forecast states, and therefore does not require implementation of the adjoint of tangent linear approximation. To evaluate the performance of the DRP-4DVar on assimilating different types of mesoscale observations, some observing system simulation experiments are conducted using MM5 and a comparison is made between adjoint-based 4DVar and DRP-4DVar using a 6-hour assimilation window.     

GRAPES伴随模式底层数据栈优化

GRAPES伴随模式是其四维变分同化系统的核心组成部分。由于其计算过程复杂,临时数据较多,实现中采用断点存储策略可以有效减少伴随模式的计算时间和存储空间。极限断点存储策略是在单积分步内以全存储策略实现为基础,将其中部分基态以计算代替的一种类断点存储策略。在该策略的支持下,需要一种新的数据管理结构,来保证程序的正确运行。文章提出了在已有栈基础上优化的新数据存储管理方式——嵌套多链栈,这种结构可有效满足使用极限断点存储技术实现GRAPES伴随模式的初态管理需求。试验表明:相比断点存储技术,在总内存增加不超过30%的情况下可使GRAPES的运行效率提高1倍。     

CMA-MESO三维变分同化系统2 m相对湿度资料同化研究

为了提高CMA-MESO （China Meteorological Administration Mesoscale model）（原GRAPES）三维变分同化系统中2 m相对湿度资料的应用效果,改善模式中相对湿度的分析和降水预报效果,分析了2015年6—8月T639（T639L60全球中期数值预报系统,0.28125°×0.28125°）分析场低层相对湿度和2 m相对湿度之差与表征稳定度的理查森数（Ri）的关系,发现二者有很好的相关,Ri＜0时,模式低层相对湿度与2 m相对湿度的差异较小,基本在同化观测误差范围内。依据该统计结果,对CMA-MESO同化系统中2 m相对湿度同化方案进行优化,Ri＜0时,将观测站地形低于模式地形的2 m相对湿度观测由观测站高度改为模式最低层高度进行同化,形成新的2 m相对湿度同化方案,旨在解决2 m相对湿度资料同化时模式地形高度与观测站高度不同的影响。2018年7月CMA-MESO三维变分同化系统（3DVar）个例和连续试验结果显示:新的2 m相对湿度同化方案同化分析资料数量明显增加,且08时多于20时（北京时）,新增观测点新息向量（背景减观测）与周围原有观测新息向量保持基本一致,分析残差偏差和均方根误差减小,降水预报效果明显改善。新2 m相对湿度同化方案通过提高观测站低于模式地形高度的观测资料合理应用,从而改善了3 km模式系统同化分析和预报效果。      

基于历史预报的四维变分资料同化（4DVar）方法中的滤波

基于历史预报的四维变分同化方法在降维的样本空间最小化代价函数,避免切线性伴随模式,是一种比较经济的算法。但是因为选取的集合样本不可能无限多,实际样本数远远小于观测资料数以及模式变量的自由度,会导致观测站点和模式格点间产生虚假的相关。介绍了在历史预报4DVar中引入的局地化滤波技术,并通过3组试验,比较了局地化前后的结果。试验结果表明:引入局地化技术后,其能有效滤去初始场中的虚假相关关系,同时由于Schur算子的作用,滤波后的分析场是光滑连续的。而且实际个例研究结果表明,局地化后能改进6 h和12 h累积降水的均方根误差,进一步提高预报效果。       

用一种新的同化方法同化降水量资料

Observations of accumulated precipitation are extremely valuable for effectively improving rainfall analysis and forecast. It is, however, difficult to use such observations directly through sequential assimilation methods, such as three-dimensional variational data assimilation or an Ensemble Kalman Filter. In this study, the authors illustrate a new approach that makes effective use of precipitation data to improve rainfall forecast. The new method directly obtains an optimal solution in a reduced space by fitting observations with historical time series generated by the model; it also avoids the implementation of tangent linear model and its adjoint. A lot of historical samples are produced as the ensemble of precipitation observations with the fully nonlinear forecast model. The results show that the new approach is capable of extracting information from precipitation observations to improve the analysis and forecast. This method provides comparable performance with the standard four- dimensional variational data assimilation at a much lower computational cost.     

面向四维变分资料同化的NSAS积云深对流参数化方案的简化及线性化研究

GRAPES全球四维变分资料同化系统需要积云深对流参数化方案的线性化与伴随方案,直接采用原始复杂参数化方案进行线性化并不可行,需要发展简化光滑方案来减缓非线性与非连续性特征。GRAPES全球模式采用NSAS积云对流参数化方案,积云深对流对环境的反馈主要通过补偿下沉来实现,研究突出补偿下沉作用,忽略降水蒸发、动量反馈等贡献,形成简化方案。采用输入温、湿度廓线加入不同幅度小扰动方法,评估参数化方案计算的温度、比湿时间倾向对输入扰动的敏感性,检验非线性与非连续特征。提出避免或减缓非连续“开关”的方法,在简化方案的基础上发展了简化光滑方案。简化光滑方案与原始积云深对流方案相比,在对流触发上一致,在对流的位温与比湿倾向、降水的时序模拟等方面相似,而在减缓非线性、避免非连续性方面显著优于原始方案。基于简化光滑方案发展的线性化方案表明,对小于2倍分析增量幅度的扰动,线性化方案可以较好地模拟非线性方案的扰动发展。发展的简化光滑方案具有合理性和实用性。      

基于奇异矢量的优化短期集合预报

在1-2 d的短期预报中,由奇异矢量构建的初始扰动主要是线性发展,为了防止在积分终止时刻,由同一奇异矢量导出的正负初始扰动的积分在集合平均时互相抵消,文中首先通过理论推导和实际计算证明了对集合成员进行优化的必要性,以及从不同奇异矢量导出的集合成员中,表现好于控制预报的一组成员相对于控制预报的离差恒大于或恒小于表现劣于控制预报的另一组成员,利用这个特征,在做集合预报时,把奇异矢量导出的正负两组预报分成离差相对大一组、离差相对小一组,就可以避免求集合平均时成员相互抵消,从而提出了一种优化基于奇异矢量的短期集合预报的方法.文中使用NCAR/PSU(美国国家大气研究中心/宾夕法尼亚大学)中尺度有限区域模式MM5第1版,及其对应切线性、伴随模式,对1999年夏季发生的两个梅雨锋低涡个例作了分析,在计算奇异矢量时采用了干能量模,分析结果表明:相对于正负两个初始扰动都入选的集合,严格按照这种方法挑选出来的优化集合可以有效地提高集合平均的精确度.在生成初始扰动的方法上,文中的计算表明,相对于用单个奇异矢量生成初始扰动,把正交的多个奇异矢量累加起来导出的初始扰动具有更大的增长率,能有效地增大集合成员间的离差,提高集合成员的预报精度.     

Conjugate Gradient Algorithm in the Four-Dimensional Variational Data Assimilation System in GRAPES

Minimization algorithms are singular components in four-dimensional variational data assimilation (4DVar). In this paper, the convergence and application of the conjugate gradient algorithm (CGA), which is based on the Lanczos iterative algorithm and the Hessian matrix derived from tangent linear and adjoint models using a non-hydrostatic framework, are investigated in the 4DVar minimization. First, the influence of the Gram-Schmidt orthogonalization of the Lanczos vector on the convergence of the Lanczos algorithm is studied. The results show that the Lanczos algorithm without orthogonalization fails to converge after the ninth iteration in the 4DVar minimization, while the orthogonalized Lanczos algorithm converges stably. Second, the convergence and computational efficiency of the CGA and quasi-Newton method in batch cycling assimilation experiments are compared on the 4DVar platform of the Global/Regional Assimilation and Prediction System (GRAPES). The CGA is 40% more computationally efficient than the quasi-Newton method, although the equivalent analysis results can be obtained by using either the CGA or the quasi-Newton method. Thus, the CGA based on Lanczos iterations is better for solving the optimization problems in the GRAPES 4DVar system.