A New Scheme for Capturing Global Conditional Nonlinear Optimal Perturbation

Abstract

Conditional nonlinear optimal perturbation (CNOP) represents the initial perturbation that satisfies a certain physical constraint condition, and leads to a maximum prediction error at the moment of prediction. The CNOP method is a useful tool in studying atmosphere and ocean predictability problems. Generally, the optimization algorithm based on the gradient of the cost function to compute CNOP requires an initial guess. The traditional scheme randomly chooses the initial guess of CNOP within the constraint range and therefore this scheme is called RIG-CNOP. However, the RIG-CNOP scheme reduces the probability of capturing the global CNOP in many cases, such as the prediction model is strongly nonlinear or long-term prediction is performed, or multiple extreme values existed in the cost function. Considering the limitations of the RIG-CNOP scheme, we propose a new initial guess selection scheme. In this scheme, we first pre-analyze a series of random initial guesses, and then, an optimal initial guess is selected. The above process replaces the initial guess selection scheme in the traditional scheme, which is called PAIG-CNOP. Numerical experiments are conducted utilizing the Lorenz-63 model. Also, to compare the performance of the PAIG-CNOP method with the RIG-CNOP method in capturing global CNOP, the CNOP and the maximum cost function value (MCFV) obtained by the filtering method (FM) are used as benchmarks (this value is called FMMCFV in brief). The experimental results show that even the prediction model is strongly nonlinear or the prediction time is long, or the cost function has multiple extreme values, the PAIG-CNOP method can capture the global CNOP with a high probability. The results show that the PAIG-CNOP method has a higher probability of capturing the global CNOP than the RIG-CNOP method. In addition, we use an ensemble-based technique in the computation of gradients, thus avoiding the use of adjoint techniques in the maximization process. Due to the attractive features of the new method, the PAIG-CNOP method is an efficient and useful method for solving CNOP, it can be more easily applied to obtain the global CNOP of operational prediction models.