A New Solar Radiation Penetration Scheme for Use in Ocean Mixed Layer Studies: An Application to the Black Sea Using a Fine-Resolution Hybrid Coordinate Ocean Model (HYCOM)

A 1/25° 1/25° cos(lat) (longitude latitude) (3.2-km resolution) eddy-resolving Hybrid Coordinate Ocean Model (HYCOM) is introduced for the Black Sea and used to examine the effects of ocean turbidity on upper-ocean circulation features including sea surface height and mixed layer depth (MLD) on annual mean climatological time scales. The model is a primitive equation model with a K-profile parameterization (KPP) mixed layer submodel. It uses a hybrid vertical coordinate that combines the advantages of isopycnal, , and z-level coordinates in optimally simulating coastal and open-ocean circulation features. This model approach is applied to the Black Sea for the first time. HYCOM uses a newly developed time-varying solar penetration scheme that treats attenuation as a continuous quantity. This scheme includes two bands of solar radiation penetration, one that is needed in the top 10 m of the water column and another that penetrates to greater depths depending on the turbidity. Thus, it is suitable for any ocean general circulation model that has fine vertical resolution near the surface. With this scheme, the optical depth–dependent attenuation of subsurface heating in HYCOM is given by monthly mean fields for the attenuation of photosynthetically active radiation (kPAR) during 1997–2001. These satellite-based climatological kPAR fields are derived from Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) data for 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 the Black Sea model applications. HYCOM simulations are driven by two sets of high-frequency climatological forcing, but no assimilation of ocean data is then used to demonstrate the importance of including spatial and temporal varying attenuation depths for the annual mean prediction of upper-ocean quantities in the Black Sea, which is very turbi d( kPAR 0.15 m 1 , in general). Results are reported from three model simulations driven by each atmospheric forcing set using different values for the kPAR. A constant solar-attenuation optical depth of 17 m (clear water assumption), as opposed to using spatially and temporally varying attenuation depths, changes the surface circulation, especially in the eastern Black Sea. Unrealistic sub–mixed layer heating in the former results in weaker stratification at the base of the mixed layer and a deeper MLD than observed. As a result, the deep MLD off Sinop (at around 42.5°N, 35.5°E) weakens the surface currents regardless of the atmospheric forcing used in the model simulations. Using the SeaWiFS-based monthly turbidity climatology gives a shallower MLD with much stronger stratification at the base and much better agreement with observations. Because of the high Black Sea turbidity, the simulation with all solar radiation absorbed at the surface case gives results similar to the simulations using turbidity from SeaWiFS in the annual means, the aspect of the results investigated in this paper.

[1]  B. Vanleer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[2]  Peter A. Rochford,et al.  The NRL Layered Global Ocean Model (NLOM) with an Embedded Mixed Layer Submodel: Formulation and Tuning , 2003 .

[3]  H. Hurlburt,et al.  Climatological SST and MLD Predictions from a Global Layered Ocean Model with an Embedded Mixed Layer , 2003 .

[4]  Melinda S. Peng,et al.  Navy Operational Global Atmospheric Prediction System (NOGAPS): Forcing for Ocean Models , 2002 .

[5]  Robert Frouin,et al.  Response of the equatorial Pacific to chlorophyll pigment in a mixed layer isopycnal ocean general circulation model , 2001 .

[6]  James C. McWilliams,et al.  Sensitivity to Surface Forcing and Boundary Layer Mixing in a Global Ocean Model: Annual-Mean Climatology , 1997 .

[7]  J. Ohlmann,et al.  Ocean Radiant Heating in Climate Models , 2003 .

[8]  E. Stanev,et al.  Oceanic response to atmospheric forcing derived from different climatic data sets. Intercomparison study for the Black Sea , 1998 .

[9]  Rainer Bleck,et al.  An oceanic general circulation model framed in hybrid isopycnic-Cartesian coordinates , 2002 .

[10]  Robert A. Arnone,et al.  Importance of solar subsurface heating in ocean general circulation models , 2001 .

[11]  P. Malanotte‐Rizzoli,et al.  Interdisciplinary Studies Integrating The Black Sea Biogeochemistry and Circulation Dynamics , 2002 .

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

[13]  J. Murray,et al.  Variations in the chemistry of the Black Sea on a time scale of decades (1960–1995) , 2001 .

[14]  D. Siegel,et al.  Radiant heating of the western equatorial Pacific during TOGA‐COARE , 1998 .

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

[16]  E. Stanev On the mechanisms of the Black Sea circulation , 1990 .

[17]  S. Maritorena,et al.  Bio-optical properties of oceanic waters: A reappraisal , 2001 .

[18]  A. Kideys,et al.  A review of zooplankton investigations of the Black Sea over the last decade , 2000 .

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

[20]  E. Özsoy,et al.  The evolution of Mediterranean water in the Black Sea: interior mixing and material transport by double diffusive intrusions , 1993 .

[21]  Peter A. Rochford,et al.  The Impact of Water Turbidity on Interannual Sea Surface Temperature Simulations in a Layered Global Ocean Model , 2004 .

[22]  H. Christian,et al.  Strength and consolidation properties of surficial sediments, Flemish Pass: effects of biological processes , 1991 .

[23]  H. Hurlburt,et al.  Air–Sea Flux Estimates And The 1997–1998 Enso Event , 2002 .

[24]  W. Timothy Liu,et al.  Evaporation and solar irradiance as regulators of sea surface temperature in annual and interannual changes , 1994 .

[25]  E. Özsoy,et al.  Coastal/deep ocean interactions in the Black Sea and their ecological/environmental impacts , 1996 .

[26]  E. Skyllingstad,et al.  Nonlocal fluxes and Stokes drift effects in the K-profile parameterization , 2002 .

[27]  Caskey,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS I . THE BASIC EXPERIMENT , 1962 .

[28]  David A. Siegel,et al.  Solar radiation, phytoplankton pigments and the radiant heating of the equatorial Pacific warm pool , 1995 .

[29]  H. I. Sur,et al.  Evolution of satellite derived mesoscale thermal patterns in the Black Sea , 1997 .

[30]  J. Beckers,et al.  Barotropic and baroclinic oscillations in strongly stratified ocean basins - Numerical study of the Black Sea , 1999 .

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

[32]  Joanna Staneva,et al.  The impact of the baroclinic eddies and basin oscillations on the transitions between different quasi-stable states of the Black Sea circulation , 2000 .

[33]  Timothy R. Parsons,et al.  Biological oceanography: an introduction. Second edition , 1997 .

[34]  R. LeVeque Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[35]  Malcolm J. Bowman,et al.  Rim current and coastal eddy mechanisms in an eddy-resolving Black Sea general circulation model , 2001 .

[36]  Temel Oguz,et al.  Circulation in the surface and intermediate layers of the Black Sea , 1993 .

[37]  A. Sobel,et al.  A role for ocean biota in tropical intraseasonal atmospheric variability , 2003 .

[38]  J. Holden,et al.  The storage and aging of continental runoff in large reservoir systems of the world , 1997 .

[39]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[40]  E. Stanev,et al.  Modelling the impact of atmospheric and terrestrial inputs on the Black Sea coastal dynamics , 2001 .

[41]  Philip B. Duffy,et al.  An extended data set of river discharges for validation of general circulation models , 1996 .

[42]  H. Hurlburt,et al.  How Does Solar Attenuation Depth Affect the Ocean Mixed Layer? Water Turbidity and Atmospheric Forcing Impacts on the Simulation of Seasonal Mixed Layer Variability in the Turbid Black Sea* , 2005 .

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

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

[45]  E. Stanev,et al.  The sensitivity of the heat exchange at sea surface to meso and sub-basin scale eddies: Model study for the Black Sea , 2001 .

[46]  Antonio J. Busalacchi,et al.  Effects of Penetrative Radiation on the Upper Tropical Ocean Circulation , 2002 .

[47]  H. Treut,et al.  A climatological field of mean optical properties of the world ocean , 1986 .

[48]  G. Halliwell,et al.  Evaluation of vertical coordinate and vertical mixing algorithms in the HYbrid-Coordinate Ocean Model (HYCOM) , 2004 .

[49]  Timothy P. Boyer,et al.  World ocean atlas 2001 : objective analyses, data statistics, and figures : CD-ROM documentation , 2002 .

[50]  Eric P. Chassignet,et al.  North Atlantic Simulations with the Hybrid Coordinate Ocean Model (HYCOM): Impact of the Vertical Coordinate Choice, Reference Pressure, and Thermobaricity , 2003 .

[51]  Rainer Bleck,et al.  A Wind-Driven Isopycnic Coordinate Model of the North and Equatorial Atlantic Ocean 1 , 1990 .

[52]  E. Özsoy,et al.  Oceanography of the Black Sea: A review of some recent results , 1997 .

[53]  A. Kostianoy,et al.  Analysis of velocity field in the eastern Black Sea from satellite data during the Black Sea '99 experiment , 2002 .

[54]  Alan J. Wallcraft,et al.  Sea Surface Temperature Sensitivity to Water Turbidity from Simulations of the Turbid Black Sea Using HYCOM , 2005 .

[55]  Mojib Latif,et al.  Warm Pool Physics in a Coupled GCM , 1996 .

[56]  T. Oguz,et al.  Observations on the Rim Current structure, CIW formation and transport in the western Black Sea , 1999 .

[57]  R. Daley Atmospheric Data Analysis , 1991 .

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

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

[60]  Charlie N. Barron,et al.  The Modular Ocean Data Assimilation System (MODAS) , 2002 .

[61]  E. Schneider,et al.  Sensitivity of the Simulated Annual Cycle of Sea Surface Temperature in the Equatorial Pacific to Sunlight Penetration , 1998 .

[62]  E. Özsoy,et al.  Hydrographic properties and ventilation of the Black Sea , 1991 .

[63]  N. Schneider,et al.  The Indonesian Throughflow and the Global Climate System , 1998 .

[64]  E. Laws,et al.  Primary production in the deep blue sea , 1990 .

[65]  D. Pierce,et al.  Anatomy of North Pacific Decadal Variability , 2002 .