Stochastic parameterizations of biogeochemical uncertainties in a 1/4° NEMO/PISCES model for probabilistic comparisons with ocean color data

In spite of recent advances, biogeochemical models are still unable to represent the full complexity of natural ecosystems. Their formulations are mainly based on empirical laws involving many parameters. Improving biogeochemical models therefore requires to properly characterize model uncertainties and their consequences. Subsequently, this paper investigates the potential of using random processes to simulate some uncertainties of the 1/4° coupled Physical–Biogeochemical NEMO/PISCES model of the North Atlantic ocean. Starting from a deterministic simulation performed with the original PISCES formulation, we propose a generic method based on AR(1) random processes to generate perturbations with temporal and spatial correlations. These perturbations are introduced into the model formulations to simulate 2 classes of uncertainties: the uncertainties on biogeochemical parameters and the uncertainties induced by unresolved scales in the presence of non-linear processes. Using these stochastic parameterizations, a probabilistic version of PISCES is designed and a 60-member ensemble simulation is performed. With respect to the simulation of chlorophyll, the relevance of the probabilistic configuration and the impacts of these stochastic parameterizations are assessed. In particular, it is shown that the ensemble simulation is in good agreement with the SeaWIFS ocean color data. Using these observations, the statistical consistency (reliability) of the ensemble is evaluated with rank histograms. Finally, the benefits expected from the probabilistic description of uncertainties (model error) are discussed in the context of future ocean color data assimilation.

[1]  Patrick M. Holligan,et al.  Basin-scale variability of phytoplankton biomass, production and growth in the Atlantic Ocean , 2000 .

[2]  L. Bopp,et al.  Globalizing results from ocean in situ iron fertilization studies , 2006 .

[3]  Momme Butenschön,et al.  A comparison of different versions of the SEEK Filter for assimilation of biogeochemical data in numerical models of marine ecosystem dynamics , 2012 .

[4]  Patrick Lehodey,et al.  Integrating biogeochemistry and ecology into ocean data assimilation systems , 2009 .

[5]  Jean-Marc Molines,et al.  Eddy compensation and controls of the enhanced sea‐to‐air CO2 flux during positive phases of the Southern Annular Mode , 2013 .

[6]  Laurent Bopp,et al.  Controlling factors of the oxygen balance in the Arabian Sea's OMZ , 2012 .

[7]  E. Woodward,et al.  Seasonal and spatial variability in plankton production and respiration in the Subtropical Gyres of the Atlantic Ocean , 2009 .

[8]  Paul G. Falkowski,et al.  Evolution of the nitrogen cycle and its influence on the biological sequestration of CO2 in the ocean , 1997, Nature.

[9]  Pierre Brasseur,et al.  The SEEK filter method for data assimilation in oceanography: a synthesis , 2006 .

[10]  Allan R. Robinson,et al.  Physical and biological modeling in the Gulf Stream region Part II. Physical and biological processes , 2001 .

[11]  Paul J. Harrison,et al.  Effects of temperature on growth rate, cell composition and nitrogen metabolism in the marine diatom Thalassiosira pseudonana (Bacillariophyceae) , 2002 .

[12]  Egil Sakshaug,et al.  Growth rates of ten diatom species from the Barents Sea at different irradiances and day lengths , 1990 .

[13]  Watson W. Gregg,et al.  Skill Assessment in Ocean Biological Data Assimilation , 2009 .

[14]  M. Kahru,et al.  Ocean Color Chlorophyll Algorithms for SEAWIFS , 1998 .

[15]  Jean-Claude Dutay,et al.  Quantifying the roles of ocean circulation and biogeochemistry in governing ocean carbon-13 and atmospheric carbon dioxide at the last glacial maximum , 2009 .

[16]  Gurvan Madec,et al.  The onset of a bloom after deep winter convection in the northwestern Mediterranean sea: mesoscale process study with a primitive equation model , 1998 .

[17]  W. Gregg,et al.  Global and regional evaluation of the SeaWiFS chlorophyll data set , 2004 .

[18]  I. Hoteit,et al.  Assimilation of ocean colour data into a Biogeochemical Flux Model of the Eastern Mediterranean Sea , 2006 .

[19]  Katja Fennel,et al.  Introduction and Assessment of Measures for Quantitative Model-Data Comparison Using Satellite Images , 2010, Remote. Sens..

[20]  Olivier Aumont,et al.  PISCES-v2: an ocean biogeochemical model for carbon and ecosystem studies , 2015 .

[21]  Jean-Michel Brankart,et al.  Impact of uncertainties in the horizontal density gradient upon low resolution global ocean modelling , 2013 .

[22]  M. Maqueda,et al.  Sensitivity of a global sea ice model to the treatment of ice thermodynamics and dynamics , 1997 .

[23]  Stéphane Blain,et al.  An ecosystem model of the global ocean including Fe, Si, P colimitations , 2003 .

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

[25]  Jacques Verron,et al.  A singular evolutive extended Kalman filter for data assimilation in oceanography , 1998 .

[26]  Jens Schröter,et al.  A comparison of error subspace Kalman filters , 2005 .

[27]  R. Geider,et al.  The role of iron in phytoplankton photosynthesis, and the potential for iron-limitation of primary productivity in the sea , 1994, Photosynthesis Research.

[28]  Olivier Aumont,et al.  Decadal variations in equatorial Pacific ecosystems and ferrocline/pycnocline decoupling , 2008 .

[29]  Christoph Heinze,et al.  Reconciling surface ocean productivity, export fluxes and sediment composition in a global biogeochemical ocean model , 2006 .

[30]  O. Talagrand,et al.  Evaluation of probabilistic prediction systems for a scalar variable , 2005 .

[31]  Pierre Brasseur,et al.  Toward a multivariate reanalysis of the North Atlantic Ocean biogeochemistry during 1998-2006 based on the assimilation of SeaWiFS chlorophyll data , 2012 .

[32]  G. Madec.  NEMO ocean engine , 2008 .

[33]  R. Buizza,et al.  A Comparison of the ECMWF, MSC, and NCEP Global Ensemble Prediction Systems , 2005 .

[34]  Jean-Michel Brankart,et al.  Stochastic estimation of biogeochemical parameters of a 3D ocean coupled physical–biogeochemical model: Twin experiments , 2011 .

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

[36]  Renate Hagedorn,et al.  Representing model uncertainty in weather and climate prediction , 2005 .

[37]  C. L. Smith,et al.  The impact of mesoscale eddies on plankton dynamics in the upper ocean , 1996 .

[38]  Thierry Penduff,et al.  An ERA40-based atmospheric forcing for global ocean circulation models , 2010 .

[39]  P. Falkowski,et al.  Photosynthetic rates derived from satellite‐based chlorophyll concentration , 1997 .

[40]  John A. Raven,et al.  The role of trace metals in photosynthetic electron transport in O2-evolving organisms , 1999, Photosynthesis Research.

[41]  John Marra.  Ocean productivity: A personal perspective since the first Liege Colloquium , 2015 .

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

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

[44]  Akinori Ito,et al.  Global modeling study of potentially bioavailable iron input from shipboard aerosol sources to the ocean , 2013 .

[45]  Cédric Jamet,et al.  Retrieving the vertical distribution of chlorophyll a concentration and phytoplankton community composition from in situ fluorescence profiles: A method based on a neural network with potential for global‐scale applications , 2015 .

[46]  Guillem Candille,et al.  Assessment of an ensemble system that assimilates Jason-1/Envisat altimeter data in a probabilistic model of the North Atlantic ocean circulation , 2014 .

[47]  Jorgen S. Frederiksen,et al.  Stochastic subgrid parameterizations for atmospheric and oceanic flows , 2012 .

[48]  Jean-Michel Brankart,et al.  Characterization of mixing errors in a coupled physical biogeochemical model of the North Atlantic: implications for nonlinear estimation using Gaussian anamorphosis , 2010 .

[49]  Stephanie Dutkiewicz,et al.  Meteorological modulation of the North Atlantic spring bloom , 2001 .

[50]  Andreas Oschlies,et al.  Eddy-induced enhancement of primary production in a model of the North Atlantic Ocean , 1998, Nature.

[51]  Gurvan Madec,et al.  Large-scale impacts of submesoscale dynamics on phytoplankton: Local and remote effects , 2012 .

[52]  Jean-Michel Brankart,et al.  Stochastic estimation of biogeochemical parameters from Globcolour ocean colour satellite data in a North Atlantic 3D ocean coupled physical–biogeochemical model , 2013 .

[53]  M. Dowd.  Estimating parameters for a stochastic dynamic marine ecological system , 2011 .

[54]  A. H. Murphy.  A New Vector Partition of the Probability Score , 1973 .

[55]  A. Longhurst.  Seasonal cycles of pelagic production and consumption , 1995 .

[56]  H. Claustre,et al.  Vertical distribution of phytoplankton communities in open ocean: An assessment based on surface chlorophyll , 2006 .

[57]  Scott C. Doney,et al.  Evaluating global ocean carbon models: The importance of realistic physics , 2004 .

[58]  T. Palmer,et al.  Stochastic representation of model uncertainties in the ECMWF ensemble prediction system , 2007 .

[59]  Olivier Aumont,et al.  Ocean biogeochemistry exhibits contrasting responses to a large scale reduction in dust deposition , 2007 .

[60]  Michael J. Follows,et al.  Nutrient streams and their induction into the mixed layer , 2006 .

[61]  Geir Evensen,et al.  A singular evolutive extended Kalman filter to assimilate ocean color data in a coupled physical–biochemical model of the North Atlantic ocean , 2001 .

[62]  Marina Lévy,et al.  The influence of mesoscale and submesoscale heterogeneity on ocean biogeochemical reactions , 2013 .

[63]  Cara Wilson,et al.  The rocky road from research to operations for satellite ocean-colour data in fishery management , 2011 .

[64]  Olivier Aumont,et al.  Response of diatoms distribution to global warming and potential implications: A global model study , 2005 .

[65]  Thierry Penduff,et al.  Impact of partial steps and momentum advection schemes in a global ocean circulation model at eddy-permitting resolution , 2006 .

[66]  Timothy J. Smyth,et al.  Assimilation of remotely-sensed optical properties to improve marine biogeochemistry modelling , 2014 .

[67]  Taro Takahashi,et al.  Redfield ratio based on chemical data from isopycnal surfaces , 1985 .

[68]  Momme Butenschön,et al.  Building the capacity for forecasting marine biogeochemistry and ecosystems: recent advances and future developments , 2015 .

[69]  Guillem Candille,et al.  A generic approach to explicit simulation of uncertainty in the NEMO ocean model , 2015 .

[70]  Lawrence M. Murray,et al.  A Bayesian approach to state and parameter estimation in a Phytoplankton-Zooplankton model , 2010 .

[71]  Cecile Cabanes,et al.  The CORA dataset: validation and diagnostics of in-situ ocean temperature and salinity measurements , 2013 .