Хакасский технический институт – филиал ФГАОУ ВПО «Сибирский федеральный университет» (г. Абакан)
Эта статья опубликована сборнике научных трудов "Фундаментальные науки и практика" с материалами Третьей Международной Телеконференции "Проблемы и перспективы современной медицины, биологии и экологии" - Том 1 - №4. - Томск - 2010.
Изменения климата Земли являются одной из наиболее обсуждаемых в настоящее время проблем. Оценить эти изменения невозможно без достоверной информации о состоянии климата в прошлом. Однако ряды прямых инструментальных измерений климатических параметров имеют, как правило, недостаточную длительность, а для многих районов такая информация вообще отсутствует. Поэтому большое значение имеют косвенные источники, позволяющие реконструировать климатические данные на основе их связей с составом и структурой природных тел. Одним из ведущих направлений в данной области является дендроклиматология. Наиболее часто при реконструкции погодных условий используется ширина годичного кольца (ШГК), что позволяет достичь временного разрешения в несколько месяцев (например, средняя летняя температура) или год. Однако годы со сходными интегральными параметрами могут значительно различаться по динамике, поэтому актуальным становится вопрос о реконструкции климата с более высоким разрешением, что возможно при использовании, например, гистометрических параметров древесины.
Ряд ранее проведенных исследований [Е.А. Ваганов, А.В. Шашкин, 2000 и др.] показал, что различные показатели годичных колец как зависят от разных климатических факторов, так и связаны между собой. Следует также учесть, что если в экстремальных условиях изменчивость ширины и деталей структуры годичных колец на 50-70% определяется основным лимитирующим фактором, то в других районах местообитания эти зависимости носят комплексный характер. К примеру, в лесостепной зоне Хакасии на формирование древесины хвойных влияют температуры и осадки всего сезона роста. Поэтому в данной работе представляло интерес провести реконструкцию ведущих климатических условий Ширинского района Республики Хакасия на основе данных по ширине годичных колец и гистометрическим параметрам древесины лиственницы сибирской (Larix sibirica Ledeb.), сосны обыкновенной (Pinus sylvestris L.) и ели обыкновенной (Picea obovata Ledeb.). Образцы древесины взяты в лесостепной зоне с двух участков, различающихся по влажности (на южном склоне – сосна PS1 и лиственница LS1, в пойме ручья – лиственница LS2 и ель PO2). В работе были использованы инструментальные климатические данные ближайшей метеостанции (Шира) – суммы осадков по месяцам, средние температуры по месяцам и декадам.
Сбор, транспортировка, первичная обработка кернов, измерение и стандартизация ширины годичных колец проводились по стандартным методикам, принятым в дендрохронологии [С.Г. Шиятов и др., 2000]. Измеренные по фотографиям срезов древесины ряды радиальных размеров клеток D и толщины клеточных стенок CWT были нормированы (количество клеток в годичном кольце приведено к 15) [E.A. Vaganov, 1990]. Для подавления взаимосвязей этих параметров с шириной годичного кольца было проведено их индексирование по методике, ранее использованной А.В. Кирдяновым с соавторами [A.V. Kirdyanov et al, 2007]. Таким образом, были получены четыре локальные индексированные хронологии радиального прироста и четыре серии локальных индексированных клеточных хронологий, охватывающие период с 1969 по 2008 годы.
Проведен корреляционный анализ взаимосвязей параметров древесины с изменениями климата, в результате определены ведущие климатические факторы (температуры и осадки мая-июля) и наиболее чувствительные к ним хронологии. На склоне повышение температуры в этот период подавляет процессы формирования древесины, увеличение осадков – стимулирует. В пойме, в связи с большей степенью увлажнения, наблюдается инверсия климатического отклика в древесно-кольцевых хронологиях. Кроме того, ель чувствительна к колебаниям температуры весной, что выражается в отрицательной реакции параметров ее древесины на температуры мая. Формирование древесины лиственницы начинается позднее и проходит интенсивнее, что также трансформирует климатический сигнал.
Для реконструкции климатических факторов использовался аппарат множественной линейной регрессии. Регрессионные модели рассчитывались методом последовательного исключения переменных, дающих в общую дисперсию наименьший вклад. В связи с большим количеством клеточных хронологий и их тесными взаимосвязями (особенно между соседними клетками), при выборе начального набора независимых переменных для построения модели использовали, во-первых, только те параметры, в которых присутствует климатический отклик на реконструируемый фактор, и во-вторых, из нескольких соседних хронологий выбирали одну, реагирующую наиболее значимо. В итоге были получены регрессионные уравнения, содержащие ШГК и гистометрические параметры либо только гистометрические параметры (для климатических факторов, слабо влияющих на радиальный прирост). Многолетняя динамика измеренных и реконструированных рядов данных показана на рисунке, статистические характеристики моделей приведены в таблице.
Модель
|
Регрессионное
уравнение
|
Критерий
Фишера, F
|
Коэфф.
корр.
с исх., R
|
Дисперсия,
объясн.
моделью
|
Ошибка
оценки,
SE
|
Тмай
|
25.65–12.94*D15iPS1+6.40*D14iLS1–9.75*D2iLS2
|
F(3,36)=14.65
p<0.00000
|
0.741
|
0.512
|
1.42
|
Тиюн
|
18.03+2.04*TRWLS2–4.00*CWT3iLS2
|
F(2,37)=16.79
p<0.00001
|
0.690
|
0.447
|
1.25
|
Тиюл
|
28.13–7.68*D2iPS1–7.09*D14iPS1+4.81*CWT6iPS1
|
F(3,36)=11.29
p<0.00002
|
0.696
|
0.442
|
1.08
|
Тмай-июл
|
24.78–4.88*D15iPS1–2.89*CWT2iLS2–
–2.40*CWT11iPO2
|
F(3,36)=22.27
p<0.00000
|
0.806
|
0.621
|
0.65
|
Рмай
|
36.3–47.0*CWT13iPS1+34.5*D8iLS2
|
F(2,37)=9.79
p<0.00039
|
0.588
|
0.311
|
9.9
|
Риюн
|
198.0+23.7*TRWPS1+17.8*TRWLS2–
–49.2*D11iLS2–51.3*TRWPO2–88.5*D8iPO2
|
F(5,34)=5.83
p<0.00054
|
0.679
|
0.383
|
17.6
|
Риюл
|
–69.6+168.2*D13iLS1–39.0*CWT8iLS2
|
F(2,37)=11.01
p<0.00018
|
0.611
|
0.339
|
30.5
|
Рмай-июл
|
389.2+87.5*TRWPS1–175.6*CWT5iPS1+
+123.2*CWT12iPS1–191.1*D5iLS1–98.0*CWT8iLS2
|
F(5,34)=5.87
p<0.00052
|
0.681
|
0.384
|
37.6
|
Коэффициенты зависимостей реконструированных рядов от исходных, показанные на врезках, близки к единице для всех моделей. Наиболее качественная модель реконструкции получена для температуры мая-июля. Качество моделей реконструкции осадков не настолько высокое, как температур, но достаточно для их дальнейшего использования. Высокая сходимость реконструированных рядов с инструментальными подтверждается сопоставлением их графиков многолетней изменчивости. При этом низкочастотные колебания идут синхронно для всех рядов, а погодичные колебания климатических факторов лучше отражают модели температур. Это связано со значительной погодичной и внутрисезонной изменчивостью осадков в районе исследования, а также с пролонгированным характером их влияния на древесно-кольцевые хронологии.
Подводя итоги, можно отметить следующее:
1. Погодичная изменчивость радиального прироста и гистометрических параметров древесины хвойных в лесостепной зоне Хакасии содержит климатический сигнал. Лимитирующими факторами для формирования древесины в данном случае являются средние температуры и суммы осадков.
2. Уровень климатического сигнала достаточно высок, чтобы реконструировать изменчивость погодных условий во время сезона роста. При этом использование гистометрических параметров структуры древесины позволяет уменьшить временное разрешение, а также реконструировать те факторы, влияние которых на ширину годичного кольца недостаточно сильно.
3. Сравнение статистических характеристик регрессионных моделей позволяет, наряду с корреляционным анализом, выявить наиболее значимые для формирования древесины хвойных периоды внутри сезона роста.
СПИСОК ЛИТЕРАТУРЫ
1. Ваганов Е.А., Шашкин А.В. Рост и структура годичных колец хвойных. Новосибирск: Наука, 2000. 232 с.
2. Шиятов С.Г., Ваганов Е.А., Кирдянов А.В., Круглов В.Б., Мазепа B.C., Наурзбаев М.М., Хантемиров P.M. Методы дендрохронологии. Ч. I. Основы дендрохронологии. Сбор и получение древесно-кольцевой информации. Красноярск: КрасГУ, 2000. 80 с.
3. Kirdyanov A.V., Vaganov E.A, Hudges M.K. Separating the climatic signal from tree-ring width and maximum latewood density records // Trees. 2007. V. 21. P. 37-44.
4. Vaganov E.A. The traheidogram method in tree-ring analysis and its application // Methods of Dendrochronology. Appl. in the Environmental Sciences. Cluwer Acad. Publ., 1990. P. 63-75.
|