An adjoint method for the assimilation of statistical characteristics into eddy-resolving ocean models

The study investigates perspectives of the parameter estimation problem with the adjoint method in eddy-resolving models. Sensitivity to initial conditions resulting from the chaotic nature of this type of model limits the direct application of the adjoint method by predictability. Prolonging the period of assimilation is accompanied by the appearance of an increasing number of secondary minima of the cost function that prevents the convergence of this method. In the framework of the Lorenz model it is shown that averaged quantities are suitable for describing invariant properties, and that secondary minima are for this type of data transformed into stochastic deviations. An adjoint method suitable for the assimilation of statistical characteristics of data and applicable on time scales beyond the predictability limit is presented. The approach assumes a greater predictability for averaged quantities. The adjoint to a prognostic model for statistical moments is employed for calculating cost function gradients that ignore the fine structure resulting from secondary minima. Coarse resolution versions of eddy-resolving models are used for this purpose. Identical twin experiments are performed with a quasigeostrophic model to evaluate the performance and limitations of this approach in improving models by estimating parameters. The wind stress curl is estimated from a simulated mean stream function. A very simple parameterization scheme for the assimilation of second-order moments is shown to permit the estimation of gradients that perform efficiently in minimizing cost functions.

[1]  T. N. Palmer,et al.  Predictability of the Atmosphere and Oceans: From Days to Decades , 1996 .

[2]  William H. Press,et al.  Numerical Recipes in Fortran 77 , 1992 .

[3]  J. Marotzke,et al.  Finding the steady state of a general circulation model through data assimilation: Application to the North Atlantic Ocean , 1993 .

[4]  Thomas Kaminski,et al.  Sensitivity of the seasonal cycle of CO2 at remote monitoring stations with respect to seasonal surface exchange fluxes determined with the adjoint of an atmospheric transport model , 1996 .

[5]  P. Courtier,et al.  Extended assimilation and forecast experiments with a four‐dimensional variational assimilation system , 1998 .

[6]  Carl Wunsch,et al.  The global ocean circulation estimated from TOPEX/POSEIDON altimetry and the MIT general circulation model , 1997 .

[7]  M. Spall,et al.  Specification of eddy transfer coefficients in coarse resolution ocean circulation models , 1997 .

[8]  P. Courtier,et al.  A strategy for operational implementation of 4D‐Var, using an incremental approach , 1994 .

[9]  Chou Jifan,et al.  Predictability of the atmosphere , 1989 .

[10]  W. C. Thacker,et al.  Relationships between Statistical and Deterministic Methods of Data Assimilation , 1986 .

[11]  Olivier Talagrand,et al.  On extending the limits of variational assimilation in nonlinear chaotic systems , 1996 .

[12]  William Carlisle Thacker,et al.  The role of the Hessian matrix in fitting models to measurements , 1989 .

[13]  William Carlisle Thacker,et al.  Fitting models to inadequate data by enforcing spatial and temporal smoothness , 1988 .

[14]  J. Eckmann.  Roads to turbulence in dissipative dynamical systems , 1981 .

[15]  A. Sy.  Investigation of large-scale circulation patterns in the central North Atlantic: the North Atlantic current, the Azores current, and the Mediterranean Water plume in the area of the Mid-Atlantic Ridge , 1988 .

[16]  D. Stammer.  Global Characteristics of Ocean Variability Estimated from Regional TOPEX/POSEIDON Altimeter Measurements , 1997 .

[17]  Robert N. Miller,et al.  Data assimilation into nonlinear stochastic models , 1999 .

[18]  Andreas Oschlies,et al.  Eddy-induced enhancement of primary production in a model of the North Atlantic Ocean , 1998, Nature.

[19]  M. Allen,et al.  Sensitivity analysis of the climate of a chaotic system , 2000 .

[20]  E. Tziperman,et al.  Finite Difference of Adjoint or Adjoint of Finite Difference , 1997 .

[21]  Jens Schröter,et al.  Variational Assimilation of Geosat Data into an Eddy-resolving Model of the Gulf Stream Extension Area , 1993 .

[22]  M. Hendershott,et al.  Stochastic closure for nonlinear Rossby waves , 1977, Journal of Fluid Mechanics.

[23]  J. Green,et al.  Transfer properties of the large‐scale eddies and the general circulation of the atmosphere , 1970 .

[24]  Andreas Schiller,et al.  The mean circulation of the Atlantic Ocean north of 30s determined with the adjoint method applied to an ocean general circulation model , 1995 .

[25]  H. Risken.  Fokker-Planck Equation , 1984 .

[26]  W. Thacker,et al.  An Optimal-Control/Adjoint-Equations Approach to Studying the Oceanic General Circulation , 1989 .

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

[28]  Jens Schröter,et al.  Assimilation of satellite altimeter data into an open ocean model , 1995 .

[29]  Isaac M. Held,et al.  Parameterization of Quasigeostrophic Eddies in Primitive Equation Ocean Models. , 1997 .

[30]  Andrew M. Moore,et al.  Data Assimilation in a Quasi-geostrophic Open-Ocean Model of the Gulf Stream Region Using the Adjoint Method , 1991 .

[31]  David J. Stensrud,et al.  Behaviors of Variational and Nudging Assimilation Techniques with a Chaotic Low-Order Model , 1992 .

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

[33]  Joseph A. C. Delaney.  Sensitivity analysis , 2018, The African Continental Free Trade Area: Economic and Distributional Effects.

[34]  E. Lorenz.  Deterministic nonperiodic flow , 1963 .

[35]  T. Palmer.  Extended-range atmospheric prediction and the Lorenz model , 1993 .

[36]  Jon M. Nese,et al.  Calculated Attractor Dimensions for Low-Order Spectral Models , 1987 .

[37]  J. Willebrand,et al.  Assimilation of altimetric data and mean sea surface height into an eddy-permitting model of the North Atlantic , 2001 .

[38]  William H. Press,et al.  Numerical recipes , 1990 .

[39]  Andreas Oschlies,et al.  Assimilation of Geosat altimeter data into an eddy-resolving primitive equation model of the North Atlantic Ocean , 1996 .

[40]  F. L. Dimet,et al.  Variational algorithms for analysis and assimilation of meteorological observations: theoretical aspects , 1986 .

[41]  Sol Hellerman,et al.  Normal Monthly Wind Stress Over the World Ocean with Error Estimates , 1983 .

[42]  Yong Li.  A note on the uniqueness problem of variational adjustment approach to four-dimensional data assimilation , 1991 .

[43]  Stephen M. Griffies,et al.  Predictability of North Atlantic Multidecadal Climate Variability , 1997, Science.

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

[45]  V. I. Oseledec.  A multiplicative ergodic theorem: Lyapunov characteristic num-bers for dynamical systems , 1968 .

[46]  Rosemary Morrow,et al.  Adjoint assimilation of altimetric, surface drifter, and hydrographic data in a quasi‐geostrophic model of the Azores Current , 1995 .

[47]  M. Farge.  Wavelet Transforms and their Applications to Turbulence , 1992 .

[48]  Guido Boffetta,et al.  An Extension of the Lyapunov Analysis for the Predictability Problem , 1998, chao-dyn/9801030.

[49]  William R. Holland,et al.  The Role of Mesoscale Eddies in the General Circulation of the Ocean—Numerical Experiments Using a Wind-Driven Quasi-Geostrophic Model , 1978 .

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

[51]  Peter H. Stone,et al.  A Simplified Radiative-Dynamical Model for the Static Stability of Rotating Atmospheres , 1972 .

[52]  B. Sonnerup,et al.  Vortex laws and field line invariants in polytropic field-aligned MHD flow , 1994 .

[53]  Robert Vautard,et al.  Four-dimensional variational assimilation and predictability in a quasi-geostrophic model , 1998 .

[54]  Pierre Gauthier,et al.  Chaos and quadri-dimensional data assimilation: a study based on the Lorenz model , 1992 .

[55]  Inverse Modeling of Seasonal Variations in the North Atlantic Ocean , 1998 .

[56]  E. Ruprecht,et al.  Determination of cloud liquid water path over the oceans from Special Sensor Microwave/Imager (SSM/I) data using neural networks , 1998 .

[57]  William H. Press,et al.  Numerical recipes in C , 2002 .

[58]  S. Meacham.  Low-frequency variability in the wind-driven circulation , 2000 .