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 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⁻³ and those depending on the atmosphere stratification at the ice geometric roughness equal to 8–10 cm.
Цель. Оценить влияние выбора вида параметризации турбулентного теплообмена на границе океан – атмосфера на основные характеристики ледового режима в вершине Таганрогского залива – цель настоящей работы.
Методы и результаты. Исследование термической сезонной динамики толщины снежно-ледяного покрова проведено с помощью нестационарной термодинамической модели морского льда. В основе алгоритма определения турбулентных потоков импульса, явного и скрытого тепла в региональной модели морского льда лежит полуэмпирическая теория Монина – Обухова. Численные эксперименты проведены для зимних сезонов 2007/08 и 2017/18 годов, метеорологические условия которых существенно различаются. В расчетах применялись как постоянные коэффициенты сопротивления, тепло- и влагообмена, так и полученные с учетом стратификации приземного слоя атмосферы и шероховатости подстилающей поверхности. При реализации устойчивой стратификации использовались три различных выражения для определения функций устойчивости импульса, тепла и влаги. Численное решение уравнений теории подобия основывается на полуэмпирической зависимости коэффициентов турбулентного обмена от объемного числа Ричардсона, что позволяет избежать итерационного процесса. Проведен анализ результатов моделирования и выявлены особенности применения вида параметризации при численном воспроизведении сезонного хода толщины снежно-ледяного покрова.
Выводы. Показано, что при условии продолжительных холодов и наличии ледяного покрова в течение почти всего сезона выбор какого-либо из рассмотренных видов параметризации турбулентного взаимодействия с подстилающей поверхностью не вносит существенных различий в определение максимальной сезонной толщины льда, а также дат замерзания и очищения. Однако при условиях крайней неустойчивости снежно-ледяного покрова зависимость результатов моделирования от способа определения коэффициентов турбулентного обмена может быть весьма заметной. Наиболее удовлетворительные результаты воспроизведения сезонных изменений толщины ледяного покрова найдены при использовании как постоянных коэффициентов турбулентного обмена СН = СЕ ≈ 1,7·10⁻³, так и коэффициентов, полученных с учетом стратификации атмосферы при значении геометрической шероховатости льда 8–10 см.