Причинами активизации оползневого процесса на территории жилой застройки в г. Воркуте являются эрозия вследствие мерзлотного и гигрогенного крипов и обводнение склона поверхностными водами. В статье дана оценка взаимосвязи динамики движения оползня, глубины сезонного промерзания грунтов и основных климатических параметров, влияющих на перечисленные показатели: мощность снежного покрова, средняя температура воздуха холодного периода, количество холодных дней, количество жидких осадков – за период 2010–2020 гг. Расчеты глубины промерзания проведены в условиях залегания техногенных грунтов и полускальных горных пород двумя методами, основанными: 1 – на термических свойствах грунтов и горных пород и некоторых климатических параметрах, 2 – на теплопроводности трехслойной среды (снежный покров, мерзлый и талый грунты) и их тепловом балансе. Среднее расхождение расчетных и некоторых наблюденных глубин промерзания техногенного грунта – 30 %. Коэффициент корреляции Ккор составил 0.86 и 0.62. По корреляционному анализу Грея средняя температура воздуха холодного периода явилась важным параметром влияния на эталонную последовательность – динамику движения оползня. Наиболее значимое влияние на глубину промерзания грунтов оказали длительность холодного периода и температура воздуха. В обеих последовательностях существенным влияющим параметром оказалось и количество жидких осадков.
глубина промерзания, оползень, климатические параметры, корреляционный анализ Грея
Введение
В регионах, где преобладает холодный климат и среднемесячная температура воздуха за холодный период ниже 0 °C, происходят циклы промерзания–оттаивания в грунтах. По некоторым данным [1, 2], примерно 50 % суши в мире подвержено вышеуказанным циклическим процессам, одновременно с которыми изменение климатических показателей по территории происходит неравномерно. В связи с этим меняются структурная целостность, физические и механические свойства горных пород. Температура воздуха, влажность, теплопроводность почвы и грунта, осадки и скорость ветра эффективно воздействуют на глубину сезонных промерзаний и оттаиваний грунта. Считается, что наиболее значимым фактором, влияющим на тепловой баланс геологической среды, является температура воздуха, в частности и на Восточно-Европейской равнине [3]. Формула Стефана – одна из первых, созданных для оценки глубины промерзания почвы и грунтов [4]. Позже различными методами многократно была доказана взаимосвязь между температурой воздуха и глубиной промерзания грунтов [5–9]. Также многолетние наблюдения за значениями глубины промерзания и оттаивания показали, что термический режим грунтов прямо пропорционально зависит от высоты снежного покрова, его плотности, теплопроводности и структуры, что оценено методами статистического анализа [3, 10–12]. В целом грунты территорий с наибольшим накоплением снега имеют более высокую температуру и промерзают на меньшую глубину. Более полувека назад В. А. Кудрявцев предложил уравнение для расчета глубины промерзания грунта, объединив попеременное отепляющее и охлаждающее действие снежного покрова на грунт в зависимости от режима снегонакопления и амплитуды годовых колебаний температуры воздуха [13]. В настоящее время на практике такие исследования широко проводятся в строительной климатологии, при проектировании, строительстве и эксплуатации инженерных сооружений [14], автомобильных и железных дорог, в сельском хозяйстве [15] и экологии [16]. Дополнительно стоит отметить, что в области распространения многолетнемерзлых толщ состояние криолитозоны отмечается как неустойчивое [17, 18], и поэтому там ускорено протекание экзогенных геологических процессов и явлений, таких как криогенное выветривание, оползни, обвалы и т. д.
В районе исследований г. Воркуты (67°30′ с. ш. 64°02′ в. д.) особого внимания заслуживает оползневый процесс на территории жилой застройки на ул. Шахтерская набережная. По результатам проведенного анализа параметров смещения точек бровки оползневого уступа за 2013–2019 гг. центральная часть оползневого склона имеет преобладающе обвально-осыпной механизм; южную и северную части можно отнести к группе оползней скольжения с нечастыми обвалами [19]. В период 2013–2018 гг. центральная и южная части проходили активную стадию оползания, далее перешли в состояние затухания, что соответствует началу пассивного периода, или нового цикла. До 2017–2018 гг. северная часть пребывала на подготовительном этапе, по окончании которого был отмечен рост трендов с невысоким ускорением. В годы наблюдений отмечалось незначительное повышение среднегодовой температуры воздуха за счет более теплых месяцев с отрицательной температурой воздуха. Наиболее тесная взаимосвязь динамики движения оползня и климатических параметров выявлена со средним значением температуры воздуха месяцев с положительными значениями, коэффициент корреляции Ккор составил 0,58. По данным Государственного казенного учреждения Республики Коми «Геокриологическая служба Республики Коми» (ГКУ РК «Геокриологическая служба»), природными причинами активизации оползневого процесса являются эрозия вследствие мерзлотного и гигрогенного крипов и обводнение склона поверхностными водами.
В связи с этим цель данной работы – оценка глубины сезонного промерзания горных пород на оползневом склоне в долине р. Воркуты и определение взаимосвязи динамики движения оползня, глубины сезонного промерзания грунтов и иных основных климатических параметров, влияющих на перечисленные показатели.
Климат и геологические условия
Город Воркута расположен в северной части Республики Коми, в южной тундре, где распространена голоценовая и голоценово-плейстоценовая криолитозона, представленная прерывистым массивно-островным распространением многолетнемерзлых пород (10–50 % покрытия) [20]. Ее подошва залегает на глубине 100 м и чуть более [21]. Климат территории – субарктический. По состоянию на 2023 г. за 76 лет среднегодовая температура воздуха составила -5.3 °C, среднегодовое количество осадков – 540 мм. Снежный покров окончательно формируется в середине октября и разрушается в конце мая. Средняя продолжительность периода со снежным покровом – 232 дня. Средняя высота покрова – 0.75 м, максимальная – 1.8 м.
Среднегодовая температура воздуха по метеостанции г. Воркуты за период исследований составила -3.88 °C, среднее значение отрицательных среднесуточных температур воздуха равно -10.99 °C, положительных – +9.16 °C, среднегодовое количество осадков – 549 мм. За период исследований в 2016 г. отмечены минимальные годовое количество осадков (424 мм) и высота снежного покрова (0.43 м) [19].
Рельеф оползневого склона – денудационный, представлющий собой эрозионный склон современной речной долины. Схема оползневого склона показана на рис. 1. Вдоль бровки оползня построен геоэлектрический разрез (рис. 2). По материалам [22], в верхней части разреза залегают техногенные образования (tH) неоднородного состава мощностью 0.4–1.5 м: щебень, дресва, галька и гравий, встречаются включения суглинка. До глубины 3.0 м состав насыпных грунтов может быть несколько иной: гравий, галька, включения песка и суглинка. Содержание крупных фракций (d > 2 мм) – около 20 %. Такая неоднородность подтверждается большим диапазоном значений удельного электрического сопротивления горизонта – 280–970 Ом·м. Далее до глубины 4.0–6.8 м следуют среднечетвертичные отложения, представленные ледниково-морскими суглинками (gmE–lrg) с увеличением содержания гравия и гальки к подошве и возможными прослоями песков и супесей, среднее содержание – от 5–10 до 20–25 % [23, 24]. Ниже (до 13–14 м) залегают полускальные пермские отложения (P1-2lv2) и представляют элювий коренных пород: переслаивание сильно выветрелых, трещиноватых, глыбо–щебенистых песчаников, алевролитов, аргиллитов, углистых аргиллитов и углей. Также они могут переслаиваться со скальными прочными и средней прочности породами. Скальный горизонт выделяется как слой удельного электрического сопротивления 275–460 Ом·м. В районе исследований находится пермский подмерзлотный водоносный комплекс терригенных отложений.
Материалы и методы
В данном исследовании глубина промерзания грунта рассчитана двумя эмпирическими методами. Затем с помощью корреляционного анализа Грея определили основные параметры, влияющие на глубину промерзания грунтов и динамику движения оползня. Метод основан на оценке геометрического сходства между главной (эталонной) последовательностью и последовательностями параметров влияния для определения тесноты связи [25]. Эталонными последовательностями в анализе выступали глубина промерзания грунтов и динамика движения оползня.
Уравнение Стефана для определения глубины промерзания горных пород при наличии снежного покрова выглядит следующим образом:
,
где λм – коэффициент теплопроводности мерзлой породы, Вт/м·°С; tn – средняя температура воздуха, °С; τ – продолжительность промерзания на конец месяца холодного периода, ч; Qф – теплота фазового перехода, Вт·ч/м3; Сом – объемная теплоемкость мерзлых пород, Вт·ч/м·°С; hm – условная мощность слоя теплоизоляции (снежного покрова), м.
Расчеты глубины промерзания осуществляли на основе задачи теплопроводности трехслойной среды (снег, мерзлые и талые породы) с фазовым переходом на границе мерзлого и талого грунтов по [26]. В расчетах принято, что температура в каждой из сред изменяется линейно. Уравнение теплового баланса, лежащее в основе расчетов, включает энергию фазового перехода, приток тепла из немерзлого грунта и отток в мерзлый грунт и атмосферу.
F1=cLV+ F2,
где F1 – отток тепла через замерзший грунт и снежный покров от фронта промерзания в атмосферу, Вт/м2; L – энергия фазового перехода, 335 кДж/кг; V – скорость движения фронта промерзания, см/с; F2 – отток тепла на охлаждение талого грунта перед фронтом промерзания, Вт/м2.
Тепловой поток выражается по закону Фурье через градиент температуры и теплопроводность и может быть рассчитан:
,
где Tвозд – средняя температура воздуха, °С; hc и hмг – толщина снежного покрова и глубина промерзания горной породы, м; λс и λмг – теплопроводность снега и мерзлого грунта, Вт/м·°С.
Отток тепла на охлаждение талого грунта перед фронтом промерзания определяется по формуле:
,
где λтг – теплопроводность талого грунта, Вт/м·°С; Т0 – температура отложений на глубине нулевых годовых колебаний. По данному методу принято считать, что глубина нулевых годовых колебаний температуры грунтов – 10 м со значением 2 °С, на начало холодного периода толщина мерзлого грунта принималась 0.005 м.
Результаты и их обсуждение
Необходимые основные климатические параметры, влияющие на глубину замерзания горных пород: средняя температура воздуха холодного периода, количество жидких осадков, количество холодных дней, высота снежного покрова – были предсталены метеостанцией г. Воркуты. Годовые данные 2010–2020 гг. – в табл. 1. В указанный период наблюдений отмечаются незначительные повышения средней температуры воздуха для холодного периода и средней высоты снежного покрова с небольшими флуктуациями (коэффициенты линейной регрессии положительны, коэффициенты детерминации R2=0.2 и 0.6 соответственно), при этом коэффициент корреляции данных параметров равен 0.6. По средним значениям отрицательных температур воздуха самым теплым оказался период 2019–2020 гг., самым холодным – 2010–2011 гг. Самому теплому периоду соответствует одна из многоснежных зим, средняя высота снежного покрова составила 0.55 м. 2020 год отмечен и высоким количеством жидких осадков – 292.8 мм. Максимальное значение – 302 мм в 2019 г. Самый холодный период был скуден на осадки: характеризуется относительно низким значением мощности снежного покрова (0.35 м) и количеством жидких осадков (201.5 мм).
По результатам анализа ежемесячных показателей климатических параметров можно сказать, что в октябре наступают условия для устойчивого сезонного промерзания верхнего слоя горных пород. Средняя температура воздуха составляет -3 °С при разбросе -0.2...-6.7 °С.
При расчетах глубины промерзания для повышения точности получаемых результатов учитывались данные о теплоемкости мерзлых и талых грунтов. По методу Стефана предполагалось, что скрытая теплота влаги в грунтах является единственной теплотой, которая отводится при замерзании. Тепловая энергия, удерживаемая в объемном виде, выделяемая при условии, когда температура почвы и грунтов падает до точки замерзания и ниже, не учитывалась. Расчетная схема, предложенная Фроловым и др. [26, 27], построена на основе перехода теплового потока в трехслойной среде: снежный покров, мерзлый и талый грунты. Для начальных условий толщина мерзлого слоя сверху была принята 0.005 м. При изменении температуры воздуха и снега в снеге варьируется содержание кристаллической, жидкой и парообразной фаз воды, поэтому в течение холодного периода плотность снега и его теплопроводность меняются. По данным ГКУ РК «Геокриологическая служба», плотность снега в районе исследований варьирует в диапазоне 180–600 кг/м3, поэтому для расчетов принята плотность 450 кг/м3. Ей соответствует коэффициент теплопроводности снега λс=0.49 Вт/м·°С. Режим промерзания грунтов определялся для верхней части оползневого склона в районе размещения термометрических скважин № 1 и 2 (рис. 2), где в верхней части разреза до глубины 3 м находятся техногенные грунты, и для склоновой части, где обнажены полускальные горные породы. Результаты расчетов и их сравнение между собой и с некоторыми наблюдаемыми данными отображены на рис. 3 и 4. Среднее расхождение расчетных и наблюденных глубин промерзания техногенного грунта составляет 30 % при значительных и сильных коэффициентах корреляции 0.67–0.92. Расчетные данные по схеме, предложенной Фроловым и др., несколько ниже фактически измеренных, но в целом они воспроизводят динамику процесса промерзания, как и значения, полученные по методу Стефана (рис. 3). В разные годы различия в глубине промерзания составляли около 50 %. Основное промерзание грунта происходит к декабрю–январю – до 90 % от всего сезонного промерзания (рис. 3). Небольшой снежный покров и низкие температуры воздуха в начале холодного периода способствуют такому стремительному промерзанию. С января при установлении большей мощности снежного покрова графики динамики промерзания становятся пологими и имеют рост от 0.02 до 0.3 м/мес. Коэффициенты корреляции мощности снежного покрова и расчетных глубин промерзания, определенные по методам Стефана и Фролова и др., составили 0.86 и 0.62 соответственно.
Результаты корреляционного анализа Грея динамики движения оползня (табл. 2), глубины промерзания грунтов на оползневом склоне и климатические параметры влияния (мощность снежного покрова, средняя температура воздуха холодного периода, количество холодных дней, количество жидких осадков) представлены в табл. 3. Коэффициент корреляции Грея рассчитан по уравнению:
,
где x0=(x0(1), x0(2),…, x0(k)) – главная последовательность, xi=(xi(1), xi(2),…, xi(k)) – последовательности сравнения и ρ – коэффициент разложения [0;1], обычно принимаемый 0.5.
Степень корреляции Грея определена по формуле:
.
Последовательности корреляции Грея выглядят следующим образом: смещение оползневых точек > средняя температура воздуха холодного периода > количество жидких осадков > мощность снежного покрова > глубина промерзания горных пород; глубина промерзания > количество холодных дней > средняя температура воздуха холодного периода > количество жидких осадков > мощность снежного покрова. Наиболее определяющее влияние на глубину промерзания грунтов оказывают длительность холодного периода и температура воздуха. При изучении динамики движения оползня по методу Грея важным параметром влияния явилась средняя температура воздуха холодного периода. В то же время в обеих последовательностях существенным влияющим параметром оказалось количество жидких осадков, иными словами, влажность грунтов – еще один необходимый наблюдаемый параметр. Для более точного анализа в аналогичных исследованиях важно было бы измерять уровень грунтовых вод и глубину проникновения поверхностных вод.
В инженерной геологии существует параметр «индекс промерзания поверхности» – искусственная величина, которая зависит от температуры воздуха и продолжительности холодного периода. На рис. 4 рассчитанные глубины промерзания сопоставлены с индексами промерзания поверхности (табл. 4) при среднегодовых температурах. Глубина промерзания в целом имеет тенденцию роста с увеличением индекса промерзания поверхности по обоим расчетным методам, что подтверждает последовательности корреляции Грея, полученные выше. При одинаковых индексах промерзания глубина промерзания более выветрелых пород выше, чем у техногенных грунтов. Это связано с тем, что полускальные породы обладают большей теплопроводностью, чем песчаные отложения с гравием и галькой, и более выветрелые полускальные породы имеют более высокую влагоемкость ввиду трещиноватости, чем техногенные грунты. Минимальные значения индексов промерзания наблюдались в самый теплый год и многоснежный холодный период 2019–2020 гг. Максимальные значения, соответственно, зафиксированы в период, когда среднегодовая температура воздуха была минимальной, однако мощность снежного покрова в холодный период 2013–2014 гг. имела среднее значение 0.45 м и максимальное – 0.65 м.
Заключение
В ходе исследований была осуществлена оценка глубины сезонного промерзания горных пород на оползневом склоне в криолитозоне в долине р. Воркуты. Для этого использовались методы Стефана и предложенный Фроловым и др. Расчетные данные частично сопоставлены с наблюденными и также сопоставлены с индексами промерзания поверхности. Кроме этого, с помощью корреляционного анализа Грея оценено геометрическое сходство между главной (эталонной) последовательностью динамики движения оползня, глубины сезонного промерзания грунтов и последовательностями параметров влияния: количества холодных дней в году, средней температуры воздуха холодного периода, количества жидких осадков и мощности снежного покрова. По результатам комплексного анализа можно сделать следующие выводы:
– в октябре наступают климатические условия для устойчивого сезонного промерзания верхнего слоя горных пород, тогда же начинает формироваться устойчивый снежный покров, средняя температура воздуха октября за период исследований составляет -3 °С;
– основное промерзание грунта и горных пород происходит к концу декабря–января – до 90 % от всего сезонного промерзания;
– глубина промерзания, определенная по методу Стефана, который основан на термических свойствах грунтов и горных пород, температуре поверхности, воздуха и других климатических параметрах, имеет значения несколько выше, чем по методу, основанному на теплопроводности трехслойной среды (снежный покров, мерзлый и талый грунты) и их тепловом балансе (по Фролову и др.), однако динамика промерзания – идентичная;
– среднее расхождение расчетных и наблюденных глубин промерзания техногенного грунта составляет 30 %, коэффициент корреляции наблюденных данных в термометрических скважинах № 1 и 2 и рассчитанных по методу Стефана составляет 0.76 и 0.92 и по методу Фролова и др. – 0.67 и 0.9 соответственно;
– в наибольшей степени определяющее влияние на глубину промерзания отложений оказывают длительность холодного периода и температура воздуха, для динамики движения оползня важное значение имеет средняя температура воздуха холодного периода, что подтверждает выводы сорокалетних наблюдений о том, что на Восточно-Европейской равнине на территории России наблюдаются значимые коэффициенты корреляции температуры почвогрунтов на глубинах от 0.9 до 3.2 м с температурой воздуха, особенно эта взаимосвязь существенна в районах с островной и прерывистой мерзлотой [3];
– режим снегонакопления в данных исследованиях имеет значительную степень корреляции Грея с динамикой движения оползня и глубиной промерзания;
– немаловажным влияющим параметром оказалось и количество жидких осадков (степень корреляции глубины промерзания и смещения оползневых точек – 0.79 и 0.86 соответственно), поэтому важно измерять уровень грунтовых вод на оползневом склоне, глубину проникновения поверхностных вод и осуществлять контроль за поверхностным стоком на ул. Шахтерская набережная;
– с целью стабилизации оползневого склона на ул. Шахтерская набережная в г. Воркуте в связи с результатами проведенных исследований предлагаются рекомендации по установке теплоизоляционных экранов на глубину сезонного промерзания грунтов и обустройству дренажных колодцев в теле оползня и выше по склону в грунтах, в настоящее время не вовлеченных в оползневый процесс.
Авторы заявляют об отсутствии конфликта интересов.
1. Hugh, A. L. Henry. Climate change and soil freezing dynamics : historical trends and projected changes / A. L. Henry Hugh // Climatic Change. – 2008. – Vol. 87. – P. 421–434. – DOIhttps://doi.org/10.1007/s10584-007-9322-8.
2. Tian, S. Cyclic behaviour of coarse-grained materials exposed to freeze-thaw cycles : experimental evidence and evolution model / S. Tian, T. Liang, L. Xianzhang, X. Kong, S. Li [et al.] // Cold Regions Science and Technology. – 2019. – Vol. 167. – 167:102815. DOIhttps://doi.org/10.1016/j.coldregions.2019.102815.
3. Шерстюков, А. Б. Корреляция температуры почвогрунтов с температурой воздуха и высотой снежного покрова на территории России / А. Б. Шерстюков // Криосфера Земли. – 2008. – Т. XII, № 1. – С. 79–87.
4. Stefan, J. Ueber die Theorie der Eisbildung, insbesondere über die Eisbildung im Polarmeere / J. Stefan // Annalen der Physik. – 1890. – Vol. 278. – Iss. 2. – 18:965.
5. Волкова, Н. Г. Климатические характеристики промерзания почвы на территории РФ / Н. Г. Волкова // Вестник МГСУ. – 2011. – № 3. – С. 235–241.
6. Jumikis, A. R. Thermal geotechnics / A. R. Jumikis. – New Brunswick : Rutgers University Press, 1977. – 375 p.
7. Steurer, P. M. Probability distributions used in 100-year return period of air-freezing index / P. M. Steurer // Journal of Cold Regions Engineering. – 1996. – Vol. 10. – Iss. 1. – P. 25–35.
8. Erlingsson, S. Correlating air freezing index and frost penetration depth – a case study for Sweden / S. Erlingsson, D. Saliko // Proceedings of the 9th International Conference on Maintenance and Rehabilitation of Pavements. – Mairepav9. – 2020. – Vol. 76. – P. 847–857.
9. Zhou, J. Practical prediction method on frost heave of soft clay in artificial ground freezing with field experiment / J. Zhou, W. Zhao, Y. Tang // Tunnelling and Underground Space Technology. – 2021. – Vol. 107 (C) : 103647. DOIhttps://doi.org/10.1016/j.tust.2020.103647.
10. Заболотник, С. И. О роли снежного покрова в формировании температурного режима многолетнемерзлых пород / С. И. Заболотник // Наука и техника в Якутии. – 2019. – № 2 (37). – С. 3–7. – DOI:https://doi.org/10.24411/1728-516Х-2019-10030.
11. Осокин, Н. И. Коэффициент теплопроводности снега и его изменчивость / Н. И. Осокин, А. В. Сосновский, Р. А. Чернов // Криосфера Земли. – 2017. – Т. XXI, № 3. – С. 60–68. – DOI:https://doi.org/10.21782/KZ1560-7496-2017-3(60-68).
12. Павлов, А. В. Теплоперенос в природных комплексах Севера России / А. В. Павлов // Криосфера Земли. – 2000. – Т. IV, № 4. – С. 22–31.
13. Кудрявцев, В. А. Температура верхних горизонтов вечномерзлой толщи в пределах СССР / В. А. Кудрявцев. – Москва : Изд-во АН СССР, 1954. – 183 с.
14. Yedek, S. A study on determination of the frost depths for building foundations in the Central Anatolia Region of Turkey / Seda Yedek, Muge Elif Orakoglu Fırat, Murat Bikce // Arabian Journal of Geosciences. – 2021. – 14 (14) : 1941. DOIhttps://doi.org/10.1007/s12517-021-08282-x.
15. Kahimba, F. Modeling soil temperature, frost depth, and soil moisture redistribution in seasonally frozen agricultural soils / Frederick Kahimba, Ramanathan Sri Ranjan, Danny D. Mann // Applied Engineering in Agriculture. – 2009. – Vol. 25. – Iss. 6. – P. 871–882. DOI :https://doi.org/10.13031/2013.29237.
16. Курочкина, В. А. Морозное пучение грунтов как фактор геоэкологического риска в условиях города Москвы / В. А. Курочкина, И. Ю. Яковлева // Геоэкология. Инженерная геология. Гидрогеология. Геокриология. – 2020. – № 1. – С. 86–92. – DOI:https://doi.org/10.31857/S0869780920010081.
17. Ананичева, М. Д. Методы оценки последствий изменения климата для физических и биологических систем / М. Д. Ананичева, Ю. А. Анохин, М. В. Болгов и др. – Москва : Планета, 2012. – 511 c.
18. Васильев, А. А. Деградация мерзлоты: результаты многолетнего геокриологического мониторинга в западном секторе Российской Арктики / А. А. Васильев, А. Г. Гравис, А. А. Губарьков и др. // Криосфера Земли. – 2020. – Т. XXIV, № 2. – С. 15–30. DOI:https://doi.org/10.21782/KZ1560-7496-2020-2(15-30).
19. Вихоть, А. Н. Исследование оползневого процесса в у??????? ?????????????? ??????? (?????????? ????? ?????? ?. ???????) / ?. ?. ?????? // ??????? ????.?? 2024. ? ?. 19, ???. 4. ? ?. 606?617. - DOI:https://doi.org/10.22227/1997-0935.2024.4.606-617.словиях сезонномерзлых грунтов (оползневый склон долины р. Воркута) / А. Н. Вихоть // Вестник МГСУ. – 2024. – Т. 19, вып. 4. – С. 606–617. - DOI:https://doi.org/10.22227/1997-0935.2024.4.606-617.
20. Westerveld, L. Arctic permafrost atlas / L. Westerveld, T. Kurvits, T. Schoolmester [et al.]. - Arendal : GRID–Arendal, 2023. – 175 p.
21. Оберман, Н. Г. Экогеология Республики Коми и восточной части Ненецкого автономного округа / Н. Г. Оберман, И. Г. Шеслер, А. И. Рубцов. – Сыктывкар : ПрологПлюс, 2004. – 256 с.
22. Государственная геологическая карта Российской Федерации. Издание второе. Масштаб: 1 : 200 000, серия: Полярно-Уральская. Q-41-V, VI (Воркута). Объяснительная записка / под ред. А. И. Водолазского. – Сыктывкар: ЗАО «ГГК МИРЕКО», 2005.
23. Андреичева, Л. Н. Плейстоцен Европейского Северо-Востока / Л. Н. Андреичева. – Екатеринбург : УрО РАН, 2002. – 323 с.
24. Андреичева, Л. Н. Природная среда неоплейстоцена и голоцена на Европейском Северо-Востоке России / Л. Н. Андреичева, Т. И. Марченко-Вагапова, М. Н. Буравская, Ю. В. Голубева. – Москва : ГЕОС, 2015. – 224 с.
25. Deng, J. L. Introduction to Grey system theory / J. L. Deng // Journal of Grey System. – 1989. – Vol. 1. – Iss. 1. – P. 1–24.
26. Frolov, D. M. Calculation of ground freezing depth under bare and covered with the snow cover ground surface on the site of the meteorological observatory of Lomonosov Moscow State University for winter seasons of 2011/12–2017/18 / D. M. Frolov // Environmental Dynamics and Global Climate Change. – 2019. – Vol. 10. – Iss. 2. – P. 86–90. DOIhttps://doi.org/10.17816/edgcc21203.
27. Голубев, В. Н. Оценка глубины промерзания грунта на севере ЕТР на основе данных о термических условиях зимних сезонов и режиме снегонакопления / В. Н. Голубев, Г. А. Ржаницын, Д. М. Фролов // Актуальные проблемы геокриологии: Сборник докладов расширенного заседания научного совета по криологии Земли РАН с участием российских и зарубежных ученых, инженеров и специалистов. – 2018. – Т. 2. – С. 193–199.