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)*
Authors: Kara A. Birol , Wallcraft Alan J. , Hurlburt Harley E. 
Published: March 2005

Abstract
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 turbid (kPAR &gt; 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.