A calibration strategy for dynamic succession models including several phytoplankton groups

A fundamental problem in water quality modeling is adequately representing the changing state of aquatic ecosystems as accurately as possible, but with appropriate mathematical relationships without creating a highly complex and overly parameterized model. A model more complex than necessary will require more input and results in unaffordable calibration times. In this work we propose and test a calibration strategy for a one-dimensional dynamic physical-ecological model (DYRESM-CAEDYM) to reproduce the seasonal changes in the functional composition of the phytoplankton community existing in El Gergal reservoir (Seville, Spain). The community is described as a succession of functional groups with different response to environmental conditions. First, we performed a sensitivity analysis to identify the parameters to include in the calibration process, and then applied a global optimization algorithm to fit the model for each algal group in a sequential fashion. Finally we simulated all the functional groups adopting parameter values established during the group-by-group calibrations. Our results show that the performance of this approach is strictly related with: (1) the level of system description (i.e. the model structure and the number of functional groups simulated); (2) the level of information included in the calibration process (i.e. the observations); and (3) the non-linear interactions among functional groups. Functional segmentation of the model should be minimized even though groups with different environmental requirements must be discriminated. Although magnitudes of biomass peaks were not always estimated correctly, the calibrated model was able to predict peak sequence and timing of dominant phytoplankton groups. Thus our study showed that: (1) model structure and nature of observations adopted have to be in agreement with the level of organization in the system; (2) integration of automatic calibration strategies is a useful approach in complex deterministic ecological models.

[1]  F. Recknagel ANNA – Artificial Neural Network model for predicting species abundance and succession of blue-green algae , 1997, Hydrobiologia.

[2]  P. Reichert,et al.  Modelling functional groups of phytoplankton in three lakes of different trophic state , 2008 .

[3]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[4]  J. R. Romero,et al.  One- and three-dimensional biogeochemical simulations of two differing reservoirs , 2004 .

[5]  Michael D. Guiry,et al.  Freshwater Algal Flora of the British Isles , 2003 .

[6]  Aysegul Aksoy,et al.  Calibration and Verification of QUAL2E Using Genetic Algorithm Optimization , 2007 .

[7]  James L. Martin,et al.  Hydrodynamics and Transport for Water Quality Modeling , 1998 .

[8]  Matthew R. Hipsey,et al.  Implementation of ecological modeling as an effective management and investigation tool: Lake Kinneret as a case study , 2009 .

[9]  Scott A. Wells,et al.  CE-QUAL-W2: A Two-dimensional, Laterally Averaged, Hydrodynamic and Water Quality Model, Version 3.5 , 2006 .

[10]  K. Lindenschmidt,et al.  The effect of water column mixing on phytoplankton succession, diversity and similarity , 1998 .

[11]  Basak Guven,et al.  A review and classification of the existing models of cyanobacteria , 2006 .

[12]  Francisco E. Werner,et al.  Calibration of the NEMURO nutrient-phytoplankton-zooplankton food web model to a coastal ecosystem: Evaluation of an automated calibration approach , 2007 .

[13]  I. Oldal,et al.  Dominant species, functional assemblages and frequency of equilibrium phases in late summer phytoplankton assemblages in Hungarian small shallow lakes , 2003, Hydrobiologia.

[14]  A. E. Irish,et al.  Exploring the potential of the PROTECH model to investigate phytoplankton community theory , 1999, Hydrobiologia.

[15]  E. Moreno-Ostos,et al.  The effects of diel changes in circulation and mixing on the longitudinal distribution of phytoplankton in a canyon-shaped Mediterranean reservoir , 2010 .

[16]  George B. Arhonditsis,et al.  Eutrophication model for Lake Washington (USA) Part II-model calibration and system dynamics analysis , 2005 .

[17]  H. Dau,et al.  A fluorometric method for the differentiation of algal populations in vivo and in situ , 2004, Photosynthesis Research.

[18]  Jason P. Antenucci,et al.  Simulating the thermal dynamics of Lake Kinneret , 2003 .

[19]  R. Spear Eutrophication in peel inlet—II. Identification of critical uncertainties via generalized sensitivity analysis , 1980 .

[20]  Dennis Trolle,et al.  Predicting the effects of reduced external nitrogen loading on the nitrogen dynamics and ecological state of deep Lake Ravn, Denmark, using the DYRESM–CAEDYM model , 2008 .

[21]  R. Hornung Numerical Modelling of Stratification in Lake Constance with the 1-D hydrodynamic model DYRESM , 2002 .

[22]  A. Mulligan,et al.  Genetic Algorithms for Calibrating Water Quality Models , 1998 .

[23]  J. Laanemets,et al.  A simple tool for the early prediction of the cyanobacteria Nodularia spumigena bloom biomass in the Gulf of Finland , 2006 .

[24]  S. Sorooshian,et al.  Effective and efficient global optimization for conceptual rainfall‐runoff models , 1992 .

[25]  Helmut Hillebrand,et al.  BIOVOLUME CALCULATION FOR PELAGIC AND BENTHIC MICROALGAE , 1999 .

[26]  Nikolaus Hansen,et al.  Completely Derandomized Self-Adaptation in Evolution Strategies , 2001, Evolutionary Computation.

[27]  M. Bonnet,et al.  DyLEM-1D: a 1D physical and biochemical model for planktonic succession, nutrients and dissolved oxygen cycling Application to a hyper-eutrophic reservoir , 2004 .

[28]  David P. Hamilton,et al.  Modelling the impact of zooplankton grazing on phytoplankton biomass during a dinoflagellate bloom in the Swan River Estuary, Western Australia , 2001 .

[29]  Soroosh Sorooshian,et al.  Toward improved calibration of hydrologic models: Multiple and noncommensurable measures of information , 1998 .

[30]  Alessandro Oggioni,et al.  A biogeochemical model of Lake Pusiano (North Italy) and its use in the predictability of phytoplankton blooms: first preliminary results , 2006 .

[31]  Francis Wilson,et al.  Total organic carbon as a predictor of biological wastewater treatment efficiency and kinetic reaction rates , 1997 .

[32]  Jasper A Vrugt,et al.  Improved evolutionary optimization from genetically adaptive multimethod search , 2007, Proceedings of the National Academy of Sciences.

[33]  J. Romero,et al.  Predicting upwelling, boundary mixing, and nutrient fluxes in lakes , 2000 .

[34]  A. E. Irish,et al.  Modelling freshwater phytoplankton communities: an exercise in validation , 2000 .

[35]  Colin S. Reynolds,et al.  The ecology of freshwater phytoplankton , 1984 .

[36]  Colin S. Reynolds,et al.  Towards a functional classification of the freshwater phytoplankton , 2002 .

[37]  David P. Hamilton,et al.  Prediction of water quality in lakes and reservoirs. Part I — Model description , 1997 .

[38]  Donald J. O'Connor,et al.  11 – Phytoplankton–Zooplankton–Nutrient Interaction Model for Western Lake Erie , 1975 .

[39]  Avi Ostfeld,et al.  A hybrid genetic—instance based learning algorithm for CE-QUAL-W2 calibration , 2005 .

[40]  C. Ancey,et al.  An exact solution for ideal dam‐break floods on steep slopes , 2008 .

[41]  D. G. George,et al.  Using a numerical model (PROTECH) to examine the impact of water transfers on phytoplankton dynamics in a Mediterranean reservoir , 2007, Limnetica.

[42]  M. B. Beck,et al.  Uncertainty, identifiability and the propagation of prediction errors: A case study of lake ontario , 1991 .

[43]  Luigi Naselli-Flores,et al.  Use and misuse in the application of the phytoplankton functional classification: a critical review with updates , 2009, Hydrobiologia.

[44]  Eligius M. T. Hendrix,et al.  A method for robust calibration of ecological models under different types of uncertainty , 1994 .

[45]  David P. Hamilton,et al.  Calibrating the Dynamic Reservoir Simulation Model (DYRESM) and filling required data gaps for one‐dimensional thermal profile predictions in a boreal lake , 2007 .

[46]  S. Jørgensen,et al.  Development of a structurally dynamic model for forecasting the effects of restoration of Lake Fure, Denmark , 2006 .

[47]  A. Blumberg,et al.  Quantifying Uncertainty in Estuarine and Coastal Ocean Circulation Modeling , 2008 .

[48]  Julian D. Olden,et al.  An artificial neural network approach for studying phytoplankton succession , 2000, Hydrobiologia.

[49]  Robert V. Thomann,et al.  Phytoplankton Modeling in the Embayments of Lakes , 1983 .

[50]  A. E. Irish,et al.  The ecological basis for simulating phytoplankton responses to environmental change (PROTECH) , 2001 .

[51]  Jeffrey G. Arnold,et al.  Automatic calibration of a distributed catchment model , 2001 .

[52]  H. Utermöhl Zur Vervollkommnung der quantitativen Phytoplankton-Methodik , 1958 .

[53]  L. Vörös,et al.  Phytoplankton biomass and chlorophyll-a in some shallow lakes in central Europe , 1991, Hydrobiologia.

[54]  Jörg Imberger,et al.  Phytoplankton patchiness and their role in the modelled productivity of a large, seasonally stratified lake , 2008 .

[55]  E. Moreno-Ostos,et al.  TEMPORAL AND SPATIAL VARIATIONS IN THE QUALITY OF WATER IN EL GERGAL RESERVOIR, SEVILLE, SPAIN , 2010 .

[56]  Blahoslav Maršálek,et al.  In situ Quantification of Phytoplankton in Reservoirs Using a Submersible Spectrofluorometer , 2005, Hydrobiologia.

[57]  Petros Koumoutsakos,et al.  Reducing the Time Complexity of the Derandomized Evolution Strategy with Covariance Matrix Adaptation (CMA-ES) , 2003, Evolutionary Computation.

[58]  J. Doherty,et al.  A hybrid regularized inversion methodology for highly parameterized environmental models , 2005 .

[59]  Mark E. Borsuk,et al.  Approaches to Evaluate Water Quality Model Parameter Uncertainty for Adaptive TMDL Implementation1 , 2007 .

[60]  J. R. Romero,et al.  Management strategies for a eutrophic water supply reservoir--San Roque, Argentina. , 2003, Water science and technology : a journal of the International Association on Water Pollution Research.

[61]  S. Sorooshian,et al.  A Shuffled Complex Evolution Metropolis algorithm for optimization and uncertainty assessment of hydrologic model parameters , 2002 .

[62]  Mingshun Jiang,et al.  A model study of the coupled biological and physical dynamics in Lake Michigan , 2002 .

[63]  Apha-Awwa-Wpcf Métodos normalizados para el análisis de aguas potables y residuales , 1992 .

[64]  R. Kalchev,et al.  Some allometric and non-allometric relationships between chlorophyll-a and abundance variables of phytoplankton , 1996, Hydrobiologia.

[65]  Ming-Der Yang,et al.  Eutrophication modelling of reservoirs in Taiwan , 2006, Environ. Model. Softw..

[66]  David P. Hamilton,et al.  Prediction of water quality in lakes and reservoirs: Part II - Model calibration, sensitivity analysis and application , 1997 .

[67]  Colin S. Reynolds,et al.  Vegetation processes in the pelagic : a model for ecosystem theory , 1997 .

[68]  E. Moreno-Ostos,et al.  The influence of external perturbations on the functional composition of phytoplankton in a Mediterranean reservoir , 2009, Hydrobiologia.

[69]  John Doherty,et al.  Efficient accommodation of local minima in watershed model calibration , 2006 .

[70]  P. Reichert,et al.  Biogeochemical model of Lake Zurich: model equations and results , 2001 .

[71]  Carlos E. Ruiz,et al.  Using spatially distributed parameters and multi‐response objective functions to solve parameterization of complex applications of semi‐distributed hydrological models , 2008 .

[72]  Donald C. Pierson,et al.  Weather driven influences on phytoplankton succession in a shallow lake during contrasting years: Application of PROTBAS , 2007 .