Multivariate assimilation in MERCATOR project: New statistical parameters from forecast error estimation

Abstract The new operational prototype of Mercator (french Global Ocean Data Assimilation Experiment contribution) is composed of a North Atlantic primitive equation ocean model OPA (Ocean Parallel Algorithm between 20°S and 70°N, [Madec, G., P. Delecluse, M. Imbard and C. Levy (1998). OPA8.1 ocean general circulation model reference manuel. Notes du pole de modelisation IPSL. n°11: 91p]) and of a multivariate and multidata assimilation scheme [De Mey, P. and M. Benkiran (2002). “A multivariate reduced-order optimal interpolation method and its application in Mediterranean basin-scale circulation.” Ocean Forecasting : Conceptual basis and application, Pinardi, N., Springer Verlag.] This system has already given some significant improvements from previous Mercator configurations (M. Benkiran, personal communication). However some biases on ocean state still remain in the tropics where the reduced-order optimal interpolation scheme is suspected to be ill-parameted in the model forecast error. Indeed the guess error covariance matrix is decomposed into an error variance value and a spatio-temporal correlation function which are assumed to have some “good” properties (spatial homogeneity of the correlation function, constant ratio between signal and error variance). This study shows how we can use ensemble methods to validate these assumptions. We can see that the correlation function can reach negative values locally, mostly in regions of high variability contradictory with the homogeneous hypothesis. The reduced space used in the operational configuration is based on the signal seasonal Empirical Orthogonal Functions (EOFs). An empirical relationship between signal and error variance has been set and the correlation function is the same on every dimension of the reduced space. By projection of the estimated guess error variance onto the reduced space, we find a repartition of this quantity quite different to what was set in the system. The error statistics is found to be inhomogeneous compared to hypothesis made in the assimilation scheme. These two new parameters tested separately in the assimilation scheme gives significant improvements of the forecast and analysis results. This is particularly obvious in the tropics. But relationship between signal and error statistics (as assumed in the optimal interpolation) is found to be complex.

[1]  M. Gavart,et al.  Isopycnal EOFs in the Azores Current Region: A Statistical Tool forDynamical Analysis and Data Assimilation , 1997 .

[2]  Franco Molteni,et al.  Ensemble prediction using dynamically conditioned perturbations , 1993 .

[3]  Mounir Benkiran,et al.  A Multivariate Reduced-order Optimal Interpolation Method and its Application to the Mediterranean Basin-scale Circulation , 2002 .

[4]  G. Caniaux,et al.  Assimilation of satellite altimeter data in a primitive-equation model of the Azores–Madeira region , 1999 .

[5]  Eric Dombrowsky,et al.  Estimation of the optimal interpolation parameters in a quasi-geostrophic model of the Northeast Atlantic using ensemble methods , 2003 .

[6]  M. Arhan,et al.  Dynamics of Eddy Motions in the Eastern North Atlantic , 1985 .

[7]  Eugenia Kalnay,et al.  Ensemble Forecasting at NMC: The Generation of Perturbations , 1993 .

[8]  J. Carton,et al.  Long waves and eddies in the tropical Atlantic Ocean: 1984–1990 , 1991 .

[9]  James M. Murphy,et al.  The impact of ensemble forecasts on predictability , 1988 .

[10]  Analysis of Errors Due to Polynomial Adjustment of Altimeter Profiles , 1991 .

[11]  Nadia Pinardi,et al.  Ocean forecasting : conceptual basis and applications , 2002 .

[12]  A. Moore.  The dynamics of error growth and predictability in a model of the Gulf stream. Part II : Ensemble prediction , 1999 .

[13]  M. Kamachi,et al.  Global statistical space-time scales of oceanic variability estimated from the TOPEX/ POSEIDON altimeter data , 2000 .

[14]  P. Delecluse,et al.  OPA 8.1 Ocean General Circulation Model reference manual , 1998 .

[15]  Michael Ghil,et al.  Dynamical Properties of Error Statistics in a Shallow-Water Model , 1993 .

[16]  B. Barnier,et al.  A new analysis of hydrographic data in the Atlantic and its application to an inverse modeling study , 1998 .

[17]  P. Poulain,et al.  Measurements of the water-following capability of holey-sock and TRISTAR drifters , 1995 .

[18]  Paul R. Julian,et al.  On Some Properties of Correlation Functions Used in Optimum Interpolation Schemes , 1975 .

[19]  P. Poulain,et al.  Quality Control and Interpolations of WOCE-TOGA Drifter Data , 1996 .

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

[21]  R. E. Kalman,et al.  A New Approach to Linear Filtering and Prediction Problems , 2002 .

[22]  Keith Haines,et al.  Altimetric assimilation with water property conservation , 1996 .

[23]  C. Maes,et al.  Sea level comparison between TOPEX/POSEIDON altimetric data and a global ocean general circulation model from an assimilation perspective , 1999 .