In this paper, a novel and simple statistical prediction model for sea surface multivariate is developed based on extended empirical orthogonal functions (referred to as the MEEOF model). This simple model embeds the temporal evolution information into the empirical orthogonal function spatial patterns, which effectively captures the spatial distribution of the sea surface variables and their evolution process over time, and can be used to improve the accuracy of intra-seasonal ocean forecasts. At the same time, it considers both the correlation between different spatial and temporal points and the dynamic balance between different sea surface variables. The performance of the MEEOF prediction model has been examined in the South China Sea (SCS) based on remote sensing satellite datasets. Experimental results demonstrate that this model has significant forecasting capability within the forecast window of 15–90 days, compared with the traditional persistence forecasts and optimal climatic normal (OCN) results. Furthermore, this model exhibits good forecast performance throughout the entire forecast window; the final prediction model (referred to as the fusion model) is established by obtaining the weighted average for the MEEOF forecasts and persistence forecast results. Numerical experimental results show that this fusion prediction model consistently outperforms the persistence model, the OCN model, and the linear regression model over the entire forecast window. A case study shows that the propagation of ocean waves and the coordination between different sea surface variables can be well predicted by this simple model, indicating that this fusion model has a potential advantage in intra-seasonal predictions of the ocean.

The joint state and parameter estimation problem is an important issue in data assimilation. An adjoint free data assimilation method, namely analytical four-dimensional ensemble variational (A-4DEnVar) data assimilation method, was developed to provide a solution for the joint estimation problem. In the algorithm, to estimate the adjoint model reasonably, the ensemble initial conditions and parameters are generated by Gaussian noise whose covariance is constructed by multiplying a very small factor by their background error covariance. The ensemble perturbations are calculated with respect to background states rather than the ensemble mean. Next, the usage of temporal cross covariances makes it possible to avoid the adjoint model and estimate the gradient in 4DVar. Furthermore, we update the solution iteratively with a linear search process to improve the stability and ensure the convergence of the algorithm. The method is tested using the three-variable Lorenz model (Lorenz-1963) to illustrate its efficiency. It is shown that A-4DEnVar results in similar performance with 4DVar. Sensitivity experiments show that A-4DEnVar is able to assimilate observations successfully with different settings. The proposed method is able to work as well as 4DVar and avoid adjoint models for the joint state and parameter estimation.

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.