Формирование ледяного покрова озер зависит от большого количества факторов, как за краткосрочный период (температура воды, температура воздуха, ветер и другие), так и задолго до наступления зимнего периода (тепловой запас озера). Поэтому лед является чувствительным индикатором изменения этих факторов [1, 2].
В связи с этим исследования ежегодно повторяющегося ледового режима озер необходимы для решения важной экологической задачи – установление тенденций глобальных и региональных изменений климата [3].
Помимо этого, формирование и разрушение льда на озерах напрямую оказывает влияние на водную экосистему озер. В частности, ледяной покров на озере уменьшает проникание солнечного света, снижает передачу тепла воды атмосфере, прекращает передачу кинетической энергии ветра воде. Эти факторы формируют экологическую систему озер в зимнее время [2].
Исследования характера протекания ледового режима озер сводится к оценке зависимости ледовитости озер (части акватории занятой льдом) от времени.
В настоящее время для получения сведений о состоянии земной поверхности, в том числе о ледяном покрове озер, широко используются спутниковые данные. Датчики спутников (MODIS, VIIRS, AIRS, MISR и многие другие) ежедневно выполняют многозональную съемку местности в различных диапазонах электромагнитной волны (видимом, инфракрасном, микроволновом) в течение нескольких последних лет (10–25 лет) [4–6].
С высокой точностью можно рассчитать ледовитость озер визуально-экспертной оценкой, исходя из спутниковых снимков датчика MODIS, выполненных в видимом диапазоне (синтезированные RGB-изображения), обладающих высоким пространственным разрешением (до 250 м в открытом доступе https://earthdata.nasa.gov/earth-observation-data).
Однако в большинстве случаев выполнить расчет ледовитости озер по спутниковым снимкам видимого диапазона невозможно из-за наличия облачности. Например, за период наблюдений 2006–2018 гг. для оценки состояния ледяного покрова Онежского озера доступно 636 снимков, Ладожского озера – 411 снимков из-за большей площади акватории (чем больше площадь акватории озера, тем меньше вероятность того, что акватория озера полностью свободна от облаков). Поэтому формировать непрерывный ряд значений ледовитости озер, используя только эти снимки, невозможно.
Однако существуют многосенсорные СВЧ-радиометры, которые могут выполнять многозональную съемку в микроволновом диапазоне (ASMU-A, ATMS, AVHHR, MODIS, VIIRS и другие), что позволяет ежедневно фиксировать сцену вне зависимости от облачности. Результатом автоматического картирования многозональных снимков этих систем являются данные о снежном и ледяном покрове планеты. Эти данные представлены Национальным центром США по снегу и льду NSDIC, а также Центром спутниковых приложений и исследований NOAA NESDIS на ftp-серверах в открытом доступе (ftp://sidads.colorado.edu/DATASETS/NOAA/G02156/, ftp://ftp.star.nesdis.noaa.gov/pub/smcd/emb/snow/binary/multisensor/global/). Несмотря на незначительное пространственное разрешение 4–6 км, эти данные пригодны для расчета ледовитости многих озер.
Однако установлено, что в некоторых случаях при автоматической интерпретации спутниковых данных возникают существенные различия при определении пространственного распределения льда на озерах. Все это приводит к ошибкам при расчете ледовитости озер. Например, по данным NSDIC и NOAA NESDIS некоторые участки акватории озер могут быть покрыты льдом (рис. 1, а, б), в то время как в действительности льда на этих участках нет (рис. 1, в).
Рис. 1. Визуализация спутниковых данных о состоянии ледяного покрова Онежского озера 25.01.2012: a) по данным NSDIC; б) по данным NOAA NESDIS; в) по данным датчика MODIS. Красным цветом показан лед, синим – вода (a, б)
В связи с этим целью данной работы является проверка адекватности использования данных NSDIC и NOAA NESDIS для определения пространственного распределения льда на озерах, а также для расчета ледовитости озер. Исследование включает следующие этапы:
1. Сбор и анализ спутниковых данных NSDIC, NOAA NESDIS о состоянии ледяного покрова озер за период 2006–2018 гг.
2. Расчет и анализ ошибок спутниковых данных NSDIC, NOAA NESDIS, возникающих при определении пространственного распределения льда на озерах, а также при расчете значений ледовитости озер в моменты времени, на которые доступны достоверные данные.
Материалы и методы исследования
Для критического анализа спутниковых данных NSDIC, NOAA NESDIS при определении пространственного распределения льда выбраны озера, расположенные в республике Карелия РФ. Поскольку в этом географическом районе одновременно расположены крупнейшие европейские озера – Онежское и Ладожское, а также небольшие (но пригодные для анализа ледяного покрова при пространственном разрешении 4–6 км) – Сегозеро и Выгозеро (рис. 2). Выбранные озера ежегодно полностью или частично покрываются льдом, что является необходимым условием для проведения данного исследования.
Рис. 2. Рассматриваемые озера в проекции Меркатора (площади акваторий): a) Онежское озеро (9720 км2), б) Ладожское озеро (17700 км2), в) Сегозеро (815 км2), г) Выгозеро (1250 км2)
Это позволило определить несоответствие спутниковых данных NSDIC, NOAA NESDIS при определении пространственного распределения льда и расчете ледовитости приблизительно равных и разных по площади акваторий озер в одинаковых климатических условиях.
Временной интервал применяемых в исследовании спутниковых данных соответствует 2006–2018 гг. В этот интервал времени доступны данные всех источников.
В данном исследовании в качестве достоверных данных использовались данные MODIS, которые имеют пространственное разрешение 500 м, в то время как данные NSDIC, NOAA NESDIS имеют пространственное разрешение 4–6 км. Поэтому сравнительный анализ пространственного распределения льда на акватории озер имеет особенности, обусловленные сравнением данных разного пространственного разрешения.
Однородные участки акватории озера 
по данным MODIS имеют соответственно площади 
. Однородные участки акватории озера 
по данным NSDIC, NOAA NESDIS имеют соответственно площади 
. При этом 
.
– множество однородных участков по данным MODIS, которые соответствуют участку 
по географическим координатам, при этом являются акваторией озера, но не являются льдом.
– множество однородных участков по данным MODIS, которые соответствуют участку 
по географическим координатам, при этом являются акваторией озера, а также являются льдом.
– множество однородных участков по данным MODIS, которые соответствуют участку 
по географическим координатам, при этом не являются акваторией озера.
Расчет ошибки спутниковых данных (NSDIC, NOAA NESDIS) при определении пространственного распределения льда на озерах для каждого момента времени осуществлялся по формуле
где 
– площадь однородного участка 
из множества 
, которые имеют несоответствие с данными MODIS; r – количество однородных участков, которые имеют несоответствие с данными MODIS; Stotal – общая площадь акватории озера, обусловленная пространственным разрешением проверяемых спутниковых данных 
.
Несоответствие каждого участка акватории озера проверяемых данных 
с данными MODIS считалось, если выполняется условие
В случае удовлетворения условия участок акватории озера 
включался во множество 
.
Расчеты δ для озер выполнялись в моменты времени, когда хотя бы по одному из источников данных (NSDIC, NOAA NESDIS, MODIS) на озере имелись ледовые образования. Затем для каждого озера вычислялось среднее значение ошибок 
за все моменты времени.
Результаты исследования и их обсуждение
За период 2006–2018 гг. количество снимков датчика MODIS без облачности в периоды ледовых явлений на озерах составило: 636 снимков Онежского озера, 411 снимков Ладожского озера, 629 снимков Сегозера, 571 снимок Выгозера. Эти снимки использовались при расчете ошибок данных NSDIC, NOAA NESDIS (для каждого момента времени, на который имеется снимок).
Средние значения ошибок при определении пространственного распределения льда на озерах 
и средние значения абсолютных отклонений ледовитости озер спутниковых данных NSDIC, NOAA NESDIS относительно фактических значений ледовитости (по данным MODIS) 
представлены в таблице.
Средние значения ошибок при определении пространственного распределения льда на озерах 
и средние абсолютные отклонения ледовитости озер 
спутниковых данных NSDIC, NOAA NESDIS
|
Озера |
|
|
||
|
NSDIC |
NOAA NESDIS |
NSDIC |
NOAA NESDIS |
|
|
Онежское озеро |
5,9 |
18,5 |
4,7 |
16,3 |
|
Ладожское озеро |
9,8 |
11,1 |
8,2 |
7,5 |
|
Сегозеро |
3,3 |
12,9 |
3,2 |
12,6 |
|
Выгозеро |
5,9 |
11,4 |
5,7 |
11,1 |
, %
, %Результаты показывают, что данные NSDIC имеют меньшие значения 
, чем данные NOAA NESDIS (
). При этом для всех озер, кроме Ладожского озера, 
.
Ледовитость характеризует только площади ледяных образований на озерах, но не характеризует их расположение, поэтому справедливо неравенство: 
. Ледовитость на озерах имеет большое практическое значение при определении сроков и длительностей фаз ледового режима. Расчет ледовитости озер с меньшим значением 
позволяют выполнить лучше данные NSDIC, чем данные NOAA NESDIS, для всех озер, кроме Ладожского озера (
).
По спутниковым данным NSDIC, NOAA NESDIS сформированы суточные ряды значений ледовитости озер: Онежское озеро, Ладожское озеро, Сегозеро и Выгозеро за период 2006–2018 гг. Графики зависимостей значений ледовитости озер от времени представлены на рис. 3.
Рис. 3. Графики зависимости ледовитости от времени: a) Онежское озеро; б) Ладожское озеро; в) Сегозеро; г) Выгозеро. Обозначения: 1 – по данным NSDIC, 2 – по данным NOAA NESDIS, 3 – по данным MODIS
Средние абсолютные отклонения значений ледовитости, полученных по данным NSDIC, от значений ледовитости, полученных по данным NOAA NESDIS, составляют: 18,9 % – для Онежского озера, 11,7 % – для Ладожского озера, 19,3 % – для Сегозера, 12,9 % – для Выгозера.
Даты начала и конца, а также длительности ледовых явлений на озерах согласно данным NSDIC иллюстрированы на рис. 4 в виде диаграмм.
Рис. 4. Даты начала и конца ледовых явлений (с указанием длительности) по данным NSDIC для озер: a) Онежское озеро; б) Ладожское озеро; в) Сегозеро; г) Выгозеро
Определение дат начала и конца ледовых явлений на рассматриваемых озерах по данным NOAA NESDIS затруднены тем, что в некоторые года ледовые явления наблюдаются и в летнее время. Например, по данным NOAA NESDIS в 2013 г. в августе и сентябре зафиксировано наличие ледовых образований на Ладожском озере, в то время как в июле и октябре их не было. Это не соответствует действительности для данных озер, которые в летнее время всегда полностью свободны ото льда. При этом невозможно определить границы начала и конца двух смежных периодов ледовых явлений. Поэтому сравнение дат начал и конца ледовых явлений по данным NSDIC и NOAA NESDIS выполнить не удалось.
Количество дней, когда на озерах наблюдались ледовые явления, за период 2006–2018 гг. по данным NSDIC и NOAA NESDIS соответственно составляет: 1846 дней и 2263 дней для Онежского озера (относительное отклонение ε = 20,3 %); 1540 дней и 2159 дней для Ладожского озера (ε = 33,5 %); 1822 дня и 2140 дней для Сегозера (ε = 16,1 %); 1941 день и 2378 дней для Выгозера (ε = 20,2 %).
Значения 
являются удовлетворительными (менее 6 %) при определении пространственного распределения льда всех озер, кроме Ладожского озера (9,8 %). Поэтому данные NSDIC представляют больший интерес при определении пространственного распределения льда, а также при расчетах ледовитости, чем данные NOAA NESDIS. Однако значение 
при определении пространственного распределения льда Ладожского озера почти в 2 раза больше, чем других озер. В связи с этим проведен детальный анализ ошибок, которые возникают при определении пространственного распределения льда на озерах по данным NSDIC, NOAA NESDIS и получена зависимость среднего значения ошибки при определении пространственного распределения льда на озерах от фактического значения ледовитости озер 
.
Также установлена зависимость 
, характеризующая распределение 
от фактических значений ледовитости (по данным MODIS).
Установлено, что зависимость 
монотонно возрастает при увеличении фактических значений ледовитости (по данным MODIS), достигая пика 
) при значении ледовитости 0,3–0,4, в дальнейшем происходит монотонное снижение ошибки до минимального значения (
при полном ледоставе), график зависимости 
изображен на рис. 5.
Рис. 5. Графики зависимости 
: a) данных NSDIC; б) данных NOAA NESDIS
При этом необходимо отметить, что на 70–80 % отобранных для анализа снимков MODIS озера на них полностью покрыты льдом, исключение составляет Ладожское озеро (лишь 30 % всех снимков). Этот факт объясняется тем, что небольшие озера: Сегозеро и Выгозеро, – а также Онежское озеро ежегодно полностью покрывается льдом на значительный интервал времени 80–120 дней. Ладожское озеро не всегда полностью покрывается льдом (рис. 3, б). В 2008–2009 гг., 2014–2017 гг. не имеется снимков MODIS, свободных от облачности, на которых Ладожское озеро полностью покрыто льдом. Все это объясняет относительно большое значение 
(9,8 %) при исследовании Ладожского озера в сравнении с другими озерами.
Также следует отметить, что данным NSDIC с небольшим средним значением 
(5,0 %) позволяют зафиксировать полное отсутствие льда на озере (рис. 5). Таким образом, с помощью данных NSDIC можно адекватно в целом определить даты начала и конца периодов, когда озеро покрыто полностью, частично и полностью не покрыто льдом, при этом большое значение 
(до 35 %) в фазы замерзания и вскрытия не влияет на расчет этих дат. Это позволяет решить важную гидрологическую задачу – определение сроков и длительностей фаз ледового режима озер. Однако формирование суточного ряда значений ледовитости озер с помощью данных NSDIC во время замерзания и вскрытия является недопустимым. Для этого необходимо использовать дополнительные более точные данные.
Зависимость 
имеет схожий характер с зависимостью 
, максимальное значение 
достигается при значении ледовитости 0,6–0,7, а минимальное значение 
при полном ледоставе (рис. 5).
Отличием зависимости 
от зависимости 
является наличие большого значения 
в моменты, когда озеро полностью свободно ото льда. Такое значение 
учитывает летние периоды, когда по данным NOAA NESDIS имеют место ледовые явления на озерах, в то время как в действительности рассматриваемые озера в летнее время всегда полностью свободны ото льда. По этой же причине за период 2006–2018 гг. количество дней, когда на озерах наблюдались ледовые явления по данным NOAA NESDIS для всех озер больше на 15–30 %, чем по данным NSDIC.
Выводы
Установлено, что данные NSDIC имеют меньшее значение 
(3–10 %), чем данные NOAA NESDIS (11–19 %) для озер (Онежского, Ладожского озер, Сегозера, Выгозера) в сравнении с данными MODIS. Поэтому при определении пространственного распределения льда на озерах целесообразно использовать данные NSDIC.
Во время замерзания и вскрытия по данным NSDIC (а также NOAA NESDIS) 
могут достигать неприемлемого значения 40 % (значения 
до 35 %). В связи с этим формирование подробного временного ряда значений ледовитости во время замерзания и вскрытия с использованием спутниковых данных NSDIC, NOAA NESDIS возможно только при корректировке их другими более точными данными, например, данными MODIS.
Однако данные NSDIC позволяют рассчитать значение ледовитости озер во время полного ледостава и отсутствия ледовых явлений на озере с достаточной точностью (
и 
). Это позволяет сделать вывод о том, что применение данных NSDIC целесообразно для определения сроков и длительностей фаз ледового режима.
Работа выполнена при финансовой поддержке гранта Президента РФ для поддержки молодых российских ученых МК-3379.2018.5.
Библиографическая ссылка
Баклагин В.Н. КРИТИЧЕСКИЙ АНАЛИЗ СПУТНИКОВЫХ ДАННЫХ О ПРОСТРАНСТВЕННОМ РАСПРЕДЕЛЕНИИ ЛЬДА НА ОЗЕРАХ // Успехи современного естествознания. – 2019. – № 2. – С. 41-48;URL: https://natural-sciences.ru/ru/article/view?id=37055 (дата обращения: 25.04.2024).