Improving the physics of a coupled physical-biogeochemical model of the North Atlantic through data assimilation: impact on the ecosystem.

Abstract Several studies on coupled physical–biogeochemical models have shown that major deficiencies in the biogeochemical fields arise from the deficiencies in the physical flow fields. This paper examines the improvement of the physics through data assimilation, and the subsequent impact on the ecosystem response in a coupled model of the North Atlantic. Sea surface temperature and sea surface height data are assimilated with a sequential method based on the SEEK filter adapted to the coupling needs. The model domain covers the Atlantic from 20°S to 70°N at eddy-permitting resolution. The biogeochemical model is a NPZD-DOM model based on the P3ZD formulation. The results of an annual assimilated simulation are compared with an annual free simulation. With assimilation, the representation of the mixed layer depth is significantly improved in mid latitudes, even though the mixed layer depth is generally overestimated compared to the observations. The representation of the mean and variance of the currents is also significantly improved. The nutrient input in the euphotic zone is used to assess the data assimilation impact on the ecosystem. Data assimilation results in a 50% reduction of the input due to vertical mixing in mid-latitudes, and in a four- to six-fold increase of the advective fluxes in mid-latitudes and subtropics. Averaged zonally, the net impact is a threefold increase for the subtropical gyre, and a moderate (20–30%) decrease at mid and high latitudes. Surface chlorophyll concentration increases along the subtropical gyre borders, but little changes are detected at mid and high latitudes. An increase of the primary production appears along the Gulf Stream path, but it represents only 12% on average for mid and high latitudes. In the subtropical gyre centre, primary production is augmented but stays underestimated (20% of observations). These experiments show the benefits of physical data assimilation in coupled physical–biogeochemical applications.

[1]  Andreas Oschlies,et al.  Improved Representation of Upper-Ocean Dynamics and Mixed Layer Depths in a Model of the North Atlantic on Switching from Eddy-Permitting to Eddy-Resolving Grid Resolution , 2002 .

[2]  N. Garcia,et al.  An estimation of annual new production and carbon fluxes in the northeast Atlantic Ocean during 2001 , 2005 .

[3]  J. Toggweiler,et al.  A seasonal three‐dimensional ecosystem model of nitrogen cycling in the North Atlantic Euphotic Zone , 1993 .

[4]  Akio Arakawa,et al.  Computational Design of the Basic Dynamical Processes of the UCLA General Circulation Model , 1977 .

[5]  S. Doney,et al.  An intermediate complexity marine ecosystem model for the global domain , 2001 .

[6]  W. Jenkins,et al.  Nitrate flux into the euphotic zone near Bermuda , 1988, Nature.

[7]  K. Haines,et al.  Altimetric assimilation in a Mediterranean general circulation model , 1997 .

[8]  A. Oschlies.  Equatorial nutrient trapping in biogeochemical ocean models: The role of advection numerics , 2000 .

[9]  O. Aumont,et al.  Dimethylsulfoniopropionate (DMSP) and dimethylsulfide (DMS) sea surface distributions simulated from a global three-dimensional ocean carbon cycle model , 2002 .

[10]  Scott C. Doney,et al.  The role of mesoscale variability on plankton dynamics in the North Atlantic , 2001 .

[11]  Gilles Reverdin,et al.  Global high-resolution mapping of ocean circulation from TOPEX/Poseidon and ERS-1 and -2 , 2000 .

[12]  D. Antoine,et al.  Oceanic primary production: 2. Estimation at global scale from satellite (Coastal Zone Color Scanner) chlorophyll , 1996 .

[13]  J. Waddington,et al.  Cutover peatlands: A persistent source of atmospheric CO2 , 2002 .

[14]  Jean-Michel Brankart,et al.  Implementation of a multivariate data assimilation scheme for isopycnic coordinate ocean models: Application to a 1993–1996 hindcast of the North Atlantic Ocean circulation , 2003 .

[15]  Daniele Iudicone,et al.  Mixed layer depth over the global ocean: An examination of profile data and a profile-based climatology , 2004 .

[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]  Charles-Emmanuel Testut.  Assimilation de données satellitales avec un filtre de Kalman de rang réduit dans un modèle aux équations primitives de l'océan Atlantique , 2000 .

[18]  Jean-Michel Brankart,et al.  Assimilation of sea-surface temperature and altimetric observations during 1992–1993 into an eddy permitting primitive equation model of the North Atlantic Ocean , 2003 .

[19]  S. Spall,et al.  A numerical model of mesoscale frontal instabilities and plankton dynamics — I. Model formulation and initial experiments , 2000 .

[20]  Richard G. Forbes,et al.  Assessment of the FOAM global data assimilation system for real-time operational ocean forecasting , 2000 .

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

[22]  T. Platt,et al.  Vertical Nitrate Fluxes in the Oligotrophic Ocean , 1986, Science.

[23]  E. Mahieu,et al.  Ground‐based infrared spectroscopic measurements of carbonyl sulfide: Free tropospheric trends from a 24‐year time series of solar absorption measurements , 2002 .

[24]  Erwin Suess,et al.  Particulate organic carbon flux in the oceans—surface productivity and oxygen utilization , 1980, Nature.

[25]  Olivier Aumont.  Etude du cycle naturel du carbone dans un modele 3d de l'ocean mondial , 1998 .

[26]  A. Oschlies.  Nutrient supply to the surface waters of the North Atlantic - a model study , 2002 .

[27]  P. Smolarkiewicz,et al.  The multidimensional positive definite advection transport algorithm: further development and applications , 1986 .

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

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

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

[31]  Corinne Le Quéré,et al.  Propagation of climatic events on ocean stratification, marine biology, and CO2: Case studies over the 1979–1999 period , 2003 .

[32]  L. A. Anderson,et al.  Physical and biological modeling in the Gulf Stream region:: I. Data assimilation methodology , 2000 .

[33]  Jean-Michel Brankart,et al.  Comparative assimilation of Topex/Poseidon and ERS altimeter data and of TAO temperature data in the tropical Pacific ocean during 1994–1998, and the mean sea-surface height issue , 2003 .

[34]  Gurvan Madec,et al.  Impacts of sub-mesoscale physics on phytoplankton production and, subduction , 2001 .

[35]  P. Delecluse,et al.  OPA 8.1 Ocean General Circulation Model reference manual , 1998 .

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

[37]  Scott C. Doney,et al.  A three‐dimensional, multinutrient, and size‐structured ecosystem model for the North Atlantic , 2004 .

[38]  H. Sverdrup,et al.  On Conditions for the Vernal Blooming of Phytoplankton , 1953 .

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

[40]  Patrice Klein,et al.  Impact of sub-mesoscale physics on production and subduction of phytoplankton in an oligotrophic regime , 2001 .

[41]  A. Oschlies.  Can eddies make ocean deserts bloom? , 2002 .

[42]  Yosuke Fujii,et al.  Application of Nonlinear Constraints in a Three-Dimensional Variational Ocean Analysis , 2005 .

[43]  Keith Haines,et al.  Altimetric assimilation with water property conservation , 1996 .

[44]  Bruno Blanke,et al.  Variability of the Tropical Atlantic Ocean Simulated by a General Circulation Model with Two Different Mixed-Layer Physics , 1993 .

[45]  Pierre Brasseur,et al.  Assimilation of altimetric data in the mid-latitude oceans using the Singular Evolutive Extended Kalman filter with an eddy-resolving, primitive equation model , 1999 .

[46]  A. Oschlies,et al.  An eddy‐permitting coupled physical‐biological model of the North Atlantic: 1. Sensitivity to advection numerics and mixed layer physics , 1999 .

[47]  J. Sarmiento,et al.  Ecosystem behavior at Bermuda Station “S” and ocean weather station “India”: A general circulation model and observational analysis , 1993 .

[48]  S. Doney.  Major challenges confronting marine biogeochemical modeling , 1999 .

[49]  B. Barnier,et al.  A new analysis of hydrographic data in the Atlantic and its application to an inverse modeling study , 1998 .

[50]  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 .