CC BY
65ЛЕСОВЕДЕНИЕ И ФИТОЦЕНОЛОГИЯ
УДК 630*561.24:58.056
С. П. Арефьев
ДЕНДРОХРОНОЛОГИЧЕСКИЙ АНАЛИЗ ФАКТОРНЫХ ПОЛЕЙ СОСТОЯНИЯ КЕДРА {РМиЭ Э1В1ШСА йи ТОиР)
В ЗАПАДНОЙ СИБИРИ
На основе анализа сходства колебаний радиального прироста кедра на обширной территории Западно-Сибирской равнины определены параметры обусловливающих их экологических (прежде всего климатических) факторных полей: напряженность и протяженность, проницаемость географической среды, их внутри-формационная модуляция, а также базисы устойчивости кедровой формации к полевым факторам. Рассмотрены особенности факторных полей для погодичных и многолетних колебаний прироста, широтно-зональные и долготные изменения их параметров.
Древесные кольца, колебания прироста, Западная Сибирь, зональность, лесные формации, кедр сибирский, устойчивость.
Известно, что ширина годичных колец (радиальный прирост древесины) является отражением равнодействующего вектора факторов внешней и внутренней среды, определяющих состояние дерева — его биопродуктивность, репродукцию, ценотический статус, устойчивость. Во временных рядах прироста запечатлена разнообразная информация об этих факторах [Ловелиус, 1979; Горчаковский, Шиятов, 1985, Грейбилл, 1990; и др.]. Уровень сходства колебаний прироста характеризует степень сходства важнейших факторов среды, силу их воздействия на состояние отдельных деревьев и леса, их пространственно-географические масштабы. Зависимость сходства от расстояния указывает на полевую природу экологических факторов, связанную с географической изменчивостью природных циклов [Константинов, Горожанкина, 1995].
Одни из первых широкомасштабных исследований географических закономерностей цикличности древесного прироста проведены С. Г. Шиятовым [1986] на Урале. Им предпринят сравнительный анализ хронологий лиственницы, сосны, ели, произрастающих по верхней границе леса вдоль всего меридионального протяжения Уральского хребта. Для анализа синхронности погодичных колебаний прироста им использован коэффициент Хубера, для сопоставления многолетних колебаний — пирсоновский коэффициент корреляции. Оба показателя с расстоянием закономерно уменьшались. Сходство рядов из ближайших пунктов увеличивалось с продвижением к северу, что связано с усилением лимитирующего климатического фактора.
Аналогичные закономерности были получены при анализе древеснокольцевых хронологий вдоль северной границы распространения древесной растительности в Урало-Сибирской Субарктике [Ваганов и др., 1996]. В работе Е. А. Ваганова с соавт. впервые предпринято широкомасштабное картирование климатических параметров, прежде всего летней температуры воздуха, рассчитанных при анализе древесно-кольцевых хронологий за период с 1611 г. По сути, авторами воссозданы ежегодные полевые картины температурного
фактора гипоарктической зоны Евразии от Урала до восточной Якутии. Вместе с тем дендроклиматические поля отдельных лет обычно довольно мозаичные и в их последовательности еще трудно разглядеть обобщающие закономерности, особенно в широтно-зональном аспекте, который этот крупнейший дендрохронологический проект практически не охватывает.
Рис. 1. Пункты исследований на территории Западно-Сибирской равнины
На большей части территории Западной Сибири подобных исследований не проводилось, хотя, в силу своей большой широтной и долготной протяженности, уникальной градиентности природных условий практически с ненарушенной широтной зональностью, она является почти идеальной моделью для анализа полевых закономерностей состояния леса [Растительный покров..., 1985]. Настоящая работа во многом базируется на материалах широкомасштабной оценки устойчивости кедровых лесов Западной Сибири, проведенной автором [Арефьев, 1997] при интерпретации погодичных колебаний ширины годичных колец кедра и их пространственно-временной динамики в
ходе XX столетия. Следует заметить, что, хотя кедр и является одним из ведущих и наиболее ценных лесообразователей Сибири, в дендрохронологиче-ских работах он используется сравнительно редко [Шиятов, 1973; Полюшкин и др., 1977; Малоквасов, 1986, Ваганов и др., 1996]. Цель настоящей работы — выявить параметры экологических факторных полей, определяющих пого-дичные и многолетние колебания радиального прироста кедра на территории Западно-Сибирской равнины, дать их оценку в связи с устойчивостью кедровой формации.
В 1990-х гг. автором настоящей статьи были взяты дендрохронологиче-ские образцы кедра из 30 пунктов, равномерно охватывающих его западносибирский ареал на протяжении 900 км в широтном и 1025 км в долготном направлении (рис. 1); представлены различные типы кедровников. При отборе и обработке образцов руководствовались общепринятыми методами [Шиятов, 1973; Methods..., 1990; Методы..., 2000], в сомкнутых древостоях брали доминирующие деревья, близкие по ценотическому статусу [Комин, 1970, 1994]. Всего составлено 52 репрезентативных обобщенных хронологий (обычно по 15 керновым образцам).
Для характеристики многолетних «низкочастотных» колебаний прироста рассчитаны индексы относительно 31-летней нормы /(31), для характеристики погодичных «высокочастотных» колебаний — относительно 3-летней нормы /(3). И хотя такой простейший способ индексирования прироста вносит некоторые искажения действительной величины прироста и конфигурации циклов, поскольку при этом не затрагиваются качественные соотношения ширины колец, мы считаем его применение целесообразным для выявления наиболее общих принципиальных закономерностей. Для оценки сходства индексированных хронологий использован параметрический пирсоновский коэффициент корреляции [Комин и др., 1973; Шиятов, 1986], для получения сравнимых значений расчет велся за 100-летний период с 1890 по 1991 г. Параллельно рассчитывали уравнения регрессии индексов прироста по годам.
Регрессионный анализ показал, что рассматриваемая полевая зависимость наиболее точно описывается ветвью параболы второго порядка, ордината основания которой близка к нулю — отсутствию сходства. Парабола обычно имеет малую кривизну, и на рассматриваемом отрезке зависимость близка к линейной, что подтверждается достоверными на высоком уровне (P = 0,001) коэффициентами линейной корреляции (табл.; рис. 2). Это позволяет описать зависимость уравнением линейной регрессии:
r(l) = r0 + bl,
где r(l) — расчетный коэффициент корреляции между хронологиями в зависимости от расстояния между пунктами l; r0 — исходный коэффициент корреляции между хронологиями при расстоянии между пунктами l = 0; b — коэффициент регрессии (имеет отрицательное значение). По аналогии с параметрами электростатического поля r0 можно интерпретировать как исходную напряженность парциального факторного поля, r(l) — как его напряженность на расстоянии l от исходного пункта, b — как проницаемость для факторного поля внешней географической среды, которая в свою очередь связана с насыщенностью этой среды различными экологическими факторами.
Параметры факторных полей радиального прироста кедра на территории Западно-Сибирской равнины
Части ареала Корреляция сходства прироста с расстоянием между пунктами и ее ошибка, ±3(1) Достоверность г Линейная регрессия г() = Г0 +100Ы Довери- тельная зона регрессии ^//Кцт Предельные параметры доверительной зоны регрессии
¿0.001 гтіп гтах /тіп /та^
Многолетние колебания ¡(31)
Широтные части:
Гипоарктические -0,45±0,12 3,71 3,50 0,47-0,027-/ 0,45 0,02 0,92 74 3407
леса
Северная тайга -0,46±0,10 4,48 3,46 0,35-0,025-/ 0,39 -0,04 0,74 0 2872
Средняя тайга -0,36±0,09 4,23 3,37 0,37-0,032/ 0,41 -0,04 0,78 0 2468
Южная тайга -0,54±0,10 5,15 3,46 0,47-0,116/ 0,45 0,02 0,92 17 791
Долготные части:
Западная -0,39±0,08 4,56 3,37 0,28-0,030/ 0,44 -0,16 0,72 0 2384
Центральная -0,48±0,05 8,69 3,29 0,38-0,045/ 0,48 -0,10 0,86 0 1906
Восточная -0,66±0,09 7,65 3,46 0,54-0,063/ 0,33 0,21 0,87 332 1375
Погодичные колебания ¡(3)
Широтные части : Гипоарктические -0,62 ±0,11 5,80 3,65 0,44-0,039/ 0,41 0,03 0,85 76 2168
леса
Северная тайга -0,79±0,07 11,36 3,46 0,49-0,053/ 0,32 0,17 0,81 321 1529
Средняя тайга -0,75±0,06 12,57 3,37 0,64-0,066/ 0,29 0,35 0,93 530 1408
Южная тайга -0,68±0,09 7,51 3,46 0,70-0,099/ 0,25 0,45 0,95 454 958
Долготные части: Западная -0,61 ±0,07 9,27 3,37 0,53-0,053/ 0,38 0,15 0,91 284 1726
Центральная -0,83±0,03 23,74 3,29 0,66-0,088/ 0,34 0,32 1,00 363 1133
Восточная -0,60±0,09 6,64 3,46 0,53-0,063/ 0,38 0,15 0,91 237 1437
Рис. 2. Параметры доверительной зоны регрессии сходства погодичных колебаний прироста кедра в южной тайге (см. табл.)
В исследованиях С. Г. Шиятова [1986] отмечено, что минимальным значениям сходства (в частности, синхронности) в пункте соответствуют сочетания контрастных пар древесно-кольцевых хронологий (например, из суходольного и заболоченного местопроизрастаний), максимальным — сочетания рядов из типологически близких местопроизрастаний. О неодинаковом влия-
нии климата на прирост кедра в разных типах леса пишет Ф. З. Глебов с со-авт. [1986]; наши данные показывают эту же закономерность. То есть, изменчивость параметров сходства анализируемых кедровых хронологий в рамках доверительной зоны уравнения г(1) ± э(г1[Лакин, 1990] определяется действием не общих для всей кедровой формации в данном пункте (при А1 ^ 0) факторов (прежде всего климатических), а диапазона типологических, популяционно-генетических, онтоценогенетических, антропогенных факторов, характерных для данной территории. С некоторой условностью эти факторы, имеющие мозаичное географическое распределение, можно отнести к внутренним факторам состояния кедровой формации, а климатические, имеющие ярко выраженную непрерывную полевую природу,— к внешним для формации факторам. Отклонение от максимального сходства хронологий на фоне общего практически одинакового климатического поля может быть интерпретировано как внутриформационная модуляция климатического факторного поля. По сути, его диапазон дает оценку базиса устойчивости кедровой формации в данных климатических условиях. При минимальном диапазоне модуляции климатические факторы могут рассматриваться как лимитирующие, базис устойчивости при этом сужается до минимального типологического и иного разнообразия кедровников (обычно на границах ареала).
По логике рассуждений, значения максимума доверительной зоны
/(Отах = Г(1) + в(г1)(
характеризуют напряженность внешнего климатического факторного поля, не модулированного внутриформационными факторами, и соответствуют кедровникам одного типа, близким по всем внутриформационным факторам развития. Значения минимума доверительной зоны
Г(1)тт = Г(1) - в(г1)(
характеризуют его напряженность при наибольшей модуляции внутриформа-ционными факторами и соответствуют контрастным типам леса. Наиболее показательна исходная напряженность поля, определяемая в точке пересечения зоны регрессии с осью ординат при I = 0, соответственно гтах и гтп (рис. 2). Соответственно по пересечению максимальной линии доверительной зоны с осью абсцисс г(1)тах = 0 определяется максимальная активная протяженность климатического поля исходной географической точки 1тах, обусловленная проницаемостью внешней среды и наличием однотипных кедровников. Аналогично по пересечению с минимальной линией доверительной зоны г(1)тт = 0 определяется его тотальная протяженность 1т|п (охватывающая все разнообразие кедровой формации), ограниченная внутриформационной модуляцией климатического поля.
Анализ расчетных параметров показывает (табл.), что для факторного поля погодичных (в том числе стрессовых) колебаний прироста по сравнению с полем многолетних колебаний исходная напряженность г0тах обычно выше, но проницаемость географической среды, напротив, ниже. Соответственно больше тотальная протяженность поля погодичных колебаний, но меньше его активная максимальная протяженность, а также и диапазон внутриформаци-онной модуляции. Здесь можно провести параллель между высокочастотными и низкочастотными электромагнитными колебаниями: первые имеют меньшую проникающую способность, но большую разрушительную силу. Ранее [Арефьев, 1997] нами была показана наибольшая опасность для кедра
резких погодичных колебаний, нежели многолетних колебаний той же амплитуды, дающих время для настройки адаптационных механизмов кедровой формации.
На широтно-зональном градиенте от предтундровых лесов до южной тайги исходная напряженность факторного поля погодичных колебаний возрастает (от 0,81-0,85 до 0,95), но проницаемость географической среды для него уменьшается (от -0,04 до -0,10 на 100 км). Соответственно увеличивается его тотальная протяженность (от 76 до 450 км), уменьшается максимальная активная протяженность (с 2170 до 960 км) и диапазон внутриформационной модуляции — базис устойчивости формации (с 0,41 до 0,25).
Для факторного поля многолетних колебаний прироста проницаемость географической среды и его максимальная активная протяженность в широтных природных хонах с севера на юг изменяется аналогично (соответственно с -0,025 до -0,116 на 100 км и от 3410 до 790 км). Но высокая исходная напряженность поля (до 0,92) отмечается в обеих крайних для кедра зонах — предтундровых лесах и южной тайге, она минимальна в средней и северной тайге (0,74-0,78).
Заметные различия параметров факторных полей прослеживаются и для долготных частей Западной Сибири, особенно для наиболее полно нами охваченных западной и центральной. Для факторного поля погодичных колебаний прироста от запада к востоку увеличиваются исходная напряженность (от 0,91 до 1,00), уменьшаются проницаемость географической среды (от -0,05 до -0,09 на 100 км) и максимальная активная протяженность поля (от 1730 до 1130 км). Подобные же тенденции видны для факторного поля многолетних колебаний прироста.
Как показали наши предыдущие исследования [Арефьев, 1997], факторы, определяющие погодичные колебания прироста (особенно с возрастанием амплитуды до стрессовой) в целом более существенны для состояния кедра, чем факторы, определяющие многолетние колебания; кроме того, они в меньшей степени модулируются внутриформационным разнообразием кедровников, составляющим базис устойчивости формации. Прибегая к показательным природным аналогиям, погодичные колебания прироста и отражаемых ими экологических факторов можно уподобить жестким высокочастотным электромагнитным излучениям, действие которых на биологические структуры намного сильнее по сравнению с низкочастотными электромагнитными колебаниями и имеет разрушительный характер. Этой аналогии соответствует и сравнение параметров проницаемости среды. Напряженность факторного поля «высокочастотных» колебаний прироста кедра особенно велика на юге Западно-Сибирской равнины, где на пределе своего ареала вид, особенно в течение последнего столетия, чрезвычайно подвержен энтомовредите-лям и пожарам. При этом последовательное уменьшение напряженности «высокочастотного» поля к северу может свидетельствовать о тенденции смещения зоны кедровой тайги (и лесной зоны в целом [Шиятов и др., 2005]) на север вследствие циклически нарастающего потепления и аридизации климата региона [Полюшкин и др., 1977].
Вместе с тем с увеличением напряженности факторных полей проницаемость их для географической среды уменьшается. На юге региона поля погодичных колебаний прироста имеют большую тотальную (450 км), но сравнительно малую максимальную активную протяженность (960 км); наиболее слабые факторные поля многолетних колебаний прироста на севере региона достигают наибольшей протяженности (до 3410 км). Ширина базиса устойчивости кедровой формации к полевым факторам проявляется в разнообразии
адаптивных путей внутри формации, модулирующих их напряженность на лесотипологическом, популяционно-генетическом, онтоценогенетическом уровнях. На крайнем юге региона, где кедр сохранился спорадически, базис его устойчивости в течение исследованного периода был узок. На севере региона базис устойчивости формации к стрессовым факторным полям, судя по величине модуляции, был широк, и они имели там небольшое тотальное распространение (76 км).
Приведенные данные позволяют заключить, что полученные результаты не только показывают полевые закономерности и географические масштабы разного рода изменчивости радиального прироста кедра на территории Западно-Сибирской равнины, но и соответствуют сделанным нами ранее выводам о низкой устойчивости кедра на его южном пределе распространения в регионе и более высокой — на северной окраине ареала [Арефьев, 1997].
ЛИТЕРАТУРА
Арефьев С. П. Оценка устойчивости кедровых лесов Западно-Сибирской равнины // Экология. 1997. № 3. С. 149-157.
Ваганов Е. А, Шиятов С. Г., Мазепа В. С. Дендроклиматические исследования в Урало-Сибирской Субарктике. Новосибирск: Наука, 1996. 246 с.
Глебов Ф. З., Черкашин В. П., Мацулева Г. Н. Влияние климата на динамику радиального прироста в двух типах кедрового леса // Дендрохронология и дендроклиматология. Новосибирск: Наука, 1986. С. 94-102.
Горчаковский П. Л., Шиятов С. Г. Фитоиндикация условий среды и природных процессов в высокогорьях. М.: Наука, 1985. 209 с.
Гоейбилл Д. А. Дендрохронологическое изучение загрязнения воздушной среды в хвойных лесах западных районов США // Лесоведение. 1990. № 2. С. 3-15.
Комин Г. Е. Изменчивость рангов деревьев по диаметру в древостое // Тр. Ин-та экологии растений и животных Уф АН СССР. Свердловск, 1970. Вып. 67: Лесообразовательные процессы на Урале. С. 252-261.
Комин Г. Е. Изменение рангов деревьев по диаметру в разновозрастных буковых древостоях // Лесоведение. 1994. № 5. С. 41-47.
Комин Г. Е., Пьянков Ю. А, Шиятов С. Г. Определение сходства между дендро-хронологическими рядами // Экология. 1973. № 4. С. 29-34.
Константинов В. Д., Горожанкина С. М. Структура зонального поля на примере лесных экосистем Приенисейской равнины // География и природные ресурсы. 1995. № 2. С. 25-31.
Лакин Г. Е. Биометрия. М.: Высш. шк., 1990. 352 с.
Ловелиус Н. В. Изменчивость прироста деревьев. Дендроиндикация природных процессов и антропогенных воздействий. Л.: Наука, 1979. 230 с.
Малоквасов Д. С. К методике дендроклиматического изучения колебаний радиального прироста в разновозрастных древостоях кедра корейского // Дендрохронология и дендроклиматология. Новосибирск: Наука, 1986. С. 39-43.
Методы дендрохронологии. Ч. 1: Сбор и получение древесно-кольцевой информации: Учеб.-метод. пособие / Шиятов С. Г., Ваганов Е. А., Кирдянов А. В. и др. Красноярск: КрасГУ, 2000. 80 с.
Полюшкин Ю. В., Ларионов Н. В., Мартынов В. Д. Периодичность засух и колебаний прироста деревьев на юге Тюменской области // Долгосрочные прогнозы природных явлений. Новосибирск: Наука, 1977. С. 83-92.
Растительный покров Западно-Сибирской равнины / Ильина И. С., Лапшина Е. И., Лавренко Н. Н. и др. Новосибирск: Наука, 1985. 251 с.
Шиятов С. Г. Дендрохронология, ее принципы и методы // Зап. Свердл. отд-ния ВБО. Свердловск, 1973. Вып. 6. С. 53-81.
Шиятов С. Г. Дендрохронология верхней границы леса на Урале. М.: Наука, 1986. 137 с.
Шиятов С. Г., Терентьев М. М., Фомин В. В. Пространственно-временная динамика лесотундровых сообществ на Полярном Урале // Экология. 2005. № 2. С. 83-90.
Methods of Dendrochronology. Applications in the Environmental Sciences / Eds. E. R. Cook, L. A. Kairiukstis. Dordrecht; Boston; L.: Kluwer Acad. Publ., 1990. 394 p.
ИПОС СО РАН, г. Тюмень
S. P. Arefyev
DENDROCHRONOLOGICAL ANALYSIS OF FACTOR FIELDS REGARDING THE STATE OF THE PINUS SIBIRICA DU TOUR IN WEST SIBERIA
Basing on the analysis of similarity in the oscillations of radial increment of the cedar through the vast territory of West Siberian plain, subject to determination were parameters of ecological (first of all, climatic) factor fields specifying these: their intensity and extension, permeability of the geographical medium, intra-formation modulation, as well as bases of stability of the Pinus sibirica formation to the field factors. The article considers distinctions of the factor fields regarding annual and perennial oscillations of the increment, as well as latitudinal-zonal and longitudinal changes of their parameters.
Rings, oscillations of increment, West Siberia, zonality, forest formations, Pinus sibirica, stability.
