Impact of improved light calculations on predicted phytoplankton growth and heating in an idealized upwelling-downwelling channel geometry

Ocean ecosystem models require accurate calculations of both hydrodynamics and biology; those calculations in turn require accurate calculation of in-water irradiance. Ecosystem models now achieve great accuracy in their hydrodynamical predictions, and the biological modules are becoming correspondingly sophisticated. The optical calculations are however often oversimplified, to the possible detriment of the physical and biological predictions. We used a recently developed, extremely fast radiative transfer code, EcoLight-S, to study differences in ecosystem and thermal development in an idealized upwelling-downwelling system when simple versus accurate irradiance calculations are used. The use of accurate irradiances gave up to 57% differences in chlorophyll concentrations after two weeks of simulated time, compared to predictions based on irradiances obtained using a simple exponential attenuation formula. Accurate irradiance calculations increased sea surface temperatures and decreased temperatures at depth, leading to increased stratification. Use of EcoLight-S couples the physical and biological calculations so that biology feeds back to physics, and vice versa. EcoLight-S outputs ancillary quantities such as remote sensing reflectance and in-water spectral irradiance, which can be used to validate ecosystem predictions using remotely sensed ocean color imagery or optical measurements from buoys or gliders, without the need to convert such measurements to chlorophyll values. After optimization, the ecosystem model total run times with EcoLight-S were less than 20% more than for the analytical irradiance models. We also found that the use of 24 h average irradiances gave factor-of-two differences in chlorophyll concentrations compared to the use of a diel irradiance pattern with the same 24 h average value.

[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]  M. Fujii,et al.  Modeling carbon and silicon cycling in the equatorial Pacific , 2007 .

[3]  A. Morel Optical modeling of the upper ocean in relation to its biogenous matter content (case I waters) , 1988 .

[4]  James J. Simpson,et al.  Alternative Parameterizations of Downward Irradiance and Their Dynamical Significance , 1981 .

[5]  James C. McWilliams,et al.  Equilibrium structure and dynamics of the California Current System , 2003 .

[6]  John S. Parslow,et al.  A model of annual plankton cycles , 2013 .

[7]  Temel Oguz,et al.  Simulation of annual plankton productivity cycle in the Black Sea by a one‐dimensional physical‐biological model , 1996 .

[8]  F. Muller‐Karger,et al.  The effects of temporal variability of mixed layer depth on primary productivity around Bermuda , 1994 .

[9]  Curtis D. Mobley,et al.  Fast light calculations for ocean ecosystem and inverse models. , 2011, Optics express.

[10]  Masahiko Fujii,et al.  The Value of Adding Optics to Ecosystem Models: A Case Study , 2007 .

[11]  Alexander F. Shchepetkin,et al.  The regional oceanic modeling system (ROMS): a split-explicit, free-surface, topography-following-coordinate oceanic model , 2005 .

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

[13]  Ronghua Ma,et al.  Impact of multiple satellite ocean color samplings in a day on assessing phytoplankton dynamics , 2012, Ocean Science Journal.

[14]  M. Fasham Variations in the seasonal cycle of biological production in subarctic oceans: A model sensitivity analysis , 1995 .

[15]  K. Fennel,et al.  Quantifying biological carbon export for the northwest North Atlantic continental shelves , 2009 .

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

[17]  K. Carder,et al.  A simple spectral solar irradiance model for cloudless maritime atmospheres , 1990 .

[18]  Kendall L. Carder,et al.  Carbon cycling in the upper waters of the Sargasso Sea: II. Numerical simulation of apparent and inherent optical properties , 1999 .

[19]  D. Antoine,et al.  Heating rate within the upper ocean in relation to its bio-optical state , 1994 .

[20]  Chuanmin Hu,et al.  Global distribution of Case-1 waters: An analysis from SeaWiFS measurements , 2006 .

[21]  Dale B. Haidvogel,et al.  Generation and unstable evolution of a density‐driven Eastern Poleward Current: The Iberian Poleward Current , 2003 .

[22]  Soo Chin Liew,et al.  Penetration of solar radiation in the upper ocean: A numerical model for oceanic and coastal waters , 2005 .

[23]  Oscar Schofield,et al.  Dynamics of turbid buoyant plumes and the feedbacks on near‐shore biogeochemistry and physics , 2008 .

[24]  P. Xiu,et al.  Connections between physical, optical and biogeochemical processes in the Pacific Ocean , 2014 .

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

[26]  M. Dinniman,et al.  Cross-shelf exchange in a model of the Ross Sea circulation and biogeochemistry , 2003 .

[27]  Richard C. Dugdale,et al.  One-dimensional ecosystem model of the equatorial Pacific upwelling system. Part I: model development and silicon and nitrogen cycle , 2002 .

[28]  James C. McWilliams,et al.  Quasi-Monotone Advection Schemes Based on Explicit Locally Adaptive Dissipation , 1998 .

[29]  G. Hurtt,et al.  A pelagic ecosystem model calibrated with BATS data , 1996 .

[30]  Curtis D. Mobley,et al.  Fast and accurate irradiance calculations for ecosystem models , 2009 .

[31]  C. Cox Statistics of the sea surface derived from sun glitter , 1954 .

[32]  P. Xiu,et al.  Spatial and temporal variability in phytoplankton carbon, chlorophyll, and nitrogen in the North Pacific , 2012 .

[33]  Curtis D. Mobley,et al.  Improved irradiances for use in ocean heating, primary production, and photo-oxidation calculations. , 2012, Applied optics.

[34]  James C. McWilliams,et al.  A method for computing horizontal pressure‐gradient force in an oceanic model with a nonaligned vertical coordinate , 2003 .

[35]  C. Mobley,et al.  Phase function effects on oceanic light fields. , 2002, Applied optics.

[36]  T. Piatt,et al.  Ocean primary production calculated by spectral and broad-band models , 1992 .

[37]  C. Roy,et al.  Modelling the shear edge eddies of the southern Agulhas Current , 2003 .

[38]  C. Paulson,et al.  Irradiance Measurements in the Upper Ocean , 1977 .

[39]  Kendall L. Carder,et al.  Carbon cycling in the upper waters of the Sargasso Sea: I. Numerical simulation of differential carbon and nitrogen fluxes , 1999 .

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

[41]  David A. Siegel,et al.  Ocean Radiant Heating. Part II: Parameterizing Solar Radiation Transmission through the Upper Ocean , 2000 .

[42]  K. Stamnes,et al.  Comparison of numerical models for computing underwater light fields. , 1993, Applied optics.

[43]  W. Munk,et al.  Measurement of the Roughness of the Sea Surface from Photographs of the Sun’s Glitter , 1954 .

[44]  Kevin E. Kohler,et al.  Influences of diurnal and intraseasonal forcing on mixed‐layer and biological variability in the central Arabian Sea , 2001 .

[45]  Oscar Schofield,et al.  Resolving the Impacts and Feedback of Ocean Optics on Upper Ocean Ecology , 2001 .