Parameter optimization and analysis of ecosystem models using simulated annealing: a case study at Station P

A simulated annealing optimization algorithm is formulated to optimize parameters of ecosystem models. The optimization is used to directly determine the model parameters required to reproduce the observed data. The optimization routine is formulated in a general manner and is easily modified to include additional information on both the desired model output and the model parameters. From the optimization routine, error analysis of the optimal parameters is provided by the error-covariance matrix which gives both the sensitivity of the model to each model parameter and the correlation coefficients between all pairs of model parameters. In addition, the optimization analysis provides a means of assessing the necessary model complexity required to model the available data. To demonstrate the technique, optimal parameters of three different ecosystem model configurations are determined from nitrate, phytoplankton, mesozooplankton and net phytoplankton productivity measurements at Station P. At Station P, error analysis of the optimal parameters indicates that the data are able to resolve up to 10 independent model parameters. This is always less than the number of unknown model parameters indicating that the optimal solutions are not unique. A simple nitrate-phosphate-zooplankton ecosystem is successful at reproducing the observations. To justify the use of a more complicated model at Station P requires additional data to constrain the optimization routine. Although there is evidence supporting the importance of the microbial loop at Station P, without additional ammonium and bacteria measurements one cannot validate a more complicated model that includes these processes.

[1]  Carl Wunsch,et al.  Oceanographic Experiment Design by Simulated Annealing , 1990 .

[2]  D. Kirchman,et al.  Carbon limitation of ammonium uptake by heterotrophic bacteria in the subarctic Pacific , 1990 .

[3]  D. Kirchman,et al.  Bacterial production and the sinking flux of particulate organic matter in the subarctic Pacific , 1992 .

[4]  T. Parsons,et al.  Comparative oceanic ecology of the plankton communities of the subarctic Atlantic and Pacific oceans , 1988 .

[5]  J. Steele The Structure of Marine Ecosystems , 1974 .

[6]  A. Vézina Construction of Flow Networks Using Inverse Methods , 1989 .

[7]  K. Baker,et al.  Relation between photosynthetically available radiation and total insolation at the ocean surface under clear skies1 , 1987 .

[8]  J. Kruger Simulated Annealing: A Tool for Data Assimilation into an Almost Steady Model State , 1993 .

[9]  Robert E. Ulanowicz,et al.  Mathematical Models in Biological Oceanography. , 1984 .

[10]  Bruce W. Frost,et al.  Life histories of large, grazing copepods in a subarctic ocean gyre: Neocalanus plumchrus, Neocalanus cristatus, and Eucalanus bungii in the Northeast Pacific , 1984 .

[11]  W. Broecker,et al.  Dynamical limitations on the Antarctic iron fertilization strategy , 1991, Nature.

[12]  Thomas D. Brock,et al.  Calculating solar radiation for ecological studies , 1981 .

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

[14]  C. D. Mcallister Aspects of Estimating Zooplankton Production from Phytoplankton Production , 1969 .

[15]  J. G. Field,et al.  The Ecological Role of Water-Column Microbes in the Sea* , 1983 .

[16]  T. Platt,et al.  Food web dynamics in the ocean. I. Best-estimates of flow networks using inverse methods , 1988 .

[17]  B. Frost,et al.  Grazing control of phytoplankton stock in the open subarctic Pacific Ocean: a model assessing the role of mesozooplankton, particularly the large calanoid copepods Neocalanus spp. , 1987 .

[18]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[19]  John H. Steele,et al.  ENVIRONMENTAL CONTROL OF PHOTOSYNTHESIS IN THE SEA , 1962 .

[20]  W. Large,et al.  Upper Ocean Thermal Response to Strong Autumnal Forcing of the Northeast Pacific , 1986 .

[21]  C. D. Mcallister Zooplankton Studies at Ocean Weather Station "P" in the Northeast Pacific Ocean , 1961 .

[22]  W. Richard,et al.  TEMPERATURE AND PHYTOPLANKTON GROWTH IN THE SEA , 1972 .

[23]  W. Thacker,et al.  An Optimal-Control/Adjoint-Equations Approach to Studying the Oceanic General Circulation , 1989 .

[24]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

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

[26]  J. Marotzke The role of integration time in determining a steady-state through data assimilation , 1992 .

[27]  B. Frost The role of grazing in nutrient-rich areas of the open sea , 1991 .