Analysis of Sentinel 1 Radar Data to Identify Frozing/Thawing Soils in the Area of Anadyr (Chukotka) and Belaya Gora (Yakutia)

Full Text

Дистанционное зондирование (ДЗ) предлагает перспективные методы для мониторинга приповерхностного состояния мерзлых и талых почв в широких географических масштабах (Khaldoune et al., 2008, 2011). Использование радиометров ограничено очень большими территориями из-за низкого пространственного разрешения. В отличие от ИК и микроволновых радиометров РСА (радар с синтезом апертуры) системы могут давать информацию о динамическом состоянии почвы и под снежным покровом (сухой снег) при высоком пространственном разрешении вне зависимости от наличия облачности и времени суток. Сигнал радара преимущественно зависит от диэлектрической проницаемости (ДП) почвы, которая напрямую связана с содержанием воды и льда.

Активные сенсоры различают мерзлую почву (МП) через изменение содержания жидкой воды в почве. Эти изменения связаны со значениями ДП почвы (Hallikainen et al., 1985). Среднее значение ДП сухой почвы порядка 2–3, в то время как жидкой воды 80. ДП почвы растет пропорционально объемному содержанию в ней жидкой воды. В критических зимних условиях большая часть воды в почве замерзает, что приводит к значительному уменьшению среднего значения ДП почвы. Значение диэлектрической проницаемости мерзлой почвы близко к значению ДП льда (3.2), что схоже с ДП сухой почвы (Boyarsky, 2003). Обратный процесс происходит весной, когда лед в почве тает, что приводит к увеличению коэффициента обратного рассеяния (КОР) на несколько децибел (Ulaby et al., 1982).

Свойства почвы (структура, текстура, дренаж) влияют на содержание воды в почве и, следовательно, на динамику процесса ее замерзания. Минерализованные почвы замерзают быстрее и на большую глубину, нежели органические почвы, из-за их меньшей возможности сохранять воду (Nyberg et al., 2001).

Не только влажность почвы влияет на КОР, но и шероховатость поверхности и растительный покров.

ХАРАКТЕРИСТИКА РАЙОНОВ ИССЛЕДОВАНИЯ

Белая Гора

Бе́лая Гора́ — поселок городского типа, административный центр Абыйского улуса Республики Саха (Якутия), находится за полярным кругом. Координаты 68.537847°, 146.184769°. Координаты метеостанции (WMO 24194) 68°32′ с. ш. 146°12' в. д.

Большую часть территории Абыйского улуса (района) занимает Абыйская низменность, расположенная в среднем течении р. Индигирка. Преобладают многолетнемерзлые породы и мерзлотные формы рельефа. Низменность сложена различного рода морскими, речными и озерными отложениями с высоким содержанием ископаемого льда в силу преобладания многолетней мерзлоты. Ландшафт в основном представлен лесотундрой, в растительности преобладает даурская лиственница. Климат варьируется от резко континентального на юге до субарктического и арктического на севере, средняя температура января –40.7⁰C, июля +14⁰C. За год выпадает в среднем 218 мм осадков в виде дождя и снега. Зимы суровые, холодные и продолжительные. Территория района является полюсом холода Северного полушария. Бывают периоды, когда температура опускается до –70⁰C.

На рис. 1 показано спутниковое изображение Белой Горы (стрелка указывает на расположение поселка) и график изменений температуры воздуха по данным метеостанции вблизи поселка (68°32' с. ш., 146°12' в. д.) за период 16.10.2014–29.10.2016, представленные на сайте rp5.ru.

Анадырь

Анадырь (Чукотка) расположен на широте 64.7337° с. ш., 177.4968° в. д. на правом берегу близ устья р. Казачка, впадающей в Анадырский залив Берингова моря, в зоне вечной мерзлоты. Центральную часть округа занимает сильно заболоченная Анадырская низменность — обширная равнина на Чукотке с высотами до 100 м над уровнем моря с множеством озер. Вечная мерзлота залегает повсеместно и начинается очень неглубоко от поверхности. Координаты метеостанции (WMO 25563) 64.78°с.ш. и 177.57° в.д.

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

На рис. 2а дано спутниковое изображение Анадыря, где в верхней части красным цветом показано положение метеостанции, а на рис. 2б представлен график изменений температуры воздуха по данным метеостанции за период 18.10.2014–12.11.2016, находящихся в архиве на сайте rp5.ru.

ИСХОДНЫЕ РАДАРНЫЕ ДАННЫЕ SENTINEL 1

В работе использованы находящиеся в открытом доступе Sentinel 1 (S1) радарные данные С‑диапазона IW (interferometric wide swath) моды с поляризацией VV и пространственным разрешением 10 м. Работа с изображениями S1 осуществлялась с помощью программы S1Toolbox и позднее SNAP (https://sentinel.esa.int/web/sentinel/toolboxes/sentinel‑1).

Предварительная обработка данных включала в себя выделение фрагмента с исследуемой областью и радиометрическую калибровку.

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

Изменение КОР для района Белой Горы с координатами 68.505° с. ш. и 146.155° в.д (метеостанция) и Анадыря с координатами 64.78° с. ш. и 177.57° в. д. (метеостанция) за указанный в таблице период съемки показано на рис. 3.

Рис. 1. (а) — спутниковое изображение Белой Горы (стрелка); (б) — температура воздуха в градусах Цельсия для Белой Горы за период 16.10.2014–29.10.2016 по данным rp5.ru.

Рис. 2. (а) - Анадырь, нахождение метеостанции отмечено красным цветом; (б) - температура воздуха в градусах Цельсия за период 18.10.2014–12.11.2016.

Рис. 3. Коэффициент обратного рассеяния в районе метеостанций Анадыря с координатами 64.78° с. ш., 177.57° в. д. за период 18.10.2014–12.11.2016 и Белой Горы с координатами 68.505° с. ш., 146.155° в. д. за период 16.10.2014–29.10.2016.

Таблица 1. Информация об использованных радарных данных Sentinel 1

Коэффициент корреляции Спирмена между температурой воздуха и КОР для Анадыря равен ρ_S = 0.36 (р = 0.04, N = 25, где N — количество сеансов съемки), и для Белой Горы ρ_S = 0.69 (р = 1*10–5, N = 30).

ОПРЕДЕЛЕНИЕ СОСТОЯНИЯ МЕРЗЛОЙПОЧВЫ ПО РАДАРНЫМ ДАННЫМ

Определение состояния МП по многовременным радарным данным

Глубина проникновения э/м волны в почву определяется по формуле

$\delta p\approx \frac{\lambda \cdot \sqrt{{\epsilon }^{\text{'}}}}{2\cdot \pi \cdot {\epsilon }^{\text{'}\text{'}}}$ ,

где λ — длина э/м волны, ε = εʹ − iεʹʹ — диэлектрическая проницаемость почвы. Для Sentinel 1 длина волны равна 5.4 см. В случае мерзлой почвы с εʹ ≈ 5.5 и 0.1 ≤ εʹʹ ≤ 0.5 получаем глубину проникновения э/м волны в почву 4 ≤ δp ≤ 20 см.

Метод идентификации талого и мерзлого состояния почвы основывается на перепаде КОР на 3–5 дБ при переходе к состоянию замерзания или оттаивания (Khaldoune et al., 2008, 2011).

Используя временной ряд S1, определим даты замерзания/оттаивания МП по значительным перепадам КОР и далее стабильно низким/высоким значениям σ⁰. На рис. 4 показан график разности значений по абсолютной величине σ⁰ в дБ для соседних дат съемки Белой Горы, где по максимумам | σ⁰ | определяем даты замерзания и оттаивания почвы. Первый максимум равен max = 3.35 дБ = | σ⁰_{20.5.2015 –} σ⁰_13.6.2015 | [дБ] — даты оттаивания почвы. Второй максимум определяет даты замерзания почвы, мах = 1.97 дБ = | σ⁰_17.9.2015 – σ⁰_11.10.2015 |. Оба максимума для 2015 г. Для 2016 г. первый максимум равен max = 2.74 дБ= | σ⁰_14.5.2016 – σ⁰_7.6.2016 | — это даты оттаивания почвы. Второй максимум определяет даты замерзания почвы, мах = 1.37 дБ = |σ⁰_11.9.2016  – – σ⁰_5.10.2016 |. Следует отметить, что точность дистанционного определения времени замерзания и оттаивания почвы определяется частотой съемки исследуемой территории.

Рис. 4. Разность значений | σ0 | в дБ для соседних дат радарной съемки Белой Горы.

Рис. 5. Разность значений | σ0 | в дБ для соседних дат радарной съемки Анадыря.

На рис. 5 показан аналогичный рис. 4 график для Анадыря. Что касается определения периода оттаивания почвы, то из-за отсутствия съемки под определенным углом для Анадыря, а именно — не было съемки с 17.04.2015 по 26.06.2015 и с 11.04.2016 по 14.07.2016, указанные на графике точки весьма отдаленно отвечают за даты оттаивания почвы. Период замерзания почвы для 2014 г., max =1.8 дБ = =| σ⁰_18.10.2014 – σ⁰_11.11.2014 |, и для 2015 г., max=2.4 дБ==| σ⁰_18.11.2015 – σ⁰_12.12.2015 |.

Следует отметить, что чем более резкий по перепадам температур климат, тем надежнее по радарным данным можно определить даты оттаивания/замерзания почвы в приповерхностном слое толщиной 5 см.

Определение состояния МП по факторусостояния поверхности и пороговому значению σ⁰

Еще один способ определения состояния МП только по радарным данным, не используя наземные измерения, приведен в работе (Mironov, Muzalevsky, 2013). Авторы предложили уравнение, характеризующее состояние поверхности почвы, аналогом для которого послужило уравнение в работе (Kim et al., 2011) для фактора, связанного с радиояркостной температурой почвы.

Рис. 6. Изменение СПФ для Белой Горы.

Рис. 7. Радарное изображение места вблизи Белой Горы (дата съемки 20.05.2015), где белый прямоугольник выделяет участок с определением замерзшей/оттаявшей почвы по σ0порог (справа), и google map данного участка (слева).

Авторы (Mironov, Muzalevsky, 2013) характеризуют состояние поверхности почвы через фактор состояния поверхности (СПФ)

$\text{СПФ}=\frac{1}{2}+\frac{{\sigma }^0\left(t\right)-{\sigma }_{\text{лето}}^0}{{\sigma }_{\text{лето}}^0-{\sigma }_{\text{зима}}^0}$, (1)

Где — средние значения КОР летом и зимой; σ⁰(t) — текущее значение КОР; t –время.

Авторы (Mironov, Muzalevsky, 2013) утверждают, что сезонные вариации СПФ сильно коррелируют с сезонными вариациями температуры почвы, принимая отрицательные значения при температуре почвы ниже 0°C. Следовательно, 5-см верхний слой почвы можно классифицировать как замерзший, если значение СПФ, измеренное радаром, становится отрицательным. На рис. 6 показан график СПФ для Белой Горы.

Для Белой Горы σ⁰ _лето и σ⁰ _зима определяем по радарным данным между точками замерзания/оттаивания почвы. По формуле (1) получается, что пересечение значения СПФ = 0 происходит для одного и того же значения σ⁰(t), но для разного времени t, т. е. значение σ⁰ одно и то же как для точки замерзания, так и для точки оттаивания почвы на глубине 5 см. И это значение можно определить для каждого места, зная средние значения σ⁰ _лето и σ⁰ _зима . Назовем это значение σ⁰ через σ⁰_порог. Для Белой Горы значения σ⁰_порог следующие:

2014–2015 гг.: σ⁰ _лето= –11.6 дБ, σ⁰ _зима = –15.2 дБ, σ⁰_порог = –13.4 дБ,

2015–2016 гг.: σ⁰ _лето = –12.3 дБ, σ⁰ _зима = –14.1 дБ, σ⁰_порог = –13.2 дБ.

Используя значения σ⁰_порог, можно построить локальные карты мест замерзания/оттаивания почвы в окрестности рассматриваемого района исследования. На рис. 7 показано радарное изображении района вблизи Белой Горы за 20.05.2015, где белым прямоугольником выделен участок размером порядка 4 × 4 км, для которого определялась замерзшая и оттаявшая почва по σ⁰_порог. На рис. 8 приведены изображения после фильтрации спеклов за 20.05.2015 и 13.6.2015 (период оттаивания почвы) и для дат 17.09.2015 и 11.10.2015 (период замерзания почвы) с порогом σ⁰_порог = –13.4 дБ. Синий цвет — значенияσ⁰ < σ⁰_порог (замерзшая почва), желтый цвет — значения σ⁰ > σ⁰_порог (оттаявшая почва). Отмечаем значительное уменьшение площади замерзшей почвы на изображении за 13.06.2015 при оттаивании и, соответственно, увеличение площади замерзшей почвы на изображении за 11.10.2015 при замерзании.

Рис. 8. Белая Гора, 20.05.2015 и 13.6.2015 (период оттаивания почвы), 17.09.2015 и 11.10.2015 (период замерзанияпочвы). Синий цвет — значения σ0 < σ0порог = 13.4 дБ (замерзшая почва), желтый цвет – значения σ0 > σ0порог (оттаявшая почва).

Определение состояния МП по текстурным признакам

Радарная съемка S1 территории Анадыря и Белой Горы осуществлялась только с одной поляризацией VV. Но и для одноканальных изображений есть возможность выделить разные типы поверхности, используя текстурные признаки Харалика (Haralick, 1979). В работе (Родионова, 2007) показана такая возможность на примере использования трех текстурных признаков — «контраста», «энтропии» и «обратного момента» — для изображений SIR-C/X-SAR. Контраст является взвешенным квадратом уровней яркости и служит естественной мерой распространения уровней яркости. Обратный момент является мерой локального сходства. Энтропия — мера беспорядка на изображении. Области с высоким значением контраста и энтропии выделяют городские застройки, области с высоким значением обратного момента — гладкие поверхности. Построив текстурное изображение, можно в качестве индикатора замерзшей/оттаявшей почвы использовать значение текстурного признака обратный момент, значение которого больше для замершей почвы. На рис. 9 показаны текстурные изображения после фильтрации спеклов для тех же дат, что и на рис. 8, в следующей RGB‑кодировке: красный цвет — контраст, зеленый — энтропия, синий — обратный момент. Сравнивая рис. 8 и 9, видим, что именно обратный момент показывает области с мерзлой почвой.

ИНДЕКС ПРОМЕРЗАНИЯ/ОТТАИВАНИЯ ПОЧВЫ

Температура воздуха, как известно (Brown, 1964), коррелирует с температурой почвы и глубиной промерзания почвы. Аккуратные предсказания глубины промерзания почвы являются важным аспектом в различных инженерных проектах. На реальную глубину промерзания влияют механический состав почвы, ее термические свойства, влажность почвы, климатические условия, такие как температура, скорость ветра, осадки, солнечная радиация. Глубину промерзания можно оценить, используя численные или аналитические методы, однако требуемые в моделях входные данные не всегда доступны или затратны для получения. Следовательно, использование численных или аналитических моделей подразумевает оценку точности при отсутствии необходимых входных данных, которые влияют на надежность результатов.

При наличии всех входных данных результаты таких моделей могут быть достаточно точными. С другой стороны, при отсутствии некоторых входных данных необходимо оценить, насколько это повлияло на точность получения результата. Из существующих моделей предсказания глубины промерзания почвы наиболее известны модель Стефана, используемая и в настоящее время, модели Berggren и Chisholm и Phang. В работе (Rajaei, Baladi, 2015) рассмотрена точность предсказания глубины промерзания почвы для этих трех аналитических и полуэмпирических моделей. Ввиду того, что ни одна из моделей не дает аккуратных результатов, авторами (Rajaei, Baladi, 2015) предложена модифицированная эмпирическая модель для разных типов почв, которая на входе требует только знания максимальных и минимальных температур воздуха за сутки.

Рис. 9. Текстурные изображения после фильтрации спеклов в RGB кодировке: красный цвет — контраст, зеленый — энтропия, синий — обратный момент.

Рис. 10. Значение AFI по годам с отрицательным трендом, характеризующее максимальную глубину промерзания почвы.

Рис. 11. Значение AFI по годам, характеризующее максимальную глубину оттаивания почвы. Отрицательный тренд для Белой Горы и положительный для Анадыря.

Авторы (Rajaei, Baladi, 2015) получили эмпирическую модель связи глубины промерзания почвы с суммарным индексом ее промерзания для двух типов почв: глинистой и песчаной, на основе измеренных значений глубин промерзания в 2010–2011 гг. в штате Мичиган и вычисленных значениях суммарного индекса промерзания CFI (the cumulative freezing index):

для глинистых почв

                            h = 4.0388 ∙ (CFI)^0.4896,

                        h = 3.3787 ∙ (CFI)^0.5423,

с коэффициентом детерминации 0.91 (число измерений 129);

$CFI=\sum _{\text{i}=\text{1}}^{\text{n}}\left(T_{\min }^{\text{воздух}}+{\text{T}}_{\text{max}}^{\text{воздух}}\right)\text{\hspace{0.33em}/2}$         

где h — глубина промерзания почвы в см, $T_{\min }^{\text{воздух}}$ и $T_{max}^{\text{воздух}}$— минимальная и максимальная температура воздуха за сутки, CFI — суммарный индекс промерзания (the cumulative freezing index).

Авторы (Rajaei, Baladi, 2015) получили обобщенное для всех типов почв эмпирическое уравнение с учетом измерения средней теплопроводности образцов почв в штате Мичиган, отмечая, что при использовании их эмпирических уравнений для других территорий необходима калибровка.

Максимальная глубина промерзания/оттаивания почвы для каждого года определятся по разнице максимального и минимального значений CFI за год. Максимум CFI соответствует точке перехода от положительных к отрицательным температурам, минимум CFI соответствует переходу от отрицательных температур к положительным. В первом случае разность между экстремумами соответствует оценке максимальной глубины оттаивания почвы, во втором — максимальной глубины промерзания почвы. Данная величина получила название AFI (air freezing index). На рис. 10 приведены значения AFI, характеризующие максимальную глубину промерзания почвы, для Белой Горы и Анадыря за 2012–2018 гг. с линиями тренда. В обоих случаях тренд отрицательный с коэффициентом детерминации R² = 0.59 и R² =0.62 соответственно, что свидетельствует об уменьшении максимальной глубины промерзания почвы для исследуемых территорий.

На рис. 11 приведены значения AFI, характеризующие максимальную глубину оттаивания почвы, для Анадыря и Белой Горы за 2012–2017 гг. с линиями тренда. Для Белой Горы тренд отрицательный, с очень малым значением R², а для Анадыря — положительный, с R² = 0.49.

Таким образом, однозначный вывод о тренде для максимальной глубины оттаивания почвы в исследуемых районах за 2012–2017 годы сделать нельзя.

ЗАКЛЮЧЕНИЕ

В работе рассматривается вопрос о различении талых и мерзлых почв в верхнем 5-см слое для наземных станций в Белой Горе (Якутия) и Анадыре (Чукотка) по радарным данным Sentinel 1 С‑диапазона за период 2014–2016 гг. Определение состояния мерзлой почвы осуществляется тремя способами: 1) по многовременным радарным данным на основе значительного на 3–5 дБ перепада значений коэффициента обратного рассеяния σ⁰ при переходе к состоянию замерзания/оттаивания почвы. Для Белой Горы максимальное значение перепада σ⁰  составило 3.35 дБ при оттаивании и 1.97 дБ при замерзании почвы и для Анадыря соответственно, 2.74 дБ и 1.37 дБ; 2) по нахождению порогового значения σ⁰_порог, при котором температура в верхнем слое почвы опускается ниже 0°С. Найденные значения σ⁰_порог позволили построить локальные карты мерзлой и талой почвы для участка вблизи Белой Горы; 3) Из-за наличия радарной съемки на исследуемые территории только на VV‑поляризации для проверки правильности определения мест мерзлой/талой почвы использовались текстурные признаки Харалика для сегментации одноканальных изображений. В результате показано, что ТП обратный момент идентифицирует области с мерзлой почвой.

Построены графики AFI, характеризующие максимальную глубину промерзания/оттаивания почвы, за период 2012–2018 гг. по архивным данным температуры воздуха для исследуемых территорий. Показано наличие отрицательного тренда AFI как для Белой Горы, так и для Анадыря c коэффициентом детерминации соответственно 0.59 и 0.62. Для оттаивания почвы однозначного ответа нет, поскольку для Белой Горы получен слабый отрицательный тренд, а для Анадыря — положительный тренд.

About the authors

N. V. Rodionova

Author for correspondence.
Email: rnv@ire.rssi.ru
Russian Federation, 11-7, Mokhovaya street, Moscow, 125009

References

Supplementary files