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

Abstract

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.