Scientific journal
Advances in current natural sciences
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,775

CRITICAL ANALYSIS OF SATELLITE DATA ON THE SPATIAL DISTRIBUTION OF ICE ON LAKES

Baklagin V.N. 1
1 Northern Water Problems Institute of the Karelian Research Centre of the Russian Academy of Sciences
The process of formation and rotting of ice on lakes is an integral part of the hydrological cycle of many lakes. The conditions of the ice regime significantly influence the ecological system of lakes. The paper includes calculation and analysis of errors in the determination of the spatial ice distribution (spatial resolution of 4-6 km) on Lake Onego, Lake Ladoga, Lake Segozero and Lake Vigozero within the period of 2006-2017 according to NSIDC, NOAA NESDIS data with regard to reliable MODIS data (spatial resolution of 500 m). It was established that within the monitoring period NSIDC data has the minimum mean values of errors in determining the spatial distribution of ice on lakes (3-10 %) comparing to NOAA NESDIS data (11-19 %) and it is also of more practical interest in estimating the ice coverage of lakes. The dependence of the mean value of errors that occur in determination of the spatial distribution of ice (according to NSIDC, NOAA and NESDIS data) on the actual value of ice coverage (according to MODIS) was revealed. The results shows that the NSIDC data allows estimating adequately the phases of the ice regime; however formation of a daily time series of ice coverage during Freeze-up and Break-up phases is possible only with a significant error (mean value of absolute deviations according to MODIS data is up to 35 %).
NSIDC
NOAA NESDIS
MODIS
spatial distribution of ice
ice coverage
lakes

Формирование ледяного покрова озер зависит от большого количества факторов, как за краткосрочный период (температура воды, температура воздуха, ветер и другие), так и задолго до наступления зимнего периода (тепловой запас озера). Поэтому лед является чувствительным индикатором изменения этих факторов [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, в).

bakl1.tif

Рис. 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). Выбранные озера ежегодно полностью или частично покрываются льдом, что является необходимым условием для проведения данного исследования.

bakl2.tif

Рис. 2. Рассматриваемые озера в проекции Меркатора (площади акваторий): a) Онежское озеро (9720 км2), б) Ладожское озеро (17700 км2), в) Сегозеро (815 км2), г) Выгозеро (1250 км2)

Это позволило определить несоответствие спутниковых данных NSDIC, NOAA NESDIS при определении пространственного распределения льда и расчете ледовитости приблизительно равных и разных по площади акваторий озер в одинаковых климатических условиях.

Временной интервал применяемых в исследовании спутниковых данных соответствует 2006–2018 гг. В этот интервал времени доступны данные всех источников.

В данном исследовании в качестве достоверных данных использовались данные MODIS, которые имеют пространственное разрешение 500 м, в то время как данные NSDIC, NOAA NESDIS имеют пространственное разрешение 4–6 км. Поэтому сравнительный анализ пространственного распределения льда на акватории озер имеет особенности, обусловленные сравнением данных разного пространственного разрешения.

Однородные участки акватории озера baklag01.wmf по данным MODIS имеют соответственно площади baklag02.wmf. Однородные участки акватории озера baklag03.wmf по данным NSDIC, NOAA NESDIS имеют соответственно площади baklag04.wmf. При этом baklag05.wmf.

baklag06.wmf – множество однородных участков по данным MODIS, которые соответствуют участку baklag07.wmf по географическим координатам, при этом являются акваторией озера, но не являются льдом.

baklag08.wmf – множество однородных участков по данным MODIS, которые соответствуют участку baklag09.wmf по географическим координатам, при этом являются акваторией озера, а также являются льдом.

baklag10.wmf – множество однородных участков по данным MODIS, которые соответствуют участку baklag11.wmf по географическим координатам, при этом не являются акваторией озера.

Расчет ошибки спутниковых данных (NSDIC, NOAA NESDIS) при определении пространственного распределения льда на озерах для каждого момента времени осуществлялся по формуле

baklag12.wmf

где baklag13.wmf – площадь однородного участка baklag14.wmf из множества baklag62a.wmfbaklag62b.wmf, которые имеют несоответствие с данными MODIS; r – количество однородных участков, которые имеют несоответствие с данными MODIS; Stotal – общая площадь акватории озера, обусловленная пространственным разрешением проверяемых спутниковых данных baklag17.wmf.

Несоответствие каждого участка акватории озера проверяемых данных baklag18.wmf с данными MODIS считалось, если выполняется условие

baklag19.wmf

В случае удовлетворения условия участок акватории озера baklag20.wmf включался во множество baklag21.wmf.

Расчеты δ для озер выполнялись в моменты времени, когда хотя бы по одному из источников данных (NSDIC, NOAA NESDIS, MODIS) на озере имелись ледовые образования. Затем для каждого озера вычислялось среднее значение ошибок baklag23.wmf за все моменты времени.

Результаты исследования и их обсуждение

За период 2006–2018 гг. количество снимков датчика MODIS без облачности в периоды ледовых явлений на озерах составило: 636 снимков Онежского озера, 411 снимков Ладожского озера, 629 снимков Сегозера, 571 снимок Выгозера. Эти снимки использовались при расчете ошибок данных NSDIC, NOAA NESDIS (для каждого момента времени, на который имеется снимок).

Средние значения ошибок при определении пространственного распределения льда на озерах baklag24.wmf и средние значения абсолютных отклонений ледовитости озер спутниковых данных NSDIC, NOAA NESDIS относительно фактических значений ледовитости (по данным MODIS) baklag25.wmf представлены в таблице.

Средние значения ошибок при определении пространственного распределения льда на озерах baklag26.wmf и средние абсолютные отклонения ледовитости озер baklag27.wmf спутниковых данных NSDIC, NOAA NESDIS

Озера

baklag28.wmf, %

baklag29.wmf, %

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 имеют меньшие значения baklag30.wmf, чем данные NOAA NESDIS (baklag31.wmf). При этом для всех озер, кроме Ладожского озера, baklag32.wmf.

Ледовитость характеризует только площади ледяных образований на озерах, но не характеризует их расположение, поэтому справедливо неравенство: baklag33.wmf. Ледовитость на озерах имеет большое практическое значение при определении сроков и длительностей фаз ледового режима. Расчет ледовитости озер с меньшим значением baklag34.wmf позволяют выполнить лучше данные NSDIC, чем данные NOAA NESDIS, для всех озер, кроме Ладожского озера (baklag35.wmf).

По спутниковым данным NSDIC, NOAA NESDIS сформированы суточные ряды значений ледовитости озер: Онежское озеро, Ладожское озеро, Сегозеро и Выгозеро за период 2006–2018 гг. Графики зависимостей значений ледовитости озер от времени представлены на рис. 3.

bakl3.tif

Рис. 3. Графики зависимости ледовитости от времени: a) Онежское озеро; б) Ладожское озеро; в) Сегозеро; г) Выгозеро. Обозначения: 1 – по данным NSDIC, 2 – по данным NOAA NESDIS, 3 – по данным MODIS

Средние абсолютные отклонения значений ледовитости, полученных по данным NSDIC, от значений ледовитости, полученных по данным NOAA NESDIS, составляют: 18,9 % – для Онежского озера, 11,7 % – для Ладожского озера, 19,3 % – для Сегозера, 12,9 % – для Выгозера.

Даты начала и конца, а также длительности ледовых явлений на озерах согласно данным NSDIC иллюстрированы на рис. 4 в виде диаграмм.

bakl4.tif

Рис. 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 %).

Значения baklag36.wmf являются удовлетворительными (менее 6 %) при определении пространственного распределения льда всех озер, кроме Ладожского озера (9,8 %). Поэтому данные NSDIC представляют больший интерес при определении пространственного распределения льда, а также при расчетах ледовитости, чем данные NOAA NESDIS. Однако значение baklag37.wmf при определении пространственного распределения льда Ладожского озера почти в 2 раза больше, чем других озер. В связи с этим проведен детальный анализ ошибок, которые возникают при определении пространственного распределения льда на озерах по данным NSDIC, NOAA NESDIS и получена зависимость среднего значения ошибки при определении пространственного распределения льда на озерах от фактического значения ледовитости озер baklag38.wmf.

Также установлена зависимость baklag39.wmf, характеризующая распределение baklag40.wmf от фактических значений ледовитости (по данным MODIS).

Установлено, что зависимость baklag41.wmf монотонно возрастает при увеличении фактических значений ледовитости (по данным MODIS), достигая пика baklag42.wmf) при значении ледовитости 0,3–0,4, в дальнейшем происходит монотонное снижение ошибки до минимального значения (baklag43.wmf при полном ледоставе), график зависимости baklag44.wmf изображен на рис. 5.

bakl5.tif

Рис. 5. Графики зависимости baklag45.wmf: a) данных NSDIC; б) данных NOAA NESDIS

При этом необходимо отметить, что на 70–80 % отобранных для анализа снимков MODIS озера на них полностью покрыты льдом, исключение составляет Ладожское озеро (лишь 30 % всех снимков). Этот факт объясняется тем, что небольшие озера: Сегозеро и Выгозеро, – а также Онежское озеро ежегодно полностью покрывается льдом на значительный интервал времени 80–120 дней. Ладожское озеро не всегда полностью покрывается льдом (рис. 3, б). В 2008–2009 гг., 2014–2017 гг. не имеется снимков MODIS, свободных от облачности, на которых Ладожское озеро полностью покрыто льдом. Все это объясняет относительно большое значение baklag46.wmf (9,8 %) при исследовании Ладожского озера в сравнении с другими озерами.

Также следует отметить, что данным NSDIC с небольшим средним значением baklag47.wmf (5,0 %) позволяют зафиксировать полное отсутствие льда на озере (рис. 5). Таким образом, с помощью данных NSDIC можно адекватно в целом определить даты начала и конца периодов, когда озеро покрыто полностью, частично и полностью не покрыто льдом, при этом большое значение baklag48.wmf (до 35 %) в фазы замерзания и вскрытия не влияет на расчет этих дат. Это позволяет решить важную гидрологическую задачу – определение сроков и длительностей фаз ледового режима озер. Однако формирование суточного ряда значений ледовитости озер с помощью данных NSDIC во время замерзания и вскрытия является недопустимым. Для этого необходимо использовать дополнительные более точные данные.

Зависимость baklag49.wmf имеет схожий характер с зависимостью baklag50.wmf, максимальное значение baklag51.wmf достигается при значении ледовитости 0,6–0,7, а минимальное значение baklag52.wmf при полном ледоставе (рис. 5).

Отличием зависимости baklag53.wmf от зависимости baklag54.wmf является наличие большого значения baklag55.wmf в моменты, когда озеро полностью свободно ото льда. Такое значение baklag56.wmf учитывает летние периоды, когда по данным NOAA NESDIS имеют место ледовые явления на озерах, в то время как в действительности рассматриваемые озера в летнее время всегда полностью свободны ото льда. По этой же причине за период 2006–2018 гг. количество дней, когда на озерах наблюдались ледовые явления по данным NOAA NESDIS для всех озер больше на 15–30 %, чем по данным NSDIC.

Выводы

Установлено, что данные NSDIC имеют меньшее значение baklag57.wmf (3–10 %), чем данные NOAA NESDIS (11–19 %) для озер (Онежского, Ладожского озер, Сегозера, Выгозера) в сравнении с данными MODIS. Поэтому при определении пространственного распределения льда на озерах целесообразно использовать данные NSDIC.

Во время замерзания и вскрытия по данным NSDIC (а также NOAA NESDIS) baklag58.wmf могут достигать неприемлемого значения 40 % (значения baklag59.wmf до 35 %). В связи с этим формирование подробного временного ряда значений ледовитости во время замерзания и вскрытия с использованием спутниковых данных NSDIC, NOAA NESDIS возможно только при корректировке их другими более точными данными, например, данными MODIS.

Однако данные NSDIC позволяют рассчитать значение ледовитости озер во время полного ледостава и отсутствия ледовых явлений на озере с достаточной точностью (baklag60.wmf и baklag61.wmf). Это позволяет сделать вывод о том, что применение данных NSDIC целесообразно для определения сроков и длительностей фаз ледового режима.

Работа выполнена при финансовой поддержке гранта Президента РФ для поддержки молодых российских ученых МК-3379.2018.5.