Climatal signal in different paramers of tree-rings of pine in the Solovetsky Archipelago.pdf Введение Древесно-кольцевые хронологии (ДКХ) - важный архив, хранящий информацию об изменении климата в прошлом. Особенно велика роль дендрохронологии для реконструкций климата последних двух тысячелетий, для которых имеется обширный материал по живым деревьям, архитектурной, археологической и погребенной древесине. Обилие ДКХ позволяет проводить пространственные реконструкции годичного разрешения на больших территориях для основных климатических параметров - температуры [Schweingruber, Shiyatov, Shishov, 2000; Briffa, Osborn, Schweingruber, 2004; Ljungqvist et al., 2019] и увлажнения [Cook et al., 2015; 2020; Anchukaitis et al., 2017; Ljungqvist et al., 2020]. Разные параметры годичных колец (например, ширина, плотность, соотношение стабильных изотопов и анатомические особенности древесины) по-разному реагируют на различные комбинации гидрометеорологических факторов, поэтому в каждом отдельном случае требуется определить наиболее подходящие для реконструкции параметры годичных колец [Briffa, Osborn, Schweingruber, 2004]. Традиционно для палеоклиматических реконструкций использовалась ширина годичных колец [Fritts, 1966; Schweingruber, 2012], но затем было показано, что во многих случаях более чувствительным индикатором климата служит максимальная плотность годичных колец [Schweingruber et al., 1978; Schweingruber, Briffa, Nogler, 1993; Briffa et al., 1998]. Максимальная плотность поздней древесины (MXD) главным образом используется для реконструкции летних температур [Frank, Esper, 2005; Esper et al., 2014, 2015, 2018; Rydval et al., 2014]. Хронологии по MXD меньше подвержены влиянию биологических особенностей деревьев, что является дополнительным преимуществом для палеоклима-тических реконструкций. Измерение плотности годичных колец обычно осуществляется методом рентгенографии, очень трудоемко. Поэтому в качестве заменителя рентгеновской плотности используют оптическую плотность - отраженный видимый свет [Yanosky, Robinove, 1986; Yanosky, Robinove, Clark, 1986; Clauson, Wilson, 1991]. McCarroll и соавт. [2002] установили, что яркость в синем диапазоне спектра имеет более высокую корреляцию с максимальной плотностью древесины, чем в зеленом, ультрафиолетовом и красном диапазонах. Именно отраженную яркость в синем диапазоне спектра, называемую Blue Intensity (BI), и используют в качестве показателя оптической плотности. Blue Intensity обратно коррелирует с плотностью древесины, т. е. древесина с более высокой плотностью отражает меньше света и выглядит темнее [Sheppard, Graumlich, Conkey, 1996; Campbell et al., 2007]. Хронологии по оптической плотности и основанные на них палеореконструкции построены в настоящее время уже для многих районов земного шара [McCarroll et al., 2013; Bjorklund et al., 2015; Linderholm et al., 2015; Fuentes et al., 2018], на северо-западе [Tene et al., 2011; Wilson et al., 2012; Rydval et al., 2017], западе [Trachsel et al., 2012] и востоке Европы [Kaczka et al., 2017, 2018; Rydval et al., 2018], в Северной Америке [Wilson et al., 2014, 2017a] и на Кавказе [Dolgova, 2016], а также в тропиках [Buckley et al., 2018]. Они использованы в полушарных реконструкциях температуры [Wilson et al., 2016; Anchukaitis et al., 2017]. В некоторых районах BI оказалась чувствительной к осадкам [Dannenberg, Wise, 2016; Seftigen et al., 2020], а также использовалась для выявления вспышек вредоносных насекомых [Arbellay et al., 2018]. Оптическая плотность нашла применение и в археологии. Как показывают некоторые работы, датировки, полученные на основе BI, часто показывают более статистически значимые результаты, чем аналогичные датировки по ширине колец [Mills et al., 2017; Wilson et al., 2017b]. Лесные насаждения Соловецкого архипелага давно являются предметом исследований [Ипатов и др., 2009; Соболев, Феклистов, 2011]. В работе Ло-велиуса и соавт. [2012] проанализированы хронологии на основе ширины годичных колец по сосне и ели продолжительностью 300 и 255 лет соответственно. Авторы установили схожие особенности изменения прироста годичных колец у двух лесообразующих пород на острове Большой Соловецкий. Была выявлена связь ширины колец с температурой воздуха, атмосферными осадками и солнечной активностью. Соломина и соавт. [2011] построили первую сводную хронологию по ширине колец продолжительностью 824 года (1185-2008 гг.) на основе живых деревьев и древесины из архитектурных памятников. В работе Мацковского [2013] было показано, что образцы сосны и ели датируются между собой, что свидетельствует о сходстве условий, на которые реагирует ширина колец этих двух пород. Долгова и соавт. [2019] установили, что большинство локальных хронологий сосны значимо коррелирует с температурой июля, причем теснота корреляции почти не зависит от типа ландшафта. В нашем исследовании мы оцениваем и сравниваем потенциал различных параметров годичных колец сосны обыкновенной - ширины кольца, ширины ранней и поздней древесины, оптической плотности поздней древесины для целей палеореконструкции на Соловецких островах. Район работ Соловецкие острова - крупнейший архипелаг Беломорского бассейна. Они находятся в сравнительно мелководной западной половине Белого моря при входе в Онежский залив, образуя западный и восточный проходы в него [Природная среда2007]. Соловецкие острова имеют ледниковое происхождение. Основа островов - принесенные ледником моренные отложения, состоящие из валунов, гравия, песка и глины. Помимо холмистоморенного типа рельефа, выделяются рельеф морских террас и заболоченные зандровые поверхности [Богославский, 1966; Колосова, 1999; Природная среда..., 2007]. На архипелаге преобладают подзолистые почвы, но встречаются и другие типы - от примитивных карликовых хрящеватых подзолов на валунах до плодородных дерново-глеевых на подветренных склонах холмов. Несмотря на то что Соловецкие острова расположены недалеко от Полярного круга, здесь сформировались климатические условия, нехарактерные для этой географической широты. Климат островов сравнительно мягкий, умерено теплый морской, переходный к континентальному и отличается от материковой части на той же широте. Различие выражается в меньшем перепаде температур по сравнению с континентом, в высокой влажности воздуха и в относительном запаздывании времен года (рис. 1). Вегетационный период со среднесуточными температурами выше +5 °С длится в среднем 128 дней, с 20 мая по 1 октября. В Архангельске он длиннее почти на декаду. Активная часть вегетационного периода с температурами выше 10 °С продолжается 74 дня, с 20 июня по 3 сентября [Ипатов и др., 2009]. Month pre tem Рис. 1. Климатические данные (осадки и температура) по данным CRU TS 4.01 [Harris et al., 2014] из узла точки (64°75' с.ш., 35°75' в.д.), ближайшего к месту наших работ Fig. 1. Climatic data (precipitation and temperature) according to CRU TS 4.01 [Harris et al., 2014] from the point node (64°75' N, 35°75' E) closest to the site of our work Для Соловецких островов характерны сильные ветровые нагрузки, которые влияют на рост деревьев. Своеобразный климат, рельеф и почвы формируют лесорастительные условия от весьма благоприятных до совершенно непригодных для произрастания леса. По лесорастительному районированию Соловецкий архипелаг относится к северной подзоне тайги [Ипатов и др., 2009]. Здесь проходит граница двух физико-географических областей - леса и тундры, между которыми существует полоса переходных ландшафтов. В работе Природная среда [2007] были выделены следующие природные комплексы: 1 - ельники; 2 -сосняки сфагновые и болота с сосной; 3 - сосняки лишайниковые; 4 - березняки и осинники; 5 - березовые криволесья в сочетании с фрагментами воро-ничных тундр и болот; 6 - вороничные тундры с фрагментами березовых криволесий и болот; 7 -необлесенные кустарничково-сфагновые болота; 8 -березово-еловые леса. Материалы и методы исследования Для проведения дендроклиматологических исследований нами была выбрана экспериментальная площадка - B75s (65°03' с.ш., 35°64' в.д.). На ней с помощью ручного приростного бура Пресслера на высоте 1-1,5 м от поверхности земли были отобраны 42 керна сосны из 21 дерева. Дальнейшая обработка образцов проводилась в дендрохронологической лаборатории Института географии РАН. Проводились измерения четырех параметров годичных колец (ширина кольца, ширина ранней и поздней древесины, оптическая плотность поздней древесины). Ранняя древесина (внутренняя часть годичного кольца) обращена к сердцевине, светлая и мягкая; поздняя (наружная часть) - к коре, темная и твердая. Различие между ранней и поздней древесиной у хвойных деревьев ясно выражено. Ранняя древесина образуется в начале лета и служит для проведения воды вверх по стволу. Она отличается наличием крупных, достаточно хорошо заметных трахеид. Поздняя древесина откладывается к концу лета и выполняет в основном механическую функцию [Шиятов и др., 2000]. Оптическая плотность - это интенсивность отраженного света в синем диапазоне, которую получают по изображениям кернов в высоком разрешении. Перед проведением измерений из древесных кернов удалялась смола. Экстракция происходила с помощью нагретого растворителя (в нашем случае этиловый спирт) в аппарате Сокслета. Прибор позволяет производить многократное повторение процесса за счет замкнутости системы. Удаление смол из древесных кернов длится не менее 48 ч. После извлечения из экстрактора аппарата Сокслета керны были высушены в вытяжном шкафу и приклеены на специальные деревянные подложки. Далее на микротоме верхний слой древесины был срезан для получения плоской поверхности. Затем на ленточно-шлифовальной машине была произведена шлифовка образцов с последовательными сменами наждачной бумаги (до 1 000 зерен) до получения гладкой поверхности, позволяющей четко наблюдать структуры древесины при увеличении. Затем поверхность кернов была отсканирована с разрешением 3 200 dpi на устройстве Epson Perfection V700, сопряженном с программным обеспечением Silverfast. Важнейшим этапом является калибровка сканера с помощью цветовых карточек «Кодак». Отсканированные керны были измерены в программе CooRecorder, которая позволяет производить измерения различных параметров годичного кольца на одном керне [Larsson, 2013]. Оператор, работающий с этой программой, в полуавтоматическом режиме устанавливает границы годичных колец. Для измерения оптической плотности важным аспектом является оптимальные размеры выбора «окна», с которого будут «считываться» данные. После измерения каждого параметра колец с помощью программного обеспечения COFECHA [Holmes, 1983] проводился контроль измерений и перекрестной датировки, поиск выпадающих и ложных колец [Fritts, 1966]. Затем для каждого параметра кольца были построены сводные хронологии путем деления годовых значений на соответствующее значение аппроксимирующей кривой (отрицательная экспонента) для полной серии измерений. Стандартные хронологии, построенные таким способом, были затем использованы для статистического анализа и определения их чувствительности со среднемесячными температурами и осадками. Вся обработка древеснокольцевых серий была проведена в дендрохронологическом пакете программ DplR в среде R [Bunn, 2008]. Качество хронологий оценивалось с помощью выраженного популяционного сигнала (EPS) [Briffa, Jones, 1990]. Значения EPS > 0,85 использовались в качестве порога, представляющего надежную и достоверную часть хронологий [Wigley, Briffa, Jones, 1984]. Климатические данные представлены рядами инструментальных наблюдений за среднемесячной температурой и осадками сеточного архива CRU TS 4.01 [Harris et al., 2014] из узла точки (64°75' с.ш., 35°75' в.д.), ближайшего к месту наших работ. Для выявления тесноты связи между приростом сосны и метеорологическими параметрами была использована климатическая функция отклика, а именно коэффициенты множественной линейной регрессии индексов древесно-кольцевой хронологии на главные компоненты месячных климатических данных [Zang, Biondi, 2015]. Статистическая значимость коэффициентов регрессии оценена бутстреп-методом [Guiot, 1991]. Этот статистический метод, состоящий в многократном расчете интересующей статистики на основе множества выборок, сформированных по исходным данным, позволяет правильно оценить значимость полученных значений статистики в условиях сложной структуры данных. Хронологии сравнивались с климатическими данными за 15 месяцев (с июня предыдущего года по сентябрь текущего года). Помимо этого был проведен анализ устойчивости сигнала во времени для каждого параметра. Для этого климатическая функция отклика рассчитывалась в плавающем 35-летнем окне с перекрытием в 2 года. В этой работе дендроклиматический анализ был проведен в дендрохронологическом пакете TreeClim [Zang, Biondi, 2015], реализованном в среде R. Также мы использовали функцию seascoorr для расчета сезонных корреляций между климатическими параметрами и хронологиями разных параметров годичных колец. Результаты На рис. 2 показаны четыре хронологии, построенные по описанной выше методике. Длина хронологий составляет 1841-2016 гг., EPS > 0,85 c 1889 г. На рис. 2 видно, что все четыре хронологии согласуются друг с другом - минимумы и максимумы на них преимущественно отмечаются в одни и те же годы, хотя есть и исключения (например, 1901 и 1936 гг.). Наибольшую изменчивость показывает ширина поздней древесины, амплитуда изменчивости ширины ранней древесины и оптической плотности существенно меньше. Закономерно высокие коэффициенты корреляции (r = 0,95, P < 0,05) наблюдаются между шириной годичного кольца и шириной ранней древесины, которая составляет существенную его часть. Все параметры ширины кольца имеют высокие корреляции между собой (r = 0,59-0,95), тогда как оптическая плотность менее тесно коррелирует с ними (r = 0,370,42) (рис. 3). Рис. 2. Стандартные хронологии ширины годичных колец (TW) сосны, ранней (EW) и поздней (LW) древесины и оптической плотности поздней древесины (BI). Общий период хронологий: 1841-2016 гг., красная вертикальная линия показывает EPS > 0,85 c 1889 г. Fig. 2. Standard chronologies of tree-ring width (TW) pine, earlywood (EW) and latewood (LW) and blue intensity (BI) the common period: 1841-2016, red vertical line shows EPS > 0.85 from 1889 Рис. 3. Коэффициенты корреляции Пирсона, рассчитанные для хронологий различных параметров годичных колец сосны за общий период 1841-2016 гг. (P < 0,05) Fig. 3. Correlations (Pearson's coefficient) between tree-ring width (TW), earlywood (EW), latewood (LW) and blue intensity (BI) chronologies of Pinus sylvestris for the common period 1841-2016 (the significance level of P < 0.05) Невысокие коэффициенты корреляции между оптической плотностью и параметрами ширины годичных колец косвенно свидетельствуют о влиянии разных климатических факторов на формирование этих параметров. На рис. 4 показана функция отклика различных параметров годичных колец сосны на среднемесячные суммы осадков. Ожидаемо связь с осадками слабая, в основном статистически незначимая, поскольку на северном пределе произрастания лесов увлажнение редко бывает недостаточным и не является фактором, лимитирующем прирост. Установлено, что связи с осадками предыдущего года статистически незначимы. Самые высокие коэффициенты обнаружены у ширины колец поздней древесины с осадками июля (r = 0,33). Ранняя древесина показала значимые положительные коэффициенты для февраля, мая и июня, а ширина кольца - для февраля, мая, июня, июля и сентября. Для оптической плотности обнаружены отрицательные коэффициенты регрессии с осадками июля и августа (-0,16 и -0,19 соответственно). На рис. 5 показана функция отклика разных параметров годичных колец на среднемесячные температуры. Для оптической плотности обнаружены значимые положительные коэффициенты с температурой теплого периода (апрель-август). У ширины годичных колец и ранней древесины фиксируются значимые отрицательные связи с температурой августа прошлого года (-0,26 и -0,23). На прирост поздней древесины отрицательно влияет температура июля и августа прошлого года (-0,16 и -0,23), а положительный отклик отмечается на температуру марта, мая и июля. Таким образом, анализ функции отклика показал, что даже в довольно суровых условиях вблизи северной границы леса у сосны на Соловецких островах наблюдается климатический сигнал не только на температуру воздуха, но и на осадки. Тем не менее устойчивость климатической функции во времени показывает, что оптическая плотность - единственный параметр из исследованных нами, для которого сигнал на температуру устойчив на всем протяжении метеонаблюдений с 1901 по 2016 г., он отражает зависимость плотности от температуры всего теплого сезона с апреля по август (см. рис. 6). Это означает, что оптическая плотность поздней древесины является наиболее подходящим параметром для реконструкции температуры на Соловецком архипелаге. Связь с осадками июля и августа у плотности также устойчива, но отрицательна. Рис. 4. Коэффициент регрессии между хронологиями по ширине колец (TW), ранней (EW), поздней (LW) древесины, оптической плотностью поздней древесины (BI) и среднемесячными осадками для общего периода 1901-2016 гг. Здесь и на рис. 5: строчными буквами обозначены месяцы предыдущего года, прописными - текущего. Красной пунктирной линией показан уровень значимости P < 0,05 Fig. 4. Regression coefficient between tree- ring width (TW), earlywood (EW), latewood (LW) and blue intensity (BI) chronologies of Pinus sylvestris with total monthly precipitation for the period 1901-2016 Here and in the fig. 5: months of the previous year are indicated in lowercase letters, and the current year in uppercase letters. The red dotted line shows the significance level of P < 0.05 Рис. 5. Коэффициент регрессии между хронологиями по ширине колец (TW), ранней (EW), поздней (LW) древесины, оптической плотностью поздней древесины (BI) и среднемесячными температурами для общего периода 1901-2016 гг. Fig. 5. Regression coefficient between tree-ring width (TW), earlywood (EW), latewood (LW) and blue intensity (BI) chronologies of Pinus sylvestris chronologies with mean monthly temperature for the period 1901-2016 dB75s_BI_R»sporise function analysisjn ■■ B75s_LW_Response function analyws_m ТМГЯ ТДОг|Ц* Рис. 6. Устойчивость климатической функции отклика во времени за период 1901-2016 гг. a - ширина колец (TW), b - ранняя (EW), c - поздняя (LW) древесина, d - оптическая плотность поздней древесины (BI). Звездочками показаны значимые коэффициенты (P < 0,05). Коэффициенты рассчитывались в плавающем 35-летнем окне с заданным перекрытием (2 года) cocf I. 00 ■ Fig. 6. Stability of the climatic response function over time for the period 1901-2016 a - ring width (TW), b - earlywood (EW), c - latewood (LW), d - blue intensity (BI) chronologies. Asterisks indicate significant coefficients (P < 0.05). The coefficients were calculated in a floating 35-year window with a given overlap (2 years) Обсуждение В рамках нашего исследования для поиска подходящего параметра также была рассчитана хронология по скорректированной поздней древесине (LWW) [Meko, Baisan, 2001]. Meko и соавт. [2001] в своем исследовании показали, что, убрав линейную зависимость ранней древесины от поздней, можно добиться усиления климатического сигнала и тем самым улучшить палеоклиматическую реконструкцию. Рассчитав климатическую функцию отклика у данного параметра, мы не обнаружили значительных отличий от показателей хронологии поздней древесины, поэтому этот параметр был исключен из дальнейшего анализа. При анализе климатической функции отклика хронологии ширины поздней древесины мы обнаружили комплексный характер сигнала. На рис. 7 показаны коэффициенты корреляции Пирсона между хронологией поздней древесины и климатическими переменными с помощью программы seascorr. Эта программа позволяет разделить первичные и вторичные климатические факторы, влияющие на прирост, что особенно важно при скоррелированно-сти климатических параметров. Значимое влияние осадков июля на ширину поздней древесины наблюдается при использовании осадков и как первичного, и как вторичного фактора. Примечателен тот факт, что если первичной переменной назначается температура, а вторичной - осадки, сигнал осадков июля становится сильнее. Таким образом, значимая частная корреляция показывает, что хронологии поздней древесины имеют потенциал для реконструкции осадков на Соловецком архипелаге, но температурный сигнал частично маскирует отклик на осадки. Таким же образом с помощью функции seascorr были рассчитаны коэффициенты корреляции за 1, 3 и 6 месяцев для хронологии по оптической плотности поздней древесины (рис. 8). На рисунке 8 видно, что коэффициент корреляции с первичным климатическим параметром (температура) доходит до значения 0,6, тем самым показывая высокий потенциал оптической плотности для реконструкции летней температуры на Соловецком архипелаге. Рис. 7. Корреляции и частные корреляции между хронологией поздней древесины и сезонными климатическими переменными a - первичная климатическая переменная - осадки (Р), частные корреляции (снизу) со вторичной климатической переменной -температурой (Т); b - первичная климатическая переменная - температура (Т), частные корреляции (снизу) со вторичной климатической переменной - осадки (Р) Fig. 7. Correlations and partial correlations latewood chronology (LW) with seasonalized climate variables a - primary climate variable - precipitation (P), partial correlations (bottom) with a secondary climatic variable - temperature (T); b -primary climatic variable - temperature (T), partial correlations (from below) with a secondary climatic variable - precipitation (P) Рис. 8. Корреляции и частные корреляции хронологии по оптической плотности поздней древесины с сезонными климатическими переменными (Вверху) простые корреляции с первичной климатической переменной, Т (температура), (снизу) частные корреляции со вторичной климатической переменной, Р (осадки) Fig. 8. Correlations and partial correlations blue intensity chronology with seasonalized climate variables (Top) Simple correlations with the primary climate variable, T (temperature), partial correlations of tree-ring index with the secondary climate variable, P (precipitation) Этот анализ (см. рис. 8) показывает статистически значимую связь с летними температурами вегетационного периода, причем наиболее сильные корреляции наблюдаются в июле и августе. Значительные сезонные (3-месячные) корреляции также последовательно выявляются для лета. Максимальная корреляция (r = 0,61) наблюдается для 6-месячного периода апрель-сентябрь. Возможно, для 4- и 5месячных периодов значения коэффициентов корреляции будут еще выше. Результаты анализа отклика разных параметров годичных колец сосны на Соловецком архипелаге не противоречат полученным ранее выводам в соседних регионах. Так, например, в работе Кононова [2018] показано, что на Кольском полуострове (Хибины) ширина годичных колец сосны реагирует на температуру июня-июля, а также на осредненные значения летних месяцев. В работе Ваганова и соавт. [1996] для обширной территории (от Полярного Урала на западе до Чукотки на востоке) был выявлен сильный климатический сигнал в клеточных хронологиях и хронологиях по ширине годичных колец с температурой летних месяцев (июнь-июль). Эти данные получены по разным видам лиственниц, ели сибирской и сосне обыкновенной. Duthorn и соавт. [Duthorn et al., 2016] исследовали климатический сигнал в хронологиях по ширине годичных колец и максимальной плотности сосны вдоль бореального градиента от 60 до 69 с.ш. Интересно, что климатический отклик хронологий MXD схож на протяжение всей трансекты и связан с летними температурами, а на ширину годичных колец сильно влияет место произрастания деревьев. Helama и соавт. [Helama et al., 2013] анализировали связь между климатическими данными и различными параметрами (ширина, ранняя, поздняя древесина и оптическая плотность) в прибрежных (влажных) и нагорных (сухих) местообитаниях в Финляндии. Авторы показали, что оптическая плотность, в отличие от других параметров, демонстрирует сильную связь с теплым периодом (апрель-сентябрь). В Фенноскандии Grudd и соавт. [2002] демонстрируют высокую корреляцию ширины годичных колец с летними месяцами (июнь-август). Важным результатом этой работы является доказательство возможности использования для перекрестного датирования субфоссированной древесины, залегающей в озерных отложениях, или плавника, лежащего на морских террасах. Это важно для удлинения хронологий годичных колец, как это сделано, например, в Фенноскандии и на севере Евразии [Grudd et al., 2002; Hellmann еt al., 2016]. Продление хронологий по живым деревьям возможно и с помощью древесины, использованной в архитектурных сооружениях. И в этом случае применим метод оптической плотности, хотя при его использовании возникает определенная специфика, связанная с различной интенсивностью синего цвета между заболонью и сердцевиной, а также между образцами живой и погребенной древесины [Bjorklund et al., 2014]. Этот факт для нас важен, так как создание длинной древесно-кольцевой хронологии на Со ловецких островах предполагает включение образцов архитектурной и археологической древесины. В работе Bjorklund и соавт. [2014] было показано, что именно ABI (А = разница (контраст) между плотностью поздней и ранней древесины) сопоставима с MXD и пригодна для использования в денд-роклиматологических реконструкциях. Выводы Результаты нашего исследования показывают: 1. Различные параметры годичных колец реагируют на различные климатические факторы. 2. Ширина ранней древесины и ширина годичного кольца положительно коррелируют с осадками февраля, мая, июня и отрицательно - с температурой летних месяцев прошлого года. 3. Ширина поздней древесины реагирует на осадки июля текущего года и температуры марта, мая и июля. Этот параметр может быть использован для реконструкции осадков июля на Соловецких островах. 4. Климатическая функция отклика оптической плотности поздней древесины существенно отличается от других параметров высокими коэффициентами корреляции с температурой, а также более продолжительным периодом: она зависит от температуры (апрель-август). Эта корреляция, кроме того, более устойчива во времени. Благодаря этому именно оптическая плотность поздней древесины сосны обыкновенной является подходящим параметром для создания палеоклиматической реконструкции на Соловецком архипелаге и в близких по физико-географически и климатическим условиям районах.

Скачать электронную версию публикации