Importance of solar subsurface heating in ocean general circulation models

The importance of subsurface heating on surface mixed layer properties in an ocean general circulation model (OGCM) is examined using attenuation of solar irradiance with depth below the ocean surface. The depth-dependent attenuation of subsurface heating is given by global monthly mean fields for the attenuation of photosynthetically available radiation (PAR), kPAR. These global fields of kPAR are derived from Sea-viewing Wide Field-of-view Sensor (SeaWiFS) data on the spectral diffuse attenuation coefficient at 490 nm (k490), and have been processed to have the smoothly varying and continuous coverage necessary for use in OGCM applications. These monthly fields provide the first complete global data sets of subsurface optical fields that can be used for OGCM applications of subsurface heating and bio-optical processes. The effect on global OGCM prediction of sea surface temperature (SST) and surface mixed layer depth (MLD) is examined when solar heating, as given by monthly mean kPAR and PAR fields, is included in the model. It is found that subsurface heating yields a marked increase in the SST predictive skill of the OGCM at low latitudes. No significant improvement in MLD predictive skill is obtained when including subsurface heating. Use of the monthly mean kPAR produces an SST decrease of up to 0.8°C and a MLD increase of up to only 4–5 m for climatological surface forcing, with this primarily confined to the equatorial regions. Remarkably, a constant kPAR value of 0.06 m−1, which is indicative of optically clear open ocean conditions, is found to serve very well for OGCM prediction of SST and MLD over most of the global ocean.

[1]  Paul R. Lowe.  An Approximating Polynomial for the Computation of Saturation Vapor Pressure. , 1977 .

[2]  C. Davis,et al.  Biological effects of Gulf Stream meandering , 1993 .

[3]  Judith L. Lean,et al.  Global upper ocean heat storage response to radiative forcing from changing solar irradiance and increasing greenhouse gas/aerosol concentrations , 1998 .

[4]  Andrew F. Loughe,et al.  A Reduced-Gravity Isopycnal Ocean Model: Hindcasts of El Niño , 1995 .

[5]  J. Zaneveld,et al.  An arc tangent model of irradiance in the sea , 1980 .

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

[7]  M. Cane,et al.  A Reduced-Gravity, Primitive Equation, Isopycnal Ocean GCM: Formulation and Simulations , 1995 .

[8]  James V. Gardner,et al.  Mapping U.S. continental shelves , 1998 .

[9]  E. Joseph Metzger,et al.  Dynamics of the Kuroshio/Oyashio current system using eddy-resolving models of the North Pacific Ocean , 1996 .

[10]  J. Carton,et al.  Annual cycle of sea surface temperature in the tropical Atlantic Ocean , 1997 .

[11]  P. Gallacher,et al.  Numerical simulations of the Arabian Sea using tracers as proxies for phytoplankton biomass , 1995 .

[12]  S. Hooker,et al.  SeaWiFS Technical Report Series. Volume 43; SeaWiFS Prelaunch Technical Report Series Final Cumulative Index , 1998 .

[13]  J. Brock,et al.  The phytoplankton bloom in the northwestern Arabian Sea during the southwest monsoon of 1979 , 1991 .

[14]  J. Turner,et al.  A one‐dimensional model of the seasonal thermocline II. The general theory and its consequences , 1967 .

[15]  Kevin E. Kohler,et al.  A four-component ecosystem model of biological activity in the Arabian Sea , 1996 .

[16]  T. Platt,et al.  Ocean primary production and available light: further algorithms for remote sensing , 1988 .

[17]  H. Hurlburt,et al.  The Nondeterministic Nature of Kuroshio Penetration and Eddy Shedding in the South China Sea , 2001 .

[18]  Pearn P. Niiler,et al.  Upper Ocean Heat Budget During the Hawaii-to-Tahiti Shuttle Experiment , 1983 .

[19]  Alexander Horch,et al.  Solar heating of the oceans—diurnal, seasonal and meridional variation , 1984 .

[20]  P. Martin.  Simulation of the mixed layer at OWS November and Papa with several models , 1985 .

[21]  James K. B. Bishop,et al.  Spatial and temporal variability of global surface solar irradiance , 1991 .

[22]  Embedding a mixed layer model into an ocean general circulation model of the Atlantic : the importance of surface mixing for heat flux and temperature , 1994 .

[23]  Gene C. Feldman,et al.  Influence of penetrating solar radiation on the heat budget of the equatorial Pacific Ocean , 1990, Nature.

[24]  L. Kantha,et al.  An improved mixed layer model for geophysical applications , 1994 .

[25]  V. Garçon,et al.  Impact of the North Equatorial Current meandering on a pelagic ecosystem: A modeling approach , 1996 .

[26]  Bifurcation of the Kuroshio Extension at the Shatsky Rise , 1998 .

[27]  Peter A. Rochford,et al.  Efficient and Accurate Bulk Parameterizations of Air-Sea Fluxes for Use in General Circulation Models , 2000 .

[28]  Peter A. Rochford,et al.  An optimal definition for ocean mixed layer depth , 2000 .

[29]  A Calculation of Ocean Heat Storage and Effective Ocean Surface Layer Depths for the Northern Hemisphere , 1984 .

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

[31]  H. Hurlburt,et al.  Coupled dynamics of the South China Sea, the Sulu Sea, and the Pacific Ocean , 1996 .

[32]  R. W. Austin,et al.  Spectral dependence of the diffuse attenuation coefficient of light in ocean waters , 1986 .

[33]  K. Arrigo,et al.  Volume 11, SeaWiFS Postlaunch Calibration and Validation Analyses, Part 3 , 2000 .

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

[35]  A. H. Murphy,et al.  Skill Scores Based on the Mean Square Error and Their Relationships to the Correlation Coefficient , 1988 .

[36]  J. Kindle,et al.  The interaction of southwest monsoon upwelling, advection and primary production in the northwest Arabian Sea , 1997 .

[37]  John C. Kindle,et al.  Hindcasting of wind-driven anomalies using a reduced-gravity global ocean model , 1992 .

[38]  Sol Hellerman,et al.  Normal Monthly Wind Stress Over the World Ocean with Error Estimates , 1983 .

[39]  W. Large,et al.  Oceanic vertical mixing: a review and a model with a nonlocal boundary layer parameterization , 1994 .

[40]  R. Lukas,et al.  Surface Buoyancy Forcing and the Mixed Layer of the Western Pacific Warm Pool: Observations and 1D Model Results , 1996 .

[41]  James C. Kitchen,et al.  Vertical structure of productivity and its vertical integration as derived from remotely sensed observations , 1993 .

[42]  Dennis K. Clark,et al.  Atmospheric effects in the remote sensing of phytoplankton pigments , 1980 .

[43]  J. Cherniawsky,et al.  An upper ocean general circulation model for climate studies: global simulation with seasonal cycle , 1992 .

[44]  J. Zaneveld,et al.  The influence of optical water type on the heating rate of a constant depth mixed layer , 1981 .

[45]  D. Houghton,et al.  Approximation of Ocean Heat Storage by Ocean–Atmosphere Energy Exchange: Implications for Seasonal Cycle Mixed Layer Ocean Formulations , 1987 .

[46]  J. Kindle,et al.  Sensitivity of the Arabian Sea mixed layer to 1994–1995 operational wind products , 2000 .

[47]  A. Lorenc.  A Global Three-Dimensional Multivariate Statistical Interpolation Scheme , 1981 .

[48]  H. Hurlburt,et al.  The contribution of the global thermohaline circulation to the Pacific to Indian Ocean Throughflow via Indonesia , 1997 .

[49]  Menghua Wang,et al.  Calibration of SeaWiFS on orbit , 2000, SPIE Optics + Photonics.

[50]  Gregory E. Terrie,et al.  Ocean Optical Database , 1992 .

[51]  P. Niiler,et al.  Effects of instability waves in the mixed layer of the equatorial Pacific , 1997 .