Due to uncertainties in initial conditions and parameters, the stability and uncertainty of grassland ecosystem simulations using ecosystem models are issues of concern. Our objective is to determine the types and patterns of initial and parameter perturbations that yield the greatest instability and uncertainty in simulated grassland ecosystems using theoretical models. We used a nonlinear optimization approach, i.e., a conditional nonlinear optimal perturbation related to initial and parameter perturbations （CNOP） approach, in our work. Numerical results indicated that the CNOP showed a special and nonlinear optimal pattern when the initial state variables and multiple parameters were considered simultaneously. A visibly different complex optimal pattern characterizing the CNOPs was obtained by choosing different combinations of initial state variables and multiple parameters in different physical processes. We propose that the grassland modeled ecosystem caused by the CNOP-type perturbation is unstable and exhibits two aspects： abrupt change and the time needed for the abrupt change from a grassland equilibrium state to a desert equilibrium state when the initial state variables and multiple parameters are considered simultaneously. We compared these findings with results affected by the CNOPs obtained by considering only uncertainties in initial state variables and in a single parameter. The numerical results imply that the nonlinear optimal pattern of initial perturbations and parameter perturbations, especially for more parameters or when special parameters are involved, plays a key role in determining stabilities and uncertainties associated with a simulated or predicted grassland ecosystem.
The response of a grassland ecosystem to climate change is discussed within the context of a theoretical model.An optimization approach,a conditional nonlinear optimal perturbation related to parameter（CNOP-P） approach,was employed in this study.The CNOP-P,a perturbation of moisture index in the theoretical model,represents a nonlinear climate perturbation.Two kinds of linear climate perturbations were also used to study the response of the grassland ecosystem to different types of climate changes.The results show that the extent of grassland ecosystem variation caused by the CNOP-P-type climate change is greater than that caused by the two linear types of climate change.In addition,the grassland ecosystem affected by the CNOP-P-type climate change evolved into a desert ecosystem,and the two linear types of climate changes failed within a specific amplitude range when the moisture index recovered to its reference state.Therefore,the grassland ecosystem response to climate change was nonlinear.This study yielded similar results for a desert ecosystem seeded with both living and wilted biomass litter.The quantitative analysis performed in this study also accounted for the role of soil moisture in the root zone and the shading effect of wilted biomass on the grassland ecosystem through nonlinear interactions between soil and vegetation.The results of this study imply that the CNOP-P approach is a potentially effective tool for assessing the impact of nonlinear climate change on grassland ecosystems.
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.
