CC BY
6
УДК 581.54
DOI: 10.24412/1728-323X-2024-3-16-21
ХАРАКТЕРИСТИКИ И РЕАКЦИЯ РАДИАЛЬНОГО ПРИРОСТА ДЕРЕВЬЕВ В ЦЕНТРАЛЬНОЙ ЯКУТИИ НА ИЗМЕНЕНИЕ КЛИМАТА
Сяохун Чжан, ведущий инженер, Северо-Восточный федеральный университет имени М. К. Аммосова, [email protected], г. Якутск, Россия
Аннотация. В работе представлено описание характеристик и исследование реакции радиального прироста деревьев (лиственница — Larix cajanderi Mayr, сосна Pinus sylvestris L.) на изменение климата (температура воздуха, осадки) в районе с. Маган (Центральная Якутия). Установлено, что в районе с. Маган за период 1929—2022 гг. средняя ширина годичных колец сосны 1,7 ± 0,829 мм выше, чем лиственницы 1 ± 0,548 мм (1907—2022 гг.). Показатель интерквартильного размаха IQR = 0,302 сосны указывает на меньшую амплитуду колебания ширины годичных колец в годовом ходе, чем лиственницы IQR = 0,340. Показатели выраженности популяционного сигнала (EPS) и соотношения сигнал/шум (SNR) у сосны также выше, чем лиственницы. Однако оба вида имеют одинаковый коэффициент чувствительности MS = 0,366, что может свидетельствовать о схожей реакции деревьев к изменениям окружающей среды. За период 1929—2022 гг. выявлены высокие коэффициенты корреляции между осадками и шириной годичных колец в июле-августе предыдущего года (для лиственницы r = 0,472, для сосны r = 0,468 при p < 0,01). Данный факт указывает на влияние осадков в эти периоды на последующий рост деревьев. В январе—марте и июне—ноябре текущего года (2023 г.) в целом наблюдаются отрицательные коэффициенты корреляции между температурой воздуха и шириной годичных колец. Положительные значения коэффициентов корреляции приходятся на декабрь предыдущего года и весенний период текущего (апрель и май) — начало вегетационного периода c максимальными значениями свыше 0,2 для лиственницы в мае при p < 0,05.
Abstract. The paper presents a description of the characteristics and study of the response of radial growth of trees (larch — Larix cajanderi Mayr, pine Pinus sylvestris L.) to climate change (air temperature, precipitation) in the area of the village of Magan (Central Yakutia). It has been established that in the area of the village of Magan for the period of 1929—2022 the average width of pine tree rings is 1,7 ± 0,829 mm higher than that of larch 1 ± 0,548 mm (1907—2022). The interquartile range indicator IQR = 0,302 for the pine indicates a smaller amplitude of fluctuation in the width of tree rings in the annual cycle than for the larch IQR = 0,340. The indicators of population signal intensity (EPS) and signal-to-noise ratio (SNR) in the pine are also higher than in the larch. However, both species have the same sensitivity coefficient MS = 0,366, which may indicate a similar response of trees to environmental changes. For the period of 1929—2022 high correlation coefficients (r) were revealed between precipitation and the width of annual rings in July—August of the previous year (for the larch r = 0,472 and pine r = 0,468 at p < 0,01). This fact indicates the influence of precipitation during these periods on the subsequent growth of trees. In January—March and June—November of the current year (2023), in general, negative correlation coefficients between air temperature and the width of tree rings are observed. Positive values of correlation coefficients occur in December of the previous year and the spring period of the current year (April and May) — the beginning of the growing season with maximum values over 0,2 for larch in May at p < 0,05.
Ключевые слова: радиальный прирост деревьев, климат, температура воздуха, осадки, Центральная Якутия.
Keywords: radial tree growth, climate, air temperature, precipitation, Central Yakutia.
Введение
Изменение климата оказывает значительное влияние на лесные экосистемы [1—3]. Республика Саха (Якутия) уникальна по географическому положению и климатическим условиям, что способствует исследованиям выявления связи между радиальным приростом деревьев и абиотическими факторами неживой природы. В работах показано [4—6], что ширина годичных колец деревьев является одним из надежных индикаторов изменения климата, влияющим на рост деревьев. Результаты анализа связи между шириной годичных колец деревьев и климатическими параметрами вносят вклад в понимание механизмов адаптации лесных экосистем Центральной Якутии на изменение климата.
Цель работы — описание характеристик и исследование реакции радиального прироста деревьев (лиственница — Larix cajanderi Mayr, сосна Pinus sylvestris L.) на изменение климатических факторов (температура воздуха, осадки) в районе с. Маган (Центральная Якутия).
Материалы и методы исследования
Исследование проведено в районе с. Маган (62°16'00.37'' с. ш., 129°34'19.73'' в. д.; рис. 1). Образцы древесных колец лиственницы и сосны собраны в 2023 г. Отбор проб проводился на открытых лесных участках, не подверженных лесным пожарам, с неглубоким сезонно-талым слоем почвы.
С помощью возрастного бура на высоте 1,3 м извлечены по два керна с разных сторон окружности ствола деревьев. Все образцы были зафиксированы, высушены на воздухе, пронумерованы и отшлифованы. Ширина годичных колец измерялась с использованием Ы№ар6.0 (разрешение 0,001 мм) и проверялась программным обеспечением TSAP-Win [7]. Качество перекрестной датировки контролировалось программным обеспечением СОБЕСНА, которое использует корреляционный анализ для проверки надежности датировки [8]. Для удаления нежелательных возрастных трендов использовалась стандартизация данных с помощью отрицательной экспоненци-
Рис. 1. Местоположение исследуемого района
альной функции [9]. После удаления тренда были получены стандартные и остаточные хронологии. Статистический анализ проводился с использованием пакета ёр1Я в программном обеспечении Я (версия 4.3.3) [10].
Среднемесячные и среднегодовые значения температуры воздуха и осадков получены из Всероссийского научно-исследовательского института гидрометеорологической информации — Мировой центр данных (http://meteo.ru).
Результаты и обсуждение
Республика Саха (Якутия) характеризуется экстремальными сезонными различиями, территория которой охватывает три климатических пояса — арктический, субарктический и умеренный.
В зимний период (декабрь—февраль) 1929— 2022 гг. по данным метеостанции «Якутск» (рис. 2, а) наблюдались экстремально низкие температуры воздуха от —35 до —41 °С с минимумом в январе —51,4 °С и количеством осадков от 7 до 9 мм, что соответствует сухим и холодным зимним условиям региона. Летом (июнь—август) температура воздуха за указанный период составляла от 16 °С до 19 °С с максимумом осадков в июле до 41 мм.
За период 1929—2022 гг. (рис. 2, в) наблюдается положительная тенденция роста температуры воздуха в Центральной Якутии. Результат
о
£
н
50 25 0
- 50
25
-25
- -50
э
о
о
1 2 3 4 5 6 7 8 9 10 11 12 Месяцы
а)
О
пТ й Й а е
к
а
е
н
150
100
50
0 20
0
-20
-40
1 2 3 4 5 6 7 8 9 10 11 12 Месяцы б)
О
£
к
а
е
н
1929 1944 1959 1974 1989 2004 2019 Годы
300 -
250
200
150
1929 1944 1959 1974 1989 2004 2019 Годы г)
0
Рис. 2. (а) Среднемесячная температура воздуха (столбцы) и осадки (сплошная линия с круглым маркером); (б) среднемесячный размах данных по температуре воздуха и осадкам; (в) среднегодовая температура воздуха; (г) годовое количество осадков. Штрихпунктирные линии на (в) и (г) — линейные тренды, серая область представляет собой 95 % доверительный интервал
s
1907 1922 1937 1952 1967 1982 1997 2012 Годы
Рис. 3. Хронология ширины годичных колец деревьев в районе с. Маган
линейного тренда указывает на скорость изменения температуры воздуха — 0,04 °С/год с коэффициентом детерминации Я2 = 0,47 (78 % влияния прочих факторов на изменчивость температуры воздуха). За период 1989—2022 гг. скорость изменения температуры воздуха повысилась — 0,07 °С/год. В динамике осадков (рис. 2, г) за период 1929—2022 гг. наблюдается положительный тренд с минимальной скоростью изменения — 0,14 мм/год.
Хронология (таблица) лиственниц включает 115 лет, а сосны — 94 г. Средняя ширина годичных колец лиственницы за период 1907—1922 гг. составила 1 ± 0,548 (рис. 3), а медиана 0,973 близка к среднему значению, что свидетельствует о равномерном распределении без значительных отклонений. Однако высокая дисперсия 0,07 и стандартное отклонение 0,3 указывают на значительную вариабельность ширины годичных колец. Такая вариабельность может отражать чувствительность лиственницы к климатическим условиям (годовые колебания температуры воздуха, осадков и т. п.). Интерквартильный размах = 0,340 подтверждает амплитуду колебаний ширины годичных колец лиственницы в годовом ходе.
Средняя ширина годичных колец для сосны 1,7 ± 0,829 мм, а медиана 1,018 мм. Низкая дисперсия 0,05 и стандартное отклонение 0,2 могут указывать на стабильный рост сосны. = 0,302 показывает на меньшую амплитуду колебания ширины годичных колец в годовом ходе, чем лиственницы.
Коэффициент автокорреляции 1-го порядка (сосна — ас1 = 0,684) указывает на более значимую непрерывность роста, чем у лиственницы (ас1 = 0,626). Межсериальный коэффициент корреляции для сосны (ЯЬаг = 0,659) выше, чем у лиственницы (ЯЬаг = 0,467), что указывает на более высокую согласованность изменений ши-
рины годичных колец. Оба вида имеют одинаковый коэффициент чувствительности MS = 0,366. Данный факт может свидетельствовать о схожей чувствительности деревьев к изменениям окружающей среды. Однако показатели выраженности популяционного сигнала (EPS) и соотношения сигнал/шум (SNR) у сосны выше, чем лиственницы.
Проведен корреляционный анализ Пирсона (r) между ш ириной г одичных колец деревьев и абиотическими факторами неживой природы (температура воздуха, осадки) за период 1929—2022 гг. — рис. 4.
В целом (рис. 4, а), за период 1929—2022 гг. наблюдается положительная корреляция свыше 0,2 (p < 0,05) между осадками и шириной годичных колец в июле и августе. Высокие коэффициенты корреляции приходятся на июль-август предыдущего года (P78 — для лиственницы r = 0,472 и сосны r = 0,468 при p < 0,01). Данный факт указывает на влияние осадков в эти периоды на последующий рост деревьев.
В январе—марте и июне—ноябре текущего года в целом наблюдаются отрицательные коэффициенты корреляции между температурой воздуха и шириной годичных колец. Низкие коэффициенты корреляции приходятся на июль-август предыдущего года (P78 — для лиственницы r = -0,335 и сосны r = -0,308 при p < 0,01). Положительные значения коэффициентов корреляции приходятся на декабрь предыдущего года
Таблица
Характеристики древостоя и статистические параметры хронологий ширины годичных колец
Характеристики древо- Лиственница Сосна
стоя и статистические 1907-2022 гг. 1929-2022 гг.
параметры (115 лет) (94 года)
Количество индивиду- 32 29
альных серий, шт.
Средняя ширина годич- 1 ± 0,548 1,7 ± 0,829
ных колец ± стандартное
отклонение, мм
Коэффициент автокор- 0,626 0,684
реляции 1-го порядка
(ас1)
Межсериальный коэф- 0,467 0,659
фициент корреляции
(Rbar)
Коэффициент чувстви- 0,366 0,366
тельности (MS)
Выраженный сигнал 0,971 0,984
популяции (EPS)
Отношение сигнала 32,951 60,438
к шуму (SNR)
РнРнРнРнРн^г^^ООООООООО ¡чгч'Г' Рн Рн Рч О О Р4
Месяцы
Сосна
Лиственница
С^ОО ^о^н гч^н ГЧ РнРнРн ^ ^ ^ ОО Рч Рн Рн
СП ^ 1Л щ
о о о о
а)
Месяцы б)
СС
Рис. 4. Коэффициенты корреляции Пирсона между осадками (а) и температурой воздуха (б) с шириной годичных колец, где P5-P12 — май—декабрь предыдущего года, Ы-11 — январь—ноябрь текущего года (2023 г.), P78 — июль-август предыдущего года. Штрихпунктирной линией указан значимый предел при p < 0,05, пунктирной линией — p < 0,01
и весенний период текущего (апрель и май) — начало вегетационного периода с максимальными значениями свыше 0,2 для лиственницы в м ае при р < 0,05. Стоит отметить, что в дальнейшем среднегодовой рост температуры воздуха [11] (рис. 2, в) может привести к увеличению испарения и уменьшению запасов почвенной влаги с негативным влиянием на приросте деревьев. Необходим дальнейший мониторинг за динамикой температуры воздуха для выявления негативных эффектов ее повышения.
В зимний период происходит снегонакопление, весной — пополнение, а летом — дефицит почвенной влаги. В целом влажность почвы является важным и лимитирующим ресурсом для лесных экосистем [12, 13]. Возможно, подобный эффект проявляется в районе с. Маган между осадками и температурой воздуха с шириной годичных колец (рис. 4) — в данном случае лимитирующим фактором может являться температура воздуха и ее положительная тенденция роста (рис. 2, в), которая в целом ослабляет влияние осадков на радиальный прирост деревьев. Видно (рис. 2, а), что в Центральной Якутии лето (июнь—август) достаточно засушливое, а пополнение почвенной влаги, скорее всего, происходит за счет инфильтрации снежного покрова весной и иных процессов.
Известно, что территория Якутии находится в зоне распространения многолетнемерзлых грунтов, в период начала (осень) промерзания которых происходит процесс миграции почвенной влаги вверх к дневной поверхности [13]. В результате аккумулируется влага в верхних и ниж-
них частях почвенно-грунтового разреза. Почвенная влага, пребывая длительное время года в твердой фазе, оттаивая в летний период (переход в жидкое состояние), становится доступной для лесной экосистемы (растительности). В настоящее время наблюдения за влажностью и температурой почвогрунтов в Центральной Якутии дали возможность оценить степень влияния гидротермического режима почв на рост и развитие древесных пород [3].
В работе [14, 15] показано, что при раннем сходе снежного покрова в апреле происходит прогревание деятельного слоя. Далее летние температуры воздуха не лимитируют радиальный прирост лиственницы, а значимая корреляционная связь с древесно-кольцевыми хронологиями в Центральной Якутии наблюдается в мае, что подтверждают исследования данной работы. Однако не исключен лимитирующий фактор летних температур воздуха, влияющий на осадки, — иссушающий эффект, что негативно отражается на радиальном приросте деревьев. Автором [12] так же показано, что для сосны значимое влияние оказывают осадки предыдущего года, особенно в летние месяцы (например, рис. 4, а). Осенние осадки предыдущего года, очевидно, накапливаясь в деятельном слое в виде почвенной влаги, используются деревьями в течение следующего вегетационного периода.
Отклик радиального прироста деревьев на изменение климата зависит от локальных и региональных геофизических условий. В дальнейшем планируется проведение дополнительных исследований в указанном районе с привлечением
иных абиотических факторов неживой природы и учетом гидротермических условий многолетне-мерзлых грунтов для выявления причинно-следственных связей.
Заключение
Таким образом, в работе показано, что в зимний период (декабрь—февраль) 1929—2022 гг. по данным метеостанции «Якутск» наблюдались экстремально низкие температуры воздуха от -35 до -41 °C с минимумом в январе -51,4 °C и количеством осадков от 7 до 9 мм, что соответствует сухим и холодным зимним условиям Якутии. Летом (июнь—август) температура воздуха за указанный период составляла от 16 до 19 °C с максимумом осадков в июле до 41 мм. Скорость изменения температуры воздуха и осадков составила 0,04 °С/год и 0,14 мм/год, соответственно.
Установлено, что в районе с. Маган за период 1929—2022 гг. средняя ширина годичных колец сосны 1,7 ± 0,829 выше, чем лиственницы 1 ± 0,548 (1907—2022 гг.). Показатель интерквар-тильного размаха IQR = 0,302 сосны указывает на меньшую амплитуду колебания ширины годичных колец в годовом ходе, чем лиственницы IQR = 0,340. Показатели выраженности попу-ляционного сигнала (EPS) и соотношения сиг-
нал/шум (ЗКЯ) у сосны также выше, чем лиственницы. Однако оба вида деревьев имеют одинаковый коэффициент чувствительности МЗ = 0,366, что может свидетельствовать о схожей реакции деревьев к изменениям окружающей среды.
За период 1929—2022 гг. выявлены высокие коэффициенты корреляции между осадками и шириной годичных колец в июле-августе предыдущего года (Р78 — для лиственницы г = 0,472 и сосны г = 0,468 при р < 0,01). Данный факт указывает на влияние осадков в эти периоды на последующий рост деревьев.
В январе—марте и июне—ноябре текущего года (2023 г.) в целом наблюдаются отрицательные коэффициенты корреляции между температурой воздуха и шириной годичных колец. Положительные значения коэффициентов корреляции приходятся на декабрь предыдущего года и весенний период текущего (апрель и май) — начало вегетационного периода с максимальными значениями свыше 0,2 для лиственницы в мае при р < 0,05.
Финансирование. Исследование выполнено в рамках государственного задания СВФУ № ЕБЯО-2023-0027.
Библиографический список
1. Lindner M. Climate change impacts, adaptive capacity, and vulnerability of European forest ecosystems / M. Lindner [et al.] // Forest ecology and management. — 2010. — V. 259. — № 4. — P. 698—709.
2. Николаев А. Н. Дендрохронологический анализ природных процессов в криолитозоне (на примере Центральной Якутии): специальность 03.02.08 «Экология (по отраслям)»: дис. ... д-ра биол. наук / Николаев Анатолий Николаевич. — Москва, 2011. — 336 с.
3. Nikolaev A. N., Fedorov P. P., Desyatkin A. R. Influence of climate and soil hydrothermal regime on radial growth of Larix cajanderi and Pinus sylvestris in central Yakutia, Russia // Scandinavian Journal of Forest Research. — 2009. — V. 24. — № 3. — P. 217—226.
4. Garcia-Valdes R. Climate change impacts on long-term forest productivity might be driven by species turnover rather than by changes in tree growth / R. Garcia-Valdes [et al.] // Global Ecology and Biogeography. — 2020. — V. 29. — № 8. — P. 1360—1372.
5. Nitschke C. R. The influence of climate and drought on urban tree growth in southeast Australia and the implications for future growth under climate change / C. R. Nitschke [et al.] // Landscape and Urban Planning. — 2017. — V. 167. — P. 275—287.
6. Wilmking M. Global assessment of relationships between climate and tree growth / M. Wilmking [et al.] // Global Change Biology. — 2020. — V. 26. — № 6. — P. 3212—3220.
7. Rinn F. TSAP—Win: Time series analysis and presentation for dendrochronology and 409 related applications // User reference. Heidelberg, 2003. — URL: https://cir.nii.ac.jp/crid/1572543024876591616 (дата обращения 24.07.2024).
8. Grissino-Mayer H. D. Evaluating crossdating accuracy: a manual and tutorial for the computer program COFECHA // Tree-ring Research. — 2001. — V. 57. — № 2. — P. 205—221.
9. Cook E. R. A time series analysis approach to tree ring standardization // Dissertation. University of Arizona. Tucson. — 1985. — URL: https://www.st-andrews.ac.uk/~rjsw/PalaeoPDFs/Cook1985-Chapter%202.pdf (дата обращения 24.07.2024).
10. Bunn A. G. Statistical and visual crossdating in R using the dplR library // Dendrochronologia. — 2010. — V. 28. — № 4. — P. 251—258.
11. De Pryck K. Intergovernmental expert consensus in the making: the case of the summary for policy makers of the IPCC 2014 Synthesis Report // Global Environmental Politics. — 2021. — V. 21. — № 1. — P. 108—129.
12. Николаев А. Н. Влияние влажности почвы на радиальный прирост лиственницы и сосны на территории Центральной Якутии // Вестник Северо-Восточного федерального университета им. М. К. Аммосова. — 2017. — № 6 (62). — С. 25—37.
13. Федоров П. П., Десяткин А. Р. Связь температурного режима мерзлотных почв и радиального прироста лиственницы в Центральной Якутии // Успехи современного естествознания. — 2016. — № 7. — С. 185—189.
14. Угаров И. С., Ефремов П. В. Влагозапасы деятельного слоя почвы бассейна реки Лена // Успехи современного естествознания. — 2022. — № 10. — С. 88—92.
15. Николаев А. Н. Влияние гидродинамического режима мерзлотных почв на радиальный прирост лиственницы и сосны в Центральной Якутии / А. Н. Николаев, П. П. Федоров, А. Р. Десяткин // Сибирский экологический журнал. — 2011. — Т. 18. — № 2. — С. 189—201.
CHARACTERISTICS AND RESPONSE OF RADIAL GROWTH OF TREES IN CENTRAL YAKUTIA TO CLIMATE CHANGE
Zhang Xiaohong, Lead Engineer, M. K. Ammosov North-Eastern Federal University, [email protected], Yakutsk, Russia References
1. Lindner M. Climate change impacts, adaptive capacity, and vulnerability of European forest ecosystems / M. Lindner [et al.]. Forest ecology and management. 2010. Vol. 259. No. 4. P. 698—709.
2. Nikolaev A. N. Dendrohronologicheskij analiz prirodnyh processov v kriolitozone (na primere Central'noj Jakutii): spe-cial'nost' 03. 02. 08 "Jekologija (po otrasljam)": dissertacija na soiskanie uchenoj stepeni doktora biologicheskih nauk [Den-drochronological analysis of natural processes in the cryolithozone (a case study of Central Yakutia): specialty 03.02.08 "Ecology (by branches)": dissertation for the degree of Doctor of Biological Sciences] / Nikolaev Anatolij Nikolaevich. Moscow. 2011. 336 p. [in Russian].
3. Nikolaev A. N., Fedorov P. P., Desyatkin A. R. Influence of climate and soil hydrothermal regime on radial growth of Larix cajanderi and Pinus sylvestris in Central Yakutia, Russia. Scandinavian Journal of Forest Research. 2009. Vol. 24. No. 3. P. 217—226.
4. Garcia-Valdes R. Climate change impacts on long-term forest productivity might be driven by species turnover rather than by changes in tree growth / R. Garcia-Valdes [et al.]. Global Ecology and Biogeography. 2020. Vol. 29. No. 8. P. 1360—1372.
5. Nitschke C. R. The influence of climate and drought on urban tree growth in Southeast Australia and the implications for future growth under climate change / C. R. Nitschke [et al.]. Landscape and Urban Planning. 2017. Vol. 167. P. 275—287.
6. Wilmking M. Global assessment of relationships between climate and tree growth / M. Wilmking [et al.]. Global Change Biology. 2020. Vol. 26. No. 6. P. 3212—3220.
7. Rinn F. TSAP—Win: Time series analysis and presentation for dendrochronology and 409 related applications. User reference. Heidelberg. 2003. Available at: https://cir.nii.ac.jp/crid/1572543024876591616.
8. Grissino—Mayer H. D. Evaluating cross-dating accuracy: a manual and tutorial for the computer program COFECHA. Tree—ring Research. 2001. Vol. 57. No. 2. P. 205—221.
9. Cook E. R. A time series analysis approach to tree ring standardization. Dissertation. University of Arizona. Tucson. 1985. Available at: https://www.st-andrews.ac.uk/~rjsw/PalaeoPDFs/Cook1985-Chapter%202.pdf.
10. Bunn A. G. Statistical and visual cross-dating in R using the dplR library. Dendrochronologia. 2010. Vol. 28. No. 4. P. 251—258.
11. De Pryck K. Intergovernmental expert consensus in the making: the case of the summary for policy makers of the IPCC 2014 Synthesis Report. Global Environmental Politics. 2021. Vol. 21. No. 1. P. 108—129.
12. Nikolaev A. N. Vlijanie vlazhnosti pochvy na radial'nyj prirost listvennicy i sosny na territorii Central'noj Jakutii [The influence of soil moisture on radial increment of larch and pine in Central Yakutia]. Vestnik Severo—Vostochnogofederal'nogo universiteta im. M. K. Ammosova. 2017. No. 6 (62). P. 25—37 [in Russian].
13. Fedorov P. P., Desyatkin A. R. Svjaz' temperaturnogo rezhima merzlotnyh pochv i radial'nogo prirosta listvennicy v Central'noj Jakutii [Dependence of larch radial increment on permafrost soil temperature regime in the Central Yakutia]. Uspehi sovremennogo estestvoznanija. 2016. No. 7. P. 185—189 [in Russian].
14. Ugarov I. S., Efremov P. V. Vlagozapasy dejatel'nogo sloja pochvy bassejna reki Lena [Active-layer soil moisture in the Lena River basin]. Uspehi sovremennogo estestvoznanija. 2022. No. 10. P. 88—92 [in Russian].
15. Nikolaev A. N. Vlijanie gidrodinamicheskogo rezhima merzlotnyh pochv na radial'nyj prirost listvennicy i sosny v Central'noj Jakutii [Effect of hydrodynamic conditions of permafrost soil on radial growth of larch and pine in Central Yakutia] / Nikolaev A. N., Fedorov P. P., Desyatkin A. R. Sibirskij jekologicheskij zhurnal. 2011. Vol. 18. No. 2. P. 189—201 [in Russian].
