Термодинамика мелководных озер умеренной зоны играет важнейшую роль в функционировании их экосистем. Температура и скорости течений определяют интенсивность биохимических процессов в водной толще и придонных слоях водоемов. В частности, температура водной толщи и режим перемешивания оказывают определяющее влияние на развитие придонной аноксии [1], на перераспределение клеток водорослей, биогенов, газов по водной толще [2]. Интенсивность динамических процессов может оказывать заметное влияние на жизненные циклы гидробионтов, создавая «благоприятные» условия для видов с разными предпочтениями по динамическому режиму [3]. Таким образом, очевидна актуальность изучения процессов переноса и перемешивания (годовой цикл) в водной толще мелководных озер и выявления их роли в функционировании озерных экосистем.
Слабая изученность зимнего динамического режима водоемов суши связана, с одной стороны, с тем, что скорости течений в период ледостава в малопроточных озерах малы [4], с другой – со сложностью и высокой стоимостью проведения прямых измерений течений в этот период. Скорости подледных течений редко превышают 1–10 мм/с [5], поэтому для их измерений требуются современные приборы с высоким разрешением, что заметно ограничивает возможности постановки натурных экспериментов. В этой ситуации численное моделирование заметно расширяет возможности изучения гидротермодинамических процессов в водоемах, покрытых льдом [6].
Цель исследования: оценить возможность применения метода трехмерного численного моделирования для расчета термогидродинамических (ТГД) процессов в небольшом озере в годовом цикле и, особенно в период ледостава.
Материалы и методы исследования
Трехмерная математическая модель гидротермодинамики внутреннего моря (МГВМ), разработанная в Институте вычислительной математики РАН (ИВМ РАН) и в Институте океанологии им. П.П. Ширшова РАН (ИО РАН) [7] и успешно адаптированная в Институте озероведения РАН (ИнОз РАН) для случая пресноводного озера [8], использовалась для воспроизведения циркуляции в небольшом мелководном озере в зимний период.
В модели для описания циркуляции водной массы и термического режима в бассейне произвольной геометрии используются трехмерные уравнения термогидродинамики. Поверхность раздела воздух – вода свободная, воспроизводятся пространственная изменчивость топографии поверхности и изменчивость среднего уровня водоема. Взаимодействие атмосферы и водной массы описывается через потоки импульса, тепла и влаги. При возникновении условий, благоприятных для формирования льда, включается модель льда, описывающая термодинамические процессы (изменение температуры, нарастание, таяние) во льду и динамические процессы (движение, торошение и т.д.). При этом потоки свойств через границу атмосфера – вода заменяются потоками через границы атмосфера – лед и лед – вода. В модели явным образом описываются потоки воды и ее свойств (минерализация, теплосодержание) через боковые границы (сток и приток рек) и поверхность раздела воздух – вода (испарение, осадки). Кроме того, в модели МГВМ для аппроксимации частных производных использован метод конечных объемов (МКО), основанный на интегральной формулировке законов сохранения массы, импульса, энергии и др. Балансовые соотношения записываются для небольшого контрольного объема; их дискретный аналог получается суммированием по всем граням выделенного объема потоков массы, импульса. Поскольку интегральная формулировка законов сохранения не накладывает ограничений на форму контрольного объема, МКО пригоден для дискретизации уравнений гидродинамики как на структурированных, так и на неструктурированных сетках с различной формой ячеек, что решает проблему сложной геометрии расчетной области.
В условиях, как правило, очень сложной геометрии расчетной области для озер Карелии – как глубоких, так и мелководных – последнее является неоспоримым и решающим преимуществом для выбора модели для расчетов.
Для выполнения расчетов МГВМ была адаптирована для использования ее в озерных условиях. Так, в частности, с учетом того, что минерализация воды в озерах Карелии и их притоках различается и ее значения много меньше морской солености, уравнение состояния морской воды было заменено на уравнение состояния слабоминерализованной воды [9]. Кроме того, в мелководных озерах важную роль в тепловом балансе играет теплообмен с дном. Поэтому МГВМ была дополнена специальным модулем для расчета данного параметра [10].
Объектом численного моделирования было выбрано озеро Вендюрское, расположенное в южной Карелии. Подробное описание озера, его батиметрии, особенностей водосбора приведено в работе [11]. Период ледостава продолжается на озере с середины ноября – начала декабря до начала – середины мая. По классификации Льюиса озеро может быть отнесено к классу полимиктических водоемов, поскольку в период открытой воды его водная масса неоднократно перемешивается. Озеро Вендюрское по морфометрическим особенностям и прозрачности, определяющим режим его перемешивания, является типичным представителем класса мелководных полимиктических водоемов, широко распространенных на территории Карелии.
С использованием данных подробной эхолокационной съемки были выполнены соответствующие работы по построению расчетной сетки для Вендюрского озера, параметры которой составили 80*80 м по горизонтали и 10 горизонтов по вертикали с шагом от 0,5 до 3,5 м. Расчеты выполнялись с шагом по времени 30 с.
В качестве исходной метеорологической информации для расчета «среднеклиматического» состояния ТГД процессов в озере использовались данные реанализа NCEP/NCAR – совместный проект Национальных центров предсказания окружающей среды (NCEP) и атмосферных исследований (NCAR), США [12].
Тестовые расчеты были проведены для одного «среднемноголетнего» года со стартовой датой 1 октября. Выбор даты для постановки начальных условий не был случайным. Дело в том, что в октябре термическая структура озера гомогенна по пространству [11]. Это практически единственный месяц в году, когда начальные условия могут быть заданы корректно.
Результаты исследования и их обсуждение
Анализ результатов численных расчетов показал, что модель адекватно воспроизводит основные особенности термического и динамического режимов озера в годовом цикле при задании среднемноголетнего для района исследований атмосферного форсинга. Даты наступления и окончания ледостава, его продолжительность, годовая амплитуда изменчивости температуры воды на разных глубинах, полученные в результате расчетов, хорошо согласуются с данными натурных исследований [11]. На рис. 1 в качестве примера приведена эволюция температуры воды на разных горизонтах водной толщи в центральной глубоководной части озера в течение одного модельного года (рис. 1, а), а также динамика нарастания и таяния льда (рис. 1, б). Как следует из результатов расчетов, с точки зрения термического режима озеро Вендюрское представляет собой типичный для Карелии и всего северо-запада РФ мелководный водоем со слабо выраженной плотностной стратификацией в период открытой воды. Отчетливо прослеживается весеннее и осеннее полное конвективное перемешивание озера от поверхности до дна, что свойственно подавляющему большинству озер умеренных широт. На этапе летнего нагревания озеро неоднократно полностью перемешивается от поверхности до дна (рис. 1, а). Таким образом, результаты моделирования подтверждают принадлежность Вендюрского озера к классу полимиктических озер.
а) б)
Рис. 1. Температура воды на разных горизонтах водной массы (сплошная линия – поверхность воды, пунктирная – придонный слой) в центральной глубоководной части оз. Вендюрского в течение года (а) и динамика нарастания и таяния льда (б) при задании среднемноголетнего климатического воздействия. Начало расчета – 1 октября
Рис. 2. Поля скоростей течений на глубинах 0,5 м (а), 1,5 м (б), 2,5 м (в) и 5,0 м (г) на 31 января
В озере отчетливо прослеживается так называемый «эффект подледного прогрева» в период ледостава. Данное явление возникает за счет перераспределения в зимний период накопленного донными отложениями тепла между ними и водной массой озера. Воспроизведение моделью подледного прогрева водной массы озера стало возможным благодаря включению в модель специального модуля для расчета теплообмена между водной массой и донными отложениями (см. выше). При этом температура в придонной области повышается до температуры выше значения, соответствующего максимальной плотности (3,98 °С). Данный эффект стал возможен благодаря тому, что в расчетах по МГВМ было использовано уравнение состояния не для пресной, а для слабоминерализованной воды, которое учитывает повышение температуры максимальной плотности за счет минерализации воды.
Также не противоречат данным наблюдений и результаты расчетов образования, нарастания и таяния льда в озере. Так, по данным моделирования первые ледовые явления наступают в озере в начале – середине ноября (рис. 1, б). Нарастание льда происходит до конца марта – середины апреля. Далее наблюдается быстрое таяние ледового покрова с полным его исчезновением в начале – середине мая.
Модельные расчеты течений показали, что вскоре после установления сплошного льда, в озере формируется циклоническая циркуляция, рис. 2, б–г. Эта циркуляция устойчива во времени и существует на протяжении большей части зимы. Скорости течений в поверхностных слоях озера максимальны и достигают 10–15 мм/с в течение первого месяца ледостава, постепенно снижаясь до 0,2–0,5 мм/с в конце зимы. С увеличением глубины скорости течений уменьшаются, достигая 1–1,5 мм/с на глубинах порядка 5 м. На рис. 2 приведены поля течений в приповерхностном слое и на глубинах 1,5; 2,5 и 5 м в середине зимнего сезона (31 января).
Сравнение модельных расчетов скоростей и направлений течений в озере, покрытом льдом, с данными натурных наблюдений [13] показало, что МГВМ удовлетворительно воспроизводит основные особенности динамики водоема в подледный период, а именно, формирование общеозерной циркуляции, существующей на протяжении большей части зимнего сезона и постепенное затухание скоростей течений в течение зимы. Появление общеозерной циркуляции в зимний период может быть обусловлено горизонтальной неоднородностью поля плотности, возникающей за счет разной интенсивности осеннего выхолаживания и теплообмена с донными отложениями зимой на прибрежных мелководьях и в центральной глубоководной части озера.
Масштабный анализ [13, 14] показал, что течения в подледный период в оз. Вендюрском носят геострофический характер. Бароклинный радиус деформации Россби (~200 м) существенно меньше средней ширины озера, а число Россби заметно меньше единицы [14], т.е. эффект вращения Земли может оказывать влияние на формирование крупномасштабной циркуляции. Результаты расчетов на МГВМ показывают формирование циклонического вихря. Натурные измерения, проведенные в разные периоды зимнего сезона, показывают наличие двух циркуляционных ячеек в начале зимнего сезона: циклонической – у северного берега и антициклонической – у южного [13]. Модельное уменьшение скоростей течений к концу зимнего сезона хорошо соответствует наблюдаемому [13] и обусловлено, скорее всего, ослаблением теплообмена водной массы озера с донными отложениями.
Заключение
С использованием трехмерной модели гидротермодинамики внутреннего моря (МГВМ), адаптированной для случая пресного озера, получены расчетные поля скоростей течений и температуры воды на разных глубинах мелководного озера Вендюрского. Результаты расчетов показали, что изменения и дополнения, внесенные в МГВМ, позволяют модели адекватно воспроизводить основные ТГД процессы в озере. Принятые в качестве граничных условий в модели потоки тепла и среднемноголетнее атмосферное воздействие позволили получить расчетные данные о полях температуры и скоростей течений, согласующиеся с данными натурных измерений.
В предстоящих численных экспериментах на адаптированной модели планируется получить схемы течений для основных этапов годового цикла – весенне-летнего нагревания, осенне-зимнего охлаждения, зимнего нагревания, а также провести модельные расчеты для разных гидрометеорологических условий, что позволит выявить основные закономерности формирования и изменчивости структуры течений в озере.
Исследование выполнено в рамках государственного задания Института водных проблем Севера КарНЦ РАН, ФИЦ «Карельский научный центр Российской академии наук».