Прогноз глубины промерзания почвы с помощью простейшей регрессионной модели | Вестн. Том. гос. ун-та. 2014. № 387. DOI: 10.17223/15617793/387/39

Прогноз глубины промерзания почвы с помощью простейшей регрессионной модели

Исследованы особенности динамики сезонно-мерзлого слоя и климатических характеристик холодного периода года по данным наблюдений на метеостанции Бакчар (Томская область) за период с 1963 по 2011 г. Анализ временного хода среднемесячных значений метеорологических характеристик позволил выявить, что статистически значимые тенденции роста температуры обнаружены для октября (0,6°C/10 лет), февраля (0,93°C/10 лет), марта (0,7°C/10 лет) и мая (0,55°C/10 лет). Значимых тенденций в суммах месячных осадков не обнаружено. Высота снежного покрова увеличивается с февраля по апрель и снижается в начале холодного периода (в октябре и ноябре). Предложенный метод прогноза глубины промерзания грунта с помощью регрессионных соотношений может использоваться для определения средних значений глубины промерзания почвогрунтов.

Soil freezing depth forecast using simple regression model.pdf Температура почвы является ключевым фактором, контролирующим многие биотические и абиотические процессы, протекающие в почвах. Температурный режим почв формируется под влиянием современных климатических условий и факторов, характеризующих условия теплового взаимодействия почвы с атмосферой [1]. Современные изменения в температуре приземного воздуха и режиме атмосферных осадков вызывают соответствующие изменения температуры почвы в верхних слоях [2-4]. Температура почвы определяет глубину развития сезонно-мерзлого или сезонно-талого слоя (в зоне залегания многолетнемерзлых пород) [5]. Реакция грунтов на потепление климата в значительной мере зависит от ландшафтных условий территории [1]. Современный рост снегонакоплений в Сибири создает частичную изоляцию почвогрунтов от температурных изменений воздуха в зимнее время [3]. Тенденция изменений температуры почвы не всегда повторяет динамику температуры воздуха, более того, они могут быть разнонаправленными во времени и пространстве [6]. Исследование связей между характеристиками сезонно-мерзлого слоя и состоянием окружающей среды является важной задачей. В данной работе исследованы особенности динамики сезонно-мерзлого слоя и климатических характеристик холодного периода года по данным наблюдений на метеостанции Бакчар (Томская область) за период с 1963 по 2011 г. Показаны тенденции изменения этих характеристик и предложен метод прогноза глубины промерзания грунта с помощью регрессионных соотношений по среднемесячным значениям температуры воздуха и высоты снежного покрова. Среднесуточные данные по температуре воздуха, сумме атмосферных осадков, высоте снежного покрова и температуре почвы на глубинах 20, 40, 80, 160 и 320 см на метеостанции получены из массива данных ВНИИГМИ-МЦД (http://www.meteo.ru/climate). Расчетная глубина промерзания определялась по глубине проникновения температуры 0°C в почву [7] путем линейной интерполяции температур почвы между двумя соседними слоями при условии, что в одном из слоев температура отрицательна. По среднесуточным данным рассчитаны среднемесячные значения, которые и использовались для дальнейшего анализа. Метеорологические характеристики. По данным метеостанции Бакчар, среднегодовая температура воздуха за 1963-2011 гг. составляет -0,2°C, амплитуда годового хода -39,8°C (рис. 1). Самый теплый месяц - июль (18,3°C), самый холодный - январь (-19,2°C). Самым теплым месяцем был июль 1969 г. (средняя температура 22,6°C). В среднем за 19632011 гг. отрицательные значения среднемесячной температуры воздуха на метеостанции наблюдаются с ноября (-9,0°C) по март (-8,9°C). В отдельные годы отрицательные температуры отмечены в октябре (-5,6°C в 1976 г.) и апреле (-6,5°C в 1983 г.). В январе температура опускается до -19,2°C. Самой холодной была зима 1968-1969 гг., когда средние месячные температуры ноября (-18,4°C), декабря (-26,3°C) и февраля (-29,7°C) были минимальными за период исследования. Средняя температура холодного периода (с ноября по март) в 1968-1969 гг. составила -23,4°C, в то время как средняя многолетняя величина равна -14,1°C. Минимальная температура января (-31,8°C) наблюдалась в 2006 г., что на один градус ниже, чем в 1968 г. Самые сильные положительные аномалии в зимние месяцы наблюдаются в конце исследуемого периода: ноябрь (-2,8°C) - 2008 г.; декабрь (-7,4°C) - 2006 г.; январь (-9,3°C); февраль (-8,1°C) - 2002 г. Сумма годовых осадков составляет 471 мм, из них 43% приходится на летние месяцы и 14% - на зимние. Наибольшее количество месячных осадков было зарегистрировано в июле 1996 г. (184 мм). Минимум осадков в годовом ходе наблюдается в феврале (15 мм). Снежный покров появляется 10 октября. Он обычно сходит под влиянием последующих оттепелей, и только после устойчивого понижения температуры устанавливается на зиму. Установление снежного покрова происходит, в среднем, к 30 октября. Максимальная высота снега за зиму наблюдается в марте (58 см) в среднем за 1963-2011 гг. Зимой 19671968 гг. наибольшая высота снега за зиму не превышала 21 см. В марте 1990 г. снежный покров достигал 98 см. Разрушение устойчивого снежного покрова происходит к 21 апреля (рис. 1). Максимальная температура почвы на глубине 20 см - в июле (17°С). Минимальная (-2,4°С) в годовом ходе температура почвы на этой глубине приходится на февраль. - Та | | Р - H^F Рис. 1. Годовой ход среднемесячных значений температуры воздуха (Та), сумм осадков (Р), высоты снежного покрова (H), глубины промерзания почвы (F) и температуры почвы (Ts) на глубинах 20, 40, 80, 160 и 320 см "1-I-I-I-I-I-I-I-I-Г "I-Г m Q. Q-(D ГО 5 в ^ < а I d о; О С глубиной амплитуда годового хода температуры уменьшается и составляет 13,4°С на глубине 80 см и 4,4°С - на глубине 320 см. Также с глубиной происходит запаздывание наступления максимальных и минимальных значений температуры почвы. Так, на глубине 320 см максимум наблюдается в октябре и составляет 7,3°С, а минимум (2,8°С) - в мае. Средняя годовая температура в слое от 20 до 320 см равна 4,9°С. Глубина промерзания определялась по глубине проникновения нулевой изотермы. Промерзание почвы начинается в конце ноября, приблизительно через месяц после установления отрицательных температур воздуха и образования устойчивого снежного покрова. Максимум глубины промерзания достигается в марте и составляет 86 см. Наибольшая глубина промерзания (160 см) была зафиксирована в феврале 1968 г. на фоне аномально низкой высоты снега (21 см) и умеренно теплой зимы (температура января -20°С). Наименьшая глубина промерзания наблюдалась зимой 1982-1983 гг., когда она составила 43 см (в марте), при этом снежный покров достигал высоты 60 см (в феврале), а температура опускалась всего до -12°С (в феврале), что на 5,1 °С выше многолетней нормы. В отдельные годы (1969-1972 гг.) мерзлый слой наблюдался в мае на значительной глубине (88110 см). Тенденции изменения метеорологических характеристик. Анализ временного хода среднемесячных значений метеорологических характеристик позволил выявить значимые тенденции их изменения. Величины линейных трендов приведены в табл. 1. Скорость увеличения среднегодовой температуры воздуха составляет 0,33°C/10 лет. Значимые тенденции роста температуры обнаружены для октября (0,6°C/10 лет), февраля (0,93°C/10 лет), марта (0,7°C/10 лет) и мая (0,55°C/10 лет). В остальные месяцы также обнаружены тенденции к потеплению, которые не являются статистически значимыми. Исключение составляют январь и июль, когда проявляются незначимые тренды уменьшения температуры. Статистически значимых тенденций в суммах месячных осадков не обнаружено (табл. 1). Т а б л и ц а 1 Коэффициенты линейного тренда температуры воздуха (aTa, °C/10 лет), сумм осадков (аР, мм/10 лет), высоты снежного покрова (аН, см/10 лет) и глубины промерзания почвы (аГ, см/10 лет) в разные месяцы Месяц аТа аР aH aF Январь -0,07 -1,06 2,39 -3,62 Февраль 0,93 0,03 4,34 -7,24 Март 0,70 1,49 5,93 -10,28 Апрель 0,44 -1,01 2,96 -13,00 Май 0,55 -3,41 - - Июнь 0,19 4,07 - - Июль -0,16 4,78 - - Август 0,16 -0,09 - - Сентябрь 0,07 2,33 - - Окттябрь 0,60 0,96 -0,46 -0,03 Ноябрь 0,37 -1,05 -1,68 1,48 Декабрь 0,22 2,24 -0,19 0,91 Год 0,33 9,28 - - Примечание. Полужирным шрис зтом выделены коэффициенты, значимые на уровне 0,05. Высота снежного покрова за 1963-2011 гг. увеличивается с февраля по апрель и снижается в начале холодного периода (в октябре и ноябре). Обнаруженные тенденции увеличения температуры воздуха в зимние месяцы и высоты снежного покрова приводят к появлению направленных изменений в температуре почвы и, следовательно, глубине промерзания почвы. Тенденция к снижению глубины промерзания получена для февраля, марта и апреля (см. табл. 1). В октябре одновременный рост температуры воздуха (что несколько уменьшает промерзание) и снижение высоты снежного покрова (что способствует увеличению глубины промерзания) не приводят к изменениям глубины промерзания почвы. Прогноз глубины промерзания. Для выявления факторов, влияющих на глубину промерзания почвы, был проведен корреляционный анализ. Значимые связи существуют между глубиной промерзания и температурой воздуха в ноябре (коэффициент корреляции R = -0,55) и марте (R = -0,38), а между глубиной промерзания и высотой снежного покрова - в декабре (R = -0,30), январе (R = -0,39), феврале (R = -0,53) и марте (R = -0,54). Увеличение температуры воздуха и высоты снежного покрова приводит к росту температуры почвы и, соответственно, уменьшению глубины промерзания. В начале развития сезонно-мерзлого слоя температура воздуха играет ведущую роль в промерзании почвы. В зимние месяцы роль управляющего фактора играет глубина снежного покрова. Развитие сезонно-мерзлого слоя в течение зимы идет постепенно. Очевидно, что толщина мерзлого слоя текущего месяца связана с толщиной мерзлоты предыдущего месяца. Коэффициенты корреляции глубины промерзания в соседние месяцы имеют высокие значения (R = 0,81 0,93). Поскольку на развитие сезонной мерзлоты в разные месяцы оказывают влияние различные факторы, нами была предложена множественная регрессионная модель для прогноза среднемесячной толщины мерзлого слоя. В общем виде в качестве независимых переменных модель включает среднемесячную температуру воздуха и высоту снежного покрова текущего и предыдущего месяца, а также глубину промерзания предыдущего месяца. Fm = a0 + a1'Fm-1 + a2'Tm-1 + a3'Tm + a4'Sm-1 + a5■Sm, (1) где a0-a5 - параметры модели; m - номер месяца (от ноября до апреля); Fm, Tm, Sm - глубина промерзания, температура воздуха и высота снежного покрова месяца m; Fm-1, Tm-1, Sm_i - те же характеристики для предыдущего месяца (m-1). Параметры модели определены c помощью программного комплекса Statistica 8.0 по имеющимся данным среднемесячных значений глубин промерзания, температур воздуха и высоты снежного покрова за 1963-2011 гг. Значения параметров и величина коэффициента детерминации модели приведены в табл. 2. Т а б л и ц а 2 Коэффициент ы регрессионной модели (1) для разных месяцев Месяц а0 а1 a2 аъ а4 а5 R2 Ноябрь 3,08 - - -2,44 - -0,85 0,51 Декабрь 31,80 0,71 -2,17 - -1,67 - 0,79 Январь 32,98 0,95 - -0,64 - -0,50 0,83 Февраль 17,30 0,86 -0,66 -0,62 0,65 -0,82 0,95 Март -14,95 0,93 -0,79 -1,75 - - 0,93 Апрель -42,47 1,10 -1,90 -3,00 - - 0,89 Набор независимых переменных определялся для каждого месяца индивидуально. Например, глубина промерзания ноября зависит только от температуры воздуха и высоты снежного покрова в ноябре. Мерзлый слой в предыдущий месяц (октябрь) наблюдается только в отдельные годы, а температура воздуха и снежный покров в октябре не имеют значимого влияния на промерзание (см. знак «-» в табл. 2). Наилучший набор параметров для прогноза глубины промерзания в декабре содержит глубину мерзлого слоя, температуру воздуха и высоту снега в ноябре и не включает температуру и высоту снега декабря. Эти характеристики декабря не входят и в регрессионное соотношение для января. Получено, что в марте и апреле температуры воздуха и глубины промерзания в предыдущий месяц достаточно для описания изменений глубины промерзания. Вклад снежного покрова в эти месяцы не является значимым. Регрессионные модели для декабря - апреля имеют высокие коэффициенты детерминации (R2 = 0,79 0,95). Несколько хуже показатели для ноября (R2 = 0,51), когда средняя за месяц глубина промерзания будет в первую очередь определяться временем установления отрицательных температур и датой начала промерзания. Как упоминалось выше, в отдельные годы мерзлый слой наблюдается на значительной глубине в мае, однако в последние годы освобождение почвы от мерзлоты сдвигается на более ранние сроки. Так, во второй половине исследуемого периода (1987-2011) мерзлота в мае была зафиксирована только в 7 случаях. Поэтому описанный выше способ прогноза толщины мерзлого слоя для мая не использовался. Для расчета многолетнего хода глубины промерзания сначала определялась глубина промерзания почвы в ноябре по известным значениям температуры воздуха и высоты снежного покрова. Затем последовательно рассчитывались глубины промерзания для каждого месяца с декабря по март (Fm) по соотношению (1). Причем в качестве независимой переменной Fm-1 использовалась расчетная глубина промерзания предыдущего месяца. Для температуры воздуха и высоты снежного покрова использовались фактические среднемесячные значения, полученные по данным измерений на метеостанции. Использование расчетной глубины промерзания для вычисления толщины мерзлого слоя в последующие месяцы приводит к последовательному накоплению ошибок. Средняя ошибка аппроксимации в течение зимы возрастает от 7 см в ноябре до 15 см в марте и 19 см в апреле. Тем не менее результаты расчетов удовлетворительно совпадают с данными фактических измерений. Корреляционная связь между рядами значений измеренной и расчетной глубины промерзания высокая (R = 0,89; R2 = 0,79). Предложенный метод позволяет оценить глубину промерзания почвы по данным широко распространенных наблюдений за температурой воздуха и снежным покровом. Для метеостанции Бакчар имеются данные измерений температуры и высоты снежного покрова с 1936 по 2013 г. Использование этих данных и регрессионного уравнения (1) дало возможность реконструировать глубину промерзания за этот период времени для ноября - апреля. На рис. 2 приведен временной ход глубины промерзания в марте, реконструированный для периода с 1936 по 2013 г. На совпадающих частях рядов согласие между фактической и расчетной глубинами промерзания 160 -| 1950 1960 1970 1980 1990 2000 Рис. 2. Фактическая (1) и реконструированная (2) глубина промерзания в марте 120 - 40 - 1940 достаточно хорошее, особенно в начале и середине интервала. Наибольшие ошибки прогноза характерны для конца исследуемого периода (см. рис. 2). По-видимому, прогноз, основанный на неизменности статистических характеристик рядов, в условиях направленных изменений независимых параметров модели работает плохо. Для качественного прогноза сезонной и многолетней динамики сезонно-мерзлой толщи необходимо использовать модели тепло-массопереноса в почвогрунтах. В работе исследованы особенности динамики се-зонно-мерзлого слоя и климатических характеристик холодного периода года по данным наблюдений на метеостанции Бакчар (Томская область) за 19632011 г. Анализ временного хода среднемесячных значений метеорологических характеристик позволил выявить, что значимые тенденции роста температуры обнаружены для октября (0,6°C/10 лет), февраля (0,93°C/ 10 лет), марта (0,7°C/10 лет) и мая (0,55 °C/10 лет). Статистически значимых тенденций в суммах месячных осадков не обнаружено. Высота снежного покрова за 1963-2011 гг. статистически значимо увеличивается с февраля по апрель и снижается в начале холодного периода (октябрь и ноябрь). Предложенный метод прогноза глубины промерзания грунта с помощью регрессионных соотношений по среднемесячным значениям температуры воздуха и высоты снежного покрова может использоваться для определения средних значений глубины промерзания за многолетние интервалы. Для качественного прогноза межгодовой динамики сезонно-мерзлой толщи необходимо использовать модели тепло-массопе-реноса в почвогрунтах.

Ключевые слова

температура почвы, сезонно-мерзлый слой, снежный покров, температура воздуха, глубина промерзания, soil temperature, seasonally-frozen soil, snow cover, air temperature, freezing depth

Авторы

ФИООрганизацияДополнительноE-mail
Дюкарев Егор АнатольевичИнститут мониторинга климатических и экологических систем Сибирского отделения Российской академии наукканд. физ.-мат. наук, вед. науч. сотр. лаборатории физики климатических системegor@imces.ru
Всего: 1

Ссылки

Павлов А.В. Мониторинг криолитозоны. Новосибирск : Гео, 2008. 229 с.
Васильев А.А., Дроздов Д.С., Москаленко Н.Г. Динамика температуры многолетнемерзлых пород Западной Сибири в связи с изменениями климата // Криосфера Земли. 2008. Т. 12, № 2. С. 10-18.
Шерстюков А.Б. Корреляция температуры почвогрунтов с температурой воздуха и высотой снежного покрова // Криосфера Земли. 2008. Т. 12, № 1. С. 79-87.
Трофимова И.Е. Современное состояние и тенденции многолетних изменений мерзлотно-термического режима почв Прибайкалья // География и природные ресурсы. 2006. № 4. С. 38-45.
Мельников Е.С., Васильев А.А., Лейбман М.О., Москаленко Н.Г. Динамика сезонноталого слоя в Западной Сибири // Криосфера Земли. 2005. Т. 9, № 2. C. 23-32.
Трофимова И.Е., Шеховцов А.И. Оценка термического режима почв котловин Прибайкалья и Северного Забайкалья // География и при родные ресурсы. 2011. № 4. С. 100-107.
Методы климатологической обработки метеорологических наблюдений / под общ. ред. О.А. Дроздова. Л. : Гидрометеоиздат, 1957. 492 с.
 Прогноз глубины промерзания почвы с помощью простейшей регрессионной модели | Вестн. Том. гос. ун-та. 2014. № 387. DOI: 10.17223/15617793/387/39

Прогноз глубины промерзания почвы с помощью простейшей регрессионной модели | Вестн. Том. гос. ун-та. 2014. № 387. DOI: 10.17223/15617793/387/39