Data assimilation and inverse methods in terms of a probabilistic formulation

Abstract The weak constraint inverse for nonlinear dynamical models is discussed and derived in term of a probabilistic formulation. The well-known result that for Gaussian error statistics the minimum of the weak constraint inverse is equal to the maximum-likelihood estimate is rederived. Then several methods based on ensemble statistics that can be used to find the smoother (as opposed to the filter) solution are introduced and compared to traditional methods. A strong point of the new methods is that they avoid the integration of adjoint equations, which is a complex task for real oceanographic or atmospheric applications. They also avoid iterative searches in a Hilbert space, and error estimates can be obtained without much additional computational effort. The feasibility of the new methods is illustrated in a two-layer quasigeostrophic ocean model.

[1]  A. Bennett,et al.  Open Ocean Modeling as an Inverse Problem: Tidal Theory , 1982 .

[2]  A. Tarantola Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[3]  A. Lorenc Optimal nonlinear objective analysis , 1988 .

[4]  W. Budgell,et al.  The Kalman smoother for a linear quasi-geostrophic model of ocean circulation , 1989 .

[5]  M. Ikeda,et al.  Eddy-Current Interactions using a Two-Layer Quasi-Geostrophic Model , 1989 .

[6]  Johnny A. Johannessen,et al.  A Process Study of Mesoscale Meanders and Eddies in the Norwegian Coastal Current , 1989 .

[7]  G. Evensen,et al.  Modeled and observed mesoscale circulation and wave‐current refraction during the 1988 Norwegian Continental Shelf Experiment , 1991 .

[8]  G. Evensen Using the Extended Kalman Filter with a Multilayer Quasi-Geostrophic Ocean Model , 1992 .

[9]  A. Bennett Inverse Methods in Physical Oceanography , 1992 .

[10]  A. Bennett,et al.  The generalized inverse of a nonlinear quasigeostrophic ocean circulation model , 1992 .

[11]  Lance M. Leslie,et al.  Tropical Cyclone Prediction Using a Barotropic Model Initialized by a Generalized Inverse Method , 1993 .

[12]  Geir Evensen,et al.  Open Boundary Conditions for the Extended Kalman Filter With a Quasi-Geostrophic Ocean Model , 1993 .

[13]  P. Courtier,et al.  Assimilation of Simulated Wind Lidar Data with a Kalman Filter , 1993 .

[14]  G. Evensen Inverse methods and data assimilation in nonlinear ocean models , 1994 .

[15]  A. Bennett,et al.  TOPEX/POSEIDON tides estimated using a global inverse model , 1994 .

[16]  G. Evensen Sequential data assimilation with a nonlinear quasi‐geostrophic model using Monte Carlo methods to forecast error statistics , 1994 .

[17]  Michael Ghil,et al.  Advanced data assimilation in strongly nonlinear dynamical systems , 1994 .

[18]  Monique Tanguay,et al.  Four‐dimensional data assimilation with a wide range of scales , 1995 .

[19]  G. Evensen,et al.  Assimilation of Geosat altimeter data for the Agulhas current using the ensemble Kalman filter with , 1996 .