Исследование природных ритмов хвойных деревьев

Введение

Актуальность проблемы

С середины XIX века средняя глобальная температура Земли увеличилась на 0,75 ± 0,2 °C градуса (IPCC, 2007).

Причина этого явления во многом до сих пор не ясна. Часть учёных полагает, что сегодняшнее потепление - это закономерное явление и оно связано с естественной цикличностью процессов. Поскольку климатические изменения подобного масштаба наблюдались и ранее, не вызывая серьезных изменений на уровне экосистем, нынешнее потепление можно рассматривать с точки зрения естественной изменчивости климата Земли.

Противники этой теории считают, что глобальное потепление - это уникальное явление, которое может принести непредсказуемые последствия и серьёзно изменить климат. Они утверждают, что потепление происходит в связи с существенным увеличением концентрации парниковых газов в атмосфере (углекислота, метан, закись углерода), которое, в свою очередь, вызвано человеческой деятельностью (антропогенный эффект). Предполагается, что эти повышение температуры будет не только превышать естественные природные колебания, но и происходить в 15-40 раз быстрее, чем аналогичные изменения в прошлом (IPCC, 2007).

Согласно разным сценариям выбросов парниковых газов, удвоение концентрации эквивалента углекислоты в атмосфере по сравнению с доиндустриальным уровнем произойдет между 2025-2050 гг., что изменит глобальную среднюю температуру воздуха от 1,5 до 4,5°С. Средняя скорость увеличения глобальной средней температуры в следующем столетии может составить примерно 0,3°С за десятилетие. Поверхность суши нагреется быстрее по сравнению с поверхностью океана.

Ожидается, что это потепление будет неравномерным: в тропических районах оно составит половину глобального среднего, а наибольшим будет в полярных районах, где превысит среднюю глобальную величину примерно в 2 раза и в основном произойдет в зимние месяцы (IPCC, 2007).

Рис.1. Тренд температурной изменчивости поверхности Земли в градусах Цельсия за десятилетие: вверху - начиная с 1870г., внизу - с 1950г (Jones et. al., 2001)

Возможности прогноза естественной и антропогенной изменчивости климата, часто ограничивается в связи с короткими рядами метеорологических данных. Как правило, их длина не превышает 30-50 лет, за редкими исключениями - 100 лет. Косвенные индикаторы климата - ежегодный прирост древесины (ширина колец), толщина снежного покрова на ледниках, кораллы и прочие - могут в какой-то мере заменить инструментальные наблюдения в прошлом и продолжить ряды. Дендрохронологический метод в этом отношении является одним из самых перспективных: эти данные имеют годовое и сезонное разрешение, поддаются калибровке с помощью инструментальных метеорологических данных и потенциально доступны во всех районах, где растут деревья.

Для прогнозирования климата необходимо иметь информацию не только и не столько о глобальном, но в большей степени - о региональном климате. Современные модели общей циркуляции атмосферы не могут быть использованы для моделирования региональных климатических изменений, а обеспеченность многих районов длительными климатическими наблюдениями явно недостаточна. Изучение климатической изменчивости на европейской части России с помощью косвенных методов индикации климата (дендрохронология) является, таким образом, одной из актуальных задач изучения климатической изменчивости. Это важно не только для реконструкции прошлых изменений, но и для оценки изменений климата в будущем.

Цель работы

Основная цель работы - при помощи вейвлет-анализа выявить структуру природных ритмов, отраженную в ширине и плотности годичных колец хвойных деревьев на территории севера европейской части России.

Для достижения цели необходимо:

1.собрать и структурировать данные по приросту древесины на севере Европейской части России и построить сводные хронологии ширины и плотности древесины

2.провести корреляционный анализ между площадками для объединения их в региональные хронологии

.провести корреляционный анализ хронологий с метеорологическими параметрами и выявить господствующий климатический сигнал, влияющий на прирост древесины

.с помощью вейвлет-анализа выявить периодичность в хронологиях и сравнить эти результаты с результатами, полученными методами преобразований Фурье и СВАН-диаграмм

.выявить сходство и различия ритмики в рядах ширины и плотности годичных колец древесины, данных метеорологических наблюдений и климатических индексов НАО

.оценить пространственную и временную изменчивость природных ритмов за последние 300 лет на субрегиональном уровне

Стратегия исследования

На схеме ниже схематично показана последовательность действий, которой мы придерживались в этой работе.

Хронологии, создаваемые до 1990х гг. на европейской части России и на прилегающих территориях за немногими исключениями (Шиятов, 1986), не подвергались хронологическому контролю (перекрестное датирование) (Fritts, 1976, Шиятов и др., 2002) и поэтому не могут использоваться для дальнейшего анализа. Хронологии, создаваемые на основе археологической древесины, не сведены в региональные, недоступны широкому пользователю, и используются лишь для датирования археологических и архитектурных объектов (Черных, 1989).

В 1990-х гг. благодаря исследованиям международной группы специалистов (Schweingruber, Briffa, 1996; Ваганов и др., 1996; Briffa et al., 2001) удалось получить большой массив дендрохронологической информации по ширине и плотности колец на северной границе леса в Евразии, эти данные доступны и хранятся в международном банке данных (www.cdc.noaa.gov). Для юга Русской равнины таких данных нет, поскольку в этом районе крайне трудно найти высоковозрастные деревья и погребенную древесину. В окружающих Русскую равнину горах построены древесно-кольцевые хронологии для Полярного и Приполярного Урала (Шиятов, 1986) и Кольского полуострова (Hiller et al., 2001, Шумилов и др., 2007).

Многочисленные древесно-кольцевые хронологии лиственницы (Larix sibirica, L.dahurica) и ели (Picea abies, P.obovata) вдоль северной границы леса между 42°E и 152°E, за последние 400-600 лет показывают хорошую взаимную согласованность (Schweingruber, Briffa, 1996; Vaganov et al., 1996; Briffa et al., 2001; 2004), поэтому они были сведены в несколько региональных хронологий. В частности, одна из них представляет европейский сектор (хронология NEUR) и включает (помимо российских) скандинавские хронологии (Briffa et al., 2001). Как ширина кольца, так и максимальная плотность древесины этих хронологий положительно коррелируют с температурой теплого периода (апрель-сентябрь r~0.7). Однако, эти реконстркукции были ограничены 1960-м годом, поскольку для более позднего периода был обнаружен эффект рассогласования дендрохронологических и климатических показателей.

Для Урала (Шиятов, 1986) и Кольского полуострова (Hiller et al., 2001) реконструкции температур теплого периода охватывают гораздо более длительный период, поскольку хронологии включают погребенную древесину. Более того, для реконструкции долгопериодных трендов в обоих районах используются также данные по динамике верхней границы леса. Работы на Кольском полуострове позволили создать несколько длительных хронологий (Hiller et al., 2001). Хронология показывает значимую положительную корреляцию с солнечной активностью и лишь невысокие положительные коэффициенты (до 0,25) с температурами с августа по ноябрь. Сходые результаты были получены и Распоповым с соавторами для сосны для этого же района (Raspopov et al. (2004).

Наиболее систематические и исчерпывающие по интересующему нас вопросу были получены К. Бриффой с соавторами (Briffa et al., 2004), Эти авторы хронологии плотности древесины осредняли по сетке 5х5 градусов и коррелировали со среднемесячной температурой осредненной для того же квадрата. Функции отклика для каждой локальной хронологии были рассчитаны с опорой на данные о среднемесячных температурах (Jones et al., 1999) и осадках (New et al., 2000) по регулярной сетке, для предшествующего формированию кольца и текущего года (за период 1881-1984 гг.). Анализ показал (рис.2) (Briffa et al., 2004), что наиболее тесные корреляции в целом обнаруживаются для максимальной плотности годичных колец и температур в период с апреля по сентябрь.

Рис.2. Функция отклика для максимальной плотности древесины хвойных на северной границе леса (Briffa et al., 2004)

Реконструкции осадков для северных районов по годичным кольцам деревьев крайне редки, поскольку лимитирующим фактором в этих условиях является в подавляющем большинстве случаев летняя температура. Оригинальную попытку реконструировать осадки на основе дендрохронологических данных Швейнгрубера на Полярном Урале предпринял Ю.М. Кононов с соавторами (Kononov, Ananicheva, Willis, 2005). Они в частности, показали, что ширина колец ранней древесины положительно коррелирует с осадками холодного периода, однако насколько устойчива эта связь пока остается неясным из-за ограниченности привлеченного Кононова с соавторами материала и неясного экологического обоснования такой связи.

Таким образом, несмотря на подробный анализ данных дендрохронологии на севере Евразии, проведенный нашими предшественниками (Ваганов и др,, Briffa et al., 2004), многие вопросы, особенно касающиеся региональной климатической изменчивости в этом районе, остались неисследованными. В частности, неизвестно насколько отличалась динамика ширины и плотности колец с запада на восток - на Кольском полуострове, в центральной части севера Русской равнины и на Урале, каковы основные ритмы этой изменчивости и насколько они устойчивы во времени за последние 300 лет, действительно ли по ширине ранней древесины возможно реконструировать осадки холодного периода в этом регионе и др. Эти проблемы будут, в частности, рассмотрены в этом исследовании.

Глава 1. Основы дендрохронологии

Дендрохронология занимается датировкой годичных слоев прироста древесины и связанных с ними событий, изучением влияния экологических факторов на величину прироста древесины, а также анализом содержащейся в годичных слоях информации для целей реконструкции условий окружающей среды. Особую значимость древесно-кольцевая информация имеет для решения глобальных, региональных и локальных проблем, связанных с многообразным воздействием человека на наземные экосистемы. Уникальность дендрохронологических методов состоит в том, что они позволяют оценивать относительный вклад различных факторов, как естественных, так и антропогенных, которые оказывают влияние на изменение и трансформацию лесных экосистем и условий окружающей среды, а так же высокое разрешение - до года или даже сезона в случае, если удается найти деревья, достаточно чувствительные к изменениям климата. Но необходимо учитывать два важных фактора: во-первых, реконструкции поддается обычно один из климатических параметров - в холодных районах это, как правило, температура теплого периода, в засушливых - осадки. Крайне редко удается получить для одного района обе характеристики. Обычно это возможно в горах, где прирост древесины на верхней границе леса обычно ограничен температурой, а на нижней - осадками. В этой работе мы предпринимаем попытки реконструировать два параметра и ниже обсуждаем разную степень достоверности реконструкции температур и осадков для севера европейской части России.

Второе обстоятельство состоит в том, что способность дендрохронологических реконструкций воспроизводить долгопериодные колебания климата ограничена в связи с тем, что при создании сводных хронологий приходится проводить индексирование серий, необходимое для вычленения возрастного тренда прироста (Шиятов и др., 2000). Это, однако, не мешает впоследствии оценить на основе древесно-кольцевых хронологий колебания климатических параметров межгодовой и внутривековой продолжительности, в частности, частоту экстремумов, периодичность колебаний, тренды и некоторые другие статистические характеристики.

Методы дендрохронологии в реконструкции климата

Древесно-кольцевая хронология представляет собой дискретный временной ряд длительностью от нескольких лет до многих тысячелетий, характеризующий тот или иной показатель годичного прироста, физико-механические свойства, анатомическую структуру и химический состав древесины.

В зависимости от целей исследования используются различные виды древесно-кольцевых хронологий. При датировке исторической и археологической древесины, а также различного рода катастрофических событий в лесу (пожары, ветровалы, вспышки массового размножения насекомых-вредителей и др.) в основном используются индивидуальные хронологии, т.е. хронологии, полученные с одного дерева или кустарника. Для реконструкции климатических и гидрологических условий используются обобщенные и генерализированные хронологии, т.е. хронологии, построенные путем использования достаточно большого числа деревьев, принадлежащих к одному или нескольким видам и произрастающих в одном или нескольких типах местообитаний.

Основные принципы дендрохронологии

Закон лимитирующих факторов: рост древесных растений не может протекать быстрее, чем это позволяется внешним или внутренним фактором, находящимся в минимуме. В случае, если этот фактор переходит в разряд оптимальных, скорость роста будет увеличиваться до тех пор, пока другой фактор не станет лимитирующим. Согласно этому закону, для дендрохронологического анализа наиболее пригодны те деревья, на прирост которых оказывает влияние тот или иной лимитирующий фактор, в предельном случае - только один.

Принцип отбора районов и местообитаний: при проведении дендрохронологических исследований образцы древесины должны браться у деревьев, которые произрастают в неблагоприятных и экстремальных климатических условиях, где наиболее полно проявляется действие лимитирующих факторов.

Принцип чувствительности: в благоприятных для роста деревьев районах и местообитаниях формируются широкие годичные кольца. При этом у таких деревьев хорошо выражены изменения прироста с возрастом, а величина прироста между соседними годами колеблется в незначительных пределах. Подобная последовательность называется «благодушной». В неблагоприятных для произрастания деревьев условиях кольца прироста узкие, их ширина значительно колеблется от года к году, возрастная кривая роста выражена слабо. Часто наблюдается выпадение колец. Такие серии колец называются «чувствительными». Они свидетельствуют о том, что на рост деревьев большое влияние оказывают факторы внешней среды.

Принцип перекрестного датирования: древесные растения, произрастающие в пределах однородного в климатическом отношении района, величиной прироста сходно реагируют на изменения лимитирующих климатических факторов. В благоприятные по климатическим условиям годы у большей части деревьев формируются широкие кольца, а в неблагоприятные - узкие. В связи с этим у таких деревьев наблюдается синхронная изменчивость величины прироста во времени. Особенно показательны узкие кольца, когда прирост в наибольшей степени лимитируется тем или иным климатическим фактором. Метод перекрестного датирования позволяет выявлять местонахождение ложных и выпавших колец и производить абсолютную и относительную датировку каждого кольца с точностью до года у сравниваемых индивидуальных древесно-кольцевых хронологий, а также продлевать хронологии далеко в глубь веков на основе использования древесины давно отмерших деревьев.

Принцип повторности: использование информации не с одного, а с определенного числа модельных деревьев, является непременным условием точной датировки колец, построения надежных древесно-кольцевых хронологий и производства более точной реконструкции условий среды в настоящем и прошлом.

Принцип униформизма (актуализма): физические и биологические процессы, обусловливающие изменения в росте дерева под воздействием факторов окружающей среды в настоящее время, вызывали подобные же изменения в прошлом. Этот принцип является обоснованием широкого использования древесно-кольцевых хронологий для реконструкции прошлых условий окружающей среды.

Достоинства и недостатки метода

Как и любой метод, дендрохронология обладает рядом достоинств и недостатков.

Достоинства дендрохронологических методов:

·Высокая разрешающая способность древесно-кольцевых хронологий (год и сезон года).

·Возможность абсолютной и относительной датировки времени формирования годичных колец, как у живых, так и давно отмерших деревьев.

·Возможность получения длительных и непрерывных хронологий (сотни и тысячи лет).

·Возможность выявлять колебания различной длительности в изменчивости различных характеристик прироста деревьев и определяющих прирост факторов внешней среды.

·Возможность получения информации для огромных территорий суши, где растут или росли в прошлом древесные растения.

Недостатки дендрохронологических методов:

·Невозможность выявления колебаний и трендов в хронологиях, длительность которых превышает возраст использованных деревьев.

·Сильная зависимость числа исследуемых деревьев на точность датировки годичных колец и надежность реконструкции параметров внешней среды.

·Различия в реакции прироста дерева на одни и те же факторы внешней среды и в случае резкого и быстрого изменения климатических и почвенно-грунтовых условий.

·Сходная реакция прироста древесины на различные факторы среды.

Процесс дендрохронологических исследований включает в себя ряд последовательных действий: в первую очередь необходимо выбрать объект исследования и отобрать древесину, затем подготовить образцы для анализа, а после провести саму датировку и измерение характеристик древесных колец.

Измерение характеристик древесины

Существует два подхода к проведению работ по окончательной датировке и измерению годичных колец. В основе обоих подходов лежит один и тот же метод - метод перекрестной датировки, основанный на использовании неповторимого во времени рисунка годичных колец древесины. Перекрестная датировка - это сравнение рисунков колец у различных деревьев и выбор точного места, где найдено соответствие в характере изменчивости показателей радиального прироста между рассматриваемыми образцами. Этот метод позволяет выявлять нарушение синхронности в пределах отдельных временных интервалов и, тем самым, выявлять точное положение ложных и выпавших колец.

При первом подходе, используемом для датировки и измерения древесных колец, сначала производится измерение колец у предварительно сдатированных и маркированных индивидуальных хронологий. На основе этих измерений строятся графики изменения прироста для каждого радиуса, которые затем визуально сопоставляются друг с другом на предмет выявления сходства и различий между ними, определения местоположения выпавших, ложных, пропущенных и лишних колец.

После того как для всех образцов будут получены графики изменения радиального прироста, производится визуальная оценка синхронности между всеми ними и окончательная абсолютная или относительная датировка колец у каждой индивидуальной древесно-кольцевой хронологии.

При втором подходе сначала производится датировка древесных колец, а затем их измерение. Датировка основана на глазомерном изучении относительной изменчивости годичного прироста древесины непосредственно на образце, в частности, на использовании экстремальных величин прироста и таких патологических образований, как морозобойные, светлые и ложные кольца и др., в качестве реперных или указательных лет.

Был разработан оригинальный метод построения так называемых скелетных графиков. Суть его состоит в том, что для каждой маркированной серии колец строится специальный график, на котором отмечаются годы с минимальным (линия вниз) и максимальным (линия вверх) приростом но отношению к соседним кольцам. Средние по ширине кольца на графике не отмечаются.

Затем скелетные графики, построенные для каждого образца древесины, сравниваются между собой. По индивидуальным хронологиям строится обобщенная (или усредненная) хронология на основе вычисления средней длины вертикальных линий за отдельные годы. Преимуществом скелетных графиков является то, что нет необходимости производить измерение величины прироста. Этот метод широко использовался для датировки исторической и археологической древесины.

Метод перекрестной датировки позволяет не только датировать кольца у ныне живущих деревьев, но и продлевать хронологии далеко в глубь веков. Длительные хронологии важны для датировки различных природных явлений и культурных памятников, а также для реконструкции динамики лесных экосистем и условий среды в далеком прошлом с высоким временным разрешением.

Рис.3. Пример построения древесно-кольцевой хронологии длительностью 1950 лет на основе использования образцов древесины, взятых с живого дерева и двух деревьев, отмерших в 1228 и 897 годах н.э. соответственно.

При сборе образцов древесины и построении длительных хронологий важно учитывать то, что перекрестно можно датировать лишь местную древесину, которая взята с остатков деревьев, произраставших в пределах однородного климатического района.

Дендроклиматические материалы и методы

В этой работе использовались результаты измерений ширины ранней древесины и плотности колец отдельных образцов из международного банка дендрохронологических данных (www.cdc.noaa.gov).

Ширина ранней древесины

Плотность древесины

Рис.4. Наглядное представление ранней ширины и плотности древесины

Процедура дендроклиматического анализа состоит в следующем. Образцы древесины (керны и спилы) отбираются на некоторой ограниченной территории с однородными орографическими, почвенными и микроклиматическими условиями. Деревья на площадках выбираются наиболее старые, отдельно стоящие, не угнетённые, не повреждённые, так как основной фактор, влияющий на прирост таких деревьев, может быть неклиматическим (Schweingruber, 1988; Cook, 1990). Для уменьшения вероятности выпавших или ложных колец в кернах, с одного ствола отбираются обычно по два керна по двум произвольным радиусам на высоте 1-1,5м от поверхности земли.

Дальнейшая работа с образцами древесины проводится в лабораторных условиях. Керны наклеиваются на специальную деревянную основу, после чего керн полируют для увеличения изображения контрастности колец (рис.5). Затем, зная точную календарную дату отбора образца, проводится разметка и предварительная датировка колец (Шиятов и др., 2000).

Рис.5. Образцы, подготовленные для дендрохронологического анализа

С помощью полуавтоматического комплекса состоящего из микроскопа и подвижного измерительного стола, измерялась ширина годичных колец с точностью до ±0,01 мм. Величина ширины колец автоматически фиксируется компьютерной программой TSAP (Time Series & Analysis Presentation - программа анализа и перекрестной датировки дендрохронологических серий и графического представления результатов обработки данных), которая помогает решить широкий спектр задач, связанных с анализом древесных кернов и спилов.

Метод перекрёстного датирования позволяет с высокой точностью определить календарную дату выпадающих или ложных колец. Кроме того, вероятность выпадения колец увеличивается с возрастом. Программа COFECHA позволяет произвести контроль качества полученных древесно-кольцевых серий. Для проверки правильности датировки программа COFECHA преобразовывает серии измерений с помощью короткого кубического сплайна, а затем сравнивает все полученные серии и выявляет проблематичные участки (Holmes, 1983). Каждая индивидуальная серия делится на сегменты по 50 лет, каждый сегмент имеет общий интервал продолжительностью 25 лет с соседним сегментом и т.д. Значения коэффициентов корреляции определяются программой для каждого сегмента серии, что наиболее точно может указать на ошибку датирования. Статистически значимым считается коэффициент корреляции, значение которого более, чем 0,3281. Кроме того, это программа позволяет определить абсолютный возраст погребённой древесины и образцов сухих деревьв относительно имеющейся хронологии.

Рис.6..Значения коэффициентов корреляции некоторых кернов со всей хронологией Кавказа

На рисунке 6 наглядно показано, как выглядят в программе результаты. Слева столбец с условными обозначениями кернов, а справа значения коэффициентов корреляции данного образца со всей хронологией с временным разбиением 50 лет.

Поскольку абсолютное значение прироста древесины зависит от различных факторов, таких как, возрастные изменения, конкурентные взаимоотношения, катастрофические явления и т.д., для сопоставления изменчивости прироста древесины необходимо перевести абсолютные значения ширины и плотности колец в относительные. Исключения или явного снижения влияния перечисленных факторов можно добиться с помощью стандартизации или индексирования (Fritts, 1976) (рис.7). Обычно для индексирования для каждого образца подбирают индивидуальную биологическую кривую роста и индексы рассчитываются путем деления или вычитания значения ширины (плотности) кольца из соответствующего значения аппроксимирующей кривой.

Рис.7. Пример индексирования - избавления от возрастного тренда

После того, как хронологии построены, требуется выяснить, от каких факторов зависит ширина или плотность кольца. Для этого проводится корреляция дендрохронологических данных со среднемесячными осадками и температурами за один-два года, предшествующих сезону прироста и строится так называемая функция отклика. Она может быть построена как для локальных хронологий, значения которых обычно коррелируют с ближайшей метеостанцией, так и для сводной - региональной. В этом случае иногда усредняют и значения метеорологических параметров. Если функции отклика для одного района имеют сходный рисунок, можно говорить о том, что хронологии отражают некий общий климатический сигнал. Для реконструкции выбирают месяцы, для которых связи дендрохронологических и метеорологических параметров являются статистически значимыми и объяснимыми с экологической точки зрения.

Глава 2. Методы обработки нестационарных сигналов

Сигналы, которые нам надо обработать имеют сложные частотно-временные характеристики, они состоят из близких по времени, короткоживущих высокочастотных компонент и долговременных, близких по частоте низкочастотных компонент. Для анализа таких сигналов нужен метод, способный обеспечить хорошее разрешение и по частоте, и по времени. Первое требуется для локализации низкочастотных составляющих, второе - для разрешения компонент высокой частоты.

Вейвлет преобразование стремительно завоевывает популярность в столь разных областях, как телекоммуникации, компьютерная графика, биология, астрофизика и медицина. Оно применяется при анализе изображений различной природы, для изучения структуры турбулентных полей, для сжатия больших объемов информации, в задачах распознавания образов, при обработке и синтезе сигналов. Благодаря хорошей приспособленности к анализу нестационарных сигналов оно стало мощной альтернативой преобразованию Фурье.

Преобразование Фурье представляет сигнал, заданный во временной области, в виде разложения по ортогональным базисным функциям (синусам и косинусам), выделяя, таким образом, частотные компоненты. Недостаток преобразования Фурье заключается в том, что частотные компоненты не могут быть локализованы во времени, что накладывает ограничения на применимость данного метода к ряду задач (например, в случае изучения динамики изменения частотных параметров сигнала на временном интервале).

Существует несколько подходов для анализа нестационарных сигналов такого типа. Первый - локальное преобразование Фурье (Short-Time Fourier Transform), оно же называется оконным преобразованием. В этом случае, мы работаем с нестационарным сигналом, как со стационарным, предварительно разбив его на сегменты (окна), статистика которых не меняется со временем (в данной работе это преобразование будем называть СВАН - спектрально-временной анализ), но у этого метода также есть ряд небольших недостатков, которые мы обсудим позже. Второй подход - вейвлет преобразование. В этом случае нестационарный сигнал анализируется путем разложения по базисным функциям, полученным из некоторого прототипа путем сжатий, растяжений и сдвигов. Функция-прототип называется материнским, или анализирующим вейвлетом.

Краткий обзор преобразования Фурье

Традиционно для анализа временных рядов используется преобразование Фурье - преобразование функции, превращающее её в совокупность частотных составляющих. Более точно, преобразование Фурье - это интегральное преобразование, которое раскладывает исходную функцию по базисным функциям, в качестве которых выступают синусоидальные (или мнимые экспоненты) функции, то есть представляет исходную функцию в виде интеграла синусоид (мнимых экспонент) различной частоты, амплитуды и фазы. В более общей форме записи:

,

где

- амплитуды гармонических колебаний соответствующей частоты, которые определяются по формуле:

Результат преобразования Фурье - амплитудно-частотный спектр, по которому можно определить присутствие некоторой частоты в исследуемом сигнале. Фурье-преобразования дают достаточно простые для расчетов формулы и прозрачную интерпретацию результатов, но не лишены и некоторых недостатков.

Преобразования:

·не отличают сигнал, являющийся суммой двух синусоид, от ситуации последовательного включения синусоид

·не дают информации о преимущественном распределении частот во времени

·могут дать неверные результаты для сигналов с участками резкого изменения.

Исследуемые ряды также далеко не всегда удовлетворяют требованию периодичности и более того, как правило, заданы на ограниченном отрезке времени.

Таким образом, когда не встает вопрос о локализации временного положения частот, метод Фурье дает хорошие результаты, но при необходимости определить временной интервал присутствия частоты приходится применять другие методы.

Одним из таких методов является обобщенный метод Фурье (локальное (оконное) преобразование Фурье). Этот метод состоит из следующих этапов:

1.в исследуемой функции создается окно - временной интервал, для которого функция f(x) 0, и f(x)=0 для остальных значений,

2.для этого окна вычисляется преобразование Фурье,

.«окно» сдвигается, и для него также вычисляется преобразование Фурье

Пройдя таким «окном» вдоль всего сигнала, получается некоторая трехмерная функция, зависящая от положения «окна» и частоты.

Но данный подход позволяет определить факт присутствия в сигнале любой частоты, и интервал ее присутствия. Это значительно расширяет возможности метода по сравнению с классическим преобразованием Фурье, но существуют и определенные недостатки. Согласно следствиям принципа неопределенности Гейзенберга в данном случае нельзя утверждать факт наличия частоты w0 в сигнале в момент времени t0 - можно лишь, что спектр частот (w1, w2) присутствует в интервале (t1, t2). Причем разрешение по частоте (по времени) остается постоянным вне зависимости от области частот (времен), в которых производится исследование. Поэтому, если, например, в сигнале существенна только высокочастотная составляющая, то увеличить разрешение можно только изменив параметры метода.

В данной работе такой метод мы будем называть СВАН-преобразованием.

В качестве метода, не обладающего подобного рода недостатками, был предложен аппарат вейвлет анализа.

Основные положения вейвлет-преобразования

Альтернатива преобразованию Фурье для исследования временных (пространственных) рядов с выраженной неоднородностью стал метод, разработанный в 80-х годах под названием вейвлет-анализ.

Вейвлеты (от англ. wavelet), всплески - это математические функции, позволяющие анализировать различные частотные компоненты данных.

Вейвлеты - это семейство функций, которые локальны во времени и по частоте («маленькие»), и в которых все функции получаются из одной посредством её сдвигов и растяжений по оси времени (так что они «идут друг за другом»).

Все вейвлет-преобразования рассматривают функцию, взятую, будучи функцией от времени, в терминах колебаний, локализованных по времени и частоте.

Вейвлет-преобразования обычно делят на дискретное вейвлет-преобразование (ДВП) и непрерывное вейвлет-преобразование (НВП). ДВП обычно используется для кодирования сигналов, в то время как НВП для анализа сигналов. В результате, ДВП широко применяется в инженерном деле и компьютерных науках, а НВП в научных исследованиях.

Вейвлет-преобразование, обладающее самонастраивающимся подвижным частотно-временным окном, одинаково хорошо выявляет как низкочастотные, так и высокочастотные характеристики сигнала на разных временных масштабах. По этой причине вейвлет-анализ часто сравнивают с "математическим микроскопом", вскрывающим внутреннюю структуру существенно неоднородных объектов. Семейства анализирующих функций, называемых вейвлетами, применяются при анализе изображений различной природы, для изучения структуры турбулентных полей, для сжатия больших объемов информации, в задачах распознавания образов, при обработке и синтезе сигналов.

Обычно, функция-вейвлет обозначается буквой ?.

Подобно тому, как в основе аппарата преобразований Фурье лежит единственная функция , порождающая ортонормированный базис пространства L²[0,2?] путем масштабного преобразования, так и вейвлет-преобразование строится на основе единственной базисной функции , имеющей солитоноподобный характер и принадлежащей пространству L²(R), т.е. всей числовой оси.

При конструировании базисной анализирующей функции должны выполняться следующие необходимые условия:

·Локализация - вейвлет должен быть локализован вблизи нуля аргумента как во временном, так и в частотном пространстве. Нулевое среднее: .

·Вейвлет должен быть знакопеременной функцией. Ограниченность: .

·Вейвлет должен достаточно быстро убывать по временной (пространственной) переменной.

Наибольшей популярностью пользуются вейвлеты:

Непрерывное вейвлет-преобразование (НВП) строится с помощью непрерывных масштабных преобразований и переносов вейвлета с произвольными значениями масштабного коэффициента a и параметра сдвига b:

,(1)

где - комплексно сопряженная функция, b - момент времени, t-ось времени, a - параметр, обратный частоте.

Вейвлет-преобразование обратимо для функций f из L2(R)

.

Итак, у нас имеется некоторая функция f(t), зависящая от времени. Результатом ее вейвлет-преобразования будет некоторая функция W(a,b), которая зависит уже от двух переменных: от времени и от частоты (обратно пропорционально). Для каждой пары a и b рецепт вычисления вейвлет преобразования следующий:

1.Функция вейвлет растягивается в a раз по горизонтали и в 1/a раз по вертикали.

2.Далее он сдвигается в точку b. Полученный вейвлет обозначается .

.Производится усреднение в окрестности точки a при помощи .

Спектр вейвлет-преобразования одномерного сигнала представляет поверхность в трехмерном пространстве. Обычно изображение спектра выполняется путем проектирования линий постоянного уровня поверхности на плоскость с переменными: параметрами сдвига (по оси абсцисс) и масштабом (по оси ординат), с градиентной заливкой оттенками серого цвета между линиями.

В результате «вырисовывается» наглядная картина, иллюстрирующая частотно-временные характеристики сигнала. По оси абсцисс откладывается время, по оси ординат - частота (иногда log(T), где T=1/a - период). А абсолютное значение вейвлет преобразования для конкретной пары a и b определяет цвет, которым данный результат будет отображен (чем в большей степени та или иная частота присутствует в сигнале в конкретный момент времени, тем более выраженный будет оттенок). На рисунке 8 показан пример вейвлет-анализа. Чётко видно, что на протяжении всего времени действия сигнала преобладают две частоты.

Рис.8. Пример вейвлет-анализа. Внизу - исходный двухчастотный сигнал, вверху - вейвлет-анализ, показывающий присутствие сигналов с частотой примерно 0,8 и 3,5. По оси абсцисс время, по оси ординат - частота

Таким образом, любая функция из L²(R) может быть представлена суперпозицией масштабных преобразований и сдвигов базисного (материнского) вейвлета с коэффициентами, зависящими от масштаба (частоты) и параметра сдвига (времени).

Двухпараметрическая функция W(a,b) дает информацию об изменении относительного вклада компонент разного масштаба во времени и называется спектром коэффициентов вейвлет-преобразования.

Вейвлеты, используемые в данной работе:

Для работы выбрали именно эти вейвлеты, так как вейвлет Морле прост в обращении и довольно чётко отражает результаты, а DOG-вейвлет хорошо показывает моменты перестройки сигнала, то есть переход от одного доминирующего периода в ряде к другому.

Получив вейвлет-спектр, можно рассчитать полную энергию сигнала:

и глобальный спектр энергии - распределение полной энергии по масштабам частоты (скейлограмму вейвлет-преобразования)

Скейлограмма соответствует спектру мощности Фурье-преобразования сигнала, сглаженному на каждом масштабе спектром Фурье анализирующего вейвлета:

,

где знак ^ обозначает Фурье-образ функции.

На практике чаще приходится иметь дело с сигналами, заданными не аналитическими функциями, а с дискретным набором данных, определенном на конечном временном интервале. В этом случае принимается, что при , а формула (1) для коэффициентов вейвлет-преобразования модифицируется следующим образом:

,

где

Дискретное вейвлет-преобразование

Базис одномерного дискретного вейвлет-преобразования (ДВП) строится на основе материнского вейвлета посредством операций сдвигов и растяжений вдоль оси t. Вводя аналог синусоидальной частоты и, принимая для простоты в качестве ее значений степени двойки, получаем для функций базиса .

Базис нормирован, если вейвлет имеет единичную норму.

Вейвлет называется ортогональным, если семейство представляет ортонормированный базис функционального пространства L2(R), т.е. . В этом случае любая функция может быть представлена в виде ряда:

,

Где

Спектрально-временной анализ (СВАН)

Как уже говорилось выше, в данной работе будет так же использоваться спектрально-временной анализ, он же - графическое представление оконного преобразования Фурье, то есть процедура спектрального анализа временных рядов в скользящем временном окне (предварительно исключается полиномиальный тренд, порядок которого соответствует той кратности, с которой заданное временное окно укладывается в полной длине ряда). Длина окна задается исходя из требуемой детальности и спектрального состава процесса (она должна быть кратна периодам исследуемых ритмов). Мы будем задавать длину 25% от длины ряда. Результат изображается в виде спектрально-временной диаграммы. На оси абсцисс откладывается календарное время (в годах), соответствующее центру скользящего временного окна. На оси ординат откладываются величины частот в циклах в единицу времени (год). Каждый вертикальный столбец представляет собой амплитудный спектр Фурье, рассчитанный в заданном скользящем временном окне. Более сильная зачерненность на диаграммах соответствует большим спектральным амплитудам. О величинах амплитуд судят по шкале уровней, помещенной около диаграммы. Повторяемость или устойчивость доминирующих ритмов выражается в виде протяженных зачерненных полос.

Но у этого анализа есть довольно существенный недостаток: в данном анализе чем шире выбираемое окно, тем больше «отрезается» значений от начального ряда после преобразования. Но, чем уже окно, тем хуже определяются частоты, доминирующие в данный момент времени. Например, если длина ряда 300 лет, а ширина окна 50 лет, то в поучившейся диаграмме будет отсутствовать по 25 лет с каждой стороны ряда.

По повторяемости или устойчивости доминирующих ритмов также можно судить о степени упорядоченности процесса.

Для оценки уровня упорядоченности процессов вычисляется т.н. «параметр хаотизации»:

,

где K = L/2+1 и L - длина скользящего временного окна анализа (количество отсчетов), W(w, t) - амплитудный спектр на частоте w и времени t. Нулевое значение параметра соответствует наличию одной спектральной линии (т.е. ряд состоит из одной синусоиды), а единичное - случаю белого шума. Таким образом, области повышенных значений параметра могут свидетельствовать о моментах повышенной «хаотизации» процесса (например, при распадной неустойчивости с образованием хаотических колебаний), а минимальные значения параметра - большей его «упорядоченности» (например, при синхронизации различных гармоник на единые предельные циклы).

Что касается достоверности определения исследуемых параметров, например, амплитуд и периодов циклических вариаций, зависит от погрешностей исходных временных рядов, дисперсии и статистического смещения рассчитываемых оценок параметров, разрешающей способности анализа. Иначе говоря, характеристика достоверности данных основывается на стандартных приемах доверительного оценивания. Как правило, мы стремимся определить спектрально-временную структуру процессов. Поэтому для нас наиболее важно знать, насколько достоверны выделяемые ритмические составляющие на СВАН-диаграммах. Интерпретируя их, мы оцениваем, во сколько раз амплитуды ритмических компонент превышают шум и сколько периодов укладывается в интервал времени, когда прослеживается та или иная выделяемая гармоника. Если амплитуда гармоники превышает амплитуду шума в несколько раз и выделяется достаточно протяженная полоса (в которую укладывается хотя бы несколько периодов), подтверждаемая наблюдениями в другие интервалы времени или на сходных объектах, то можно уверенно говорить о существовании в этом процессе данного ритма.

Сравнение методов анализа

Для того, чтобы понять, как все эти методы работают, наглядно представим результат обработки синусоидного сигнала каждым из них (рис.9)

Рис.9. Сравнение методов анализа.

На рисунке 9 представлен синусоидный сигнал, частота которого увеличивается со временем, пусть для простоты единицы по оси x будут года. Как видно из рисунка, Фурье-преобразование даёт информацию о наличии той или иной частоты в данном сигнале, но при этом, когда именно эти частоты присутствовали нам не известно. В то же время СВАН- и вейвлет- анализы дают чёткую картину динамики изменения частотных характеристик во времени. На обеих рисунках видно, что доминирующий сигнал периодом примерно 120 лет, затем сменился на период длиной 35 лет, а в конце - примерно 5-ти летний период. Так же в СВАН-диаграмме по степени яркости всплесков видно, что амплитуда всех трёх составляющих сигналов одинакова. Но, как говорилось выше, в СВАН-диаграмме есть недостатки - потеря части ряда при преобразовании: это видно на рисунке 9 - в начале и в конце ряда отрезается несколько значений. Таким образом, можно сделать вывод, что вейвлет-преобразование и СВАН-анализ более информативны, по сравнению с преобразованием Фурье. Однако, для более детального анализа ряда всё же требуется Фурье-спектр, чтобы видеть, какие именно частоты доминируют в хронологии. Например, по данным Фурье-анализа видно, что две доминирующие частоты находятся в интервале от 0 до 0,05 [циклов/год]. В этом случае СВАН-анализ логичнее проводить в диапазоне именно этих частот (как и было сделано), чтобы более «крупно» видеть получившуюся картинку.

Для оценки степени упорядоченности сигнала используется параметр хаотизации. Он вычисляется для каждого момента времени. Если он равен 1, это означает, что в исходном ряду практически невозможно выделить доминирующие частоты («состояние хаоса»). Чем ниже этот коэффициент, тем более однозначен и ярко выражен сигнал в исследуемом ряду. В нашем примере этот коэффициент равен практически нулю, что свидетельствует об однозначности сигнала (некоторое его ослабление заметно только во время перехода от одной частоты к другой).

Итак, вейвлет анализ лучше всего демонстрирует частотно-временной спектр, однако спектрально-временной анализ в сочетании с преобразованием Фурье могут составить важное дополнение к вейвлет-преобразованию. График «хаотизации» может подсказать, какие временные интервалы требуют особенного внимания при анализе структуры ряда.

Глава 3. Ритмика природных процессов, отраженная в ширине и плотности колец на севере Русской равнины

Построение локальных хронологий и их выбор для включения в региональные; преобразование локальных хронологий в региональные

На рисунке 10 представлена северная часть европейской территории России и положение всех площадки, по которым в международном банке данных (www.cdc.noaa.gov) представлены первичные материалы по ширине и плотности колец. Площадки 1,3,5,6,8,22,27 сразу были отсеяны, так как по ним есть только данные по ширине колец, а данные по максимальной плотности отсутствуют.

Рис.10. Положение дендрохронологических площадок на севере Русской равнины (www.cdc.noaa.gov)

Из набора измерений ширины и плотности колец отдельных образцов необходимо было составить сводные хронологии. Для этого мы использовали программу ARSTAN, которая исключает возрастной тренд путём индексирования ряда, а также исключает эффект эндогенности (внутренние причины) и помех со стороны других деревьев (Cook, Holmes, 1986). Пример такого индексирования представлен на рисунке 7. Таким образом, мы получили по 20 хронологий (всего 27 площадок минус 7 - те, по которым нет необходимых нам данных) максимальной плотности и ранней ширины древесины. После этого, с помощью коэффициента корреляции Пирсона мы рассчитывали взаимную корреляцию между площадками, чтобы выяснить, насколько сходны полученные хронологии между собой. Площадки, коэффициенты корреляции которых с остальными оказались незначимыми были исключены из дальнейшего анализа. По ширине колец были исключены площадки №№ 10,11,15,20, по максимальной плотности - №№ 10,13,15,20,24 (коэффициент корреляция ниже 0,3). В результате этой процедуры осталось 15 и 16 площадок, соответственно, которые затем Були объединены в несколько региональных хронологий.

Таким образом, по этому критерию площадки максимальной плотности разбились на 3 группы, из которых было составлено соответственно 3 региональных хронологии, а для ранней ширины - на 4. Объединение проведено путем погодичного осреднения индексов локальных хронологий. На рисунках 11 и 12 показано, как именно локальные площадки были объединены в региональные хронологии, под рисунком показаны коэффициенты корреляции и стрелками указано, к какой именно части территории относятся те или иные коэффициенты.

Коэффициент корреляции четырёх региональных хронологий по ширине равен от 0,27 до 0,56, а для трёх рядов по максимальной плотности - 0,37-0,67. Все коэффициенты корреляции значимы на 95% уровне. Это означает, что и ширина и плотность колец отражают некий общий климатический сигнал, причем каждый параметр - свой сигнал (ширина и плотность также скорелированы между собой, но с более низким коэфициентом корреляции).

Максимальная плотность несколько более тесно скоррелированый параметр и это дает основание предполагать, что ее связь с климатическими параметрами будет более заметной. Из работ наших предшественников (Ваганов и др., 1996, Briffa et al., 2004) известно, что максимальная плотность коррелирует с температурой теплого периода, а температура сама по себе - характеристика менее изменчивая в пространстве, чем осадки. Наш анализ косвенно подтвержадет этот вывод.

Таким образом у нас получилось всего 7 региональных хронологий (см.рис.13), которые далее исследуем с помощью вейвлет-, СВАН- и Фурье-анализа. Выявим доминирующие частоты в данных рядах. Кроме того, изучим связь метео данных с нашими рядами.

Рис.11. Разбиение площадок на региональные хронологии по ширине ранней древесины.

Региональные хронологии соответсвенно относятся к Кольскому полуострову, Карелии, северной части Русской равнины и Уралу. Внизу - соответствующие коэффициенты корреляции

Рис.12. Разбиение площадок на региональные хронологии по максимальной плотности древесины.

Региональные хронологии соответственно относятся к Кольскому полуострову и Карелии, северной части Русской равнины и Уралу. Внизу - соответствующие коэффициенты корреляции

Рис.13. хронологии ранней ширины и максимальной плотности колец хвойных на северной границе леса на Русской равнине. По оси ординат - индексы ширины и плотности колец

Корреляция хронологий с метеорологическими данными. Выявление доминирующего климатического сигнала

Для того чтобы выяснить, от каких факторов зависит ширина или плотность кольца, необходимо провести корреляционный анализ этих параметров со среднемесячными осадками и температурами ближайших метеостанций или данными реанализа. Для этого мы использовали программный пакет DendroClim2002 и построили так называемую функцию отклика для нескольких локальных и всех семи региональных хронологий и метеостанций, показанных на рис.14 и перечисленных в таблице 2. Если функции отклика для одного района имеют сходный рисунок, можно говорить о том, что хронологии отражают некий общий климатический сигнал.

Рис.14. Расположение метео станций использванных для анализа (синие точки)

Таблица 2. Список метеостанций и региональных и локальных хронологий, использованных для расчета функций отклика ширины и плотности колец деревьев на климатический сигнал

Ширина ранней древесиныРегиональная хронологияМетео станция№ локальной хронологииМетео станцияКольский п-овМурманск, Кандалакша2Мурманск12АрхангельскКарелияРеболыСеверная часть равниныАрхангельск, Шенкурск, Нарьян-Мар, Усть-Цильма, Печора16Нарьян-Мар18Печора24ЧердыньУралТулпан, Чердынь, Ивдель25-26ИвдельМаксимальная плотностьРегиональная хронологияМетео станция№ локальной хронологииМетео станцияКольский п-ов + КарелияМурманск, Кандалакша, Реболы, Санкт-Петербург2Мурманск11Санкт-ПетербургСеверная часть равниныАрхангельск, Шенкурск, Нарьян-Мар, Усть-Цильма, Печора12Архангельск16Нарьян-Мар18ПечораУралТулпан, Ивдель25-26ИвдельРезультаты корреляционного анализа представлены на рисунках 15,16,17. Согласно этим результатам, максимальная плотность древесины как на уровне локальных хронологий, так и для региональных, устойчиво коррелирует с температурой летнего периода (май-август) (рис.15). Коэффициент корреляции трёх региональных хронологий с данными летней температуры ближайших метеостанций и коэффициент отдельных локальных с этими же станциями составляет от 0,22 до 0,78) с уровнем значимости 95%. Таким образом, наши выводы согласуются с выводами К.Бриффы о корреляции максимальной плотности региональных хронологий для Сибири и севера Европы (включая Скандинавию) с температурой лета (см.рис.2).

Рис.15. Функция отклика максимальной плотности региональных хронологий и температурных данных метеостанций, находящихся внутри данного региона (по оси х - месяцы, по оси у - значении функции отклика, по оси z - метеорологические станции)

Рис.16. Функция отклика максимальной плотности локальных хронологий и температуры отдельных станций (подписи см.рис.15)

Согласно нашей рабочей гипотезе, основанной на выводах Кононова и др. (200,), ширина колец ранней древесины в этом регионе отражает количество осадков холодного периода. Однако наш анализ показывает, что функции отклика, как для осадков, так и для температур в этом случае хаотичны, и, следовательно, не дают представления о том, какой климатический фактор оказывает господствующее влияние на параметр ширины ранней древесины годичных колец ни для локальных, ни для региональных хронологий.

Рис.17. Функция отклика ранней ширины региональных хронологий с количеством осадков (подписи см рис.15)

Рис.18. Функция отклика ранней ширины региональных хронологий с температурными данными метеостанций, находящиеся внутри региона (подписи см рис.15)

Таким образом, в северной части Европейской России максимальную плотность колец можно использовать как индикатор, отражающий летнюю температуру в прошлом, в то время как ширина колец ранней древесины непригодна для реконструкции температур и осадков.

Вейвлет-анализ хронологий и метеорологических данных

На рис.19 и 20 показаны результаты вейвлет-анализа каждой из региональных хронологий ширины колец ранней древесины и максимальной плотности годичного кольца. Анализ проводился с помощью математического пакета MatLab 6.0. Кроме того, мы провели вейвлет-анализ летних температур и осадков холодного периода, измеренных на самых длиннорядных метеостанций, расположенных в этом регионе и соотносимые с исследуемыми региональными хронологиями: Санкт-Петербург (Карельская хронология), Архангельск (северная часть), Усть-Цильма (северная часть), Чердень (Урал). На получившихся диаграммах (рис.19,20,21,23) по осям абсцисс и ординат откладываются время (в годах) и период соответственно. Шкала справа показывает амплитуду сигнала (относительная величина).

Рис.19. Вейвлет-анализ рядов ширины ранней древесины с помощью DOG-вейвлета (по осям абсцисс - время в годах, по осям ординат - длина периода в годах, справа колонка - спектр мощности)

Рис.20. Вейвлет-анализ рядов максимальной плотности древесины с помощью DOG-вейвлета (подписи см.рис.19)

Рис.21. Вейвлет-анализ метео данных по осадкам с октября по апрель следующего года на различных станциях (подписи см.рис.19)

Рис.22. Упрощённое представление циклов, полученные с помощью вейвлет-анализа. По оси абсцисс - время (в годах), по оси ординат - длина периода (в годах).

Рис.23. Вейвлет-анализ метео данных по температуре с мая по август на различных станциях (подписи см.рис.19)

В результате этого исследования замечены следующие особенности региональных хронологий ширины и плотности колец (рис. 20-21). Для большей наглядности на рис.22 выделенные периоды показаны в упрощенном виде.

1.В первую очередь, хочется обратить внимание на то, что периодичность в рядах ширины колец ранней древесины в целом ярче выражены, чем в рядах плотности. И при этом степень хаотизации для всех хронологий по ширине примерно 0,25, а по плотности около 0,4. Это связано с большей амплитудой изначального ряда ширины колец (примерно 1), по сравнению с максимальной плотностью примерно 0,4 (см.рис.13). То же самое показывает и преобразование Фурье - по ранней ширине хорошо отслеживаются доминирующие частоты в сигнале.

2.На диаграммах ширины и плотности колец видно, что во всех хронологиях присутствует в первую очередь примерно столетний период. Несмотря на то, что анализируемые ряды (около 300 лет) слишком коротки для достоверного выявления периода такой продолжительности, полностью игнорировать присутствие данного ритма нельзя. Как будет показано ниже, этот период, так как он часто фигурирует в других дендроклиматических рядах большей длительности в субарктике.

3.На всех рисунках (кроме Урала) виден сигнал с периодом в 55-60 лет, который обрывается примерно в 1860 году. На Урале по ширине древесины он, напротив, появляется в 19 в. (чётко начинает проявляться с 1830 года).

4.Во всех хронологиях отслеживается сигнал, близкий к 20-ти годам, который по ширине ранней древесины существует примерно до 1800 года, а по максимальной плотности - постоянно (кроме Урала, там - до 1900 года).

.Рассмотрим, как изменяются хронологии ширины колец ранней древесины на всей Европейской части России с запада на восток. Кольский полуостров и Карелия очень схожи во временной и частотной локализации циклов. 100 и 60-ти летние присутствуют одновременно в обеих хронологиях, на Кольском п-ове 20-ти летний цикл сменился 30-ти летним примерно в 1830 году, а в Карелии на несколько лет раньше, но потом он перешёл в 35-ти летний. Таким образом, хотя корреляционный анализ показал различие локальных хронологий в этих двух районах, в них отмечается сходная ритмическая структура климатического сигнала, и, скорее всего, климатические изменения к этом регионе в последние 300 лет имеют сходный характер.

На севере Русской равнины присутствуют периоды той же длины, что и на западе, причём, 60-тилений период так же сменяется 30-32-хлетним примерно в 1860 году. Так же как и на западе здесь существует период длиной 18-20 лет до 1850 года. Отличие только в том, что на севере после 1900 года появляется циклы длиной 17 и 10-11 лет. Несмотря на то, что ширина ранней древесины плохо коррелирует с метеоданными, осадки холодного периода (рис.23, станция Архангельск) показывают наличие аналогичных периодов.

На Урале картина несколько иная. Здесь (как было сказано выше) 60-ти летний период появляется только в 1830 год; заметен так же цикл длиной 35-37, просуществовавший до 1900 года. Но при этом так же, как на западе и на севере Русской равнины, отмечается периодичность длиной 20 лет, присутствовавшая до 1830 года. Кроме того, здесь (как и на западе) заметен 10-11-ти летний период, который появился примерно с 1850 года, но с 1900 он перешёл в 6-ти летний. Диаграмма метеоданных (осадков) со стации Чердынь так же начиная с 1900 года показывает присутствие 10-11-ти и 6-ти летних периодов.

Таким образом, можно заключить, что ширина колец ранней древесины в районах Кольского полуострова и Карелии с одной стороны и в северной части Русской равнины, с другой имеют сходную структуру ритмов в последние 300 лет. в то время как региональная хронология Урала существенно отличается от первых трех по своей ритмике. Это различие подтверждается и низкой корреляцией этой хронологи с остальными. Сходство ритмической структуры у разных региональных хронологий и их совпадение со структурой соответствующих метеорологических рядов указывает на общность климатического сигнала ими управляющего. Показательно, что в этом смысле запад и север рассматриваемого региона существенно отличается от восточной, уральской хронологии.

6.Теперь рассмотрим изменения ритмов в рядах максимальной плотности на трех региональных площадках. В первую очередь видно, что во всех трёх хронологиях практически на протяжении всего времени присутствует 10-11-ти летний период (по ширине они наблюдались только эпизодически). Кроме того, в западном и северном регионе в течение всех 300-х лет виден постоянный сигнал длиной 18-20 лет, а на Урале схожий сигнал (15-17 лет) просуществовал до 1900 года, а в 1800 он изменился и стал более долгопериодным - 22-25 лет. Если рассматривать цикл длиной 60 лет, то на западе он так же, как и для ранней ширины прослеживается до 1850 года, на Севере сходный сигнал проявляется тогда же, но длиной примерно 50 лет, а на Урале очень слабо проявляется, начиная с 1930 года. Хронология запада равнины схожа с северной частью, в то время как север довольно тесно коррелирует с Уралом (корреляция 0,42 при значимости 95%). Но так как во всех трёх рядах в основном доминируют короткопериодные колебания, можно сделать вывод, что хронологии по максимальной плотности хорошо реагируют на климатические колебания, периоды которых меньше 30-ти лет. Это видно и на самих рядах - малые амплитуды сигнала свидетельствуют о доминирующем высокочастотном сигнале.

7.Сравнивая хронологии максимальной плотности за последние 100 лет с летними температурами за этот же период, видно, что, за исключением периода длиной 60 лет, сходство в диаграммах очень заметное. В Карелии присутствуют 40, 20 и 11-ти летние периоды, на севере 15-18 и 8-11-ти летние, а на Урале 25-28 и 8-11-ти летние периоды. То есть, подтверждается предположение, что максимальная плотность хорошо отражает температуру летнего периода.

.Если рассмотреть отдельно северную часть Русской равнины, то в первую очередь необходимо отметить самый низкий показатель степень хаотизации по ширине ранней древесины. Так же нигде кроме севера нет 80-ти летнего цикла, который появился в хронологии по ширине примерно с 1870, а по плотности с 1770 года. Сравнив хронологию ранней ширины с метеоданными по осадкам (Архангельск), видно, что с 1900 г. в обоих рядах в одно и то же время присутствуют периоды длиной 25-30, 17 и 10-11 лет, которые обрываются так же одновременно. А если сравнить с данными станции Усть-Цильма, то так же видны одновременно существующие периоды длиной 17 и 8-11 лет. То есть, на севере Русской равнины короткие циклы (до 30-ти лет) одновременно присутствуют в метео и дендроклиматических данных.

.Для хронологий Урала по плотности и по ранней ширине степень хаотизации максимальна (рис.24). Кроме того, как было сказано выше, это единственное место, где по ширине древесины 60-ти летний цикл не исчезает в 1860, а наоборот - появляется в 1840 году. По максимальной плотности период этой длины начиная с 1950 года проявляется слабо, что видно так же и по ряду летних температур метеостанции Чердынь. Так же на диаграмме по плотности видно, что в первые примерно 100 лет не были чётко ясны доминирующие периоды, потом начали появляться циклы длиной 25 и 40 лет. И только за последние 100 лет структура ряда начала выстраиваться менее хаотично (коэффициент хаотизации так же это подтверждает, с 1900 года он опустился примерно с 0,6 до 0,4) - появились ярко выраженные 100, 25 и 10-11 летние периоды. Но, несмотря на высокую степень хаотизации, при этом, как было сказано выше, Уральская региональная хронология максимальной плотности хорошо коррелирует с севером Русской равнины.

Рис.24. Графический анализ максимальной плотности хронологии Урала (подписи см.рис.9)

Обсуждение возможных причин периодов разной продолжительности, выявленных в древесно-кольцевых хронологиях

Ритмы разной продолжительности могут иметь разное происхождение. Более того, иногда и ритмы одинаковой или близкой продолжительности также бывают вызваны разными причинами (Атлас временных вариаций). Понятно, что явления с одинаковой ритмической структурой могут быть никак не связанными причинно-следственными связями, поэтому к интерпретации ритмических сигналов в рядах разного просхождения следует подходить с осторожностью и прежде всего, учитывать физические процессы, которые могут объяснить ритмическое сходство разных рядов. Имея ввиду эти положения, мы лишь гипотетически можем объяснить происхождение ритмов, обнаруженных нами в дендрохронологических данных.

Самый длинный период - около 100 лет, обычно связывают с длиной солнечного цикла. Хотя в рядах длиной 300 лет этот цикл, возможно появляться как краевой эффект, присущий вейвлет-анализу, и таким образом является артефактом. Однако, этот период часто обнаруживается в дендросериях большей продолжительности. Так например, Кляшторин и Любушин (2005), изучая периодические флуктуации температуры по годичным кольцам сосны в Скандинавии за 1400 лет так же выявили 108-летний доминирующий период (рис.25). Скандинавия - наиболее близкий к исследуемому регион не только в пространственном, но и в климатическом отношении. Поэтому вероятность того, что 100-летний период - это не краевой эффект, а истинная периодичность в анализируемых рядах очень велика.

Рис.25. Спектр периодических флуктуаций температуры по кольцам арктической сосны за 1400 лет

-ти летний цикл, который ярче всего проявился в хронологиях максимальной плотности, т.е. в рядах, отражающих температуру летнего времени года, с большой долей вероятности отражает именно влияние Солнечной активности. Доминирование этого сигнала в дендрохронологических рядах из самых разных районов хорошо известно и описано в соответствующей литературе (см. например, Fritts, 1976) На рис.26 показан Фурье и вейвлет-анализ ряда солнечной активности. На этих рисунках также виден доминирующий сигнал длиной 100 и 11 лет. Кроме того, здесь виден 60-ти летний цикл солнечной активности, который так же проявляется во всех наших хронологиях до 1850 года, а на Урале начиная с 1840 года.

Рис.26. Фурье и вейвлет-преобразование ряда солнечной активности

Известно, что на климат рассматриваемой территории существенное влияние оказывает так называемый North Atlantic Oscillation Index (индекс Северо-Атлантического колебания - разница высоты поверхности 500 гПа между пунктами в Исландии (64° с.ш., 24° з.д.) и вблизи Азорских островов (39° с.ш., 24° з.д.) за зимне-весенний период). Этот показатель является мерой активности западного переноса воздушных масс и поэтому крайне важен для формирования климата Европы и северной части России. Чтобы понять, какие периоды доминируют в этом индексе, приведём вейвлет-анализ ежегодного и сезонных (май-август, октябрь-апрель) индексов НАО (www.cpc.ncep.noaa.gov) (рис.27). Из диаграмм видно, что его основные периоды составляют примерно 5, 10-15, 20 и 40 лет.

Рис.27. Вейвлет-анализ ежегодного и сезонных индексов NAO (подписи см.рис.19)

Самое сильное сходство видно между индексом NAO за летний период (май-август) и ряда плотности по Карелии и Кольскому полуострову. В обоих рядах присутствуют периоды длиной 40, 20-22 и 5 лет. Схожесть индекса с хронологией севера и Урала проявляется только в цикле длиной порядка 20-ти лет. Следовательно, можно предположить, что этот индекс за летний период большее влияние оказывает именно на западную часть равнины, в то время как на север и на Урал в меньшей степени. Что касается индекса в зимний период, то здесь видно сходство со всеми тремя рядами плотности - везде присутствуют периоды 5, 10 и 20 лет. Так же замечено сходство индекса NAO в зимний период со всеми тремя хронологиями максимальной плотности в периодах длиной 20 и 10 лет, а на западе равнины ещё схожий и 5-ти летний цикл. Таким образом, если предположить, что периоды длиной 20 и 5 лет обусловлены влиянием Северо-Атлантического колебания, то этот индекс влияет именно на максимальную плотность, так как сильно выраженного сходства индекса за летний и зимний период с хронологиями ранней ширины не наблюдается.

Заключение

В ходе дипломной работы было:

1.Изучен метод дендрохронологии в обработке данных по ежегодному приросту древесины.

2.Освоены на практике лабораторные методы обработки образцов древесины.

.Знакомство с международной дендрохронологической базой данных (National Climatic Data Centre).

.Создан и структурирован массив данных, необходимых для дальнейшего анализа климата Русской равнины.

.Проведён корреляционный анализ хронологий

.Изучена связь метеоданных и хронологический рядов

.Изучены основы метода вейвлет-анализа.

.С помощью Фурье-, СВАН и вейвлет-анализа изучена ритмика региональных хронологий европейской части России.

В результате наших исследований мы установили, что по индикатору ранней ширины кольца все площадки лучше всего разбить на 4 региональные хронологии - Кольский полуостров, Карелия, Северная часть Русской равнины и Урал, а по максимальной плотности - на 3: Кольский полуостров с Карелией вместе, северная часть и Урал. Длина всех рядов - примерно 300 лет. Таким образом, мы получили 7 хронологий. Хронологии по максимальной плотности несколько лучше коррелируют между собой, чем по ранней ширине (средний коэффициент корреляции 0,51 против 0,41 соответственно при уровне значимости 95%). Более высокая корреляция максимальной плотности в пространстве скорее всего связана с тем, что после летних температур, от которой в основном зависит этот параметр, обычно мало изменчиво в пределах европейского севера России.

Корреляционный анализ показал, что максимальная плотность колец довольно сильно коррелирует с температурой летнего периода (май-август), коэффициент корреляции при этом составляет от 0,22 до 0,78. К сожалению, сказать то же самое о связи между ранней шириной и осадками в зимний период (октябрь-апрель следующего года) нельзя. Корреляция практически отсутствует, и нет никакого выраженного сезона, где бы она была бы постоянной. Таким образом, максимальная плотность отражает летнюю температуру, тогда как климатические параметры, определяющие прирост ранней древесины остались невыясненными. Гипотеза Кононова и др. о связи этого параметра с зимними осадками не подтвердилась. Однако, скоррелированность ширины ранней древесины на площадках, удаленных друг от друга на тысячу и более километров показывает, что такой общий внешний управляющий параметр(ы) существует. Сходство ритмической структуры рядов ширины ранней древесины и осадков зимы косвенно свидетельствует, что это в какой-то мере могут быть именно зимние осадки. Отсутствие непосредственной корреляции с метеоданными может быть вызвано большой пространственной изменчивостью снегонакопления и удаленностью метеостанций от мест отбора образцов.

Сопряженный анализ региональных хронологий, метеорологических данных, солнечной активности и Индекса Северо-Атлантического колебания был проведен с помощью вейвлет-преобразований, простого и оконного анализа Фурье (в работе он фигурирует как Спектрально-временной анализ - СВАН). Мы установили, что

·Хронологии ранней ширины колец хорошо отражают циклы только высокой продолжительности (40-100 лет), а хронологии максимальной плотности - ещё и более короткие циклы (5-20 лет).

·Во всех семи рядах присутствуют 100-летний цикл и 11-летний циклы которые, скорее всего, связаны с Солнечной активностью.

·Во всех рядах кроме Урала присутствует и сигнал 60-ти летней продолжительности примерно до 1850 годы, а на Урале он появился с 1840 года.

·Более высокочастотные сигналы 40, 20-22 и 5, обнаруженные в древесно-кольцевых хронологиях, скорее всего связаны с климатическими флуктуациями, вызванными колебаниями Северо-Атлантического Индекса, имеющего сходную ритмическую структуру.

·В целом, хронологии Кольского полуострова, Карелии и северной части Русской равнины демонстрируют большое сходство ритмического рисунка, Уральская - существенно от них отличается. Это согласуется с современными представлениями о климате севера европейской части России. В то время как на Урале периоды (не считая 100-летнего цикла), выявленные в хронологии ширины отличаются и по частоте и по временной локализации.

Список литературы

древесина вейвлет прирост

1.Шиятов С.Г. Ваганов Е.А., Кирдянов А.В., Круглов В.Б., Мазепа В.С., Наурзбаев М.М., Хантемиров Р.М. Методы дендрохронологии. Часть I. Основы дендрохронологии. Сбор и получение древесно-кольцевой информации, 2000.

2.Ваганов Е.А., Шиятов С.Г., Мазепа В.С. Дендроклиматические исследования в Урало-Сибирской Субарктике. - Новосибирск: Наука,1996.

.Ваганов Е.А., Шишкин А.В. Рост и структура годичных колец хвойных. - Новосибирск: Наука, 2000.

.Кляшторин Л.Б., Любушин А.А. Циклические изменения климата и рыбопродуктивности. - М.: Издательство ВНИРО, 2005

5.Briffa K.R., Osborn T.J., Schweingruber F.H. Large-scale temperature influence from tree rings. Global and planetary change 40 (2004).

.Briffa K.F., Osborn T.J., Schweingruber F.H., Jones P.D., Shiyatov S.G., Vaganov E.A. Tree-ring width and density data around the Northern Hemisphere: Part 1, local and reginal climate signals.//The Holocene, 2002.

.Briffa K.F., Schweingruber F.H., Jones F.N., Osborn T.J., Shiyatov S.G., Vaganov E.A. Reduced sensitivity of recent tree-growth to temperature at high northern latitudes.// Nature, 1998.

.IPCC, 2007: Climate change: The Scientific Basis. Contribution of Working group 1 to the Third Assessment Report of the Intergovernmental Panel on Climate Change. Houghton, J.T., Ding,Y., Griggs D.J et al. (eds).Cambridge, United Kindom and New York, USA, Cambridge University Press.

.Cook E.R., A time series analysis approach to tree-ring standardization. The University of Arizona, 1985

.Raspopov, O.M., Dergachev, V.A., Kolstrom, T. Periodicity of climate conditions and solar variability derived from dendrochronological and other paleoclimatic data in high latitudes. Palaeogeography, Palaeoclimatology, Palaeoecology, 2004

11.Смоленцев Н.К. Основы теории вейвлетов. Вейвлеты в MATLAB. - М.: ДМК Пресс, 2005.

.Воробьёв В.И., Грибунин В.Г. Теория и практика вейвлет-преобразования. - СПб: ВУС, 1999.

.Астафьева Н.М. Вейвлет-анализ: основы теории и примеры применения. УФН, 1996. - Т.166. - №11.

.Дремин И.М, Иванов О.В., Нечитайло В.А. Вейвлеты и их использование. УФН, 2001. - Т.177. - № 5

.Дьяконов В.П., Абраменкова И.В. MATLAB. Обработка сигналов и изображений. Специальный справочник. - СПб.: Питер, 2002.

16.http://matlab.exponenta.ru/wavelet/index.php

17.http://www.oceanographers.ru/index.php?option=com_content&task=view&id=324&Itemid=53