A Variant Constrained Genetic Algorithm for Solving Conditional Nonlinear Optimal Perturbations简

A variant constrained genetic algorithm （VCGA） for effective tracking of conditional nonlinear optimal perturbations （CNOPs） is presented. Compared with traditional constraint handling methods, the treatment of the constraint condition in VCGA is relatively easy to implement. Moreover, it does not require adjustments to indefinite pararneters. Using a hybrid crossover operator and the newly developed multi-ply mutation operator, VCGA improves the performance of GAs. To demonstrate the capability of VCGA to catch CNOPS in non-smooth cases, a partial differential equation, which has ＂on off＂ switches in its forcing term, is employed as the nonlinear model. To search global CNOPs of the nonlinear model, numerical experiments using VCGA, the traditional gradient descent algorithm based on the adjoint method （ADJ）, and a GA using tournament selection operation and the niching technique （GA-DEB） were performed. The results with various initial reference states showed that, in smooth cases, all three optimization methods are able to catch global CNOPs. Nevertheless, in non-smooth situations, a large proportion of CNOPs captured by the ADJ are local. Compared with ADJ, the performance of GA-DEB shows considerable improvement, but it is far below VCGA. Further, the impacts of population sizes on both VCGA and GA-DEB were investigated. The results were used to estimate the computation time of ~CGA and GA-DEB in obtaining CNOPs. The computational costs for VCGA, GA-DEB and ADJ to catch CNOPs of the nonlinear model are also compared.     

Impact of Different Guidances on Sensitive Areas of Targeting Observations Based on the CNOP Method

The conditional nonlinear optimal perturbations(CNOPs) obtained by a fast algorithm are applied to determining the sensitive area for the targeting observation of Typhoon Matsa in 2005 using an operational regional prediction model-the Global/Regional Assimilation and PrEdiction System(GRAPES).Through a series of sensitivity experiments,several issues on targeting strategy design are discussed,including the effectivity of different guidances to determine the sensitive area(or targeting area) and the impa...     

Application of the Conditional Nonlinear Optimal Perturbation Method to the Predictability Study of the Kuroshio Large Meander

A reduced-gravity barotropic shallow-water model was used to simulate the Kuroshio path variations.The results show that the model was able to capture the essential features of these path variations.We used one simulation of the model as the reference state and investigated the effects of errors in model parameters on the prediction of the transition to the Kuroshio large meander (KLM) state using the conditional nonlinear optimal parameter perturbation (CNOP-P) method.Because of their relatively large uncertainties,three model parameters were considered:the interfacial friction coefficient,the wind-stress amplitude,and the lateral friction coefficient.We determined the CNOP-Ps optimized for each of these three parameters independently,and we optimized all three parameters simultaneously using the Spectral Projected Gradient 2 (SPG2) algorithm.Similarly,the impacts caused by errors in initial conditions were examined using the conditional nonlinear optimal initial perturbation (CNOP-I) method.Both the CNOP-I and CNOP-Ps can result in significant prediction errors of the KLM over a lead time of 240 days.But the prediction error caused by CNOP-I is greater than that caused by CNOP-P.The results of this study indicate not only that initial condition errors have greater effects on the prediction of the KLM than errors in model parameters but also that the latter cannot be ignored.Hence,to enhance the forecast skill of the KLM in this model,the initial conditions should first be improved,the model parameters should use the best possible estimates.     

A Preliminary Application of the Differential Evolution Algorithm to Calculate the CNOP

A projected skill is adopted by use of the differential evolution （DE） algorithm to calculate a conditional nonlinear optimal perturbation （CNOP）. The CNOP is the maximal value of a constrained optimization problem with a constraint condition, such as a ball constraint. The success of the DE algorithm lies in its ability to handle a non-differentiable and nonlinear cost function. In this study, the DE algorithm and the traditional optimization algorithms used to obtain the CNOPs are compared by analyzing a theoretical grassland ecosystem model and a dynamic global vegetation model. This study shows that the CNOPs generated by the DE algorithm are similar to those by the sequential quadratic programming （SQP） algorithm and the spectral projected gradients （SPG2） algorithm. If the cost function is non-differentiable, the CNOPs could also be caught with the DE algorithm. The numerical results suggest the DE algorithm can be employed to calculate the CNOP, especially when the cost function is non-differentiable.     

遗传算法求解最大预报误差的可行性研究

利用条件非线性最优扰动（conditional nonlinear optimal perturbation,CNOP）可以实现最大预报误差的上界估计。CNOP通常由基于梯度信息的约束优化算法进行求解,且其中的梯度信息由伴随模式提供。然而当非线性模式中含不连续＂开关＂时,传统伴随方法不能为优化过程提供正确的梯度方向,从而导致优化失败。为此,采用自适应变异和混合交叉的遗传算法,联赛选择机制和小生境技术的约束处理方法来求解最大预报误差上界。为检验新方法的有效性,以修改的Lorenz模型作为预报模式,对3个初始态分别用新方法和传统伴随方法进行比较,数值试验结果显示新方法求解出的最大预报误差的上界更加精确。     

The Roles of Spatial Locations and Patterns of Initial Errors in the Uncertainties of Tropical Cyclone Forecasts

In this study,a series of sensitivity experiments were performed for two tropical cyclones (TCs),TC Longwang (2005) and TC Sinlaku (2008),to explore the roles of locations and patterns of initial errors in uncertainties of TC forecasts.Specifically,three types of initial errors were generated and three types of sensitive areas were determined using conditional nonlinear optimal perturbation (CNOP),first singular vector (FSV),and composite singular vector (CSV) methods.Additionally,random initial errors in randomly selected areas were considered.Based on these four types of initial errors and areas,we designed and performed 16 experiments to investigate the impacts of locations and patterns of initial errors on the nonlinear developments of the errors,and to determine which type of initial errors and areas has the greatest impact on TC forecasts.Overall,results from the experiments indicate the following:(1) The impact of random errors introduced into the sensitive areas was greater than that of errors themselves fixed in the randomly selected areas.From the perspective of statistical analysis,and by comparison,the impact of random errors introduced into the CNOP target area was greatest.(2) The initial errors with CNOP,CSV,or FSV patterns were likely to grow faster than random errors.(3) The initial errors with CNOP patterns in the CNOP target areas had the greatest impacts on the final verification forecasts.     

Effect of Stochastic MJO Forcing on ENSO Predictability

Within the frame of the Zebiak-Cane model,the impact of the uncertainties of the Madden-Julian Oscillation(MJO) on ENSO predictability was studied using a parameterized stochastic representation of intraseasonal forcing.The results show that the uncertainties of MJO have little effect on the maximum prediction error for ENSO events caused by conditional nonlinear optimal perturbation(CNOP);compared to CNOP-type initial error,the model error caused by the uncertainties of MJO led to a smaller prediction uncertainty of ENSO,and its influence over the ENSO predictability was not significant.This result suggests that the initial error might be the main error source that produces uncertainty in ENSO prediction,which could provide a theoretical foundation for the data assimilation of the ENSO forecast.     

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.     

Using the Lund-Potsdam-Jena Model to Understand the Different Responses of Three Woody Plants to Land Use in China

In this study, the approach of conditional nonlinear optimal perturbation related to initial perturbation (CNOP-I) was employed to investigate the maximum variations in plant amount for three main woody plants (a temperate broadleaved evergreen, a temperate broadleaved summergreen, and a boreal needleleaved evergreen) in China. The investigation was conducted within a certain range of land use intensity using a state-of-the-art Lund-Potsdam-Jena dynamic global vegetation model (LPJ DGVM). CNOP-I represents a class of deforestation and can be considered a type of land use with respect to the initial perturbation. When deforestation denoted by the CNOP-I has the same intensity for all three plants, the variation in plant amount of the boreal needleleaved evergreen in northern China is greater than the variation in plant amount of both the temperate broadleaved evergreen and temperate broadleaved summergreen in southern China. As deforestation intensity increases, the plant amount variation in the three woody plant functional types carbon changes, in a nonlinear fashion. The impact of land use on plant functional types is minor because the interaction between climate condition and land use is not considered in the LPJ model. Finally, the different impacts of deforestation on net primary production of the three plant functional types were analyzed by modeling gross primary production and autotrophic respiration. Our results suggest that the CNOP-I approach is a useful tool for exploring the nonlinear and different responses of terrestrial ecosystems to land use.     

二阶精度半隐式半拉格朗日轨迹计算和时间积分方案在GRAPES区域模式中的应用

将两时间层稳定外插方案(Stable Extrapolation Two-Time-Level Scheme,SETTLS)引入GRAPES区域模式,并将其用于上游点和非线性项的时间外插计算。对线性项采用二阶精度的非中央权重时间平均,并取等温参考大气的温度大于实际大气平均温度,以保证半隐式积分方案的稳定性。原GRAPES时间积分方案对线性项做一阶非中央权重时间平均,对参考大气的选择并无限制,而为保证稳定性,须取较大的非中央权重系数,但非中央权重系数会对低波数波动产生较强的衰减作用。理想试验结果表明,相比原GRAPES半隐式半拉格朗日(SISL)时间积分方案,新SISL时间积分方案计算稳定且对波动的衰减作用较弱。