热带气旋奇异向量在GRAPES全球集合预报中的初步应用

基于副热带奇异向量的初值扰动方法已应用于GRAPES （Global and Regional Assimilation PrEdiction System）全球集合预报系统,但存在热带气旋预报路径离散度不足的问题。通过分析发现,热带气旋附近区域初值扰动结构不合理导致预报集合不能较好地估计热带气旋预报的不确定性,是路径集合离散度不足的可能原因之一。通过建立热带气旋奇异向量求解方案,将热带气旋奇异向量和副热带奇异向量共同线性组合生成初值扰动,以弥补热带气旋区域初值扰动结构不合理这一缺陷,进而改进热带气旋集合预报效果。利用GRAPES全球奇异向量计算方案,以台风中心10个经纬度区域为目标区构建热带气旋奇异向量求解方案,针对台风“榕树”个例进行集合预报试验,并开展批量试验,利用中国中央气象台最优台风路径和中国国家气象信息中心的降水观测资料进行检验,对比分析热带气旋奇异向量结构特征和初值扰动特征,评估热带气旋奇异向量对热带气旋路径集合预报和中国区域24 h累计降水概率预报技巧的影响。结果表明,热带气旋奇异向量具有局地化特征,使用热带气旋奇异向量之后,热带气旋路径离散度增加,路径集合平均预报误差和离散度的关系得到改善,路径集合平均预报误差有所减小,集合成员更好地描述了热带气旋路径的预报不确定性;中国台风降水的小雨、中雨、大雨、暴雨各量级24 h累计降水概率预报技巧均有一定提高。总之,当在初值扰动的生成中考虑热带气旋奇异向量后,可改进热带气旋初值扰动结果,并有助于改善热带气旋路径集合预报效果。      

基于条件非线性最优扰动的目标观测中瞄准区 不同引导性变量的影响试验研究

基于GRAPES区域业务预报模式,采用一种快速算法计算出来的条件非线性最优扰动对实际台风个例麦莎(No.0509)开展了目标观测研究,应用数值模式,进行一系列的敏感性试验,讨论了与目标观测设计相关的一些问题,包括确定瞄准区时使用不同的引导性变量对目标观测效果的影响、及瞄准区范围变化对预报效果的影响。文中分别以提高麦莎在检验区(20.125°—35.3125°N,116.8125°—129.75°E)内的24 h海平面气压预报和24 h累积降水量预报为目的,基于条件非线性最优扰动使用了3种不同的引导性变量寻找敏感区(又称瞄准区),对这些敏感区的分布特点和有效性进行了比较和讨论。试验结果表明,在使用的3种引导性变量中,用不同的引导性变量识别的敏感区是有差别的,总体上说,文中使用的3种引导性变量识别的瞄准区对提高预报都是有效的,特别是第2和第3种的效果更好些,且两者识别的瞄准区常显示出类似的特点。文中进一步针对检验区内24 h累积降水量预报误差问题,将前面确定的瞄准区范围扩大相同的幅度,讨论瞄准区范围变化对改进预报的影响。试验结果表明,增加瞄准格点数,有可能使预报效果得到改善,但是试验结果同时也暗示了单纯靠扩大瞄准...     

条件非线性最优扰动方法在适应性观测研究中的初步应用

针对适应性观测中敏感性区域的确定问题，考虑初始误差对预报结果的影响, 比较了条件非线性最优扰动(CNOP)与第一线性奇异向量(FSV)在两个降水个例中的空间结构的差异，考察了它们总能量范数随时间发展演变的异同。结合敏感性试验的分析，揭示了预报结果对CNOP类型的初始误差的敏感性要大于对FSV类型的初始误差的敏感性，因而减少初值中CNOP类型误差的振幅比减少FSV类型的收益要大。这一结果表明可以把CNOP方法应用于适应性观测来识别大气的敏感区。关于将CNOP方法有效地应用于适应性观测所面临的挑战及需要采取的对策等也进行了讨论。     

条件非线性最优扰动在可预报性问题研究中的应用

本文总结了近年来条件非线性最优扰动方法的应用研究的主要进展.主要包括四个方面:(1)将条件非线性最优扰动(CNOP)方法拓展到既考虑初始扰动又考虑模式参数扰动,形成了拓展的CNOP方法.拓展的CNOP方法不仅能够应用于研究分别由初始误差和模式参数误差导致的可预报性问题,而且能够用于探讨初始误差和模式参数误差同时存在的情形;(2)将拓展的CNOP方法分别应用于ENSO和黑潮可预报性研究,考察了初始误差和模式参数误差对其可预报性的影响,揭示了初始误差在导致ENSO和黑潮大弯曲路径预报不确定性中的重要作用;(3)考察了阻塞事件发生的最优前期征兆(OPR)以及导致其预报不确定性的最优增长初始误差(OGR),揭示了OPR和OGR空间模态及其演变机制的相似性及其局地性特征;(4)研究了台风预报的目标观测问题,用CNOP方法确定了台风预报的目标观测敏感区,通过观测系统模拟试验(OSSEs)和/或观测系统试验(OSEs),表明了CNOP敏感区在改进台风预报中的有效性.具体地,台风OGR的主要误差分布在某一特定区域,空间分布具有明显的局地性特征,在台风OGR的局地性区域增加观测,有效改进了台风的预报技巧,该区域代表了台风预报的初值敏感区.事实上,上述El Ni(n)o事件、黑潮路径变异以及阻塞事件的OGR的空间分布也具有明显的局地性特征,这些事件的OGR刻画的局地性区域可能也代表了初值敏感区.     

条件非线性最优扰动(CNOP):简介与数值求解

介绍了条件非线性最优扰动(Conditional Nonlinear Optimal Perturbation,CNOP)的定义及其在大气和海洋等可预报性研究中的应用。根据研究对象不同,CNOP分为与初始扰动有关的CNOP(CNOP-I)方法、与模式参数扰动有关的CNOP(CNOP-P)方法和同时考虑初始扰动和模式参数扰动的CNOP方法。目前,CNOP-I方法已经应用于ENSO、黑潮和阻塞可预报性以及热盐环流和草原生态系统稳定性的研究。此外,CNOP-I方法也被应用于探讨台风目标观测的研究,利用CNOP-I方法能够识别出台风预报的初值敏感区,通过观测系统模拟试验表明在初值敏感区增加观测能够有效改进台风的预报技巧。CNOP-P方法也在ENSO和黑潮可预报性以及热盐环流和草原生态系统稳定性研究中得到了应用。为了将CNOP方法应用于更多的领域,本文利用一个简单的Burgers方程,介绍了如何通过建立Burgers方程的切线性模式和伴随模式,从而利用非线性最优化算法计算获得CNOP。这一数值试验为将CNOP方法应用于更多的领域提供了借鉴。     

GRAPES全球模式静力平衡奇异向量改进及应用试验

采用线性化物理过程方案的GRAPES全球模式奇异向量在进行非线性模式积分时会有部分奇异向量出现崩溃问题,这说明奇异向量结构可能存在扰动变量之间不协调之处,需要对奇异向量扰动的计算方法优化,进而改进基于奇异向量的集合预报初值扰动,提高GRAPES全球集合预报效果。基于原有的GRAEPS全球奇异向量计算方法,在求解奇异向量时,对气压扰动的处理进行改进,将初始时刻的气压扰动分量通过位温扰动根据静力平衡关系导出获得,其他保持一致,发展了静力平衡奇异向量改进方法。基于有两个台风过程的个例（2019年8月8日12时（世界时））,分别采用原奇异向量方法和静力平衡奇异向量改进方法进行热带气旋目标区奇异向量的计算求解,并进行相应奇异向量的非线性模式积分,对比分析奇异向量非线性积分的稳定性。进而,对比分析奇异向量求解方法改进前、后热带气旋奇异向量的结构特征和初值扰动特征,开展了集合预报试验,评估改进后的奇异向量求解方法对GRAPES全球集合预报系统预报性能的影响。试验结果表明,静力平衡奇异向量改进方法通过产生协调的气压扰动和位温扰动场,解决了奇异向量非线性积分崩溃的问题,消除了原来不利于积分稳定性的气压扰动过于局地化的小尺度结构。静力平衡奇异向量改进方法对奇异向量中位温扰动分量和纬向风扰动分量结构影响较小,使得气压扰动分量的大值区位于台风附近,更好地描述热带气旋初值不确定性,与位温扰动分量的分布更加协调。采用静力平衡奇异向量改进方法,可以提高GRAPES全球集合预报在北半球和南半球等压面要素集合预报技巧和中国地区24 h累计降水概率预报技巧,增大台风路径集合离散度。      

条件非线性最优扰动在长江中下游地区冬季暴雨中的应用研究

为了提高长江中下游地区高影响天气的数值预报,利用条件非线性最优扰动(CNOP)方法,对一次长江中下游地区冬季降水个例(高影响天气事件)进行目标观测研究,并通过观测系统模拟试验(OSSE)检验了该方法确定敏感区的有效性和可行性。试验结果表明,CNOP方法可有效识别对应于高影响天气事件的敏感区。通过对敏感区进行初始场修正后,可明显改善验证区内24 h累积降水预报误差和总能量预报误差。进一步分析发现,通过改善敏感区内的初始场信息(如水汽通量场和低层冷空气活动等),使得数值模式不仅能更真实刻画该天气系统的初始结构,还能更好模拟出该天气系统随时间的演变特征,因而减少了验证区内对该天气系统的预报误差。这一结果表明可以把CNOP方法应用于长江中下游地区高影响天气事件的目标观测研究或实践中。      

模式含有开关时遗传算法求解条件非线性最优扰动的有效性研究

在基于条件非线性最优扰动(CNOP)的台风适应性观测研究中,针对预报模式的湿物理参数化产生的“on-off”开关导致传统伴随方法不能为最优化过程提供正确梯度这一现象,将模式含有“on-off”开关时求解CNOP的问题视为非光滑最优化问题,引入遗传算法,在给出详细的算法流程后,以一个在强迫项中含“on-off”开关的理想模式,分析了“on-off”开关对求解CNOP的影响,三个数值试验检验了模式含有“on-off”开关时遗传算法求解CNOP的有效性,并分析了不同初始种群对最优化结果的影响。结果显示,所采用的含有“on-off”开关的理想模式下,遗传算法能有效求解CNOP,最后对遗传算法求解CNOP的优缺点进行了详细讨论。     

伴随敏感性方法、第一奇异向量方法以及条件非线性最优扰动方法在台风目标观测敏感区识别中的比较研究

本文通过深入分析伴随敏感性（ADS）方法、第一奇异向量（LSV）方法、以及条件非线性最优扰动（CNOP）方法在目标观测敏感区识别方面的原理,提出了非线性程度的概念和计算方法,考察了转向型和直线型台风的非线性程度,分析了上述三种方法在不同非线性程度下识别的敏感区的异同,同时对比了转向型和直线型台风的敏感区的差异,并通过敏感性试验探讨了在不同非线性程度下以及在转向型与直线型台风中,预报对敏感区内初值的敏感性程度,进而探讨台风目标观测在不同情况下的有效性。结果表明,转向型台风的非线性程度差别比较大,或者特别强,或者特别弱;而直线型台风非线性程度居中,不同台风个例之间的非线性程度差别较小。对于非线性较弱的台风,三种方法识别的敏感区较为相似,而对于非线性较强的台风,LSV方法与ADS方法识别的敏感区较为相似,但是与CNOP方法识别的敏感区具有较大的差别。对于转向型台风,敏感区主要位于行进路径的右前方,而对于直线型台风,敏感区主要位于初始台风位置的后方。敏感性试验表明,不论台风非线性强弱,转向还是直行,CNOP敏感区内的随机扰动发展最大,而LSV敏感区内叠加的随机扰动发展次之,ADS敏感区内叠加的扰动发展最小;此外,非线性弱的台风,扰动的发展大于非线性强的台风的扰动的发展,表明非线性弱的台风预报受初值影响更大,目标观测的效果可能会更明显。     

ETKF初值扰动方法中真实观测及扰动调节因子研究

目前国家气象中心业务GRAPES区域集合预报系统中集合变换卡尔曼滤波（ETKF）方法采用的是模拟观测信息，为进一步完善ETKF方法，拟对ETKF初值扰动通过引入真实探空观测资料，使扰动场能够代表真实观测的不确定信息，改善区域集合预报技巧。真实观测资料的引入会使得每日的观测数目和分布发生变化，这对ETKF方法而言可能会引起扰动振幅的不稳定，因此在引入真实观测资料的基础上设计了新的扰动振幅调节因子，通过格点空间中离散度和均方根误差关系来对初值扰动振幅进行自适应调整。从初值扰动结构、概率预报技巧以及降水预报效果等方面对比分析了基于模拟观测、真实观测以及真实观测结合新型调节因子的ETKF方案的差异，结果表明:真实探空资料能够有效应用于GRAPES区域集合预报系统中，真实观测资料与模拟观测资料相比较为稀疏，可以获得更大量级的初值扰动振幅；真实观测资料有助于提高区域集合的离散度，但对集合预报准确度以及概率预报结果的提高有限，对于降水预报效果提高也有限；新型的扰动振幅调节因子可以有效获得稳定的初值扰动振幅，并保持ETKF扰动结构，真实观测资料与扰动振幅自适应调节因子相结合，可以有效提高区域集合的概率预报结果，并有效提高降水预报效果。     

A Fast Algorithm for Solving CNOP and Associated Target Observation Tests

Conditional Nonlinear Optimal Perturbation （CNOP） is a new method proposed by Mu et al. in 2003, which generalizes the linear singular vector （LSV） to include nonlinearity. It has become a powerful tool for studying predictability and sensitivity among other issues in nonlinear systems. This is because the CNOP is able to represent, while the LSV is unable to deal with, the fastest developing perturbation in a nonlinear system. The wide application of this new method, however, has been limited due to its large computational cost related to the use of an adjoint technique. In order to greatly reduce the computational cost, we hereby propose a fast algorithm for solving the CNOP based on the empirical orthogonal function （EOF）. The algorithm is tested in target observation experiments of Typhoon Matsa using the Global/Regional Assimilation and PrEdiction System （GRAPES）, an operational regional forecast model of China. The effectivity and feasibility of the algorithm to determine the sensitivity （target） area is evaluated through two observing system simulation experiments （OSSEs）. The results, as expected, show that the energy of the CNOP solved by the new algorithm develops quickly and nonlinearly. The sensitivity area is effectively identified with the CNOP from the new algorithm, using 24 h as the prediction time window. The 24-h accumulated rainfall prediction errors （ARPEs） in the verification region are reduced significantly compared with the ＂true state,＂ when the initial conditions （ICs） in the sensitivity area are replaced with the ＂observations.＂ The decrease of the ARPEs can be achieved for even longer prediction times （e.g., 72 h）. Further analyses reveal that the decrease of the 24-h ARPEs in the verification region is attributable to improved simulations of the typhoon＇s initial warm-core, upper level relative vorticity, water vapor conditions, etc., as a result of the updated ICs in the sensitivity area.     

Applications of Conditional Nonlinear Optimal Perturbation in Predictability Study and Sensitivity Analysis of Weather and Climate

Considering the limitation of the linear theory of singular vector (SV), the authors and their collaborators proposed conditional nonlinear optimal perturbation (CNOP) and then applied it in the predictability study and the sensitivity analysis of weather and climate system. To celebrate the 20th anniversary of Chinese National Committee for World Climate Research Programme (WCRP), this paper is devoted to reviewing the main results of these studies. First, CNOP represents the initial perturbation that has largest nonlinear evolution at prediction time, which is different from linear singular vector (LSV) for the large magnitude of initial perturbation or/and the long optimization time interval. Second, CNOP, rather than linear singular vector (LSV), represents the initial anomaly that evolves into ENSO events most probably. It is also the CNOP that induces the most prominent seasonal variation of error growth for ENSO predictability; furthermore, CNOP was applied to investigate the decadal variability of ENSO asymmetry. It is demonstrated that the changing nonlinearity causes the change of ENSO asymmetry. Third, in the studies of the sensitivity and stability of ocean’s thermohaline circulation (THC), the nonlinear asymmetric response of THC to finite amplitude of initial perturbations was revealed by CNOP. Through this approach the passive mechanism of decadal variation of THC was demonstrated; Also the authors studies the instability and sensitivity analysis of grassland ecosystem by using CNOP and show the mechanism of the transitions between the grassland and desert states. Finally, a detailed discussion on the results obtained by CNOP suggests the applicability of CNOP in predictability studies and sensitivity analysis.     

Algorithm Studies on How to Obtain a Conditional Nonlinear Optimal Perturbation (CNOP)

The conditional nonlinear optimal perturbation (CNOP), which is a nonlinear generalization of the linear singular vector (LSV), is applied in important problems of atmospheric and oceanic sciences, including ENSO predictability, targeted observations, and ensemble forecast. In this study, we investigate the computational cost of obtaining the CNOP by several methods. Differences and similarities, in terms of the computational error and cost in obtaining the CNOP, are compared among the sequential quadratic programming (SQP) algorithm, the limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm, and the spectral projected gradients (SPG2) algorithm. A theoretical grassland ecosystem model and the classical Lorenz model are used as examples. Numerical results demonstrate that the computational error is acceptable with all three algorithms. The computational cost to obtain the CNOP is reduced by using the SQP algorithm. The experimental results also reveal that the L-BFGS algorithm is the most effective algorithm among the three optimization algorithms for obtaining the CNOP. The numerical results suggest a new approach and algorithm for obtaining the CNOP for a large-scale optimization problem.     

What Kind of Initial Errors Cause the Severest Prediction Uncertainty of EI Nino in Zebiak-Cane Model

With the Zebiak-Cane (ZC) model, the initial error that has the largest effect on ENSO prediction is explored by conditional nonlinear optimal perturbation (CNOP). The results demonstrate that CNOP-type errors cause the largest prediction error of ENSO in the ZC model. By analyzing the behavior of CNOP- type errors, we find that for the normal states and the relatively weak EI Nino events in the ZC model, the predictions tend to yield false alarms due to the uncertainties caused by CNOP. For the relatively strong EI Nino events, the ZC model largely underestimates their intensities. Also, our results suggest that the error growth of EI Nino in the ZC model depends on the phases of both the annual cycle and ENSO. The condition during northern spring and summer is most favorable for the error growth. The ENSO prediction bestriding these two seasons may be the most difficult. A linear singular vector (LSV) approach is also used to estimate the error growth of ENSO, but it underestimates the prediction uncertainties of ENSO in the ZC model. This result indicates that the different initial errors cause different amplitudes of prediction errors though they have same magnitudes. CNOP yields the severest prediction uncertainty. That is to say, the prediction skill of ENSO is closely related to the types of initial error. This finding illustrates a theoretical basis of data assimilation. It is expected that a data assimilation method can filter the initial errors related to CNOP and improve the ENSO forecast skill.     

What Kind of Initial Errors Cause the Severest Prediction Uncertainty of El Nino in Zebiak-Cane Model

With the Zebiak-Cane （ZC） model, the initial error that has the largest effect on ENSO prediction is explored by conditional nonlinear optimal perturbation （CNOP）. The results demonstrate that CNOP-type errors cause the largest prediction error of ENSO in the ZC model. By analyzing the behavior of CNOPtype errors, we find that for the normal states and the relatively weak E1 Nifio events in the ZC model, the predictions tend to yield false alarms due to the uncertainties caused by CNOP. For the relatively strong E1 Nino events, the ZC model largely underestimates their intensities. Also, our results suggest that the error growth of E1 Nifio in the ZC model depends on the phases of both the annual cycle and ENSO. The condition during northern spring and summer is most favorable for the error growth. The ENSO prediction bestriding these two seasons may be the most difficult. A linear singular vector （LSV） approach is also used to estimate the error growth of ENSO, but it underestimates the prediction uncertainties of ENSO in the ZC model. This result indicates that the different initial errors cause different amplitudes of prediction errors though they have same magnitudes. CNOP yields the severest prediction uncertainty. That is to say, the prediction skill of ENSO is closely related to the types of initial error. This finding illustrates a theoretical basis of data assimilation. It is expected that a data assimilation method can filter the initial errors related to CNOP and improve the ENSO forecast skill.     

Applications of Conditional Nonlinear Optimal Perturbation to the Study of the Stability and Sensitivity of the Jovian Atmosphere

A two-layer quasi-geostrophic model is used to study the stability and sensitivity of motions on smallscale vortices in Jupiter’s atmosphere. Conditional nonlinear optimal perturbations (CNOPs) and linear singular vectors (LSVs) are both obtained numerically and compared in this paper. The results show that CNOPs can capture the nonlinear characteristics of motions in small-scale vortices in Jupiter’s atmosphere and show great difference from LSVs under the condition that the initial constraint condition is large or the optimization time is not very short or both. Besides, in some basic states, local CNOPs are found. The pattern of LSV is more similar to local CNOP than global CNOP in some cases. The elementary application of the method of CNOP to the Jovian atmosphere helps us to explore the stability of variousscale motions of Jupiter’s atmosphere and to compare the stability of motions in Jupiter’s atmosphere and Earth’s atmosphere further.     

Nonlinearity and finite-time instability in a T21L3 quasigeostrophic model

In this paper, a nonlinear optimization method is used to explore the finite-time instability of the atmospheric circulation with a three-level quasigeostrophic model under the framework of the conditional nonlinear optimal perturbation (CNOP). As a natural generalization of linear singular vector (SV), CNOP is defined as an initial perturbation that makes the cost function the maximum at a prescribed forecast time under certain physical constraint conditions. Special attentions are paid to the different structures and energy evolutions of the optimal perturbations.     

Optimal Initial Error Growth in the Prediction of the Kuroshio Large Meander Based on a High-resolution Regional Ocean Model

Based on the high-resolution Regional Ocean Modeling System(ROMS) and the conditional nonlinear optimal perturbation(CNOP) method, this study explored the effects of optimal initial errors on the prediction of the Kuroshio large meander(LM) path, and the growth mechanism of optimal initial errors was revealed. For each LM event, two types of initial error(denoted as CNOP1 and CNOP2) were obtained. Their large amplitudes were found located mainly in the upper 2500 m in the upstream region of the LM, i.e., southeast of Kyushu. Furthermore, we analyzed the patterns and nonlinear evolution of the two types of CNOP. We found CNOP1 tends to strengthen the LM path through southwestward extension. Conversely,CNOP2 has almost the opposite pattern to CNOP1, and it tends to weaken the LM path through northeastward contraction.The growth mechanism of optimal initial errors was clarified through eddy-energetics analysis. The results indicated that energy from the background field is transferred to the error field because of barotropic and baroclinic instabilities. Thus, it is inferred that both barotropic and baroclinic processes play important roles in the growth of CNOP-type optimal initial errors.     

An Adjoint-Free CNOP–4DVar Hybrid Method for Identifying Sensitive Areas in Targeted Observations: Method Formulation and Preliminary Evaluation

This paper proposes a hybrid method, called CNOP–4 DVar, for the identification of sensitive areas in targeted observations, which takes the advantages of both the conditional nonlinear optimal perturbation(CNOP) and four-dimensional variational assimilation(4 DVar) methods. The proposed CNOP–4 DVar method is capable of capturing the most sensitive initial perturbation(IP), which causes the greatest perturbation growth at the time of verification; it can also identify sensitive areas by evaluating their assimilation effects for eliminating the most sensitive IP. To alleviate the dependence of the CNOP–4 DVar method on the adjoint model, which is inherited from the adjoint-based approach, we utilized two adjointfree methods, NLS-CNOP and NLS-4 DVar, to solve the CNOP and 4 DVar sub-problems, respectively. A comprehensive performance evaluation for the proposed CNOP–4 DVar method and its comparison with the CNOP and CNOP–ensemble transform Kalman filter(ETKF) methods based on 10 000 observing system simulation experiments on the shallow-water equation model are also provided. The experimental results show that the proposed CNOP–4 DVar method performs better than the CNOP–ETKF method and substantially better than the CNOP method.