The extended Kalman filter for forecast of algal bloom dynamics

Abstract

A deterministic ecosystem model is combined with an extended Kalman filter (EKF) to produce short term forecasts of algal bloom and dissolved oxygen dynamics in a marine fish culture zone (FCZ). The weakly flushed FCZ is modelled as a well-mixed system; the tidal exchange with the outer bay is lumped into a flushing rate that is numerically determined from a three-dimensional hydrodynamic model. The ecosystem model incorporates phytoplankton growth kinetics, nutrient uptake, photosynthetic production, nutrient sources from organic fish farm loads, and nutrient exchange with a sediment bed layer. High frequency field observations of chlorophyll, dissolved oxygen (DO) and hydro-meteorological parameters (sampling interval Δt = 1 day, 2 h, 1 h, respectively) and bi-weekly nutrient data are assimilated into the model to produce the combined state estimate accounting for the uncertainties. In addition to the water quality state variables, the EKF incorporates dynamic estimation of algal growth rate and settling velocity. The effectiveness of the EKF data assimilation is studied for a wide range of sampling intervals and prediction lead-times. The chlorophyll and dissolved oxygen estimated by the EKF are compared with field data of seven algal bloom events observed at Lamma Island, Hong Kong. The results show that the EKF estimate well captures the nonlinear error evolution in time; the chlorophyll level can be satisfactorily predicted by the filtered model estimate with a mean absolute error of around 1-2 μg/L. Predictions with 1-2 day lead-time are highly correlated with the observations (r = 0.7-0.9); the correlation stays at a high level for a lead-time of 3 days (r = 0.6-0.7). Estimated algal growth and settling rates are in accord with field observations; the more frequent DO data can compensate for less frequent algal biomass measurements. The present study is the first time the EKF is successfully applied to forecast an entire algal bloom cycle, suggesting the possibility of using EKF for real time forecast of algal bloom dynamics. © 2009 Elsevier Ltd. All rights reserved.