Sequential data assimilation applied to a physical-biological model for the Bermuda Atlantic time series station

Abstract In this study, we investigate sequential data assimilation approaches for state estimation and prediction in a coupled physical–biological model for the Bermuda Atlantic Time Series (BATS) site. The model is 1-dimensional (vertical) in space and based on the General Ocean Turbulence Model (GOTM). Coupled to GOTM is a biological model that includes phytoplankton, detritus, dissolved inorganic nitrogen, chlorophyll and oxygen. We performed model ensemble runs by introducing variations in the biological parameters, each of which was assigned a probability distribution. We compare and contrast here 2 sequential data assimilation methods: the ensemble Kalman filter (EnKF) and sequential importance resampling (SIR). We assimilated different types of BATS observations, including particulate organic nitrogen, nitrate + nitrite, chlorophyll a and oxygen for the 2-year period from January 1990 to December 1991, and quantified the impact of the data assimilation on the model's predictive skill. By applying a cross-validation to the data-assimilative and deterministic simulations we found that the predictive skill was improved for 2-week forecasts. In our experiments the EnKF, which exhibited a stronger effect on the ensemble during the assimilation step, showed slightly higher improvements in the predictive skill than the SIR, which preserves dynamical model consistency in our implementation. Our numerical experiments show that statistical properties stabilize for ensemble sizes of 20 or greater with little improvement for larger ensembles.

[1]  Dale B. Haidvogel,et al.  Nitrogen cycling in the Middle Atlantic Bight: Results from a three‐dimensional model and implications for the North Atlantic nitrogen budget , 2006 .

[2]  Branko Ristic,et al.  Beyond the Kalman Filter: Particle Filters for Tracking Applications , 2004 .

[3]  G. Evensen,et al.  An Ensemble Kalman Filter with a complex marine ecosystem model: hindcasting phytoplankton in the Cretan Sea , 2003 .

[4]  Mark R. Abbott,et al.  Configuring an ecosystem model using data from the Bermuda Atlantic Time Series (BATS) , 2001 .

[5]  T. Kana,et al.  Dynamic model of phytoplankton growth and acclimation: responses of the balanced growth rate and the chlorophyll a:carbon ratio to light, nutrient-limitation and temperature , 1997 .

[6]  Yasuhiro Yamanaka,et al.  NEMURO—a lower trophic level model for the North Pacific marine ecosystem , 2007 .

[7]  Jens Schröter,et al.  Sequential weak constraint parameter estimation in an ecosystem model , 2003 .

[8]  Hugh L. MacIntyre,et al.  A dynamic model of photoadaptation in phytoplankton , 1996 .

[9]  K. Lindsay,et al.  Marine Biogeochemical Modeling: Recent Advances and Future Challenges , 2001 .

[10]  Wolfgang Rodi,et al.  Examples of calculation methods for flow and mixing in stratified fluids , 1987 .

[11]  Eileen E. Hofmann,et al.  Biogeochemical Data Assimilation , 2001 .

[12]  G. Kitagawa.  A self-organizing state-space model , 1998 .

[13]  R. Murtugudde,et al.  Annual ecosystem variability in the tropical Indian Ocean: Results of a coupled bio-physical ocean general circulation model , 2006 .

[14]  Michael Dowd,et al.  Bayesian statistical data assimilation for ecosystem models using Markov Chain Monte Carlo , 2007 .

[15]  Scott C. Doney,et al.  Modeling Ocean Ecosystems: The PARADIGM Program , 2006 .

[16]  Katja Fennel,et al.  Subsurface maxima of phytoplankton and chlorophyll: Steady‐state solutions from a simple model , 2003 .

[17]  K. Brusdala,et al.  A demonstration of ensemble-based assimilation methods with a layered OGCM from the perspective of operational ocean forecasting systems , 2003 .

[18]  Neil J. Gordon,et al.  A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking , 2002, IEEE Trans. Signal Process..

[19]  G. Evensen.  Data Assimilation: The Ensemble Kalman Filter , 2006 .

[20]  T. D. Dickey,et al.  Influence of mesoscale eddies on new production in the Sargasso Sea , 1998, Nature.

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

[22]  Peter R. Oke,et al.  Ensemble member generation for sequential data assimilation , 2008 .

[23]  G. Kitagawa.  Monte Carlo Filter and Smoother for Non-Gaussian Nonlinear State Space Models , 1996 .

[24]  Barbara A. Adams-Vanharn,et al.  Evaluation of the current state of mechanistic aquatic biogeochemical modeling: citation analysis and future perspectives. , 2006, Environmental science & technology.

[25]  R. Wanninkhof.  Relationship between wind speed and gas exchange over the ocean , 1992 .

[26]  Louis I. Gordon,et al.  Oxygen solubility in seawater : better fitting equations , 1992 .

[27]  Sovan Lek.  Uncertainty in ecological models , 2007 .

[28]  K. Fennel,et al.  Denitrification effects on air‐sea CO2 flux in the coastal ocean: Simulations for the northwest North Atlantic , 2008 .

[29]  B. Martin.  PARAMETER ESTIMATION , 2012, Statistical Methods for Biomedical Research.

[30]  Michael Dowd,et al.  A Bayesian approach to the ecosystem inverse problem , 2003 .

[31]  N. Gordon,et al.  Novel approach to nonlinear/non-Gaussian Bayesian state estimation , 1993 .

[32]  Jens Schröter,et al.  Testing a marine ecosystem model: Sensitivity analysis and parameter optimization , 2001 .

[33]  Karline Soetaert,et al.  Application of an Ensemble Kalman filter to a 1-D coupled hydrodynamic-ecosystem model of the Ligurian Sea , 2007 .

[34]  G. Evensen,et al.  Sequential Data Assimilation Techniques in Oceanography , 2003 .

[35]  E L Ionides,et al.  Inference for nonlinear dynamical systems , 2006, Proceedings of the National Academy of Sciences.

[36]  James D. Annan,et al.  Parameter estimation in an intermediate complexity earth system model using an ensemble Kalman filter , 2005 .

[37]  Adrian E. Raftery,et al.  Weather Forecasting with Ensemble Methods , 2005, Science.

[38]  Soroosh Sorooshian,et al.  Dual state-parameter estimation of hydrological models using ensemble Kalman filter , 2005 .

[39]  Scott C. Doney,et al.  Eddy‐driven sources and sinks of nutrients in the upper ocean: Results from a 0.1° resolution model of the North Atlantic , 2003 .

[40]  Andreas Oschlies,et al.  Parameter estimates of a zero-dimensional ecosystem model applying the adjoint method , 2001 .

[41]  Nicholas R. Bates,et al.  Overview of the US JGOFS Bermuda Atlantic Time-series Study (BATS): a decade-scale look at ocean biology and biogeochemistry , 2001 .

[42]  G. Evensen,et al.  Assimilation of ocean colour data into a biochemical model of the North Atlantic: Part 1. Data assimilation experiments , 2003 .

[43]  Marjorie A. M. Friedrichs.  A data assimilative marine ecosystem model of the central equatorial Pacific: Numerical twin experiments , 2001 .

[44]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[45]  R. Reynolds,et al.  The NCEP/NCAR 40-Year Reanalysis Project , 1996, Renewable Energy.

[46]  Lakshmi Kantha,et al.  A general ecosystem model for applications to primary productivity and carbon cycle studies in the global oceans , 2004 .

[47]  Scott C. Doney,et al.  Assessment of skill and portability in regional marine biogeochemical models : Role of multiple planktonic groups , 2007 .

[48]  G. Kitagawa.  Theory and Methods , 1998 .