дендроклиматическии анализ годичных колец каменной березы (betula ermanп cham.)
на южной камчатке
Е.С. ПУГАЧЕВА, асп. Института географии РАН
Дендроклиматология - это научное направление, которое занимается реконструкцией климатических условий прошлого по годичным кольцам деревьев. Древесные растения реагируют на изменение лимитирующих климатических факторов, таких как температура воздуха, количество осадков и др. В благоприятные по климатическим условиям годы у большей части деревьев формируются широкие кольца, а в неблагоприятные - узкие. Такое чередование колец неповторимо во времени, и чем сильнее годичная изменчивость величины прироста деревьев, тем более надежным индикатором изменений условий среды она является [8].
Многие районы мира, такие как Альпы, Аляска, юго-запад США, исследованы в ден-дроклиматическом отношении очень подробно. Камчатка, которой посвящена эта работа, пока к ним не относится, тем не менее в последние
pugachevaelena@gmail.com годы был выполнен ряд исследований в этой области. Существуют реконструкции температур мая-июня по центральной и восточной Камчатке, выполненные на основе ширины колец лиственницы [12, 16, 17]. Поскольку на южной Камчатке, за исключением кедрового стланика, не растут хвойные деревья, здесь этот подход неприменим. Первые попытки использования каменной березы в дендроклиматологических целях на Камчатке [16] носили экспериментальный характер. Они показали, что использовать каменную березу можно, но она крайне трудно поддается перекрестному датированию из-за неясных границ между годичными кольцами. В международной базе дендрохронологических данных (NOAA, www.noaa.gov), где собраны хронологии со всего мира, нет хронологий каменной березы.
В этой работе приводятся результаты дендроклиматического исследования образ-
ЛЕСНОЙ ВЕСТНИК 1/2009
67
цов каменной березы, отобранных с четырех площадок на Южной Камчатке, которые позволяют судить о связи приростов березы с климатическими параметрами.
Полуостров Камчатка располагается в восточной части Евразии. Камчатка - это горная территория, с развитым оледенением (446 ледников, площадью 900 км2), с активной вулканической деятельностью, где 29 активных и около 300 потухших вулканов. Климат Камчатки находится преимущественно под влиянием процессов циркуляции атмосферы, развивающихся над северной частью Тихого океана. Главным климатооб-разующим фактором является интенсивная циклоническая деятельность. По количеству выпадающих осадков Камчатка относится к зоне достаточного увлажнения, а по гидротермическому режиму преобладают территории с гумидным климатом. Этот характерный и общий для всего ареала каменной березы признак соответствует биологическим требованиям породы к условиям среды [7]. Среднегодовые температуры воздуха - около 0 °С в южной части полуострова в районе города Петропавловск-Камчатский. Зимы холодные и снежные, а лето влажное и облачное. Летние температуры, по данным метеостанции Петропавловск-Маяк, составляют +10 °С, зимние температуры около - 6 °С.
Мощный снежный покров является важной особенностью климата Камчатки, он оказывает влияние на многие природные процессы. Количество выпадающих и отложенных твердых осадков весьма неравномерно, и их распределение зависит от расчлененности рельефа и ветрового режима [5]. Наибольшее количество осадков (2000 мм) выпадает на наветренных склонах хребтов южного вулканического района и в центральной части Кро-ноцкого полуострова. Наименьшее (300 мм) - на севере области и в центральной части долины реки Камчатки. На остальной территории суммы осадков изменяются от 1600 мм в восточном горно-вулканическом районе до 400 мм на северо-западном побережье [3]. Твердые осадки выпадают с октября по май и составляют около 45 % годовой суммы [2].
Согласно [6] на Камчатском полуострове выделяются следующие растительные
пояса (снизу вверх): 1) горно-таежный пояс светлохвойных и темнохвойных лесов; 2) пояс каменноберезовых лесов; 3) пояс ольхового и кедрового стлаников; 4)пояс горных кустарничково-лишайниковых тундр.
Каменная береза является наиболее распространенной древесной породой. Роль лесообразующей породы она играет на меньшей территории. Отдельные ее рощи, а затем большие массивы появляются севернее мыса Лопаткина, занимая почти весь юг полуострова. Далее на север, разделяясь на две полосы, леса каменной березы тянутся до 58-60° с.ш. (рис. 1). Наибольшая высота распространения каменноберезняков отмечается на внутренних склонах Центральной депрессии, где граница сплошных массивов проходит на высоте 700-800м н.у.м. На противоположных склонах, обращенных к Охотскому морю и Тихому океану, граница снижается до 600 м н.у.м. Также происходит снижение верхней границы березняков к северу и к югу полуострова по мере сужения суши и усиления морского влияния [7].
Рис. 1. Ареалы распространения растительности на Камчатке: 1 - сфагновые болота; 2 - береза каменная; 3 - кедровый стланик; 4 - горные тундры; 5 - лиственничные леса; 6 - горнотаежные темнохвойные леса (Атлас СССР, 1962 г.)
Таблица 1
Характеристика площадок отбора образцов каменной березы, южная Камчатка
Площадка Расположение N E Н, м Кол-во кернов Средняя взаимная корреляция образцов Средняя чувствительность
NACH Начикский пер. 53 06 441 157 50 459 424 21 0,642 0,396
ZAV Завойко 52 57 324 158 40 114 120 22 - -
PAR Паратунка 52 58 229 158 15 308 60 18 0,624 0,376
VIL Мутновский пер. 52 41 473 158 09 699 350 20 0,612 0,329
Ствол каменной березы, как правило, искривленный, рано теряет преобладание роста над боковыми ветвями, поэтому крона начинается невысоко (4-7 метров) над поверхностью почвы. Под влиянием внешних условий каменная береза на Камчатке образует в условиях влажного и туманного климата широкую раскидистую крону, что повышает ее устойчивость к неблагоприятным факторам. Каменноберезняки в верхней части пояса задерживают кронами и конденсируют низкую облачность и увеличивают снегонакопление, улучшая тем самым теплообмен с воздухом и снижая промерзание почвы. Деревья достигают возраста 30-400 лет, однако они очень часто поражаются сердцевинной гнилью, так что реальный набор сохранившихся годичных колец у них редко превышает 100-150 лет. Большое влияние на распространение и рост данной породы оказывают такие факторы, как теплообеспеченность, инсоляция, высота снежного покрова, степень промерзания почвы, влажность почвы и воздуха. Каменная береза хорошо переносит содержание в воздухе сернистых газов и погребение почв вулканическими пеплами, что дает ей возможность расти вблизи от действующих вулканов и горячих сернистых источников [7].
В 2006 г. на южной части полуострова Камчатка был произведен отбор кернов березы каменной (Betula Ermani Cham.) для дендрохронологического анализа. Были отобраны образцы на четырех площадках: две на высотах 350 и 424 м.н.у.м., одна - на высоте 60 м.н.у.м. Эти три площадки располагаются на территории, в основном свободной от антропогенного воздействия. Четвертая пробная площадка была выбрана в пределах города Петропавловск-Камчатский на высоте 120 м.н.у.м., (табл. 1). Для отбора и датиро-
вания образцов использованы стандартные дендрохронологические методы [8].
С каждого дерева отбиралось по два керна. Если дерево оказывалось молодым или поврежденным сердцевинной гнилью, так что набор сохранившихся колец был небольшим, второй керн не отбирали. Всего был отобран 81 образец. Ширина годичных колец измерена с точностью около 0,01 мм с помощью полуавтоматической установки для измерения ширины годичных колец LINTAB. Наименьшая полученная ширина кольца составила 0,06 мм, а наибольшая - 5 мм. Выпадающих колец ни в одном из образцов не обнаружено. Перекрестное датирование проводилось при помощи программы COFECHA [12], которая позволила отбраковать так называемые «ленивые керны» с низкой межгодовой изменчивостью ширины колец, не поддающиеся перекрестному датированию. Это часто случается в тех районах, где деревья растут в условиях, близких к оптимальным [15]. Мы не включили в хронологию образцы, которые коррелировали с основной хронологией с коэффициентом корреляции менее 0,40, полагая, что доминирующим в их росте является неклиматический сигнал. Интересно, что образцы, отобранные на площадке 2ЛУ, которая располагается непосредственно на территории города Петропавловск-Камчатский, также не коррелируют с общей хронологией; более того, они даже не датируются между собой (рис.2). Это скорее всего связано с активной антропогенной деятельностью в районе их отбора, в том числе и по причине сведения леса в связи с хозяйственными нуждами населения в XX в.
Ширина колец образцов каменной березы, отобранных на расстоянии в несколько десятков километров друг от друга, коррелируют между собой, и эта корреляция достаточно тесная (табл. 2).
Таблица 2
Коэффициенты корреляции ширины колец образцов каменной березы,
прошедшие тест контроля качества перекрестного датирования, полученные при помощи программы COFECHA
Seq Series Time span 1800 1825 1850 1875 1900 1925 1950 1975
1849 1874 1899 1924 1949 1974 1999 2024
1 parl3b 1900 2005 0,48 0,44 0,79 0,84
2 parl4a 1890 2005 0,49 0,46 0,56 0,70 0,72
3 parl4b 1910 2005 0,39 0,49 0,53 0,61
4 par5a 1864 2005 0,42 0,63 0,49 0,33 0,51 0,60
5 par5b 1870 2005 0,55 0,64 0,60 0,56 0,39 0,47
6 par9a 1880 2005 0,42 0,58 0,44 0,43 0,48
7 par9b 1878 2005 0,58 0,60 0,59 0,66 0,69
8 par3a 1940 2005 0,57 0,59 0,66
У par3b 1911 2005 0,41 0,62 0,64 0,69
10 par8a 1920 2005 0,51 0,50 0,45 0,52
11 par8b 1880 2005 0,51 0,62 0,43 0,45 0,52
12 par2b 1890 2005 0,43 0,58 0,54 0,44 0,56
13 parl6a 1890 2005 0,47 0,57 0,67 0,40 0,45
14 parl6b 1890 2005 0,36 0,48 0,58 0,64 0,66
15 par4a 1880 2005 0,46 0,63 0,67 0,66 0,71
16 nach21b 1920 2005 0,56 0,53 0,63 0,66
17 Nach21a 1910 2005 0,57 0,45 0,58 0,59
18 nachl9a 1910 2005 0,55 0,45 0,54 0,55
19 nachl9b 1900 2005 0,59 0,48 0,53 0,59
20 nachl5a 1848 2005 0,52 0,48 0,50 0,38 0,41 0,60 0,64
21 nachl5b 1856 2005 0,59 0,70 0,69 0,52 0,56 0,58
22 nach7a 1840 2005 0,56 0,59 0,74 0,61 0,50 0,71 0,72
23 nach7b 1842 2005 0,46 0,60 0,73 0,58 0,55 0,60 0,63
24 nach8a 1870 2005 0,77 0,76 0,64 0,68 0,80 0,85
25 nach8b 1824 2005 0,51 0,50 0,69 0,74 0,68 0,65 0,78 0,78
26 nach20a 1880 2005 0,49 0,51 0,60 0,63 0,65
27 nach20b 1850 2005 0,53 0,67 0,59 0,57 0,49 0,54
28 nachl2a 1875 2005 0,57 0,77 0,66 0,69 0,71
29 nachl2b 1870 2005 0,58 0,60 0,65 0,61 0,70 0,70
30 nachl4a 1916 2005 0,55 0,52 0,69 0,71
31 nachl4b 1870 2005 0,42 0,52 0,71 0,63 0,65 0,66
32 nach6a 1890 2005 0,53 0,51 0,57 0,59 0,58
33 nach6b 1890 2005 0,48 0,49 0,46 0,59 0,65
34 nachl3a 1870 2005 0,53 0,63 0,75 0,67 0,67 0,65
35 nachllb 1909 2005 0,55 0,62 0,60 0,63
36 nachlla 1910 2005 0,63 0,63 0,58 0,63
37 vil3b 1891 2005 0,56 0,68 0,45 0,60 0,57
38 vil6b 1841 2005 0,64 0,61 0,59 0,59 0,56 0,69 0,77
39 vil6a 1822 1999 0,62 0,63 0,70 0,70 0,59 0,60 0,79
40 vil2a 1855 2005 0,52 0,41 0,43 0,56 0,46 0,48
41 vil2b 1920 2005 0,46 0,42 0,39 0,45
42 vill8a 1861 2004 0,68 0,59 0,48 0,37 0,55 0,43
43 vill8b 1841 2005 0,66 0,67 0,52 0,36 0,42 0,48 0,52
44 vill4a 1870 2005 0,46 0,44 0,36 0,40 0,69 0,69
45 vill4b 1890 2005 0,45 0,42 0,56 0,65 0,71
46 villa 1890 2005 0,50 0,55 0,54 0,59 0,58
47 villb 1887 2005 0,46 0,52 0,45 0,51 0,55
48 vill3a 1852 2005 0,71 0,76 0,45 0,47 0,65 0,61
49 vill3b 1872 2005 0,74 0,74 0,64 0,52 0,58 0,58
50 villla 1864 2005 0,69 0,66 0,61 0,57 0,58 0,63
51 villlb 1820 2005 0,42 0,40 0,64 0,51 0,45 0,55 0,48 0,49
52 vill5a 1852 2005 0,56 0,63 0,53 0,55 0,69 0,70
53 vill5b 1840 2005 0,53 0,75 0,66 0,50 0,60 0,80 0,80
54 vill9a 1900 2005 0,56 0,61 0,60 0,59
55 vill9b 1878 2005 0,52 0,55 0,47 0,51 0,58
56 vill2a 1850 2005 0,63 0,50 0,41 0,58 0,76 0,76
57 vill2b 1945 2005 0,62 0,69 0,69
Av segment correlation 0,52 0,54 0,60 0,57 0,55 0,54 0,60 0,63
Таблица 3 Коэффициенты взаимной корреляции хронологий ширины колец каменной березы
VIL NACH
VIL
NACH 0,63
PAR 0,33 0,38
По-видимому, это означает, что прирост березы в большой степени определяется общим внешним фактором, и это, скорее всего, фактор климатический. Различий в годовом приросте колец березы на площадках, отличающихся по высоте на 300-400 м, не выявлено. Это позволяет нам объединить образцы в сводную общую хронологию. Кроме того, мы построили три локальных хронологии (PAR, NACH, VIL), используя для этой цели программу ARSTAN [9, 13] (рис. 2). Эта программа рассчитывает индивидуальные индексы прироста каждого образца, аппроксимируя его возрастной тренд линейной или криволинейной зависимостью (в зависимости от выбора оператора). Это необходимо для
удаления возрастного тренда и достижения сопоставимости ширины колец отдельных образцов между собой и приведения их к безразмерному виду. Хронологии представляют собой осредненные погодично значения индексов индивидуальных образцов. На рис. 2 видно, что наши хронологии имеют явное сходство (табл. 3), особенно заметное для последних 100 лет, что связано с более низкой обеспеченностью образцами ранней части хронологий.
Поскольку все площадки располагаются недалеко друг от друга и коррелируют между собой, они были объединены в одну общую хронологию, построенную RCS методом (рис. 3). Этот метод позволяет максимально сохранить долгопериодный сигнал в древесно-кольцевых хронологиях и, следовательно, более адекватно реконструировать и прогнозировать изменения климата на этой основе [11]. Всего в хронологию вошли 57 образцов с трех площадок (1822-2004 гг.). Мы использовали как минимум три керна для построения хронологии для каждого года (рис. 4).
2 1,8 1,6 1,4 1,2 1 0,8 0,6 0,4 0,2 0
i
1 t л AM1 , /1
11 tí í A M ш L 1 l\ U.i/\. Uli N.'i M .я/w к
i ,д A ívvsp, ж . ,* J к m\\J , ]Ч> A i i w 1 1 Ail V \ -M1
41 »VI A Л И A " » I/»l/' r'j\ WWi' W Si VF1 n ЦЯХ iv 'wvv шш * II \í 1
( vJ 1 \ ™ 1 , У i vvfl \ i/I ti
v м/ ' V ' v ^ v - ' 1
V I V 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
о W
(D
ч и s
1840
1880
1920
1960
2000
_ _ _ VIL-NACH-PAR Годы
Рис. 2. Стандартная ARSTAN хронология для площадок VIL, NACH, PAR
1822 1842 1862 1882 1902 1922 1942 1962 1982 2002
Годы
Рис. 3. Сводная RCS хронология для всех образцов
01........."........"....................................у...................................................У...................................................У.....
1822 1842 1862 1882 1902 1922 1942 1962 1982 2002
Годы
Рис. 4. Количество образцов для каждого года в сводной стандартной хронологии (не менее трех)
h I
<D
s э s
-е -е
£
0,4 0,3 0,2 0,1 0 -0,1 -0,2 -0,3 -0,4
Температура
1 I |
1 I ■ 1 1 1
■ III 11.. II II
■ J A S О ■ D j f m a m j j a s
Месяц
Рис.5. Коэффициенты корреляции сводной стандартной хронологии каменной березы со среднемесячной температурой
н
к s
а —
я а
а й
s R
•& g
й о
О ¡5
0,4 0,2 0 -0,2 -0,4
Осадки
M.A.S.9.N.PM.I
■ i1 i
a ■ J J
■ ■ i i ■
Месяц
Рис.6. Коэффициенты корреляции сводной стандартной хронологии каменной березы со среднемесячными осадками
Для выявления лимитирующих факторов произрастания каменной березы были использованы данные по осадкам и температуре, полученные на метеостанции Петропавловск-Маяк, которая расположена в районе пробных площадей. Непрерывный ряд метеоданных с выбранной метеостанции, которые были использованы для получения корреляционных связей: осадки с 1895 по 1994 гг., температура воздуха с 1891 по 1994 г.
Для получения функции отклика были рассчитаны коэффициенты корреляции прироста древесины со средней температурой и осадками для каждого месяца с июня предшествующего года по сентябрь текущего года (всего 16 месяцев) для периода с 1891 по 1994 г. для температуры (рис. 5) и с 1895 по 1994 г. для осадков (рис.6). В результате была получена функция отклика как для отде-
льных локальных хронологии, так и для сводной стандартной хронологии, рассчитанной с помощью программы ARSTAN.
На рис. 5 и 6 видно, что имеется статистически значимая корреляция индексов прироста со среднемесячными значениями температуры воздуха с июня по сентябрь текущего года, которые составили: для июня 0,23, июля 0,38, августа 0,33, сентября 0,26 (рис. 5). Температуры предыдущих лет, видимо, не оказывают существенного влияния на прирост, хотя предшествующие теплые сезоны также имеют тенденцию положительно влиять на ширину кольца. Корреляции с осадками не являются статистически значимыми и составляют не более - 0,2 (рис. 6), но являются отрицательными в летние месяцы текущего года. В целом береза скорее отрицательно реагирует на избыток осадков, что и понятно, принимая во внимание высокую увлажненность на Камчатке.
Полученные нами данные интересно сравнить с аналогичными исследованиями, проводившимися в центральной Японии на верхней границе леса горы Норикура (Mount Norikura) (1600-2500 м н.у.м.) [14]. В работе получены статистически значимые положительные корреляции прироста каменной березы с летними температурами текущего года (коэффициент корреляции 0,2-0,3) и отрицательная корреляция со среднемесячными осадками июня-сентября текущего года (коэффициент корреляции -0,25/-0,35). Эти результаты согласуются с нашими данными для южной Камчатки. Существенно иные статистические связи получены для каменной березы в северной части Китая, в горах Чангбаи (Changbai Mountain) [10] на высоте 17002000 м н.у.м. Для этого района обнаружены высокие коэффициенты корреляции ширины колец каменной березы со среднемесячными значениями температуры октября предыду-
щего года (+0,38) и среднемесячными значениями осадков ноября и декабря предыдущего года (+0,45). Положительная связь с осадками в данном случае объясняется авторами тем, что маломощные с большим включением вулканических пород почвы на этой территории имеют незначительную влагоемкость, поэтому растения испытывают недостаток в почвенной влаге. На рост каменной березы в этом районе оказывают влияние как жидкие осадки, которые выпадают преимущественно во время сезона роста (максимум в июле-августе), так и снегозапасы предыдущего зимнего периода, которые положительно влияют на ранний прирост древесины [10].
Помимо климатических факторов на прирост березы оказывают влияние и другие процессы, в частности, упоминавшиеся здесь ранее антропогенные воздействия и термальные воды, повышающие температуру почвы и локальную температуру воздуха. Влияние вулканической деятельности на растительность очень существенно [4]. В частности, деревья страдают от термических и химических ожогов, физических повреждений во время извержений и сопутствующих им катастрофических событий. Это вызывает уменьшение годичного прироста или ведет к выпадающим кольцам [18]. С другой стороны, пепел, покрывающий снег, может вызвать более раннее снеготаяние и продлить, таким образом, вегетационный период, что, скорее всего, положительно скажется на ширине годичного кольца. Тонкий слой пепла, богатый разными минеральными примесями, также обычно вызывает увеличение годового прироста. С другой стороны, если этот слой превосходит 20-30 см, деревья могут погибнуть или испытать заметное угнетение. Все это, однако, вряд ли применимо к рассмотренным здесь случаям, поскольку выбранные нами площадки расположены далеко от действующих вулканов.
По полученной нами хронологии каменной березы можно четко выделить годы с наименьшей шириной кольца (рис. 2) - 2004, 1985, 1935, 1918, 1887, 1831, 1828 гг. Первые четыре года покрываются интервалом инструментальных исследований. Два из них (как видно на рис. 6), совпадают с резким
снижением среднелетних температур (1985 и 1918 гг.). В эти же годы наблюдался пониженный прирост лиственницы [16]. Мы затрудняемся на этом этапе исследования объяснить минимумы приростов в 2004 и 1935 гг.; в хронологии лиственницы эти годы ничем не примечательны, и прирост близок к среднему многолетнему значению. 1887 г. был, по-видимому, холодным на всей Камчатке - средняя ширина кольца в сводной хронологии лиственницы меньше двух стандартных отклонений. В хронологию березы в период до середины XIX в. входит небольшое число образцов; таким образом, ее пониженные приросты в 1831 и 1828 гг. нуждаются в дополнительном обосновании.
Помимо межгодовой изменчивости, хронология каменной березы отражает и более долгопериодные колебания климата. В частности, на рис 3. видно, что за последние полтора столетия выделяется два периода, неблагоприятных в целом для ее роста в 1880-е и в 1970-е гг. Эти периоды были холодными и по метеоданным. По данным о ширине колец лиственницы на Камчатке в 1860-1880-е гг. наблюдалось самое значительное похолодание за последние 350 лет [16]. Температуры начала летнего сезона были на 1,5 °С ниже по сравнению с периодами потепления в середине XX в. [12]. В последний этап этого периода (1880 гг.) наблюдалась аномально высокая аккумуляция снега, что привело к наступа-нию ледников. Многие ледники стабилизировались в отступании в период похолодания 1970-1980-е гг. [16]. В динамике ширины колец березы на Южной Камчатке отсутствует долгопериодный тренд. Не наблюдается он и за последнее столетие по инструментальным данным (рис. 7).
Каменная береза - широко распространенная древесная порода, которая может быть с успехом использована для ден-дрохронологических целей. В этой работе мы приводим первую сводную региональную хронологию каменной березы для южной Камчатки для периода XIX и XX вв. Прирост березы в этих условиях зависит от температуры вегетационного периода (май-сентябрь) текущего года. Влияния осадков не обнаружено.
13
с 12
й 11
й 10
& 9 е § 8
£ 7
Средняя температура, июль-август
■A iwvvaa. . а<чкл
f-^—--———
1891 1901 1911 1921 1931 1941 1951 1961 1971 1981 1991
Год
Рис. 7. Среднелетняя температура воздуха (м/с Петропавловск-Маяк)
Сходная реакция на колебания температуры обнаружена ранее в Японии. Полученные результаты дают надежду на продолжение исследований в этом направлении, в том числе в плане использования остатков древесины, погребенных в вулканических отложениях, для датирования событий, связанных с вулканической активностью, и для продолжения дендрохроно-логического ряда.
Автор выражает благодарность за помощь в проведении полевых исследований и ценные советы О.Н. Соломиной, В.Н. Михаленко и Я.Д. Муравьеву Частичное финансирование работы было осуществлено в рамках программы РАН «Природные катастрофы».
Библиографический список
1. Атлас СССР. Карта растительности СССР. - М.: ГУГК, 1962.
2. Виноградов, В.Н. Современное оледенение районов активного вулканизма. Результаты исследований по международным геофизическим проектам / В.Н. Виноградов. - М.: Наука, 1975.
3. Глазырин, Г.Е. О климатическом фоне оледенения Камчатки. Гляциологические исследования № 27 / Г.Е. Глазырин, Я.Д. Муравьев, В.Н. Виноградов // Результаты исследования по международным геофизическим проектам АН СССР. - М., 1985. -С. 51-66.
4. Гришин, С.Ю. Влияние на растительность извержения группы Туйла на Ключевской сопке, Камчатка / С.Ю. Гришин // Вестник КАРУНЦ, 2007. - Т. 2 - № 10. - С. 9-16.
5. Муравьев, Я.Д. Снежный покров горных районов Камчатки / Я.Д. Муравьев. // Вопросы географии Камчатки. - Вып. 9. - Петропавловск-Камчатский: АН СССР, 1985. - С. 30-40.
6. Турков, В.Г. Лесоводственно-таксационная характеристика каменноберезовых древостоев Камчатки / В.Г. Турков, В.А. Шамшин // Леса Камчатки и их лесохозяйственное значение. - М.: АН СССР, 1963.
7. Шамшин, В.А Каменноберезовые леса Камчатки / В.А. Шамшин. - М.: ГЕОС, 1999.
8. Шиятов, С.Г. Методы дендрохронологии / С.Г. Ши-ятов, Е.А. Ваганов, А.В. Кирдянов и др.: уч.-мет. пос. - Красноярск: КГУ, 2001.
9. Cook E.R. A time-series approach to tree-ring standartization. Ph.D. Thesis. University of Arizona, Tucson, Arizona, 1985.
10. Dapao Yu, Geoff G. Wang, Limin Dai, Qingli Wang. Dendroclimatic analysis of Betula Ermanii forest at their upper limit of distribution in Changbai Mountain, Northeast China. Forest Ecology and Management 2007, 240: 105-113.
11. Esper, J., Cook, E.R., Schweingruber, F.H. (2002). «Low-frequency signals in long tree-ring chronologies for reconstructing past temperature variability.» Science 295: 2250-2253.
12. Gostev M., Wiles G., D'Arrigo R., Jacoby G., and Khomentovsky M. Early summer temperature since 1670 A.D. for Central Kamchatka reconstructed based on Siberian larch tree-ring width chronology, Can. J. For. Res., 26, 2048-2052, 1996.
13. Holmes, R.L. Dendrochronology Program Library, Laboratory of Tree-Ring Research, University of Arizona, Tucson, Arizona, 1992.
14. Kochi Takahashi, Yohei Tokumitsu, Koh Yasue. Climatic factors affecting the tree-ring width of Betula Ermanii at the timberline on Mount Norikura, central Japan. Ecol. Res. 2005, 20: 445-451
15. Schweingruber, F. H. Tree rings: Basics and applications of dendrochronology. Dordrecht, Riedel Publ., 1988
16. Solomina O., Wiles G., Shiraiwa T., D'Arrigo R., Multiproxy records of climate variability for Kamchatka for the past 400 years. Clim. Past, 3, 119128, 2007
17. Solomina O.N., Muravyev Y.D., Braeuning A., Kravchenko G.N. Two new ring width chronologies of Larch and Birch from the Kamchatka Peninsula (Russia) and their relationship to climate and volcanic activities. Cryospheric studies in Kamchatka II, Institute of low temperature science, Hokkaido University, 111-124, 1999.
18. Yamaguchi, D.K., Hoblitt, R.P. 1995. Tree-ring dating of pre-1980 volcanic flowage deposits at Mount St. Helens, Washington. Geological Society of America Bulletin 107(9): 1077-1093.