(Copernicus Publications on behalf of the European Geosciences Union & the American Geophysical Union, 2017-07-19)
Gebbie, Geoffrey A.; Hsieh, Tsung-Lin
The Lagrange multiplier method for combining observations and models (i.e., the adjoint method or "4D-VAR") has been avoided or approximated when the numerical model is highly nonlinear or chaotic. This approach has been adopted primarily due to difficulties in the initialization of low-dimensional chaotic models, where the search for optimal initial conditions by gradient-descent algorithms is hampered by multiple local minima. Although initialization is an important task for numerical weather prediction, ocean state estimation usually demands an additional task – a solution of the time-dependent surface boundary conditions that result from atmosphere–ocean interaction. Here, we apply the Lagrange multiplier method to an analogous boundary control problem, tracking the trajectory of the forced chaotic pendulum. Contrary to previous assertions, it is demonstrated that the Lagrange multiplier method can track multiple chaotic transitions through time, so long as the boundary conditions render the system controllable. Thus, the nonlinear timescale poses no limit to the time interval for successful Lagrange multiplier-based estimation. That the key criterion is controllability, not a pure measure of dynamical stability or chaos, illustrates the similarities between the Lagrange multiplier method and other state estimation methods. The results with the chaotic pendulum suggest that nonlinearity should not be a fundamental obstacle to ocean state estimation with eddy-resolving models, especially when using an improved first-guess trajectory.
