Implementation of a multivariate data assimilation scheme for isopycnic coordinate ocean models: Application to a 1993–1996 hindcast of the North Atlantic Ocean circulation

[1] A 1993–1996 hindcast of the North Atlantic circulation is presented and analyzed. Multivariate satellite observations (sea surface height and sea surface temperature) are assimilated into an isopycnic coordinate ocean model (Miami Isopycnic Coordinate Ocean Model [MICOM]) using a reduced order Kalman filter (the SEEK filter). The objective is to validate the assimilation system in hindcast mode before running a near real-time exercise as part of the European DIADEM project. The analysis scheme is upgraded to make the statistical parameterization as realistic as possible, and this is shown to be particularly important as soon as a multivariate data set is used. Second, the paper describes the modifications to the SEEK filter that are necessary to meet the requirements of a layered model and to benefit from the advantages of the isopycnic coordinate. In particular, the need for an initialization procedure is explained, and this dynamical adjustment is shown to be more natural and robust in the context of the isopycnic coordinate. The hindcast experiment is validated by comparisons with observations, by studying innovation statistics, and by comparing the results with independent observations (expendable bathythermograph [XBT] profiles).

[1]  Rainer Bleck,et al.  Salinity-driven Thermocline Transients in a Wind- and Thermohaline-forced Isopycnic Coordinate Model of the North Atlantic , 1992 .

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

[3]  Alexey Kaplan,et al.  Mapping tropical Pacific sea level : Data assimilation via a reduced state space Kalman filter , 1996 .

[4]  P. L. Traon,et al.  AN IMPROVED MAPPING METHOD OF MULTISATELLITE ALTIMETER DATA , 1998 .

[5]  Jacques Verron,et al.  A singular evolutive extended Kalman filter for data assimilation in oceanography , 1998 .

[6]  Frank O. Bryan,et al.  A Model Comparison: Numerical Simulations of the North and Equatorial Atlantic Oceanic Circulation in Depth and Isopycnic Coordinates , 1996 .

[7]  Jacques Verron,et al.  An extended Kalman filter to assimilate satellite altimeter data into a nonlinear numerical model of the tropical Pacific Ocean: Method and validation , 1999 .

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

[9]  G. Evensen,et al.  Data assimilation and inverse methods in terms of a probabilistic formulation , 1996 .

[10]  George L. Mellor,et al.  Continuous assimilation of Geosat altimeter data into a three-dimensional primitive equation Gulf Stream model , 1994 .

[11]  K. Haines,et al.  Satellite altimetry data assimilation in the OCCAM global ocean model , 1999 .

[12]  P. Malanotte‐Rizzoli,et al.  An approximate Kaiman filter for ocean data assimilation: An example with an idealized Gulf Stream model , 1995 .

[13]  Robert T. Guza,et al.  Observations of nearshore circulation: Alongshore uniformity , 2003 .

[14]  Rainer Bleck,et al.  Wind‐driven spin‐up in eddy‐resolving ocean models formulated in isopycnic and isobaric coordinates , 1986 .

[15]  Pierre Brasseur,et al.  Assimilation of altimetric data in the mid-latitude oceans using the Singular Evolutive Extended Kalman filter with an eddy-resolving, primitive equation model , 1999 .

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

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

[18]  Ichiro Fukumori,et al.  Chapter 5 Data Assimilation by Models , 2001 .

[19]  D. Chelton,et al.  Chapter 2 Large-Scale Ocean Circulation , 2001 .

[20]  K. Haines,et al.  A Direct Method for Assimilating Sea Surface Height Data into Ocean Models with Adjustments to the Deep Circulation , 1991 .

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

[22]  Howard P. Hanson,et al.  Mixed Layer-Thermocline Interaction in a Three-Dimensional Isopycnic Coordinate Model , 1989 .

[23]  Rosemary Morrow,et al.  Chapter 3 Ocean Currents and Eddies , 2001 .

[24]  Richard G. Forbes,et al.  Assessment of the FOAM global data assimilation system for real-time operational ocean forecasting , 2000 .

[25]  P. Houtekamer,et al.  An Adaptive Ensemble Kalman Filter , 2000 .

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

[27]  Pierre Brasseur,et al.  Dynamical evolution of the error statistics with the SEEK filter to assimilate altimetric data in eddy‐resolving ocean models , 2001 .

[28]  D. Dee.  Simplification of the Kalman filter for meteorological data assimilation , 1991 .

[29]  G. Veronis.  Large Scale Ocean Circulation , 1973 .

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