Initial ensemble generation and validation for ocean data assimilation using HYCOM in the Pacific

A method to initialize an ensemble, introduced by Evensen (Physica, D 77:108–129, 1994a; J Geophys Res 99(C5):10143–10162, 1994b; Ocean Dynamics 53:343–367, 2003), was applied to the Ocean General Circulation Model (OGCM) HYbrid Coordinate Ocean Model (HYCOM) for the Pacific Ocean. Taking advantage of the hybrid coordinates, an initial ensemble is created by first perturbing the layer interfaces and then running the model for a spin-up period of 1 month forced by randomly perturbed atmospheric forcing fields. In addition to the perturbations of layer interfaces, we implemented perturbations of the mixed layer temperatures. In this paper, we investigate the quality of the initial ensemble generated by this scheme and the influence of the horizontal decorrelation scale and vertical correlation on the statistics of the resulting ensemble. We performed six ensemble generation experiments with different combinations of horizontal decorrelation scales and with/without perturbations in the mixed layer. The resulting six sets of initial ensembles are then analyzed in terms of sustainability of the ensemble spread and realism of the correlation patterns. The ensemble spreads are validated against the difference between model and observations after 20 years of free run. The correlation patterns of six sets of ensemble are compared to each other. This study shows that the ensemble generation scheme can effectively generate an initial ensemble whose spread is consistent with the observed errors. The correlation pattern of the ensemble also exhibits realistic features. The addition of mixed layer perturbations improves both the spread and correlation. Some limitations of the ensemble generation scheme are also discussed. We found that the vertical shift of isopycnal coordinates provokes unrealistically large deviations in shallow layers near the islands of the West Pacific. A simple correction circumvents the problem.

[1]  G. Evensen,et al.  Analysis Scheme in the Ensemble Kalman Filter , 1998 .

[2]  Ocean data assimilation with background error covariance derived from OGCM outputs , 2004 .

[3]  M. Buehner,et al.  Atmospheric Data Assimilation with an Ensemble Kalman Filter: Results with Real Observations , 2005 .

[4]  Christian L. Keppenne,et al.  Assimilation of temperature into an isopycnal ocean general circulation model using a parallel ensemble Kalman filter , 2003 .

[5]  S. Ricci,et al.  Ensembles of global ocean analyses for seasonal climate prediction: impact of temperature assimilation , 2005 .

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

[7]  D. Legates,et al.  Mean seasonal and spatial variability in gauge‐corrected, global precipitation , 1990 .

[8]  Arlindo da Silva,et al.  Data assimilation in the presence of forecast bias , 1998 .

[9]  Andrew P. Morse,et al.  DEVELOPMENT OF A EUROPEAN MULTIMODEL ENSEMBLE SYSTEM FOR SEASONAL-TO-INTERANNUAL PREDICTION (DEMETER) , 2004 .

[10]  W. Teague,et al.  A Comparison Between the Generalized Digital Environmental Model and Levitus climatologies , 1990 .

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

[12]  Geir Evensen,et al.  Assimilation of SLA and SST data into an OGCM for the Indian Ocean , 2002 .

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

[14]  Geir Evensen,et al.  Coordinate Transformation on a Sphere Using Conformal Mapping , 1999 .

[15]  J. Whitaker,et al.  Ensemble Data Assimilation without Perturbed Observations , 2002 .

[16]  G. Evensen,et al.  Assimilation of ocean colour data into a biochemical model of the North Atlantic: Part 1. Data assimilation experiments , 2003 .

[17]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[18]  M. Rienecker,et al.  Initial testing of a massively parallel ensemble Kalman filter with the Poseidon isopycnal ocean general circulation model , 2002 .

[19]  E. Kalnay,et al.  C ○ 2007 The Authors , 2006 .

[20]  J. Whitaker,et al.  Accounting for the Error due to Unresolved Scales in Ensemble Data Assimilation: A Comparison of Different Approaches , 2005 .

[21]  Olwijn Leeuwenburgh,et al.  Assimilation of along‐track altimeter data in the tropical Pacific region of a global OGCM ensemble , 2005 .

[22]  Olwijn Leeuwenburgh,et al.  Validation of an EnKF System for OGCM Initialization Assimilating Temperature, Salinity, and Surface Height Measurements , 2007 .

[23]  A. O'Neill.  Atmospheric Data Assimilation , 2000 .

[24]  G. Evensen,et al.  Sequential Data Assimilation Techniques in Oceanography , 2003 .

[25]  D. Legates,et al.  Mean seasonal and spatial variability in global surface air temperature , 1990 .

[26]  Frederic Vitart,et al.  An Ensemble Generation Method for Seasonal Forecasting with an Ocean-Atmosphere Coupled Model , 2005 .

[27]  Thomas M. Smith,et al.  An Improved In Situ and Satellite SST Analysis for Climate , 2002 .

[28]  Rainer Bleck,et al.  A Wind-Driven Isopycnic Coordinate Model of the North and Equatorial Atlantic Ocean 1 , 1990 .

[29]  O. Alves,et al.  Tropical Pacific Ocean model error covariances from Monte Carlo simulations , 2005 .

[30]  J. Whitaker,et al.  Distance-dependent filtering of background error covariance estimates in an ensemble Kalman filter , 2001 .

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

[32]  K. Brusdala,et al.  A demonstration of ensemble-based assimilation methods with a layered OGCM from the perspective of operational ocean forecasting systems , 2003 .

[33]  Geir Evensen,et al.  Assimilation of ice concentration in a coupled ice–ocean model, using the Ensemble Kalman filter , 2003 .

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

[35]  P. Houtekamer,et al.  A Sequential Ensemble Kalman Filter for Atmospheric Data Assimilation , 2001 .

[36]  H. Hurlburt,et al.  Black Sea Mixed Layer Sensitivity to Various Wind and Thermal Forcing Products on Climatological Time Scales , 2005 .

[37]  Rainer Bleck,et al.  An oceanic general circulation model framed in hybrid isopycnic-Cartesian coordinates , 2002 .

[38]  D. Pham.  Stochastic Methods for Sequential Data Assimilation in Strongly Nonlinear Systems , 2001 .