Бесплатный автореферат и диссертация по наукам о земле на тему

Географические закономерности термического режима разнотипных озер Северо-Запада России

ВАК РФ 25.00.27, Гидрология суши, водные ресурсы, гидрохимия

Автореферат диссертации по теме "Географические закономерности термического режима разнотипных озер Северо-Запада России"

российская академия наук

ИНСТИТУТ ОЗЕРОВЕДЕНИЯ

На правах рукописи

Ефремова Татьяна Владимировна

ГЕОГРАФИЧЕСКИЕ ЗАКОНОМЕРНОСТИ ТЕРМИЧЕСКОГО РЕЖИМА РАЗНОТИПНЫХ ОЗЕР СЕВЕРО-ЗАПАДА РОССИИ

25.00.27 - гидрология суши, водные ресурсы, гидрохимия

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата географических наук *

Санкт-Петербург 2005

Работа выполнена в Институте водных проблем Севера Карельского научного центра РАН

Научный руководитель: доктор географических наук

Шумейко Михаил Арсеньевич

Официальные оппоненты: доктор географических наук,

профессор

Догановский Аркадий Михайлович,

кандидат географических наук Потахин Сергей Борисович

Ведущая организация: Лимнологический институт

СО РАН (Иркутск)

Защита диссертации состоится на заседании

диссертационного совета Д 002.064.01 при Институте озероведения РАН по адресу: 196105 Санкт-Петербург, ул. Севастьянова, д.9

С диссертацией можно ознакомиться в библиотеке Института озероведения РАН.

Автореферат разослан.....апреля 2005 г.

Ученый секретарь Диссертационного совета кандидат биологических наук

В.П.Беляков

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы. Термический режим озер в значительной мере определяет развитие физических, химических и биологических процессов в водоемах, что, в конечном счёте, сказывается на уровне их трофии. В свою очередь, пространственно-временная изменчивость распределения температуры воды в озерах зависит одновременно от зональных климатических особенностей региона, определяющих результирующий тепловой поток в обмене энергией между водоемом и атмосферой, и от озерной морфометрии, влияющей на процессы перемешивания и перераспределения тепла в водной толще.

В настоящее время влияние различных географических факторов на термический режим пресноводных водоемов проанализировано во многих работах. При глобальном рассмотрении озер Земли в основном принимается во внимание зональность инсоляции и обычно исследуется влияние географической широты и высоты над уровнем моря Z на элементы теплового баланса водоемов (Адаменко, 1985), кратность полного перемешивания водных масс в годовом термическом цикле (Hutchinson, Loftier, 1956), средние и экстремальные значения температуры поверхности воды (Рянжин, 1991) и ее годовой ход (Страшкраба, Гнаук, 1989).

При исследовании озер одной климатической зоны главным образом учитывается влияние геометрических размеров и формы строения озерных котловин на термическую стратификацию и перемешивание (Хатчинсон, 1969; Хомскис, 1969), глубины летнего эпилимниона и термоклина (Arai, 1981; Gorfiam, Boyce, 1989; Patalas, 1984), тепловой бюджет (Birge, 1915; Goriiam, 1964; Форш, 1974;Якушко, 1971).

Разделение на глобальный и региональный подходы упрощает анализ зависимости элементов термического режима озер от отдельных географических факторов при выделении приоритетных параметров. На территории Северо-Запада России и Финляндии, промышленно развитых и интенсивно осваиваемых районов, находится более 275 000 озер. Несмотря на огромное количество озер в регионе, регулярные наблюдения за температурой воды проводятся менее чем на 100 водоемах. За последние годы озерная сеть гидрологических постов страны сократилась на 25%, причем в основном за счет средних и малых водоемов озерных районов (Румянцев и др., 2004). Поэтому возникает необходимость оценки и анализа сезонного хода температуры воды на различных глубинах в зависимости от морфометрических и географических особенностей водоёмов. Особую важность решение этой проблемы приобретает для озёр с полным отсутствием гидрометеорологической информации.

Таким образом, актуальной проблемой лимнологии является выявление генетически обусловленных связей и закономерностей, определяющих формирование и развитие термической структуры в озерах, разнообразных по своим морфометрическим характеристикам и расположенных в разных климатических условиях, на основе данных регулярных наблюдений.

Цели и задачи исследования. Цель исследования - установить закономерности изменчивости термического режима разнотипных водоемов Северо-Запада России в зависимости от физико-географических факторов и разработать методы их оценки на основе анализа натурных данных.

Для достижения поставленной цели необходимо решить следующие задачи:

- сформировать ряды данных по температуре воды разнотипных озер Северо-Запада России и их физико-географическим характеристикам;

- получить количественные выражения для условий возникновения термической стратификации в разнотипных озерах в зависимости от их геометрических размеров;

- оценить влияние различных географических факторов на вертикальную термическую структуру разнотипных водоемов и установить регрессионные зависимости между температурой воды на стандартных горизонтах и характеристиками озер;

- оценить влияние географической широты и морфометрических особенностей водоемов на сезонный ход температуры поверхности воды;

- исследовать применимость известных функций для аппроксимации натурных данных годового хода температуры воды на различных горизонтах в разнотипных водоемах;

- исследовать пространственно-временную изменчивость среднемноголетней температуры воды в крупном водоеме региона на примере Онежского озера;

- оценить влияние изменчивости температуры воздуха, как одного из основных климатических факторов для зоны избыточного увлажнения, на вертикальную термическую структуру Петрозаводской губы Онежского озера.

Объектом исследования являются озера Северо-Запада России и Финляндии. Среди крупных озер особое внимание уделено Онежскому озеру, снабжающему водой крупные населенные пункты Карелии, в том числе г. Петрозаводск.

Предметом исследования являются среднемноголетние распределения температуры воды разнотипных водоемов.

В соответствии с задачами определены основные методы исследований -методы корреляционного и регрессионного анализа данных и метод стохастического моделирования.

Научная новизна и основные научные положения, выносимые на защиту:

- на примере разнотипных озер Северо-Запада России и Финляндии впервые разработаны доступные широкому кругу специалистов методы оценки влияния физико-географических закономерностей и озерной морфометрии на годовой ход температуры поверхности воды и на вертикальную термическую структуру водоемов;

- установлено, что скорость нагревания водоемов умеренной зоны в весенне-летний период зависит одновременно от их географического положения и

геометрических размеров, в то время как охлаждение в осенний период — преимущественно от глубины водоёма. Среднемноголетняя температура верхнего пятиметрового слоя озер в период максимального прогрева определяется главным образом зональными климатическими факторами, а на температуру градиентного слоя и гиполимниона значительное влияние оказывают морфометриче-ские характеристики водных объектов;

- предложена аппроксимирующая функция для описания среднемноголет-них ежедневных значений температуры поверхности воды 54 озер Европы и ежедневного хода температуры воды на стандартных горизонтах в глубоководных озерах (Онежском, Инари, Каллавеси);

- благодаря разработанной методике впервые получены ежедневные средние пространственные распределения температуры воды Онежского озера, построены карты-схемы распределения температуры воды на стандартных горизонтах на 1-ое число каждого месяца, количественно оценена интенсивность нагревания и охлаждения различных слоев воды;

- установлены регрессионные зависимости между эмпирическими параметрами модели для температуры воздуха и воды, которые позволяют прогнозировать сезонные изменения вертикальной термической структуры Петрозаводской губы Онежского озера при различных сценариях изменения регионального климата.

Практическая значимость работы. Полученные результаты могут быть использованы для вычисления характерной температуры воды в озерах, на которых нет стационарных наблюдений, а также для совершенствования термических классификаций озер и калибровки боксовых экологических моделей. Разработанные стохастические модели позволяют оперативно проводить оценки термического режима водоемов (по вертикальной термической структуре, сезонным изменениям температуры поверхности воды, толщине эпилимниона), используя при этом ограниченную географическую информацию о водных объектах. В силу своей простоты эмпирические модели доступны широкому кругу специалистов для решения многих научных и практических задач, например, при строительстве водохранилищ или при изменении уровенного и водного режима озер.

Основные результаты получены автором в рамках научно-исследовательских тем, выполнявшихся в Институте водных проблем Севера Карельского научного центра РАН (ИВПС КарНЦ РАН), а также при проведении исследований по проекту РФФИ (грант №02-05-97508р).

Апробация работы. Результаты исследований по теме диссертации были представлены на 5 международных научных конференциях: «Биологические ресурсы водоемов бассейна Балтийского моря» (Петрозаводск, 1991), на Ш международном Ладожском симпозиуме «Мониторинг и устойчивое управление Ладожским озером и другими крупными озерами» (Петрозаводск, 1999); «Экология северных территорий России. Прогноз ситуации, пути развития, решения»

(Архангельск, 2002), на IV международном Ладожском симпозиуме (Новгород, 2002); на VI конференции «Динамика и термика рек, водохранилищ и прибрежной зоны морей» (Москва, 2004).

Работа полностью заслушана на совместном заседании Лаборатории гидрологии и Лаборатории географии и природопользования Института озероведения РАН в ноябре 2004 г.

Публикации. Результаты диссертации представлены в 20 научных публикациях. Диссертация написана на основе материалов многолетних исследований автора. Все работы, результаты которых приведены в диссертации, выполнены лично автором или вместе с соавторами.

Структура и объем работы. Диссертация состоит из введения, четырех глав и заключения, изложена на 136 страницах, включает 24 рисунка, 5 таблиц, список использованной литературы из 117 наименований, в том числе 36 на иностранных языках.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Во введении обоснована актуальность работы, сформулированы цели и задачи исследования, показана научная новизна и практическая значимость полученных результатов.

Глава 1. Влияние географических факторов и морфометрических характеристик на термическую структуру озер

Состояние озера как географического объекта в первую очередь определяется географическим положением (географические координаты, высота над уровнем моря и т.п.). На влияние физико-географической зональности на термический режим пресноводных водоемов впервые обратил внимание основоположник лимнологии Ф. Форель и предложил термическую классификацию (Forel, 1892). Позже термические классификации, основанные на принципе физико-географической зональности, уточняли, дополняли и развивали многие исследователи (Runner, 1931; Yoshimura, 1936; Hutchinson, Lofiler, 1956; Zafar, 1959; Lewis, 1983; Рянжин, 1989 и др.).

В пределах одной климатической зоны определяющую роль в формировании термического режима водных масс имеют морфометрические особенности водоема (Gorham, 1964; Хомскис, 1969; Форш, 1974; Драбкова и др., 1979; Изотова, 1982), что нашло отражение во многих региональных термических классификациях. Для умеренного климатического пояса это классификации П.Ф. Домрачева, А.А. Браславского и З.А. Викулиной, В. Хомскиса, И.С. Захаренкова, О.Ф. Якушко, А.И. Пехович и В.М. Жидких, Л.В. Несиной и Т.А Огневой, СП. Китаева и др. Наибольший интерес представляет класси-

фикация А.И. Тихомирова, в которой термический режим водоемов рассматривается в годовом цикле в зависимости от средней глубины.

Рассматриваемая территория расположена на северо-западе Европейской части России между 55°-69° с. ш. (протяженность около 1450 км) и 27°- 42° в. д. (около 800 км). В административном отношении она включает республику Карелия, Мурманскую, Ленинградскую, Новгородскую, Псковскую области, а также часть Архангельской, Вологодской и Калининской областей (рис.1).

20 30 40

Рис.1. Карта-схема расположения исследованных озер

Климатические особенности исследуемого региона определяются географическим положением территории в субарктических и умеренных широтах. Географическое положение региона и связанное с ним количество поступающей солнечной радиации, близость Баренцева и Балтийского морей, преобладание западного переноса воздушных масс атлантического и арктического происхождения обусловливают интенсивную циклоническую деятельность. Территория относится к зоне избыточного увлажнения и имеет хорошо развитую гидрографическую сеть.

Для изучения сезонной изменчивости среднемноголетнего хода температуры поверхности воды использованы данные по 52 озерам Европейской территории России, собранные за многолетний период (1945-1980 гг.) на сети станций и постов Росгидромета и данные по двум озерам Финляндии, обобщенные государственной системой гидрометеорологических наблюдений Финляндии (1961-1991 гг.). Для изучения среднемноголетней вертикальной структуры использованы данные о температуре воды, собранные гидрологическими постами с 1958 по 1989 г. на рейдовых вертикалях 60 водоемов России и 8 озер Финляндии (1961-1975 гг.)

Выборка рассматриваемых в работе водоемов представляет большое их разнообразие по морфометрическим характеристикам, формам строения озерных котловин и проточности. По геометрическим размерам в нее включены все типы водоемов: от крупнейших - Ладожского и Онежского (с площадью зеркала 17872 км2 и 9777 км2) до самых малых с площадью менее 1 км2. В выборку входят крупные и глубоководные финские озера (Инари, Каллавеси, Паяние, Пиели-нен) и озера Северо-Запада России (Имандра, Пяозеро, Топозеро, Сегозеро); крупные, но мелководные (Выгозеро, Чудское, Ильмень). Больше половины водоемов по площади относятся к разряду «средних» или «малых», т. е. их площади меньше 100 км2 и 10 км2, соответственно. Диапазон некоторых характеристик исследованных озер приведен в таблице 1.

Таблица 1

Характеристики исследованных озер

Для исследуемых озер получена зависимость стратификации водоемов от максимальной глубины ( Н^^ м) и площади (5, км2) (Ефремова, Палыпин, 2003). Для проведения границ между гипо-, эпи- и метатермическими типами водоемов (рис. 2) согласно терминологии из классификации А. И. Тихомирова (1982) использовалась безразмерная температура 0:

Глубокие озера (гипотермические - © < 0,3) - типичные представители димиктических озер с наличием по вертикали трех слоев: эпи-, мета- и ги-полимниона. Водные массы мелких озер (эпитермические - 0 > 0,7) летом однородны, или из-за неустойчивой термической стратификации периодически перемешиваются до дна под воздействием ветров или из-за охлаждения воды. Озера со средними глубинами (метатермические - 0,3 < 0 < 0,7) обычно имеют два слоя: верхний квазиоднородный и нижний стратифицированный. Придонный квазиоднородный слой в них, как правило, отсутствует или имеет небольшую толщину и наблюдается не каждый год.

Из рис. 2 следует, что для определенной площади озера должна быть минимальная глубина, при которой оно будет стратифицировано; если площадь озера увеличивается, то минимальная глубина также увеличивается.

Количественные соотношения между длиной разгона ветра и глубиной, при которых озера должны быть стратифицированы, получены Р. Латропом и Р. Лили (Lathrop, LШie, 1980), Е. Горхам и Ф. Бойсом (Gorham, Boyce, 1989). В отличие от предложенного выше разделения озер на три класса Латроп и Лили, Горхам и Бойс делили озера на две группы: стратифицированные и

нестратифицированные, поэтому границы разделения озер не совпадают. Е. Гор-хам и Ф. Бойс дали теоретическое обоснование различий в режимах ветрового перемешивания малых (Ь < 5 км) и больших ^ > 5 км) озер, основываясь главным образом на работе Р. Шпигеля и Й. Имбергера (Spigel, bnberger, 1980).

К сожалению, для озер Северо-Запада России не представляется возможным оценить глубину перемешанного слоя, так как наблюдения за температурой воды проводились с разной дискретностью по глубине: в глубоких озерах обычно через 5 м, в средних - через 2 м. Поэтому была рассмотрена возможность применения формул К Паталаса (РаЫаз, 1984) к нашему региону. По его данным, собранным для разных климатических зон (170 озер Северной Америки, Польши и Японии), была рассчитана глубина залегания эпилимниона (Ие) не только в зависимости от длины разгона ветра (Ь), но еще от <р. Поправки к Ие от (р для его формулы оказались не большими, но значимыми. В озерах разных геометрических размеров преобладают различные механизмы ветрового перемешивания. Для того чтобы оценить, как это отражается на глубине эпилимниона в зависимости от длины разгона ветра, выборка озер была разбита на три группы: малые озера ^ < 5,5 км), малые и средние озера (К 15 км) и все 170 озер (включая крупные с Ь < 33 км и площадью до 1000 км2).

Для различных групп озер получены близкие значения регрессионных коэффициентов и свободных членов уравнений, что свидетельствует об устойчивости полученных зависимостей и позволяет использовать при расчетах единое уравнение для всей выборки озер:

Не = 4,83 + 0,119 <Р 1§ Ь (±1,56; Л= 0,88;)

В скобках приведено стандартное отклонение и коэффициент корреляции. Таким образом, эту формулу можно использовать для расчетов Ие одновременно как для малых, так и для крупных озер, расположенных в диапазоне от 45° до 66° с. ш. (в том числе для нашего региона).

Глава 2. Средняя многолетняя вертикальная термическая структура озер в период максимального прогрева

Вертикальную термическую структуру любого озера на конкретную дату можно представить в виде детерминированной и стохастической составляющих:

Т(г,с1)= Т(г,(1) + АТ(х> ¿0,

где сГ) - средняя температура в озере на горизонте z за конкретную дату й, АТ(г, с1) - сумма поправок на пространственную и временную изменчивость температуры:

Средняя температура в озере на конкретную дату Т{2,<Х) является климатической характеристикой и зависит от зональных факторов (географическая широта, высота над уровнем моря) и морфометрических характеристик самого озера.

Для исследования вертикальной термической структуры озер Северо-Запада выбран период с прямой термической стратификацией (20 июля - 10 августа) при максимальном накоплении тепла в озерах. Вертикальная термическая структура озер оценивалась на основании измерения температуры воды на стандарт-

ных горизонтах (0.1,2,5,10,15,20,25,30,40 м). Ряды измеренных в озерах значений Т(т) неоднородны и имеют длину от 5 до 27 лет, для большинства озер они составляют около 15 лет. Первоначальные данные Дг) осреднялись по каждому горизонту и озеру с учетом всех измерений. Результаты предварительного статистического анализа свидетельствуют об отсутствии зависимости между геометрическими размерами озер и ф, показателем условного водообмена (коэффициенты парной корреляции г< 0Д5). Между отдельными морфометрическими характеристиками озер существуют довольно тесные связи. Например, между средними Н и максимальными Нтах глубинами озер, или между их площадью 5 и объемом Г наблюдаются явные линейные зависимости. Однако, связь между 5 и Н или для водоемов данной выборки невысокая.

Для статистической обработки данных были применены методы нелинейного оценивания Ньютона и множественного регрессионного анализа. Задача заключалась в том, чтобы выбрать основные географические факторы, оказывающие влияние на температуру воды в озерах, и получить наиболее простое выражение с ограниченным количеством предикторов, которые широко используются в лимнологии и их количественные характеристики известны для многих озер.

В результате проведенного анализа в качестве приоритетных характеристик были выбраны: географическая широта, высота над уровнем моря, площадь озера, средняя глубина, коэффициент условного водообмена (Палыпин, Ефремова, 2002; Ефремова, Пальшин, 2003). Для крупных глубоководных озер в связи с неоднородностью поля температуры по пространству была добавлена глубина в точке измерения, а для учета влияния разгона ветра на вертикальную структуру - коэффициент открытости.

Следующий этап исследований заключался в нахождении оптимального регрессионного уравнения, описывающего связи температуры воды с выбранными географическими характеристиками. Было апробировано около 50 уравнений с различными комбинациями выбранных предикторов и проверкой регрессионных

моделей на устойчивость и_адркватность. Найдено, что осредненные за многолетний период значения Т(г) на различных горизонтах для периода летней стратификации в озерах наилучшим образом оцениваются с помощью следующего регрессионного уравнения:

Г(г) = а0 + а,(72-<р)+аг2 + а3 +а4(1§я)2 + а5й, + а6 ^К + ауЕ (1)

где Z - высота над уровнем моря; , м - глубина озера в точке измерения; ЯГ, год - коэффициент условного водообмена;. , км2/м - коэффициент открытости.

Таблица 2

Коэффициенты множественной корреляции R, стандартные отклонения температуры воды Е И коэффициенты а0 — £К7 для различных горизонтов для

уравнения(1)

м Л £°С ао, °С Ль °С/° с.ш. аь °С/м °С 04, °С сч, °С/м «6, "С Л7, /м)

0 0.94 0.78 14.16 0.564 -0.0048 -0.12 -0.31 •0.0038 0.095 -0.0058

2 0.93 0.78 13.52 0.554 -0.0048 -0.04 -0.31 -0.0040 0.103 -0.0059

5 0.86 1.09 11.69 0.553 -0.0033 0.52 -066 -0.0118 0.364 -0.0054

10 0.83 1.36 8.22 0391 0.0025 2.14 -125 -0.0349 1.663 -0.0049

15 0.86 1.09 5.70 0.227 0.0085 2.45 -1.45 ■«0292 2.253 -0.0039

20 0.89 0.86 4.94 0.214 0.0097 225 -1.50 -0.0202 2.141 -0.0049

25 0.92 0.70 2.78 0.261 0.0112 2.00 -139 •0.0062 1.992 -0.0058

30 0.95 0.63 2.22 0.282 0.0115 1.72 -1.32 -0.0002 1.934 -0.0059

40 0.97 0.53 1.73 0.280 0.0115 1.52 -1.11 -0.0012 1.754 -00058

Коэффициенты а0 — Оу для различных горизонтов и их размерность приведены в табл. 2, где также показаны коэффициенты множественной корреляции Я и стандартные отклонения температуры е. Из табл. 2 видно, что значения Л достаточно высокие, а е < 10С для эпи- и гиполимнионов и немного выше для градиентного слоя.

Коэффициенты Л0 — Л7 плавно изменяются с глубиной. Это свидетельствует об устойчивости подобранной регрессионной модели при описании вертикальной термической структуры. Изменение свободного члена а0 с глубиной по форме представляет собой осредненный для озер вертикальный профиль температуры воды, но абсолютные значения его несколько ниже. Параметры с коэффициентами следует рассматривать как поправки к нему для уче-

та зональных факторов и озерной морфометрии.

Анализ адекватности модели показал, что коэффициенты множественной корреляции достаточно высокие (0,86-0,97), соответственно высокие и коэффициенты детерминации, характеризующие объясненную долю дисперсии и свидетельствующие, насколько хорошо результаты моделирования согласуются с натурными данными. Стандартные отклонения < 1°С для эпи- и гипо-лимнионов и немного выше для металимниона.

С помощью табл. 2 и формулы (1) можно оценить влияние каждого фактора на температуру воды на различных горизонтах в озерах. Это особенно важно при сравнительном анализе их термических режимов. На горизонтах 0-5 м температура воды в первую очередь зависит от широты. С увеличением глубины на горизонтах зависимость от уменьшается, а от морфо-метрических показателей водоемов и их проточности возрастает. Увеличение 5 приводит к слабому уменьшению Т(т) на горизонтах 0-2 м и к существенному росту Т(г) на глубине > 10 м ( СК3 ). При увеличении н в озерах уменьшаются значения Т(г) на всех горизонтах ( ), причем отрицательная поправка к Т\г) больше для глубинных слоев. При глубине на рейдовой вертикали 100 м поправка к Т(г) на горизонтах 0 и 2 м составляет минус 0.4°С, а в слое скачка (10-20 м) увеличивается по абсолютной величине до 2-3 °С ( ¿3 ). Для проточных озер наблюдается небольшое повышение температуры воды в верхнем слое (0-2 м) и значительное увеличение ее в слое скачка и гиполимнионе ( а6 ).

Для всех исследуемых озер построены вертикальные профили температуры воды по осредненным за многолетний период измеренным значениям на различных горизонтах и рассчитанные по уравнению (1). Пример для трех озер показан на рис. 3.

Рис. 3. Рассчитанные по (1) среднемноголетние профили температуры воды в озерах (штриховые кривые) и данные наблюдений (сплошные)

а - Топозеро, б — Онежское озеро, в - Валдайское.

Для проверки регрессионной модели (1) выполнены расчеты вертикальной термической структуры в озерах, данные по которым не входили в выборку, а именно Панаярви и Плещеево. Предсказанные по регрессионной модели значения температуры воды на различных горизонтах озер достаточно близко соответствуют из-меренным_величинам. Расхождения между наблюденными и рассчитанными значениями T(z) до 1-1,5°С отмечены для придонных слоев этих озер (в оз. Плещеево — на горизонтах 15-20 м, а в оз. Панаярви - на горизонтах 25-40 м). Они могут быть связаны как с тем, что наблюдения в этих озерах характеризуют не средне-многолетние значения, а относятся к конкретным годам, так и с неучтенными в регрессионной модели особенностями строения форм озерных котловин.

Глава 3. Средний многолетний годовой ход температуры поверхности воды

Первое описание изменения температуры воды в озере с помощью непрерывной функции времени было сделано на примере оз. Мендота (Birge, 1929). Позже М. Страшкраба (Strashkraba, 1980) показал, что в глобальном масштабе общие закономерности зональной изменчивости среднегодовых вариаций температуры поверхности водоемов «средних размеров» хорошо описываются гармонической функцией в зависимости от географической широты и высоты над уровнем моря. Нагревание озер умеренной зоны обычно происходит быстрее, чем их охлаждение, поэтому одним из важных условий для аппроксимирующей функции является асимметричность возрастающей и убывающей ветвей. Такой годовой ход температуры водоемов можно представить в виде склейки двух косинусоид (Румянцев и др., 2001) или композиционными гармоническими функциями (Страшкраба, Гнаук, 1989). Моделирование среднего многолетнего цикла температуры воды непериодическими функциями приводится для Великих озер Америки (Lesht, Brandner, 1992) и Европы (Ефремова, Петров, 1992,2000; Науменко и др., 2000; Науменко, Каретников, 2002).

Для моделирования годового цикла температуры поверхности исследуемых водоемов были апробированы различные функции. Наилучшим образом наблюденные данные температуры поверхности замерзающих озер Северо-Запада аппроксимируются непериодическими временными функциями (Lesht, Brandner, 1992) и (Ефремова, Петров, 1992). В расчетах использовалась формула (Ефремова, Петров, 1992), в которую дополнительно введен аддитивный член (b0), что существенно улучшило приближение к натурным данным при низких значениях температуры воды ранней весной и поздней осенью:

0 '{ l+cxrKtf-iJbJJ [ l+exrKrf-A.felr'

где Т(ф - температура поверхности озера в данный момент времени (°С); й- дата в днях от 1 января; Ьо — Ь}- эмпирические параметры.

Эмпирические параметры модели для каждого конкретного озера подбирались итерационным квази-ньютоновским методом нелинейного оценивания. Параметры 2>о и Ь1 связаны с минимальными и максимальными значениями температуры, но численно не равны им. Безразмерные параметры Ь; И Ъз определяют форму аппроксимирующей кривой, т.е. крутизну возрастающей и убывающей ветвей функции, а их численные различия между собой приводят к асимметрии. Два последних параметра связаны с точками перегиба функции и соответствуют датам максимальных скоростей роста температуры воды весной ( ) и понижения осенью ( ).

Параметры модели (2), полученные для всех водоемов, исследовались в зависимости от географических координат, озерной морфометрии и показателя условного водообмена. Наилучшим образом связи оцениваются единым для всех параметров регрессионным уравнением:

Ь, ~ал + ал(72-ф)+аа2+ав +аа(1ёя)2 +ал\%К+алР, (3) /=0,1,...5

где - коэффициент «емкости». Регрессионные коэффициенты

для расчета модельных параметров даны в табл. 3, где также приведены коэффициенты множественной корреляции.

Таблица 3

Регрессионные коэффициенты аЛ — а1б к формуле (3) и коэффициенты

корреляции Я

Ь, Л «Л ал в.2 «.3 О,А «.5 а,6

Ьо,°С 0,96 -4,5963 -0,0311 -0,00144 0,984 -0,507 -0,242 -1Д7

Ьь°С 0,99 53717 0,134 -0,00149 -0,347 0,647 0,129 0,112

Ььсуг 0,84 271,68 0,842 -0,0287 -3,863 0,138 -1,130 -5,09

Ьз 0,94 -0,0429 0,00071 -0,00002 -0,00191 0,00922 0,00119 -0,0007

Ь«,суг 0,95 148,54 -1,718 •0,0000 2,662 11,02 1,607 -1,282

Ь5 0,87 -0,0588 0,00139 -0,00003 -0,0026 0,00683 -0,00008 -0,0128

Примечание Размерность регрессионных коэффициентов и й^ —соответствует размерности параметров Ъ,, а в размерности О,, и Оп дополнительное с ш)"1 ил*"1

В результате регрессионного анализа установлено, что наибольшее влияние на средний годовой ход температуры поверхности воды оказывают географическая широта, площадь озер и их средняя глубина. При этом минимальная температура (параметр ) зависит главным образом от площади озер, максимальная

температура (параметр Ы) - от географической широты, а, например, скорость понижения температуры в осенний период (параметр Ьз) определяется преимущественно глубиной озер. В меньшей мере проявляются зависимости параметров модели от условного водообмена и форм строения озерных котловин (коэффициент «емкости»). Однако, не учет водообмена озер приводит к погрешности температуры воды в период весенне-летнего прогревания до ± 0.9°С, а коэффициента «емкости» ранней весной и в период осеннего охлаждения - до ± 0.6°С. Влияние высоты над уровнем моря для данной выборки озер, расположенных от О до 200 м также небольшое, но значимое.

По уравнению (3) вычислены количественные значения параметров Ь( для всех озер, по которым были построены модельные кривые годовых циклов температуры воды. Наилучшее соответствие между ними и измеренными значениями температуры характерно для средних и малых водоемов. Для Ладожского и Онежского озер получены модельные кривые и осредненные натурные данные по акваториям в целом и трем разным зонам (центральной, промежуточной и прибрежной). Это вызвано тем, что поля температуры в крупных озерах характеризуются значительной пространственной неоднородностью (Бояринов и др., 1999; Науменко, Каретников, 2002). Для крупных озер (Выгозеро, Сегозеро, Топозеро, Имандра, Инари и др.) различия между модельными кривыми и натурными данными могут достигать 1°С (Пальшин, Ефремова, 2003). Для примера на рис. 4 показаны модельные кривые крупного Ладожского озера и небольшого оз. Глубокого.

2« - --- -------

5-20 •

^ 18

(О О.

а. 2 ' г

£ 4

о

60 120 180 240 300 в0 120 180 240 300 360

Дни

Рис. 4, Ежедекадные среднемноголетние данные температуры поверхности воды (точки), модельные кривые с наилучшей подгонкой уравнения (2) — жирные линии и модельные кривые, полученные при расчете параметров Ь( по уравнению (3) — тонкие линии для а) - Ладожского озера (1 - прибрежная; 2 - промежуточная; 3 - центральная зоны; 4 — акватория всего озера), б) - для озера Глубокого

Максимум температуры поверхности воды во всех озерах наблюдается приблизительно в одно время (середина-конец июля). Форма ветвей кривых температуры поверхности воды озер, расположенных южнее 60° с. ш., близка

к симметричной. В приполярных озерах наблюдается большая асимметричность кривых, возрастающая ветвь более крутая. Это вызвано тем, что процесс прогревания в северных широтах происходит более интенсивно в связи с круглосуточным поступлением солнечной радиации в июне-июле. Наибольшая асимметрия характерна для глубоководных районов крупных озер (Ладожское, Онежское, Сегозеро), в которых период охлаждения на один-два месяца продолжительнее периода нагревания.

о

Наблюденная температура (°С)

Рис. 5. Диаграмма рассеяния натурных данных температуры поверхности воды и значений, (а) - рассчитанных при наилучшей подгонке уравнения (2), (б) - при расчете параметров Ь, по уравнению (3) для 47 озер.

Для 47 озер по параметрам Ь,, вычисленным по уравнению (3), с помощью модели (2) рассчитаны декадные значения температуры воды. На диаграмме рассеяния наблюдаемых и предсказанных значений все точки расположены вдоль прямой линии (рис. 5). Это свидетельствует о том, что модель учитывает асимметрию ветвей роста и понижения температуры воды и одинаково хорошо описывает годовой цикл, как при высоких, так и при низких значениях температуры. Стандартные отклонения между измеренными значениями температуры воды и полученными в результате расчетов при наилучшей подгонке уравнения (2) составляет 0.27°С, а при вычислении параметров Ь1 по географическим характеристикам озер 0.7°С. Такая точность делает предложенную модель полезной для описания "типичного" годового температурного цикла.

Верификация стохастической модели осуществлялась по годовому ходу температуры воды в семи разнотипных озерах, которые расположены на рассматриваемой территории от 56° до 67° с. ш. Эти озера (Ондозеро, Кимасозеро, Нюк, Совдо-зеро, Ковдозеро, Торбино, Охват) не были включены в основную выборку. Для них по уравнению (3) были определены параметры Ь0 в соответствии с их географическими характеристиками, а по модели (2) построены годовые термические

циклы. Результаты верификации показали хорошее соответствие между модельными расчетами и натурными данными, особенно по малым (Совдозеро, Торбино, Охват) и средним (Ондозеро, Кимасозеро) озерам.

Глава 4. Пространственно-временная изменчивость среднемноголетней температуры воды глубоководных озер

Глава посвящена стохастическому моделированию годового хода вертикальной термической структуры глубоководных озер. Для Онежского озера проанализирована пространственная и временная изменчивость вертикальной термической структуры за период открытой воды.

Для озер Инари и Каллавеси подробно рассматривается проверка модели на устойчивость и адекватность, а также проанализировано поведение эмпирических параметров модели в зависимости от глубины. Выбор этих озер связан с наличием для них подробных данных по вертикальной термической структуре - до глубины 20 метров через один метр, а глубже - через два метра до глубины 44-46 м. Наблюдения проводились регулярно один раз в декаду в течение 1981-2001 гг. Среднемноголетний годовой ход температуры воды в Инари и Каллавеси существенно различается, что связано с географическим положением и морфомет-рическими особенностями: оз. Инари - более крупный и глубокий водоем, расположенный севернее на 6° с.ш. (Korhonen J., 2002).

Для вычисления среднемноголетнего хода температуры воды на различных горизонтах использовалась формула (2). Эмпирические параметры модели для каждого конкретного горизонта подбирались итерационным квази-ньютонов-ским методом. На рис. 6 представлены модельные кривые для нескольких горизонтов. Модель описывает зимнюю стратификацию, переход температуры воды через 4°С, весенний рост и осеннее понижение температуры. Модельные кривые позволяют определить на различных горизонтах среднемноголетние даты начала весеннего нагревания, осеннего охлаждения, прохождения весеннего и осеннего термобара, даты наступления осенней гомотермии, максимальной температуры, продолжительность зимней и летней стратификации.

Для обоих озер проанализированы значения эмпирических параметров модели, отмечено, что они в основном плавно изменяются по глубине. Так как

Т'С

Дни

Рис. 6. Средняя многолетняя температура воды на горизонтах 1, 5, 10, 15, 20, 40 м а) оз. Инари; б) оз. Каллавеси.

расчеты по горизонтам выполнялись по независимым выборкам, это может свидетельствовать об устойчивости применяемой модели. Числовые значения Ьо практически соответствуют зимнему профилю температуры. Вертикальное распределение параметра Ь1 по форме очень похоже на летний профиль температуры, но числовые его значения существенно меньше. Параметры Ьо и Ь^ характеризуют диапазон изменения температуры воды на различных горизонтах от минимальных значений зимой до максимальных в период летнего прогрева. Параметры определяют формы аппроксимирующих кривых, т.е. ост-

рую или пологую вершину и асимметрию ветвей функции. В работе анализируется распределение по глубине каждого из параметров.

Стандартные отклонения вычислялись для каждой декадной группы наблюдений по всем горизонтам для обоих озер. Коэффициенты корреляции при аппроксимации годового хода для всех горизонтов больше 0,99. Коэффициенты детерминации, также близки к единице и свидетельствуют о хорошем приближении модельных кривых к наблюдаемым данным и о качестве построения модели.

Во втором разделе четвертой главы основное внимание уделено сезонной пространственной изменчивости вертикальной термической структуры Онежского озера.

Пространственная неоднородность полей температуры по акватории озера в значительной степени может быть охарактеризована температурой приповерхностного слоя, поэтому на первом этапе основное внимание было уделено определению сезонного хода среднемноголетней температуры поверхности воды (ТПВ). Для характеристики этого показателя использовались данные ИК-съемок за 1972-1984 г. В зависимости от условий теплонакопления и теплоотдачи в озере было выделено шесть районов (Ефремова, 1990).

С учетом районирования по ТПВ для характеристики среднемноголетнего вертикального термического режима Онежского озера были выбраны 14 наиболее репрезентативных рейдовых вертикалей. Исходными материалами послужили данные измерений за период с 1958 по 1989 г. Измерения проводились в навигационный период на рейдовых вертикалях с борта судна на постоянной сетке станций на стандартных горизонтах (0,1; 2; 5; 10; 15; 20; 30; 50 м; дно) с дискретностью 10-30 дней.

Для вычисления среднемноголетнего хода температуры воды на различных горизонтах использовалась непрерывная временная параметрическая функция (2), которая была впервые получена и применена для характеристики термического цикла Онежского озера (Ефремова, Петров, 1992; 2000). Данная аппроксимирующая функция позволяет получить ежедневные средние значения температуры воды на стандартных горизонтах в различных районах Онежского озера. По полученным модельным кривым построены карты-схемы распределения температуры воды на стандартных горизонтах на 1-ое число каждого месяца с июня

по ноябрь. Основное достоинство этих карт заключается в возможности наглядно проследить изменчивость поля температуры в различных районах Онежского озера на разных глубинах для "типичного" года.

Разрушение льда в Онежском озере обычно начинается в первых числах мая. Переход температуры воды через 4°С сопровождается формированием термического бара. По мере нагревания воды фронт термобара постепенно перемещается в область больших глубин. 1 июня изотерма 4°С оконтуривает озеро и располагается над глубинами 20-25 м. ТПВ в это время характеризуется наибольшей пространственной неоднородностью.

Термический бар в Онежском озере по среднемноголетним данным исчезает в третьей декаде июня. В центральной глубоководной части озера средняя скорость повышения температуры воды в это время в поверхностном слое 0-5 м возрастает в 3,5-4 раза и составляет приблизительно 0,35-0,40°С в сутки. На горизонте 20 м она близка весеннему значению - около 0,1 °С в сутки, а на глубинах более 40 м составляет менее 0,03 °С в сутки. На прибрежных станциях во второй половине июня, когда в центре озера еще только происходит переход температуры воды через 4°С, скорости роста температуры уже уменьшаются до 0,15-0,25°С в сутки.

1 июля минимальная температура в приповерхностном слое наблюдается в центральной части озера — 9°С, максимальная - в заливах и губах (15°С в Заонежском заливе). В это время сезонный слой скачка в глубоководной части озера располагается на глубине 5-10 м. Температура воды в нем изменяется от 8,5 до 5,5°С.

В июле отмечается постепенное выравнивание температуры поверхности воды по всей акватории озера. В это время происходит переход температуры воды в поверхностном слое через 10°С на всей акватории озера, наступает "биологическое лето". Несмотря на более северное положение, в первую очередь "биологическое лето" наступает в мелководном Заонежском заливе (в первых числах июня). Запаздывание наступления "биологического лета" для центральной глубоководной части озера составляет в среднем 30 дней, а, например, для Ладожского озера - 50 дней (Науменко, Каретников, 2002). Продолжительность "биологического лета" для приповерхностного слоя воды варьирует от 118 суток в Заонежском заливе до 83 суток в центральной глубоководной части озера, верхнего 10 метрового слоя - от 105 суток в Заонежском заливе до 67 в глубоководном Большом Онего.

Конец июля - начало августа характеризуются наибольшим прогревом верхнего 5-метрового слоя воды. Средняя температура поверхности воды достигает своих максимальных значений - 16-17°С, а ее распределение по акватории становится довольно однородным. К 1 августа сезонный термоклин опускается на глубину 15-25 м в зависимости от района озера. На горизонте

50 м температура воды в центральном плесе составляет приблизительно 4,5°С,авгубах-5,5°С.

Следует отметить, что продолжительность процесса нагревания и охлаждения верхнего слоя воды (до 15 м) для прибрежных районов приблизительно одинакова. Так, например, в Петрозаводской губе продолжительность нагревания от 4°С до максимума составляет 90 суток, а охлаждения до 4°С - 95 суток; в центральной части озера процесс нагревания поверхностных слоев происходит приблизительно в 1,5 раза быстрее, чем охлаждение. В слое воды глубже 15 м, ниже залегания сезонного термоклина, процесс нагревания длится в 2-2,5 раза дольше, чем охлаждение.

В первой декаде августа начинается период осеннего охлаждения. Температура квазиоднородного слоя воды уменьшается в среднем на 0,10-0,13°С в сутки. Скорость заглубления перемешанного слоя с августа по октябрь постепенно возрастает с 0,3 до 3,0 м-сут Толщина квазиоднородного слоя с температурой 13°С на 1 сентября составляет 10 м. К 1 октября этот слой увеличивается до 20 м, а температура воды уменьшается до 9,5°С. В южной части центрального плеса 1 октября на глубинах 20-30 м по-прежнему отмечаются более низкие температуры.

К концу октября в озере устанавливается гомотермия с температурой 6-8°С. Охлаждение озера происходит с постепенным понижением температуры во всей толще вод, приблизительно со скоростью 0,10С-сут -1. В двадцатых числах октября на горизонте 80 м наблюдается максимальная температура воды - 6,6°С. Таким образом, запаздывание наступления максимума температуры на нижних горизонтах по сравнению с верхним слоем в Онежском озере достигает трех месяцев. 1 ноября в центральном плесе наблюдается гомотер-мия с температурой выше 6°С. В прибрежных районах, с глубиной до 20 м температура воды, как правило, ниже 5°С. У берега вода в это время охлаждается до температуры максимальной плотности 4°С, и возникает осенний термический бар.

В заливах и губах термобар обычно проходит в первой половине ноября, в центральном плесе в районах с глубиной до 30 м - в середине ноября, 50 м -приблизительно 20 ноября. Исчезает термобар 22-25 ноября. 1 декабря в прибрежных районах температура воды составляет 1,9-2,2°С, а в центральном плесе - более 3°С.

Разработанная методика позволила впервые оценить пространственное и вертикальное распределение среднемноголетней температуры воды Онежского озера на любую дату. Полученные многолетние характеристики отражают многообразие температурных условий в крупном глубоководном озере и позволяют выявить особенности формирования температурного режима в разных частях его акватории, количественно оценить интенсивность нагревания и охлаждения отдельных слоев,

определить сроки, их максимального прогревания. Нивелируя межгодовую и синоптическую изменчивость, среднемноголетняя температура воды полезна для понимания процесса формирования термической структуры озера, а также для моделирования "типичного" года.

В третьем разделе оценивается зависимость температуры воды на различных глубинах от температуры воздуха (на примере Петрозаводской губы Онежского озера).

В стохастических моделях температура воды в озерах наилучшим образом связана с температурой воздуха (Robertson, Ragotzkie, 1990; Hongzo, Stefan, 1993). Среднемесячные ее значения достаточно хорошо воспроизводятся моделями глобального изменения климата в отличие от других метеорологических элементов. Поэтому в работе рассмотрена зависимость температуры воды на различных горизонтах от температуры воздуха (одной из основных характеристик климата) на примере Петрозаводской губы Онежского озера (Ефремова, Палынин, 2004).

Для этого использовались характерные для разных типов лет среднемно-голетние кривые годового хода температуры воды на различных горизонтах (по данным в центральной части губы) и соответствующие им значения температуры воздуха по метеостанции Петрозаводск за этот период. Выборка лет с разнообразными периодами (жаркими, теплыми, средними, прохладными, холодными) для года в целом, весны, лета и осени осуществлялась по температуре воздуха. Таким образом, удалось получить 20 осредненных наиболее характерных годовых циклов по температуре воздуха и воды на различных горизонтах.

Для разных характерных типов лет получены результаты моделирования годовых циклов температуры воздуха и воды с использованием единого уравнения (2). Стандартные отклонения между наблюденными значениями температуры воздуха и полученными в результате расчетов по модели составляют около 0,3°С.

На рис. 7 показан совмещенный график модельных кривых по температуре воздуха и воды на различных горизонтах для среднего года, полученных с использованием осредненных данных за весь период наблюдений. Были исследованы зависимости эмпирических параметров модели Ьо — Ь5 от аналогичных параметров для температуры воздуха. Установленные эмпирические зависимости позволяют рассчитывать температуру воды в Петрозаводской губе при наличии данных по годовому ходу температуры воздуха, что очень важно, например, при прогнозировании сезонных изменений вертикальной термической структуры воды в губе при различных сценариях изменения регионального климата, или при восстановлении отсутствующих данных. Установлено, что отклонения температуры воздуха в

период максимального прогрева на 3-4°С приводят к близким изменениям температуры воды в верхнем квазиоднородном слое (до 5 м), однако, в слое скачка и гиполимнионе зависимость температуры воды от температуры воздуха существенно уменьшается, но прослеживается связь в сроках и скоростях изменения температур воздуха и воды при весеннем нагревании и осеннем охлаждении.

Основные выводы, полученные в результате выполненных исследований, сводятся к следующему.

1. В результате обобщения многолетних данных инструментальных наблюдений за температурой воды озер Северо-Запада России и Финляндии установлено, что среднемноголетняя температура верхнего пятиметрового слоя определяется преимущественно географической широтой, а в слое скачка и гиполимнионе наибольшее влияние на нее оказывают морфометрические характеристики водных объектов. В годовом цикле температуры поверхности озер влияние географической широты отражается, прежде всего, на процессах их нагревания и на значениях температуры воды в период максимального прогрева, а скорости выхолаживания водоемов связаны преимущественно с их глубиной и площадью.

2. Показано, что для возникновения устойчивой термической стратификации в озерах должно существовать определенное сочетание геометрических размеров водоемов. Для данного региона получены количественные выражения для границ между эпи-, мета- и гипотермическими типами озер в зависимости от их площади и максимальной глубины.

3. Для описания среднемноголетних ежедневных значений температуры воды на различных горизонтах предложена 6-параметрическая аппроксимирующая

Т°С

20 г—

Рис. 7. Совмещенный график модельных кривых для температуры воздуха (жирная линия) и температуры воды на различных горизонтах.

IV V VI упуш к х XI хц Месяцы

функция. Для реконструкции термического состояния малоизученных или вновь создаваемых водных объектов разработана стохастическая модель, которая позволяет оперативно рассчитывать годовой ход температуры поверхности воды и вертикальную термическую структуру в период максимального прогрева для любого водоема в регионе на основе доступной географической информации (широта, высота над уровнем моря, геометрические размеры, коэффициент условного водообмена).

4. Впервые получены среднемноголетние пространственные распределения полей температуры воды в Онежском озере на любую дату для периода открытой воды. Они позволили выявить характерные особенности в формировании температурного режима в различных районах озера на разных горизонтах, количественно оценить интенсивность нагревания и охлаждения отдельных слоев, определить сроки максимального прогревания различных горизонтов, время прохождения термобара, сроки наступления «биологического лета».

5. На примере Петрозаводской губы Онежского озера установлены эмпирические зависимости между годовым ходом температуры воздуха и температуры воды в губе на различных горизонтах, что позволяет прогнозировать сезонные изменения вертикальной термической структуры губы при различных сценариях изменения регионального климата.

Основные публикации по теме диссертации:

1. Физико-географическая характеристика озера. Экосистема Онежского озера. Л.: Наука, 1990. С. 5-11.

2. Температура поверхности воды по данным ИК-съемок. Экосистема Онежского озера. Л.: Наука, 1990. С. 37-46.

3. Моделирование термического режима глубоководной зоны Большой губы Повенецкого залива. // Повенецкий залив Онежского озера. Петрозаводск, 1992. С. 12-18. (соавтор М. П. Петров).

4. Термика и динамика вод озера // Онежское озеро. Экологические проблемы. Петрозаводск: Изд-во КарНЦ РАН, 1999. С. 22-48 (соавторы П. М. Бояринов, Н. И. Пальшин, М. П. Петров, Н. К Филатов).

5. Характеристики термического режима Онежского озера. // Ладожское озеро. Петрозаводск, 2000. С. 330-335 (соавтор М. П. Петров).

6. Географические закономерности термического режима озер. // Экология северных территорий России. Проблемы, прогноз ситуации, пути развития, решения. Архангельск, 2002. С. 366—371. (Материалы международной конференции; Т. 1) (соавтор Н. И. Пальшин).

7. Формирование вертикальной термической структуры озер Северо-Запада России и Финляндии. // Вод. ресурсы.2003, том 30, № 6. С. 696-706 (соавтор Н. И. Пальшин).

8. Стохастическая модель вертикальной термической структуры Петрозаводской губы Онежского озера. Труды VI конференции «Динамика и термика рек, водохранилищ и прибрежной зоны морей». М., 2004. С. 140-143 (соавтор Н. И, Пальшин).

9. Effect ofmorphometry and climate on the thermal regime oflarge lakes. Proceedings ofthe Fourth International Lake Ladoga Symposium 2002. Joensuu. 2003. P. 26-30 (соавтор Н. И. Пальшин).

10. Characteristic of thermal regime in lake Onego. Proceedings of the third international lake Ladoga symposium. 2000. P. 202-205 (соавтор М. П. Петров).

11. Annual cycle of the surface temperature in large lakes. Proceedings of the Fourth International Lake Ladoga Symposium 2002. Joensuu. 2003. P. 96-101 (соавтор Н. И. Пальшин).

Изд. лиц. № 00041 от 30.08.99. Подписано в печать 15.03.05. Формат 60x84 '/16. Бумага офсетная. Гарнитура «Times». Печать офсетная. Уч.-изд. л. 1,3. Усл. печ. л. 1,4. Тираж 100 экз. Изд. № 13. Заказ № 480

Карельский научный центр РАН 185003, Петрозаводск, пр. А. Невского, 50 Редакционно-издательский отдел

Содержание диссертации, кандидата географических наук, Ефремова, Татьяна Владимировна

ВВЕДЕНИЕ.

Глава 1. ВЛИЯНИЕ ГЕОГРАФИЧЕСКИХ ФАКТОРОВ И МОРФОМЕТРИЧЕСКИХ ХАРАКТЕРИСТИК НА ТЕРМИЧЕСКУЮ СТРУКТУРУ ОЗЕР.

1.1. Термический режим водоемов как следствие физико- > географической зональности.

1.2. Значение морфометрических характеристик для исследований термического режима озер.

1.3. Влияние глубины и площади озер на глубину расположения термоклина.

1.4. Физико-географическая характеристика Северо-Запада России.

1.5. Характеристика исследуемых озер.

1.6. Зависимость термической стратификации озер от их геометрических размеров.

Глава 2. СРЕДНЯЯ МНОГОЛЕТНЯЯ ВЕРТИКАЛЬНАЯ ТЕРМИЧЕСКАЯ СТРУКТУРА ОЗЕР В ПЕРИОД МАКСИМАЛЬНОГО ПРОГРЕВА.

2.1. Анализ данных и выбор приоритетных характеристик.

2.2. Зависимость вертикального распределения температуры воды в озерах от географических факторов.

2.3. Анализ полученных результатов.

Глава 3. СРЕДНИЙ МНОГОЛЕТНИЙ ГОДОВОЙ ХОД ТЕМПЕРАТУРЫ

ПОВЕРХНОСТИ ВОДЫ.

3. 1. Аппроксимация годового хода температуры воды функциями от времени.

3. 2. Сезонная изменчивость годовых циклов температуры воды озер Северо-Запада.

3.3. Зависимость температуры поверхности воды в озерах от географических факторов.

Глава 4. ПРОСТРАНСТВЕННО-ВРЕМЕННАЯ ИЗМЕНЧИВОСТЬ СРЕДНЕМНОГОЛЕТНЕЙ ТЕМПЕРАТУРЫ ВОДЫ ГЛУБОКОВОДНЫХ ОЗЕР.

4.1. Сезонная изменчивость вертикальной термической структуры озер Инари и Каллавеси.

4.2. Пространственно-временная изменчивость температуры воды Онежского озера.

4.2.1. Краткая физико-географическая характеристика озера.

4.2.2. Температура поверхности воды по данным ИК-съемок.

4.2.3. Материалы и методика.

4.2.4. Термический режим Онежского озера.

4.3.3ависимость температуры воды на различных горизонтах от температуры воздуха (на примере Петрозаводской губы).

Заключение Диссертация по теме "Гидрология суши, водные ресурсы, гидрохимия", Ефремова, Татьяна Владимировна

ЗАКЛЮЧЕНИЕ

1. В результате обобщения многолетних данных инструментальных наблюдений за температурой воды озер Северо-Запада России и Финляндии установлено, что среднемноголетняя температура верхнего пятиметрового слоя определяется преимущественно географической широтой, а в слое скачка и гиполимнионе наибольшее влияние на нее оказывают морфометрические характеристики водных объектов. В годовом цикле температуры поверхности озер влияние географической широты отражается, прежде всего, на процессах их нагревания и на значениях температуры воды в период максимального прогрева, а скорости выхолаживания водоемов связаны преимущественно с их глубиной и площадью.

2. Показано, что для возникновения устойчивой термической стратификации в озерах должно существовать определенное сочетание геометрических размеров водоемов. Для данного региона получены количественные выражения для границ между эпи-, мета- и гипотермическими типами озер в зависимости от их площади и максимальной глубины.

3. Для описания среднемноголетних ежедневных значений температуры воды на различных горизонтах предложена 6-параметрическая аппроксимирующая функция. Для реконструкции термического состояния малоизученных или вновь создаваемых водных объектов разработана стохастическая модель, которая позволяет оперативно рассчитывать годовой ход температуры поверхности воды и вертикальную термическую структуру в период максимального прогрева для любого водоема в регионе на основе доступной географической информации широта, высота над уровнем моря, геометрические размеры, коэффициент условного водообмена). Модели в силу своей простоты, могут служить как рабочий инструмент для широкого круга специалистов при решении научных и практических задач.

4. Впервые получены среднемноголетние пространственные распределения полей температуры воды в Онежском озере на любую дату для периода открытой воды. Они позволили выявить характерные особенности в формировании температурного режима в различных районах озера на разных горизонтах, количественно оценить интенсивность нагревания и охлаждения отдельных слоев, определить сроки максимального прогревания различных горизонтов, время прохождения термобара, сроки наступления «биологического лета».

5. На примере Петрозаводской губы Онежского озера установлены эмпирические зависимости между годовым ходом температуры воздуха и температуры воды в губе на различных горизонтах, что позволяет прогнозировать сезонные изменения вертикальной термической структуры воды в губе при различных сценариях изменения регионального климата. Установлено, что отклонения температуры воздуха в период максимального прогрева на 3-4°С приводят к близким изменениям температуры воды в верхнем квазиоднородном слое (до 5 м), однако, в слое скачка и гиполимнионе зависимость температуры воды от температуры воздуха существенно уменьшается, но прослеживается связь в сроках и скоростях изменения температур воздуха и воды при весеннем нагревании и осеннем охлаждении.

Библиография Диссертация по наукам о земле, кандидата географических наук, Ефремова, Татьяна Владимировна, Петрозаводск

1. Адаменко В. Н. Климат и озера. Л.: Гидрометеоиздат, 1985. 263 с.

2. Алисов Б. П. Климат СССР. М.: «Высшая школа», 1969. 131 с.

3. Астраханцев Г.П., Мальгин А.Н., Минина Т.Р., Петрова Н.А., Полосков В.Н., Руховец Л.А., Белкина Н.А., Ефремова Т. В., Сабылина А.В., Тержевик А.Ю., Тимакова Т.М., Филатов Н.Н. Моделирование экосистемы Онежского озера. Препринт ЭМИ. 2004. 31 с.

4. Багрова Г. М. Оценка радиационного баланса в зависимости от широты и характеристик водоема.//Водные ресурсы. № 1, 1983. С. 94-101.

5. Бискэ Г. С., Лак Г. Ц., Лукашов А. Ц., Горюнова П. Н., Ильин В. А. Строение и история котловины Онежского озера. Петрозаводск, 1971. 73 с.

6. Богословский Б. Б. Озероведение. М.: Изд-во Московского университета, 1960.335 с.

7. Большие озера Кольского полуострова. Л.: Наука, 1975. 350 с.

8. Боровиков В. П. STATISTICA: искусство анализа данных на компьютере. Для профессионалов. СПб., Питер, 2001. 656 с.

9. Бояринов П. М., Ефремова Т. В., Пальшин Н. И., Петров М. П., Филатов Н.Н. Термика и динамика вод озера // Онежское озеро. Экологические проблемы. Петрозаводск: Изд-во КарНЦ РАН, 1999. С. 22-48.

10. Бояринов П. М., Петров М. П. Процессы формирования термического режима глубоких пресноводных водоемов. Л.: Наука. 1991. 175 с.

11. Бояринов П. М., Руднев С. Ф. Инструментальные исследования течений. Экосистема Онежского озера. Л.: Наука, 1990. С. 53-71.

12. Браславский А. П., Викулина 3. А. Нормы испарения с поверхности водохранилищ. J1., 1954. 212 с.

13. Верещагин Г. Ю. Методы морфологической характеристики озер. // Труды Олонецкой научной экспедиции. J1., 1930. Ч. 2, вып. 1. 106 с.

14. Гидрология озер Воже и Лача. Л.: Наука, 1979. 288 с.

15. Государственный водный кадастр. Многолетние данные о режиме и ресурсах поверхностных вод суши. Л.: Гидрометеоиздат, Т. 1, вып. 5, 1986. 688 с.

16. Государственный водный кадастр. Многолетние данные о режиме и ресурсах поверхностных вод суши. Л.: Гидрометеоиздат, Т. 1, вып. 7, 1987, 220 с.

17. Государственный водный кадастр. Основные гидрологические характеристики. — Л.: Гидрометеоиздат, Т. 2, 1978, 670 с.

18. ДоманицкиЙ А. П., Дубровина Р. Г., Исаева А. И. Реки и озера Советского Союза. Л.: Гидрометеоиздат, 1971. 104 с.

19. Домрачев П. Ф. К вопросу о классификации озер Сев. Зап. Края. // Изв. ГГИ. Вып. 4, №53. 1922. С. 1-43.

20. Драбкова В. Г., Сорокин И. Н. Озеро и его водосбор единая природная система. Л.: Наука, 1979. 195 с.

21. Драйпер Н., Смит Г. Прикладной регрессионный анализ. М., 1986.Т.2. 350 с.

22. Ефремова Т. В. Физико-географическая характеристика озера. Экосистема Онежского озера. Л.: Наука, 1990. С. 5-11.

23. Ефремова Т. В. Температура поверхности воды по данным ИК-съемок. Экосистема Онежского озера. Л.: Наука, 1990. С. 37-46.

24. Ефремова Т. В. Термическая структура вод по данным многолетних наблюдений. // Биологические ресурсы водоемов бассейна Балтийского моря. Тез. докл. XXIII науч. конф. по изуч. Водоемов Прибалтики, Петрозаводск, 1991. С. 4-5.

25. Ефремова Т. В. Термическая структура вод по данным многолетних наблюдений.//Деп. ВИНИТИ 25.11.91. №4400-В91. 1991.23 с.

26. Ефремова Т. В., Пальшин Н. И. Формирование вертикальной термической структуры озер Северо-Запада России и Финляндии. // Вод. ресурсы.2003, том 30, № 6. С. 696-706.

27. Ефремова Т. В., Пальшин Н. И. Стохастическая модель вертикальной термической структуры Петрозаводской губы Онежского озера. Труды VI конференции «Динамика и термика рек, водохранилищ и прибрежной зоны морей». М., 2004. С. 140-143.

28. Ефремова Т. В., Петров М. П. Моделирование термического режима глубоководной зоны Большой губы Повенецкого залива. // Повенецкий залив Онежского озера. Петрозаводск, 1992. С. 12—18.

29. Ефремова Т. В., Петров М. П. Многолетние характеристики термического режима Онежского озера. // Деп. ВИНИТИ 14.04.95. № 1033-В95, 1995.30 с.

30. Изотова А. Ф. Турбулентный тепло- и влагообмен больших озер. Л.: Наука, 1982. 144 с.

31. Кириллова Т. В. Радиационный режим озер и водохранилищ. Л.: Гидрометеоиздат, 1970. 253 с.

32. Китаев С. П. Термические классификации озер мира. // Водные ресурсы. 1978. № 1.С. 97-103.

33. Китаев С. П. Экологические основы биопродуктивности озер разных природных зон. М.: Наука, 1984. 207 с.

34. Лифшиц В. X. Краткая физико-географическая характеристика и некоторые элементы гидрологического режима Петрозаводской губы. Гидробиология Петрозаводской губы Онежского озера. Петрозаводск, 1980. С.5-10.

35. Молчанов И. В., Былинкина В. Н., Викулина 3. А., Горшунова Т. А. Онежское озеро. Л., 1946. 207 с.

36. Муравейский С. Д. Очерки по теории и методам морфометрии озер. // Реки и озера. М., 1960. С. 91-125.

37. Науменко М. А. О термических различиях между центральной и прибрежной частями крупного озера в период развитой стратификации. Водные ресурсы, 1994, т. 21, №6. С. 590-594.

38. Науменко М. А. Закономерности пространственно-временной изменчивости термических процессов в крупных димиктических озерах. Автореф. дис. док. геог. наук. СПб, 1998.38 с.

39. Науменко М. А., Гузиватый В. В., Каретников С. Г. Ежедневные средние пространственные распределения температуры поверхности воды Ладожского озера с мая по ноябрь // Ладожское озеро. Петрозаводск, 2000. С. 335 -345.

40. Науменко М. А., Каретников С. Г. Морфометрия и особенности гидрологического режима Ладожского озера. // Ладожское озеро: прошлое, настоящее, будущее. СПб., Наука, 2002. С. 16-49.

41. Науменко М. А., Каретников С. Г. Сезонная эволюция пространственного распределения температуры поверхности воды Ладожского озера в связи с его морфометрией. // Доклады Академии Наук, 2002. Т.386, N 2. С. 250 -253.

42. Несина Л. В., Огнева Т. А. Солнечная радиация и тепловой баланс водоемов. // Круговорот вещества и энергии в озерных водоемах. Новосибирск, 1975. С. 308-313.

43. Николаев И. И. Современное состояние Петрозаводской губы Онежского озера как экологической системы. Гидробиология Петрозаводской губы Онежского озера. Петрозаводск, 1980. С. 169-174.

44. Пальшин Н. И., Ефремова Т. В. Пространственная и временная изменчивость термохалинной структуры и динамики вод Выгозерского водохранилища. Деп. ВИНИТИ 13.03.92. № 869-В92, 1992.45 с.

45. Пальшин Н. И., Ефремова Т. В. Гидрологическая структура вод Петрозаводской губы. Гидрологические исследования на водоемах Карелии. Петрозаводск, 1994. С. 3-9.

46. Пальшин Н. И., Ефремова Т. В Стохастическая модель годового хода температуры поверхности воды в озерах. // Метеорология и гидрология. 2005, № 3, (в печати).

47. Пановский Г. А., Брайер Г. В. Статистические методы в метеорологии. Л.: Гидрометеоиздат, 1972. 209 с.

48. Пехович А. И., Жидких В. М. Расчет термического режима водохранилищ и озер. //Тр. IV Всесоюз. гидрол. съезда. Т. 6. Л., 1976. С. 172-185.

49. Пивоварова 3. И. Радиационные характеристики климата СССР. Л.: Гидрометеоиздат, 1977. 336 с.

50. Ресурсы поверхностных вод СССР. Том 1. Кольский полуостров. Л., 1970. 316с.

51. Ресурсы поверхностных вод СССР. Том 2. Карелия и Северо-Запад. Часть 1. Л., 1972. 527 с.

52. Ресурсы поверхностных вод СССР. Том 2. Карелия и Северо-Запад. Часть 2. Л., 1972. 278 с.

53. Румянцев В. А., Тройская Т. П., Догановский А. М. Озера и водохранилища: проблемы изучения, охраны и рационального использования. VI Всероссийский гидрологический съезд. Тезисы докладов, пленарное заседание. 2004. С. 39-42.

54. Румянцев В. А., Трапезников Ю. А., Григорьев А. С. Стохастические модели влияния климата на гидрологию озер. СПб.: НИИ Химии СПбГУ, 2001. 156 с.

55. Рянжин С. В. Закономерности температурного режима пресноводных озер мира. Л., 1989. 70 с.

56. Рянжин С. В. Температура поверхности пресноводных озер северного полушария в зависимости от географической широты и высоты озера над уровнем моря. //Докл. АН СССР, 1990. Т. 312, № 1. С. 209-214.

57. Рянжин С. В. Зональные закономерности элементов термического режима пресноводных озер северного полушария // Вод. ресурсы. 1991. № 4. С. 1529.

58. Смирнова Н. П. Классификация озер по температуре воды и элементам теплового баланса. // Теоретические вопросы классификации озер. СПб.: Наука, 1993. С. 35-51.

59. Советов С. Л. Онежское озеро. Пгр., 1917. 164 с.

60. Сорокин И. Н. Морфологические, морфометрические и гидрологические показатели и их роль в комплексной классификации озер и районировании. // Теоретические вопросы классификации озер. СПб.: Наука, 1993. С. 24-35.

61. Состояние окружающей среды Северо-западного и Северного регионов России. СПб.: Наука. 1995. 370 с.

62. Страшкраба М., Гнаук А. Пресноводные экосистемы. Математическое моделирование. М.: Мир. 1989. 373 с.

63. Тихомиров А. И. Некоторые результаты исследований термического режима Онежского озера// Изв. ВГО. Т.98, вып. 6, 1966. С. 513-520.

64. Тихомиров А. И. Температура воды, теплозапасы, тепловой баланс и термический режим Онежского озера // Тепловой режим Онежского озера. Л.: Наука, 1973. С. 202 323.

65. Тихомиров А. И. Термика крупных озер. Л.: Наука, 1982. 232 с.

66. Филатов Н. Н. Зона проявления, повторяемость прибрежного апвеллинга и циркуляция вод озера // Экосистема Онежского озера. Петрозаводск, 1990. С. 47-53.

67. Форш Л. Ф. Термический и тепловой режим озер Кольского полуострова. // Озера различных ландшафтов Кольского полуострова. 4.1, Л.: Наука,, 1974. С.156-194.

68. Хатчинсон Дж. Лимнология. М.: Прогресс, 1969. 592 с.

69. Хендерсон-Селлерс Б. Инженерная лимнология. Л.: Гидрометеоиздат, 1987. 335 с.

70. Химмельблау Д. Анализ процессов статистическими методами. М.: Мир, 1973.957 с.

71. Хомскис В. Динамика и термика малых озер. Вильнюс: Минтис, 1969. 220 с.

72. Хуббатуллин В. J1. Особенности температурного режима малых озер (на примере Запада и Северо-Запада ETC). Автореф. дисс. канд. геогр. наук, Л.,1989. 20 с.

73. Чудско-Псковское озеро. Под ред. А. А. Соколова. Л.: Гидрометеоиздат, 1983. 162 с.

74. Экосистема озера Плещеево. Л.: Наука, 1989. 263 с

75. Якушко О. Ф. Белорусское Поозерье. Минск, 1971. 336 с.

76. Arai Т. Climatic and geomorphological influences on lake temperature // Int. Ver. Theor. Angew. Limnol. Verh. 1981. V. 21. P. 130-134.

77. Bennett E. B. Characteristics of the thermal regime of Lake Superior. J. Gr. Lakes Res., 1978. V. 4, №3-4. P. 310-319.

78. Birge E. A. Plankton studies on Lake Mendota. Ttans. Wis. Acad. Sci. Arts. Lett., II., 1897. P. 274-448.

79. Birge E. A. On the evidence for temperature seiches. Ttans. Wis. Acad. Sci. Arts. Lett., 16, 1910. P. 1005-1016.

80. Birge E. A. The heat budget of American and Europian lakes. Trans. Wis. Acad. Sci. Arts. Lett. 18, 1915. P. 166-391.

81. Birge E. A., Juday C., March M.W. The temperature of the bottom deposits of Lake Mendota. — A chapter in the heat exchange of the lake. Trans. Wise. Acad. 1927, vol. 23. P. 187—232.

82. Blumberg Л. F., Di Того D.M. Effect of Climate Wanning on Dissolved Oxygen Concentrations in Lake Erie. 210-223.

83. Carmack E. C., Farmer D. M. Cooling processes in deep, temperate lakes. A review with examples from two lakes in British Columbia // J. Mar. Res. 1982. 40. Suppl. P. 85-111.

84. Efremova Т. V., Palshin N. I. Effect of morphometry and climate on the thermal regime of large lakes. Proceedings of the Fourth International Lake Ladoga Symposium 2002. Joensuu. 2003. P. 26-30.

85. Efremova Т. V., Petrov M. P. Characteristic of thermal regime in lake Onego. Proceedings of the third international lake Ladoga symposium. 2000. P. 202-205.

86. Forel F. A. La thermique des lacs Lean douce. Verh. Schweiz. Naturforsch. Ges. 1882. 75. P. 5-8.

87. Geistbeck A. Die Seen der Deutschen Alpen mit 128 Figuren geologischen and geographischen profilen, Tiffenschichtenkarten und Diagrammen. Leipzig, 1885. 47 p.

88. Gorham E. Morphometric control of annual heat budget in temperate lakes. Limnol. Oceanogr., 1964, № 9, P. 525-529.

89. Gorham E., Boyce F. M. Influence of lake surface area and depth upon thermal stratification and the depth of the summer thermocline // J. Gr. Lakes Res. 1989. V. 15(2). P. 233-245.

90. Hondzo M., Stefan H.G. Regional water temperature characteristics of lakes subjected to climate change. Climate change, 1993. 24. P 187-211.

91. Korhonen J. Suomen vesistojen lampotilaolot 1900-luvulla. Suomen ymparisto 566, Helsinki, 2002. P. 116. In Finish..

92. Kuusisto E. Suomenvesistojen lampotilat kaudella 1961-1975. Vesihallitus -National Board of Waters, Finland. Helsinki, 1981. 40 p. In Finish..

93. Lathrop R. C., Lillie R. A. Thermal stratification of Wisconsin lakes. J. Wise. Acad. Sci. 68, 1980. P. 90-96.

94. Lesht В. M., Brandner D. J. Functional representation of Great Lakes Surface temperatures —J. Great Lakes Res. 1992. Vol. 18, N 1. P. 98—107.

95. Lewis W. M. A revised classification of lakes based on mixing. Can. J. Fish. Aquat. Sci. 1983. 40, N 10. P 1779-1787.

96. McCormic M. J. Potential Changes in Thermal Structure and Cycle of Lake Michigan Due to Global Wanning. Transact. Amer. Fisher. Soc. 1990, 119(2). P. 183-194.

97. Mortimer С. H. Lake hydrodynamics. Int. Ver. Theor. Angew. Limnol. Verh., 20. 1974. P. 124-197.

98. Rogers G.K. Time and onset of full thermal stratification in Lake Ontario in relation to lake temperatures in winter. Canadian Journal of Fisheries and Aquatic Sciences, 1987. 44(12). P. 2025-2229.

99. Ruttner F. Hydrographische und hydrochemische Beobachtungen auf Java, Sumatra und Bali. Arch. Hydrobiol., Suppl., 1931, 8. P. 197-454.

100. Yoshimura S. A. A contribution to the knowledge of deep water temperatures of Japanese lakes. Part I, Summer temperatures. Japan. J. Astr. Geophs., V. 13, 1936. P. 20-61.

101. Zalar A. R. Taxonomy oflakes. Hydrobiologia. 1959. Vol. 13, N 3. P. 287-299.