Features of Parameterizing Turbulent Interaction with Underlying Surface in the Regional Thermodynamic Model of Sea Ice

D. D. Zavyalov, T. A. Solomakha

Marine Hydrophysical Institute of RAS, Sevastopol, Russian Federation

e-mail: zavyalov.dd@mhi-ras.ru


Purpose. The study is aimed at assessing the impact of choice of parameterizing the turbulent heat transfer at the ocean – atmosphere boundary upon the basic characteristics of ice regime in the Taganrog Bay apex.

Methods and Results. Thermal seasonal dynamics of the snow-ice cover thickness was studied using the non-stationary thermodynamic model of sea ice. The algorithm for defining the turbulent fluxes of momentum, sensible and latent heat in the sea ice regional model is based on the Monin – Obukhov similarity theory. The numerical experiments were performed for the winter seasons of 2007–2008 and 2017–2018, the meteorological conditions of which differed significantly. The transfer coefficients were determined both as the constant values, and as those depending on the atmosphere stratification and the aerodynamic roughness of underlying surface. Implementing stable stratification implied application of three different expressions for determining the stability functions. To avoid the iteration process required for numerical solving of the equations of the Monin – Obukhov similarity theory, the transfer coefficient parameterizations based on the approach relating these coefficients to the bulk Richardson number, were used in the model. Having been analyzed, the results of simulating the evolution of seasonal snow-ice thickness permitted to reveal the features in applying the parameterization of turbulent fluxes.

Conclusions. It is shown that the type of parameterizing the turbulent fluxes for the winters characterized by stable frosty weather and ice cover, does not impact significantly the basic elements of ice regime in the Taganrog Bay apex. However, in case the snow-ice cover is highly unstable during a season, the simulation results significantly depend on the method of determining the turbulent transfer coefficients. The best results in reconstructing the seasonal changes in ice cover thickness were obtained when using both the constant coefficients of turbulent transfer CH = CE equal to ≈ 1.7·10-3 and those depending on the atmosphere stratification at the ice geometric roughness equal to 8–10 cm.


Monin – Obukhov theory, parameterization, turbulent fluxes, sea ice


The study was carried out within the framework of the state assignment on theme FNNN-2021-0004 “Oceanological processes”

Original russian text

Original Russian Text © D. D. Zavyalov, T. A. Solomakha, 2023, published in MORSKOY GIDROFIZICHESKIY ZHURNAL, Vol. 39, Iss. 4 (2023)

For citation

Zavyalov, D.D. and Solomakha, T.A., 2023. Features of Parameterizing Turbulent Interaction with Underlying Surface in the Regional Thermodynamic Model of Sea Ice. Physical Oceanography, 30(4), pp. 385-397.


  1. Monin, A.S. and Obukhov, A.M., 1954. Basic Laws of Turbulent Mixing in the Surface Layer of the Atmosphere. Trudy Akademii Nauk SSSR Geophizichesky Institut, 24(151), pp. 163-187. (in Russian).
  2. Zav’yalov, D.D. and Solomakha, T.A., 2021. Influence of Thermodynamic Model Resolution on the Simulation of Ice Thickness Evolution in the Sea of Azov. Russian Meteorology and Hydrology, 46, pp. 474-482. doi:10.3103/S1068373921070062
  3. Zavyalov, D.D. and Solomakha, T.A., 2021. Parameterization of Solar Radiation Absorption by Snow-Ice Cover in the Thermodynamic Sea Ice Model of the Sea of Azov. Physical Oceanography, 28(5), pp. 499-513. doi:10.22449/1573-160X-2021-5-499-513
  4. Smith, S.D., 1980. Wind Stress and Heat Flux over the Ocean in Gale Force Winds. Journal of Physical Oceanography, 10(5), pp. 709-726. doi:10.1175/1520- 0485(1980)010<0709:WSAHFO>2.0.CO;2
  5. Banke, E.G., Smith, S.D. and Anderson, R.J., 1980. Drag Coefficient at AIDJEX from Sonic Anemometer Measurements. In: R. S. Pritchard, ed., 1980. Sea Ice Processes and Models. Seattle: University of Washington Press, pp. 430-442.
  6. Zilitinkevich, S.S., Grachev, A.A. and Fairall, C.W., 2001. Scaling Reasoning and Field Data on the Sea Surface Roughness Lengths for Scalars. Journal of the Atmospheric Sciences, 58(3), pp. 320-325. doi:10.1175/1520-0469(2001)058<0320:NACRAF>2.0.CO;2
  7. Andreas, E.L., 1987. A Theory for the Scalar Roughness and the Scalar Transfer Coefficients over Snow and Sea Ice. Boundary-Layer Meteorology, 38(1-2), pp. 159-184. doi:10.1007/BF00121562
  8. Launiainen, J., 1995. Derivation of the Relationship between the Obukhov Stability Parameter and the Bulk Richardson Number for Flux-Profile Studies. Boundary-Layer Meteorology, 76(1- 2), pp. 165-179. doi:10.1007/BF00710895
  9. Li, Y., Gao, Z., Lenschow, D.H. and Chen, F., 2010. An Improved Approach for Parameterizing Surface-Layer Turbulent Transfer Coefficients in Numerical Models. Boundary-Layer Meteorology, 137(1), pp. 153-165. doi:10.1007/s10546-010-9523-y
  10. Högström, U., 1988. Non-Dimensional Wind and Temperature Profiles in the Atmospheric Surface Layer: A Re-Evaluation. Boundary-Layer Meteorology, 42(1-2), pp. 55-78. doi:10.1007/BF00119875
  11. Beljaars, A.C.M. and Holtslag, A.A.M., 1991. Flux Parameterization over Land Surfaces for Atmospheric Models. Journal of Applied Meteorology and Climatology, 30(3), pp. 327-341. doi:10.1175/1520-0450(1991)030<0327:FPOLSF>2.0.CO;2
  12. Chenge, Y. and Brutsaert, W., 2005. Flux-Profile Relationships for Wind Speed and Temperature in the Stable Atmospheric Boundary Layer. Boundary-Layer Meteorology, 114(3), pp. 519-538. doi:10.1007/s10546-004-1425-4
  13. Grachev, A.A., Andreas, E.L., Fairall, C.W., Guest, P.S. and Persson, P.O.G., 2007. SHEBA Flux-Profile Relationships in the Stable Atmospheric Boundary Layer. Boundary-Layer Meteorology, 124(3), pp. 315-333. doi:10.1007/s10546-007-9177-6

Download the article (PDF)