An ensemble conditional nonlinear optimal perturbation approach: Formulation and applications to parameter calibration

Abstract

The conditional nonlinear optimal perturbation (CNOP) method proposed by Mu et al. (2003) has been a useful tool for studying predictability dynamics. Its further applications are, however, hampered by the need of an adjoint model. Aiming to solve this problem, an ensemble conditional nonlinear optimal perturbation (EnCNOP) approach is proposed in this paper by merging the Monte Carlo method and the proper orthogonal decomposition (POD) technique into the CNOP to transform an implicit optimization problem into an explicit one. This approach is further formulated for parameter calibration. Numerical experiments with a 1-D model of the soil water equation show that the EnCNOP approach outperforms a dual-pass optimization framework on the basis of the ensemble Kalman filter (EnKF) and a shuffled complex evolution approach (SCE-UA) in terms of both increasing the calibration precision and reducing the computational costs. In particular, the EnCNOP approach outperforms the SCE-UA when only an approximate value range of the input parameters is known. Our EnCNOP method is further implemented within a complicated highly nonlinear land surface model (namely, NCAR CLM3) and evaluated through a case study which illustrates its applicability to complex, highly nonlinear problems. Copyright 2010 by the American Geophysical Union.