Application of a SEEK filter to a 1D biogeochemical model of the Ligurian Sea: Twin experiments and real in-situ data assimilation

The Singular Evolutive Extended Kalman (SEEK) filter has been implemented to assimilate in-situ data in a 1D coupled physical-ecosystem model of the Ligurian Sea. The biogeochemical model describes the partly decoupled nitrogen and carbon cycles of the pelagic food web. The GHER hydrodynamic model (1D version) is used to represent the physical forcings. The data assimilation scheme (SEEK filter) parameterizes the error statistics by means of a set of empirical orthogonal functions (EOFs). Twin experiments are first performed with the aim to choose the suitable experimental protocol (observation and estimation vectors, number of EOFs, frequency of the assimilation,…) and to assess the SEEK filter performances. This protocol is then applied to perform real data assimilation experiments using the DYFAMED data base. By assimilating phytoplankton observations, the method has allowed to improve not only the representation of the phytoplankton community, but also of other variables such as zooplankton and bacteria that evolve with model dynamics and that are not corrected by the data assimilation scheme. The validation of the assimilation method and the improvement of model results are studied by means of suitable error measurements.

[1]  I. Charpentier,et al.  Can biogeochemical fluxes be recovered from nitrate and chlorophyll data? A case study assimilating data in the Northwestern Mediterranean Sea at the JGOFS-DYFAMED station , 2003 .

[2]  Ibrahim Hoteit,et al.  A singular evolutive extended Kalman filter to assimilate real in situ data in a 1-D marine ecosystem model , 2003 .

[3]  Hervé Claustre,et al.  Phytoplankton pigment distribution in relation to upper thermocline circulation in the eastern Mediterranean Sea during winter , 2001 .

[4]  J. Marty.  The DYFAMED time-series program (French-JGOFS) , 2002 .

[5]  Ibrahim Hoteit,et al.  A singular evolutive interpolated Kalman filter for efficient data assimilation in a 3-D complex physical biogeochemical model of the Cretan Sea , 2003 .

[6]  Marilaure Grégoire,et al.  Revisited ecosystem model (MODECOGeL) of the Ligurian Sea: seasonal and interannual variability due to atmospheric forcing , 2002 .

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

[8]  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 .

[9]  O. Connan,et al.  Nutrients in the Mediterranean Sea, mass balance and statistical analysis of concentrations with respect to environmental change , 1998 .

[10]  Triantafyllou George,et al.  Towards a data assimilation system for the Cretan Sea ecosystem using a simplified Kalman filter , 2004 .

[11]  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 .

[12]  Dinh-Tuan Pham,et al.  An adaptively reduced-order extended Kalman filter for data assimilation in the tropical Pacific , 2004 .

[13]  L. Merlivat,et al.  The relevant time scales in estimating the air-sea CO2 exchange in a mid-latitude region , 2002 .

[14]  R. Todling,et al.  Data Assimilation in the Presence of Forecast Bias: The GEOS Moisture Analysis , 2000 .

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

[16]  Arthur Gelb,et al.  Applied Optimal Estimation , 1974 .

[17]  J. Marty,et al.  Seasonal and interannual dynamics of nutrients and phytoplankton pigments in the western Mediterranean Sea at the DYFAMED time-series station (1991 1999) , 2002 .

[18]  F. Rassoulzadegan,et al.  Full-depth profile (0–2000 m) of bacteria, heterotrophic nanoflagellates and ciliates in the NW Mediterranean Sea: Vertical partitioning of microbial trophic structures , 2002 .

[19]  Jean-Claude Marty,et al.  Seasonal and interannual variations in phytoplankton production at DYFAMED time-series station, northwestern Mediterranean Sea , 2002 .

[20]  Dinh-Tuan Pham,et al.  A simplified reduced order Kalman filtering and application to altimetric data assimilation in Tropical Pacific , 2002 .

[21]  Geneviève Lacroix,et al.  Influence of meteorological variability on primary production dynamics in the Ligurian Sea (NW Mediterranean Sea) with a 1D hydrodynamic/biological model , 1998 .

[22]  Karline Soetaert,et al.  FEMME, a flexible environment for mathematically modelling the environment , 2002 .

[23]  K. Taylor.  Summarizing multiple aspects of model performance in a single diagram , 2001 .

[24]  Salim Djenidi,et al.  Perspective in Three-Dimensional Modelling of the Marine System , 1987 .

[25]  P. Morin,et al.  Temporal trends in nutrient ratios: chemical evidence of Mediterranean ecosystem changes driven by human activity , 2002 .

[26]  G. Evans,et al.  A vertically resolved model for phytoplankton aggregation , 2000 .

[27]  C. Tamburini,et al.  Biopolymer hydrolysis and bacterial production under ambient hydrostatic pressure through a 2000 m water column in the NW Mediterranean , 2002 .

[28]  Geir Evensen,et al.  A singular evolutive extended Kalman filter to assimilate ocean color data in a coupled physical–biochemical model of the North Atlantic ocean , 2001 .

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

[30]  A. Jazwinski.  Stochastic Processes and Filtering Theory , 1970 .

[31]  J. Nihoul,et al.  Three-dimensional models of marine and estuarine dynamics , 1987 .

[32]  Jean-Claude Marty,et al.  Phytoplankton pigment variations during the transition from spring bloom to oligotrophy in the northwestern Mediterranean sea , 2000 .

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

[34]  K. Soetaert,et al.  Study of the seasonal cycle of the biogeochemical processes in the Ligurian Sea using a 1D interdisciplinary model , 2005 .

[35]  Jean-Michel Brankart,et al.  Assimilation of sea-surface temperature and altimetric observations during 1992–1993 into an eddy permitting primitive equation model of the North Atlantic Ocean , 2003 .

[36]  N. Mahowald,et al.  Inverse methods in global biogeochemical cycles , 2000 .

[37]  Jean-Michel Brankart,et al.  Comparative assimilation of Topex/Poseidon and ERS altimeter data and of TAO temperature data in the tropical Pacific ocean during 1994–1998, and the mean sea-surface height issue , 2003 .

[38]  J. Beckers,et al.  Dissection of the GHER turbulence closure scheme , 1999 .

[39]  Philippe Courtier,et al.  Unified Notation for Data Assimilation : Operational, Sequential and Variational , 1997 .