An adjoint data assimilation approach for estimating parameters in a three-dimensional ecosystem model

In this paper an ecosystem model, including phytoplankton, zooplankton, nitrate, ammonium, phosphate and detritus, is described. The model is driven by physical fields derived from a three-dimensional physical transport model. Simulation includes nitrate input from a river. Simulated results are then sampled and the sampled data are used in sequential numerical experiments to assess the ability of using an adjoint data assimilation approach for estimating the poorly known parameters of the ecosystem model, such as growth and death rate, half-saturation constant of nutrients, etc. Data with different spatial and temporal resolution over 1 week are assimilated into the ecosystem model. Assimilation of data at 30 grid stations with a sampling interval of 6 h is proved to be adequate for recovering all the parameters of the ecosystem model. Both the spatial and temporal resolution of the data are mutually complementary in the assimilative model. Thus, improvement of either of them can result in improvement of model parameter recoveries. The assimilation of phytoplankton data is essential to recover the model parameters. Phytoplankton is the core of the food web and without the information on phytoplankton, the structure of the ecosystem cannot be constructed correctly. The adjoint method can work well with the noisy data. In the twin experiments with noisy data, the parameters can be recovered but the error is increased. The results of the model and parameter recovery are sensitive to the initial conditions of state variables, so the determination of the initial condition is as important as that of the model parameter. The spatial and temporal resolution and the data type of the observations in Analysis and Modelling Research of the Ecosystem in the Bohai Sea (AMREB) are suitable for the recovery of the model parameters used in this study.

[1]  J. Vallino.  Improving marine ecosystem models: Use of data assimilation and mesocosm experiments , 2000 .

[2]  A. Watson,et al.  A modelling investigation of the role of phytoplankton in the balance of carbon at the surface of the North Atlantic , 1991 .

[3]  Lisan Yu,et al.  Variational Estimation of the Wind Stress Drag Coefficient and the Oceanic Eddy Viscosity Profile , 1991 .

[4]  Daniel R. Lynch,et al.  An adjoint data assimilation approach to diagnosis of physical and biological controls on Pseudocalanus spp. in the Gulf of Maine-Georges bank region , 1998 .

[5]  Changsheng Chen,et al.  Plankton production in tidal fronts : A model of Georges Bank in summer , 1996 .

[6]  Richard J. Matear,et al.  Parameter optimization and analysis of ecosystem models using simulated annealing: a case study at Station P , 1995 .

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

[8]  H. Ducklow,et al.  A nitrogen-based model of plankton dynamics in the oceanic mixed layer , 1990 .

[9]  Liu Dongyan,et al.  Estimating biomass of phytoplankton in the Jiaozhou bay. I . Phytoplankton biomass estimated from cell volume and plasma volume , 2000 .

[10]  R. Hoffman,et al.  A four-dimensional analysis exactly satisfying equations of motion. [for atmospheric circulation , 1986 .

[11]  E. Hofmann,et al.  Time series sampling and data assimilation in a simple marine ecosystem model , 1996 .

[12]  Allan R. Robinson,et al.  Coupled physical and biological modeling of the spring bloom in the North Atlantic (I): model formulation and one dimensional bloom processes , 1995 .

[13]  Jan O. Backhaus,et al.  A finite difference general circulation model for shelf seas and its application to low frequency variability on the North European shelf , 1987 .

[14]  Liang Zhao,et al.  Phytoplankton dynamics in the Bohai Sea—observations and modelling , 2004 .

[15]  A. Moll.  Regional distribution of primary production in the North Sea simulated by a three-dimensional model , 1998 .

[16]  Ionel M. Navon,et al.  Conjugate-Gradient Methods for Large-Scale Minimization in Meteorology , 1987 .

[17]  B. Frost.  A modelling study of processes regulating plankton standing stock and production in the open Subarctic Pacific Ocean , 1993 .

[18]  Glenn R. Flierl,et al.  Behavior of a simple plankton model with food-level acclimation by herbivores , 1986 .

[19]  D. Wiesenburg,et al.  Influences of river discharge on biological production in the inner shelf : A coupled biological and physical model of the Louisiana-Texas Shelf , 1997 .

[20]  Y. Sasaki.  SOME BASIC FORMALISMS IN NUMERICAL VARIATIONAL ANALYSIS , 1970 .

[21]  Liang Zhao,et al.  Annual cycle and budgets of nutrients in the Bohai Sea , 2002 .

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

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

[24]  James G. Richman,et al.  Data assimilation and a pelagic ecosystem model: parameterization using time series observations , 1998 .

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

[26]  R. W. Lardner,et al.  VARIATIONAL PARAMETER ESTIMATION FOR A TWO-DIMENSIONAL NUMERICAL TIDAL MODEL , 1992 .

[27]  Allan R. Robinson,et al.  Coupled physical and biological modelling of the spring bloom in the North Atlantic (II): three dimensional bloom and post-bloom processes , 1995 .

[28]  Ionel Michael Navon,et al.  Practical and theoretical aspects of adjoint parameter estimation and identifiability in meteorology and oceanography , 1998 .

[29]  J. Toggweiler,et al.  A seasonal three‐dimensional ecosystem model of nitrogen cycling in the North Atlantic Euphotic Zone , 1993 .

[30]  J. O'Brien,et al.  Variational data assimilation and parameter estimation in an equatorial Pacific ocean model , 1991 .

[31]  Peter G. Challenor,et al.  A Markov chain Monte Carlo method for estimation and assimilation into models , 1997 .

[32]  Jan O. Backhaus,et al.  A three-dimensional model for the simulation of shelf sea dynamics , 1985 .

[33]  M. J. D. Powell,et al.  Restart procedures for the conjugate gradient method , 1977, Math. Program..

[34]  G. Radach,et al.  Estimation of the variability of production by simulating annual cycles of phytoplankton in the central North Sea , 1993 .

[35]  Michio J. Kishi,et al.  A data assimilation technique applied to estimate parameters for the NEMURO marine ecosystem model , 2004 .

[36]  E. Hofmann,et al.  Plankton dynamics on the outer southeastern U.S. continental shelf. Part II: A time-dependent biological model , 1988 .

[37]  J. Wroblewski,et al.  A model of phytoplankton plume formation during variable Oregon upwelling , 1977 .

[38]  F. A. Lootsma,et al.  Numerical methods for non-linear optimization , 1974 .