Бесплатный автореферат и диссертация по наукам о земле на тему
Почвенно-климатические сценарии в моделях круговорота углерода и азота в лесных экосистемах умеренного пояса
ВАК РФ 25.00.30, Метеорология, климатология, агрометеорология

Автореферат диссертации по теме "Почвенно-климатические сценарии в моделях круговорота углерода и азота в лесных экосистемах умеренного пояса"

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ имени М.В. ЛОМОНОСОВА Географический факультет

На правах рукописи

ООЭ47В32Э

БЫХОВЕЦ Сергей Станиславович

ПОЧВЕ1ШО-КЛИМАТИЧЕСКИЕ СЦЕНАРИИ В МОДЕЛЯХ КРУГОВОРОТА УГЛЕРОДА И АЗОТА В ЛЕСНЫХ ЭКОСИСТЕМАХ УМЕРЕННОГО ПОЯСА

Специальность 25.00.30 - Метеорология, климатология, агрометеорология

- 8 ОПТ

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата географических наук

Москва 2009

003479329

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

Научный руководитель:

Официальные оппоненты:

доктор биологических наук, профессор

Комаров Александр Сергеевич

доктор биологических наук, профессор

Шульгин Игорь Александрович

кандидат географических наук, Ольчев Александр Валентинович

Ведущая организация: Институт географии РАН

Защита состоится «29» октября 2009 г. в 15-00 часов на заседании диссертационного совета Д 501.001.68 при Московском государственном университете имени М.В.Ломоносова по адресу:

119991, Москва, ГСП-1, Ленинские горы, МГУ, географический факультет, 18 этаж, аудитория 1801.

С диссертацией можно ознакомиться в научной библиотеке географического факультета МГУ на 21 этаже.

Автореферат разослан «25»_сентября_2009 г.

Ученый секретарь диссертационного совета, доктор геолого-минералогических наук

В.С. Савенко

Общая характеристика работы

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

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

При этом возникает ряд задач, связанных с:

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

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

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

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

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

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

Для достижения поставленной цели решались следующие задачи:

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

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

3. Разработка субмодели водно-теплового режима почвы.

4. Определение поправок для оценки температурного режима лесных почв в зависимости от типа леса.

5. Верификация моделей и их применение в рамках системы моделей биологического круговорота углерода и азота ЕР1МСЮ для сравнительного

Л

анализа последствии сильных внешних воздействии разных климатических ■.коёаЫ- и типах леса.

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

Практическая ценность результатов. Полученные результаты используются в научных исследованиях для разработки прогноза развития лесных почв при различных сценариях изменений климата и усиления антропогенной нагрузки. Разработанное программное обеспечение также используется в практикумах для студентов Пущинского государственного университета и Филиала МГУ имени М.В.Ломоносова в Пущино,

Личный вклад автора. Автором проанализированы возможности известных в мировой науке статистических генераторов метеорологических данных, разработаны простые статистические модели водно-теплового режима почв, учитывающих тип леса, принято активное участие в создании системы моделей биологического круговорота углерода и азота в лесных экосистемах ЕНМСЮ, проведены вычислительные эксперименты и анализ результатов моделирования по теме диссертации.

Связь темы диссертации с плановыми исследованиями. Тема диссертации связана с основным планом научно-исследовательских работ

Института физико-химических и биологических проблем почвоведения РАН по теме «Математическое моделирование циклов элементов и сукцессионных процессов в системе почва - растительность - атмосфера (№ 01.2 006 073924.1.16.11)».

Основные положения диссертации, выносимые на защиту:

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

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

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

Апробация и внедрение результатов. Результаты работы были представлены на: I Всероссийской научной школе "Математические методы в экологии" (Петрозаводск, 2001), 2й Международной конференции «Эмиссия и сток парниковых газов на территории Северной Евразии» (Пущино, 2003), Всероссийской конференции «Биосферные функции почвенного покрова» поев. 100-летию со дня рождения чл.-корр. АН СССР В.А. Ковды. (Пущино, 2005), 5й Европейской конференции по математическому моделированию в экологии (ЕСЕМ-05, Pushchino, 2005), а также на рабочих совещаниях по проектам ИНТАС (97-20355 "FORMOD", 01-0512 "PODZOL", 01-0633 "SILVICS"), EU FP6 INCO (013388 "OMRISK") в Финляндии, Великобритании и Италии, семинарах Лаборатории моделирования экосистем ИФХиБПП РАН и Учебного

Центра почвоведения, экологии и природопользования Путинского государственного университета.

Модель SCLISS [Быховец, Комаров, 2002; Быховец, 2007] является составной частью системы моделей EFIMOD-2 [Komarov et al., 2003], разработанной в Лаборатории моделирования экосистем ИФХиБПП РАН и Лаборатории биохимии почв Биологического НИИ СПбГУ; полученные на ее основе сценарии использовались и используются при моделировании динамики органического вещества в лесных экосистемах Европейской части России, ряда стран Северной и Центральной Европы (Финляндии, Швеции, Германии, Чешской Республики) и Канады.

Публикации. Результаты работы изложены в 24 публикациях, в числе которых 9 статей в рецензируемых журналах, рекомендованных ВАК, 4 главы в монографиях, учебно-методическое пособие, 7 статей в сборниках научных трудов и 3 тезисов докладов.

Структура и объем диссертации. Диссертация состоит из введения, б глав, заключения и библиографии. Она включает 148 страниц текста, 24 рисунка, 17 таблиц. Список литературы состоит из 220 наименований на русском и английском языках.

Автор искренне благодарен своему руководителю проф., д.б.н. А.С. Комарову, а также проф., д.б.н. О.Г. Чертову за помощь в постановке задачи и постоянное внимание к работе и коллегам - сотрудникам Лаборатории моделирования экосистем ИФХиБПП РАН и Учебного Центра почвоведения, экологии и природопользования Пущинского государственного университета за помощь в проведении исследований и моральную поддержку.

Работа выполнена при поддержке РФФИ (проект № 05-04-49284), Европейского института леса (Йоэнсу, Финляндия; проекты № 619 и 649), INTAS (97-20355 "FORMOD", 01-0512 "PODZOL", 01-0633 "SILVICS"), EU FP6 INCO (013388 "OMRISK"), a также Программы фундаментальных

исследований РАН № 16 «Окружающая среда в условиях изменяющегося климата: экстремальные природные явления и катастрофы».

Краткое содержание работы

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

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

Для генерировании искусственных рядов метеорологических переменных с заданными статистическими свойствами служат имитаторы, или генераторы погоды, разработка которых интенсивно ведется в течение последних 20-25 лет; в основном - применительно к агрометеорологическим моделям [Richardson, 1981; Сиротенко, Павлова, 1986; Бельченко, 1989; Топаж, 1992; Racsko et al„ 1991], и пока в меньшей степени - к моделям лесных экосистем [Strandman et al, 1993]. Практически все известные имитаторы имеют суточное разрешение что, безусловно, необходимо для большинства современных агрометеорологических моделей. При этом одной из основных проблем является моделирование достаточно сложной временной структуры рядов

суточных сумм осадков. Для ее решения чаще всего используется один из двух подходов: представление последовательности сухих и дождливых дней в виде марковского процесса [Richardson, 1981], или оценка распределения вероятностей продолжительности дождливых / сухих периодов [Racsko et al., 1991]. Лишь изредка встречаются генераторы с декадным [Бельченко, 1989], или месячным [.Голубятников и др., 2002] разрешением. При моделировании же лесных экосистем обычно не требуется столь высокая степень временной дискретизации, как в агрометеорологических задачах; во многих случаях бывает достаточно средних месячных значений.

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

Модели, ориентированные на лесные условия (либо достаточно универсальные, чтобы позволить учесть характерные для леса особенности почвенного тепловлагообмена) сравнительно немногочисленны как в отечественной литературе [Мамихин, 1997; Гусев, Насонова, 1998, 2002, 2007; Кучмент, Демидов, 2004] так и мировой [Halldin, 1980; Janssen, 1980, 1991; Arp, Yin, 1992; Yin, Arp, 1993; Tiktak, van Grinsven, 1995]. В то же время, в нашей стране выполнен целый ряд серьезных экспериментальных исследований [Молчанов, I960, 1961; Федоров, 1977; Галенко, 1983; Еруков, Власкова, 1986; Абражко, 1998, и др.], результаты которых могут быть успешно использованы в подобной работе.

Статистические модели требуют обширного экспериментального материала, но при этом недостаточно гибки и универсальны. Они могут успешно применяться для решения определенных прикладных задач в условиях, достаточно близких к тем, для которых были получены. Хотя имеется опыт применения подобного подхода к температурному режиму лесных почв [Bocock et al. 1977; Paul et al, 2004], довольно сложно найти достаточное

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

Модели водного режима более многочисленны, чем модели теплового режима или сопряженные модели тепло- и влагообмена. Собственно статистический подход применяется редко, наиболее широко распространены как простые воднобалансовые модели (в англоязычной литературе часто называемые «bucket models»), так и собственно динамические модели, основанные на решении уравнения влагопереноса, а также ряд промежуточных вариантов: многослойные воднобалансовые и «квазидинамические» (semi-dynamic) [de Jong, Bootsma, 1996; Tiktak, van Grinsven, 1995]. Динамические модели, как и в предыдущем случае, наиболее универсальны, но требуют знания большого количества трудноопределяемых параметров.

В качестве примеров рассмотрены «климатические» компоненты наиболее известных моделей динамики органического вещества почвы, при этом показано, что поставленная задача может решаться на разном уровне сложности, от использования весьма детализированной динамической модели водно-теплового режима (система моделей CoupModel [Jansson, Karlberg, 2002]) до полного отсутствия почвенно-гидротермического блока, и использования регрессионных зависимостей скоростей почвенно-биогеохимических процессов непосредственно от «внешних» климатических переменных (модель Yasso [Liski et al, 2005]).

Кратко рассмотрена система моделей биологического круговорота углерода и азота в лесных экосистемах EFIMOD, проанализированы методы

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

Вторая глава посвящена формированию «входных» рядов климатических данных (в базовом варианте - температуры воздуха и осадков). При этом могут использоваться как ряды фактических наблюдений (данные близлежащих метеорологических станций), так и ряды, сгенерированные с помощью модели, основанной на ряде предположений, вытекающих из анализа климатологической литературы [Жуковский и др., 1976; Кобышева, Наровлянский, 1978; Исаев, 1988, 2002; Бельченко, 1989; Швер, 1976, 1984; Климат ..., 2001; Heino, 1994; и др.]. Месячные суммы осадков гт представляются независимой случайной величиной с логнормальным распределением:

rm =exp(/^m +Ут 'ni) (1)

где т - номер месяца, рт - lnrm и vm = а (hl гт) - параметры логнормального распределения, которые могут быть выражены через традиционно используемые климатологами характеристики: многолетнюю среднюю гт и коэффициент вариации сv(rm):

Мт = In fc/Jl + C v2(rj) vm = J\n{\ + Cv2(rm))

«, - нормально распределенная случайная величина N(0,1).

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

Та,т =та,г, +Ваа,т -Та„_)+Вагт^пгт-ЫгпУпг-БТат

(2)

где: Т„ шА - температура предыдущего месяца, гт - сумма осадков текущего месяца, ТаП1, 7;лМ и 1пг„ - средние значения соответствующих величин, п2-нормально распределенная случайная величина N(0,1);

) - Р{Та,т . 1ч г» ) ■ Гт ' Т„.т-\ )

в р(7')~')•р(1пг"')

- соответствующие коэффициенты регрессии,

^ -1-рг(1пг..Г^)-

остаточное отклонение; - коэффициент корреляции температур ти-

го и (т-1)-го месяцев (значение автокорреляционной функции со сдвигом 1), р{Т^аМгт) и р(1пг„,Тйт.,) - коэффициенты корреляции логарифма месячной суммы осадков т-то месяца со средними температурами т-то и (т-1)-го месяцев (значения взаимной корреляционной функции со сдвигом 0 и 1).

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

региональном масштабе можно считать относительно однородньми в пространстве).

В качестве модельных климатических сценариев, учитывающих ожидаемые изменения климата в XXI веке, могут использоваться как сгенерированные, так и «экстраполированные» ряды аномалий температуры и осадков, в сочетании с трендами, прогнозируемыми на основании моделирования «будущего» климата с использованием моделей общей циркуляции атмосферы и модельных сценариев роста содержания С02 в атмосфере [IPCC, 2000, 2003; Mitchell etal., 2004].

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

-тепловое состояние корнеобитаемого слоя минеральной почвы может быть охарактеризовано ее температурой на глубине 20 см [Димо, 1972]\

- температура подстилки считается равной температуре воздуха в теплое время года (Та >0), и температуре почвы на глубине 20 см - в холодное.

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

т,.т-7Ü + 2U .{т.м-TZ)+n,.STsjl (3)

где Tsm и Г5„.] - температура почвы за текущий и предыдущий месяцы, Тат-температура воздуха (текущая), Т,„, и Tajn - соответствующие

многолетние средние, В„„ и В„т - коэффициенты регрессии, 5П„ - остаточное отклонение (определяются подобно параметрам предыдущей субмодели).

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

станций). А затем, с учетом разностей температур «лес -метеорологическая площадка», полученных для разных типов леса из анализа литературных данных {Молчанов, 1961; Елагин, Изотов, 1968; Павлов, 1975; Казимиров и др., 1978; Воронков и др., 1979; Галенко, 1983; Еруков, Власкова, 1986 и др.], оценивается температура почвы в лесу.

Температура почвы измеряется на более редкой сети станций, чем «основные» метеорологические элементы, а для ряда интересовавших нас стран Европы сведений по температуре почвы найти вообще не удалось. Поэтому нами, на основе многолетних данных Справочника по климату СССР [19641968], с использованием идей В.Н. Димо [1967, 1982] и В.Р. Волобуева [7983], были получены регрессионные зависимости средней температуры почвы на глубине 20 см от температуры воздуха. Это позволяет оценить многолетние нормы средней месячной температуры почвы для каждого месяца в случае отсутствия фактических данных:

( = [а<>+% при '„2 0

1 \а0+а21о при 1а <0' (4)

Значения коэффициентов а0, а\ и а2 приведены в таблице 1.

Месяц

Рисунок 1. Годовой ход разностей температур почвы (на глубине 0.2 м) под пологом леса и под косимой травой (на метеорологической станции) в зависимости от доминанты древостоя

Таблица 1.

Значения коэффициентов в уравнении регрессии (4)

Месяц ao±s.e. а/± s.e fl2±s.e R2 SEE

I 2.53 ± 0.17 0.57± 0.13 0.45 ± 0.01 0.753 2.48

II 2.55 ± 0.17 0.69 ± 0.09 0.51 ± 0.01 0.791 2.35

III 1.66± 0.12 0.86± 0.03 0.59 ± 0.01 0.894 1.68

IV -0.69 ± 0.08 1.10± 0.01 0.69± 0.02 0.951 1.28

V -2.95 ± 0.18 1.21± 0.01 - 0.925 1.65

VI -4.26 ± 0.34 1.26± 0.02 - 0.861 2.02

VII -4.34± 0.38 1.25± 0.02 - 0.868 1.85

VIII -2.25 ± 0.26 1.20 ± 0.01 - 0.919 1.47

IX 0.41 ± 0.13 1.14± 0.01 - 0.951 1.15

X 2.80± 0.07 0.97 ± 0.01 0.43 ± 0.02 0.958 1.04

XI 2.77± 0.08 0.90± 0.02 0.41 ± 0.01 0.914 1.37

хп 2.26± 0.13 0.76± 0.06 0.40± 0.01 0.802 2.04

Примечание: s.e. - стандартная ошибка соответствующего коэффициента, R -коэффициент детерминации, SEE - ошибка уравнения регрессии.

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

В Главе 4 описана субмодель влажности почвы. Построение для влажности почвы статистической модели, аналогичной описанной в предыдущей главе, в принципе возможно, но мы не располагаем для этого достаточным количеством данных. Поэтому, мы использовали простую воднобалансовую модель. За основу такой модели принят широко известный комплексный метод расчета испарения [Будыко, 1956, 1971; Зубенок, 1976], основанный на сопряженном расчете составляющих водного баланса и их зависимости от влажности почвы. При этом суммарная испаряемость оценивалась упрощенным способом по температуре воздуха, поскольку ее определение в соответствии с комплексным методом требует дополнительных входных данных.

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

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

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

- Детализация учета структуры эвапотранспирации. Суммарное испарение (эвапотранспирация) лесной экосистемы может быть представлено в виде суммы трех компонент: испарения осадков, задержанных (перехваченных) кронами деревьев Ег; транспирации древесного яруса ЕТ и эвапотранспирации напочвенного покрова Е0 (без разделения последней на собственно транспирацию живого напочвенного покрова и физическое испарение):

Е ~ Ер + Ет + ,

причем, очевидно, ЕР, Ет и Еа являются функциями как метеорологических условий, так и характеристик растительности. Для параметризации этих зависимостей использованы эмпирические обобщения [Федоров, 1977; Крестовский, 1986; Книзе, Крестовский, 1993; Бондарик, Карпечко, 1999; и др.] основанные на экспериментальном материале для лесов Европейской России.

Для оценки параметров субмодели влажности почвы (объемной плотности, влажности устойчивого завядания, наименьшей полевой и полной влагоемкости, а также отношения влажностей подстилки и нижележащего минерального горизонта) были собраны и проанализированы литературные данные о водно-физических свойствах лесных почв [Абражко и др., 1973; Алифанов, Лошакова, 1981; Васильев, 1950, 1959; Верхоланцева, 1972, 1972а, 1977; Верхоланцева, Бобкова, 1972; Воронина, 1975; Газизуллин, Хасаншин,

1980; Еруков, Власкова; 1986; Забоева, 1975; Кашанский, 1974; Кононенко, 1986; Подзолистые почвы ..., 1977, 1979, 1980, 1981; Почвы зоны переброски..., 1983; Руднева, 1973, 1983; Структура ..., 1973; Структурно-функциональная .... 2001 и др.]. В настоящее время входные параметры модели задаются на основе подбора аналогов (по типу леса, типу и гранулометрическому составу почвы). Представляется перспективной параметризация, основанная на их зависимости от гранулометрического состава почвы и содержания органического вещества [Полуэктов, Терлеев, 2005; Rawls et al, 2003; Wösten et al, 2001].

Глава 5 посвящена верификации модели.

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

В таких случаях верификация модели может сводиться к:

- исследованию "географической" и "почвенно-микроклиматической" изменчивости параметров или определению пределов применимости их конкретного набора;

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

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

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

Что же касается влажности почвы: во-первых, в основу положена известная модель, многократно верифицированная как на зональном, так и (пусть в меньшей степени), на локальном уровне; во-вторых, достаточно показательным представляется непосредственное сравнение результатов расчетов с измеренными величинами. Такое сравнение произведено по материалам многолетних наблюдений Валдайского филиала ГГИ [Материалы ..., 1961-1989; Федоров, 1977]. Результаты расчетов приведены в таблице 2 и на рисунке 2.

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

20 t 5 1 D

W. .„..,%

25 20 1 5

1 0 -

Рисунок 2. Многолетняя динамика измеренных (1) и рассчитанных по модели (2) значений температуры на глубине 0.20 м (а) и весовой влажности почвы, средней в слое 0-1 м (б). Валдай, Новгородская обл., ельник-кисличник, почва подзолистая, суглинистая

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

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

Таблица 2.

30-летние средние значения и стандартные отклонения (с.к.о.) фактических и рассчитанных по модели значений температуры на глубине 0.20 м и весовой влажности

почвы, средней в слое 0-1 м (Валдай, Новгородская обл., ельник-кисличник, почва _подзолистая, суглинистая)_

Характер данных Величина Месяц

1 2 3 4 5 6 7 8 9 10 11 12

Средняя месячная температура почвы на глубине 0.20 м (°С)

Фактические данные 1957-1986 средняя с.к.о. - - - - 6.2* 10.1 12.5 12.6 9.7 6.0* - - - - - - 2.4* 1.1 0.9 0.7 0.8 1.3* - -

Расчет по модели (с фактической погодой) средняя с.к.о. -1.1 -1.0 -0.8 0.7 5.7 10.8 13.4 12.6 9.6 5.2 1.5 -0.3 0.6 0.3 0.3 0.9 1.2 1.0 1.2 1.0 1.0 1.4 0.7 0.5

Расчет по модели (с «генерированной» погодой) средняя С.К.О. -1.3 -1.2 -0.9 0.4 5.0 10.8 13.4 12.6 9.6 4.8 1.4 -0,3 1.2 0.7 0.5 1.3 1.6 1.4 1.4 0.9 0.8 1.4 0.7 0.7

Весовая влажность почвы в конце месяца, средняя в слое 0-1 м (%)

Фактические данные 1957-1986 средняя с.к.о. 18.3 17.9 18.8 19.7 16.7 15.2 14.3 13.7 14.8 16.1 17.3 17.9 1.9 1.7 1.6 2.0 1.8 1.8 2.2 3.2 2.9 1.7 1.8 1.2

Расчет по модели (с фактической погодой) средняя с.к.о. 17.6 17.6 17.7 21.1 18.3 15.8 14.6 14.2 14.9 17.2 17.7 17.7 1.8 1.8 1.9 1.1 1.6 2.1 1.9 1.7 1.7 1.6 1.9 1.9

Расчет по модели (с «генерированной» погодой) средняя с.к.о. 17.8 17.8 18.1 20.9 18.2 15.6 14.2 14.5 15.4 17.5 17.8 17.9 2.1 2.1 2.1 1.3 2.4 2.1 1.6 1.7 1.9 1.8 2.1 2.1

Примечания:

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

Глава 6: посвящена применению модели БСЫББ как составной части системы моделей ЕР1МОО [Котагоу е1 а!., 2003; Моделирование динамики..., 2007]. Модель неоднократно использовалась для исследования динамики органического вещества и азота в лесных экосистемах Европейской части России (в частности, Московской, Ленинградской, Тверской, Костромской областей, Республики Коми, Республики Марий Эл), ряда стран Северной и Центральной Европы (Финляндия, Швеция, Германия, Нидерланды, Чешская

Республика), а также центральной Канады (провинции Саскачеван и Манитоба). Многие результаты опубликованы [Моделирование динамики..., 2007; Chertov et al, 2004, 2009; Shaw et al, 2007; Надпорожская и др., 2009].

Один из наиболее показательных примеров использования модели - анализ влияния предполагаемых изменений климата в сочетании с различными уровнями выпадений соединений азота и сценариями лесопользования на динамику запасов органического вещества на примере сосновых лесов ОЛХ «Русский лес» (юг Московской области, почвы дерново-подзолистые, супесчаные и песчаные). Рассматривались 4 сценария лесопользования [Михайлов и др., 2007]: естественное развитие без рубок (А), выборочные рубки (Б), нормативные рубки (В) и имитация элементов самовольной рубки (Г); 4 варианта выпадения соединений азота (0, 6, 12 и 20кгК/га/год) и 2 варианта климатических сценариев: стационарный климат и «экстремальный» вариант глобального потепления (сценарий AlFi+HadCM3 [Mitchell et al, 2004\).

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

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

ГОДЫ

годы

2100

2000

годы

2000

годы

2050

2100

Рисунок 3. Разность между накоплением углерода в экосистеме при «экстремальном» потеплении по сравнению со стационарным климатом. Слева -при выпадениях азота 6 кг/га в год; справа - 12 кг/га в год. Вверху приведены данные для древостоев, внизу - для почвы. Л, Б, В, Г - разные сценарии лесопользования. Начало отсчета -1950 г

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

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

углерода и азота в лесных экосистемах БИМСЮ, должен позволить сравнение влияния изменений климата и внешних воздействий, таких как сплошные рубки и выпадения соединений азота из атмосферы.

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

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

4. Разработанная модель БСЫЗБ в составе системы моделей ЕР1МСЮ как источник сценариев климата почвы для моделировании круговорота углерода и азота в лесных экосистемах позволила промоделировать сукцессионные процессы при изменениях климата и антропогенных воздействиях на лесные экосистемы.

5. Вычислительные эксперименты с системой моделей продуктивности БИМС®, показали, что:

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

• «сильные воздействия» - рубки и пожары - приводят к значительным потерям органического вещества из экосистемы (как из надземного, так, в конечном счете, и из почвенного пулов), компенсация которых требует длительного времени; при повторяющихся «сильных воздействиях» запасы углерода в экосистеме не восстанавливаются полностью, и остаются существенно ниже, чем в ненарушенных лесах;

• влияние рубок и пожаров более существенно, чем эффект потепления климата; причем, чем интенсивнее рубки, тем меньше влияние как потепления климата, так и увеличения уровней поступления атмосферного азота (но при этом, больше деградация почв и меньше продуктивность древостоев).

Основные публикации по теме диссертации Статьи в рецензируемых журналах, рекомендованных ВАК:

1. Chertov О.G. Komarov A.S., Nadporozhskaya M.A., Bykhovets S.S., Zudin S.L. ROMUL - a model of forest soil organic matter dynamics as a substantial tool for forest ecosystem modelling // Ecological Modelling. 2001. Vol. 138. No. 1-3, p. 289308.

2. Быховец C.C, Комаров A.C. Простой статистический имитатор климата почвы с месячным шагом // Почвоведение. 2002. № 4, с. 443-452.

3. Chertov O.G., Komarov A.S., Bykhovets S.S., Kobak K.I. Simulated soil organic matter dynamics in forests of the Leningrad administrative area, northwestern Russia // Forest Ecology and Management. 2002. Vol. 169, No. 1-2. P. 29-44.

4. Komarov A., Chertov O., Zudin S., Nadporozhskaya M., Mikhailov A., Bykhovets S., Zudina E., Zoubkova E. EFIMOD 2 - a model of growth and cycling of elements in boreal forest ecosystems // Ecological Modelling. 2003. Vol. 170. No. 2-3, P. 373392.

5. Shaw C., Chertov O., Komarov A., Bhatti J., Nadporozskaya M., Apps M., Bykhovets S., Mikhailov A. Application of the forest ecosystem model EFIMOD 2 to Jack pine along the Boreal Forest Transect Case Study // Canadian Journal of Soil

Science. 2006. Vol. 86, No. 2, p. 171-185.

6. Быховец C.C., Сороковиков B.A., Мартуганов P.A, Мамыкин В.Г., Гиличинский Д.А. История наблюдений за температурой почвы на сети метеорологических станций России // Криосфера Земли. 2007. Т. 11, JV® 1, с. 720.

7. Larionova А.А., Yevdokimov I.V., Bykhovets S.S. Temperature response of soil respiration is dependent on readily decomposable C. Biogeosciences, 2007, Vol. 4, No. 6, 1073-1081.

8. Chertov O., Bhatti J.S., Komarov A., Mikhailov A., Bykhovets S. Influence of climate change, fire and harvest on the carbon dynamics of black spruce in Central Canada // Forest Ecology and Management. 2009. Vol. 257, No. 3, p. 941-950.

9. Надпорожская M.A., Цудлин П., Новак Ф., Быховец С.С., Чертов О.Г., Комаров А.С., Михайлов А.В. Применение математической модели ROMUL для анализа устойчивости почв ельников Крконоше в Чехии // Почвоведение. 2009, №6, с. 708-718.

Прочие публикации:

10. Быховец С.С. Статистический имитатор климата почвы // Новые подходы и методы в изучении природных и природно-хозяйственных систем. Доклады к Международной конференции. Алматы: Казак Университету 2000. С. 109-113. П.Комаров А.С., Чертов О.Г., Михайлов А.В., Надпорожская М.А., Быховец С.С., Зудин С.Л., Зудина Е.В., Зубкова Е.В. EFIMOD - система имитационных моделей циклов элементов в бореальных лесных экосистемах и ее область применения // Эмиссия и сток парниковых газов на территории Северной Евразии / Под ред. Н.П. Лаверова. Пущино: ОНТИ ПНЦ РАН, 2004. С. 167-172.

12. Быховец С.С. Климатологическая информация для моделирования баланса углерода в лесных экосистемах // там же. С. 190-195.

13. Bhatti J., Shaw S., Chertov О., Komarov A., Nadporozhskaya M., Apps M., Bykhovets S. Mikhailov A. Influence of climate change, fire and harvest on the soil С dynamics for Jack Pine in Central Canada: simulation approach with the EFIMOD

model // The 5~ European Conference on Ecological Modelling - ECEM 2005. Proceedings. Pushchino, 2005, P. 27-28.

14. Bykhovets S. Statistical modeling of soil climate in models of elements cycles in forest ecosystems // там же, P. 31-32.

15. Чертов О.Г., Комаров A.C., Надпорожская M.A., Михайлов A.B., БыховецС.С., Зудин С.Л., Зубкова Е.В. Динамическое моделирование процессов трансформации органического вещества почв. Имитационная модель ROMUL: Учебно-методич. пособие / Науч. ред. Б.Ф. Апарин. СПб.: СПбГУ, 2007.96 с.

16. Чертов О.Г., Быховец С.С., Надпорожская М.А., Комаров A.C., Ларионова A.A. Оценка скоростей трансформации органического вещества почвы в модели ROMUL // Моделирование динамики органического вещества в лесных экосистемах / Отв. ред. В.Н. Кудеяров. М,: Наука, 2007, с. 83-99.

17. Быховец С.С. Статистическое моделирование почвенного климата // там же, с. 105-131.

18. Михайлов A.B., Лукьянов А.М., Быховец С.С., Припутана И.В. Влияние уровней выпадения азота и климатических изменений на различные способы ведения лесного хозяйства // там же, с. 297-305.

19.Бхатти Дж., Шоу С., Чертов О.Г., Комаров A.C., Anne М., Михайлов A.B., Быховец С.С., Надпорожская М.А. Калибрация и применение модели EFIMOD для оценки совместного влияния изменения климата, рубок, пожаров и инвазий вредителей в бореальных лесах центральной Канады // там же, с. 314-323.

Подписано в печать: 24.09.2009

Заказ № 2582 Тираж -100 экз. Печать графарешая. Типография «11-й ФОРМАТ» ИНН 7726330900 115230, Москва, Варшавское ш., 36 (499) 788-78-56 www.autoreferat.ru

Содержание диссертации, кандидата географических наук, Быховец, Сергей Станиславович

ВВЕДЕНИЕ.

Глава 1. ОБЗОР СОВРЕМЕННОГО СОСТОЯНИЯ ПРОБЛЕМЫ.

1.1. К постановке задачи.

1.2. Фактические данные.

1.2.1. Материалы метеорологических наблюдений.

1.2.2. Экспериментальные исследования водно-теплового режима лесных почв.

1.3. Сценарии предполагаемых изменений климата.

1.4. Статистические имитаторы («генераторы») погоды.

1.5. Моделирование климата почвы.

1.5.1. Температурный режим почвы.

1.5.2. Режим влажности почвы.

1.5.3. Климатические субмодели в моделях динамики органического вещества почвы.

1.6. Краткое описание системы моделей ЕИМОО.

Глава 2. "ВХОДНЫЕ" МЕТЕОРОЛОГИЧЕСКИЕ РЯДЫ (ИЛИ

СОБСТВЕННО КЛИМАТИЧЕСКИЕ СЦЕНАРИИ).

2.1. Сценарии, основанные на фактических данных.

2.2. Статистическая имитация квазистационарных рядов температуры воздуха и осадков (имитатор "входных" метеорологических данных).

2.2.1. Обоснование выбора допущений и описание модели.

2.2.2. Оценка параметров.

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

Глава 3. СТАТИСТИЧЕСКАЯ МОДЕЛЬ ТЕМПЕРАТУРЫ ПОЧВЫ

3.1. Обоснование и описание модели.

3.2. Оценка параметров модели по данным метеорологических станций.

3.3. Разности температуры «лес - метеорологическая площадка»

Глава 4. ПРОСТАЯ БАЛАНСОВАЯ МОДЕЛЬ ВЛАЖНОСТИ ПОЧВЫ

4.1. Описание модели.

4.1.1. «Базовый» вариант модели.

4.1.2. Суммарная испаряемость (потенциальная эвапотранспирация).

4.1.3. Детализация описания лесных почв.

4.1.4. Детализация учета структуры эвапотранспирации лесных биогеоценозов.

4.2. Оценка параметров.

Глава 5. ВЕРИФИКАЦИЯ МОДЕЛЕЙ.

Глава 6. МОДЕЛЬНЫЕ ИССЛЕДОВАНИЯ КРУГОВОРОТА

УГЛЕРОДА И АЗОТА В ЛЕСНЫХ ЭКОСИСТЕМАХ.

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

6.2. Анализ динамики запасов углерода и азота в лесных экосистемах Костромской области при различных лесохозяйственных и климатических сценариях.

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

Введение Диссертация по наукам о земле, на тему "Почвенно-климатические сценарии в моделях круговорота углерода и азота в лесных экосистемах умеренного пояса"

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

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

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

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

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

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

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

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

Для достижения поставленной цели решались следующие задачи:

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

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

3. Разработка субмодели водно-теплового режима почвы.

4. Определение поправок для оценки температурного режима лесных почв в зависимости от типа леса.

5. Верификация моделей и их применение в рамках системы моделей биологического круговорота углерода и азота ЕРI МО 13 для сравнительного анализа последствий сильных внешних воздействий в разных климатических условиях и типах леса.

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

Практическая ценность результатов. Полученные результаты используются в научных исследованиях для разработки прогноза развития лесных почв при различных сценариях изменений климата и усиления антропогенной нагрузки. Разработанное программное обеспечение также используется в практикумах для студентов Пущинского государственного университета и Филиала МГУ имени М.В. Ломоносова в Пущино.

Личный вклад автора. Автором проанализированы возможности известных в мировой науке статистических генераторов метеорологических данных, разработаны простые статистические модели водно-теплового режима почв, учитывающих тип леса, принято активное участие в создании системы моделей биологического круговорота углерода и азота в лесных экосистемах ЕР1МСЮ, проведены вычислительные эксперименты и анализ результатов моделирования по теме диссертации.

Связь темы диссертации с плановыми исследованиями. Тема диссертации связана с основным планом научно-исследовательских работ Института физико-химических и биологических проблем почвоведения РАН по теме «Математическое моделирование циклов элементов и сукцессионных процессов в системе почва - растительность — атмосфера (№ 01.2 006 073924.1.16.11)».

Основные положения диссертации, выносимые на защиту:

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

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

3. Верифицированная модель динамики углерода и азота в лесных экосистемах ЕЕ1МСЮ, включающая в себя модель динамики органического вещества почвы 1ЮМиЬ, позволяют, используя разработанные модели водно-теплового режима почв, провести сравнительное исследование роли климатических изменений и сильных внешних воздействий (сплошные рубки, выпадения соединений азота из атмосферы) на лесные экосистемы.

Апробация и внедрение результатов. Результаты работы были представлены на: I Всероссийской научной школе "Математические методы в экологии" (Петрозаводск, 2001), 2~ Международной конференции «Эмиссия и сток парниковых газов на территории Северной Евразии» (Пущино, 2003), Всероссийской конференции «Биосферные функции почвенного покрова» поев. 100-летию со дня рождения чл.-корр. АН СССР В.А. Ковды. (Пущино, 2005), 5й Европейской конференции по математическому моделированию в экологии (ЕСЕМ-05, Pushchino, 2005), а также на рабочих совещаниях по проектам ИНТАС (97-20355 "FORMOD", 01-0512 "PODZOL", 01-0633 "SILVICS"), EU FP6 INCO (013388 "OMRISK") в Финляндии, Великобритании и Италии, семинарах Лаборатории моделирования экосистем ИФХиБПП РАН и Учебного Центра почвоведения, экологии и природопользования Пущинского государственного университета.

Модель SCLISS [Быховец, Комаров, 2002; Быховец, 2007] является составной частью системы моделей EFIMOD-2 [Komarov et ai, 2003], разработанной в Лаборатории моделирования экосистем ИФХиБПП РАН и Лаборатории биохимии почв Биологического НИИ СПбГУ; полученные на ее основе сценарии использовались и используются при моделировании динамики органического вещества в лесных экосистемах Европейской части России, ряда стран Северной и Центральной Европы (Финляндии, Швеции, Германии, Чешской Республики) и Канады.

Публикации. Результаты работы изложены в 24 публикациях, в числе которых 9 статей в рецензируемых журналах, рекомендованных ВАК, 4 главы в монографиях, учебно-методическое пособие, 7 статей в сборниках научных трудов и 3 тезисов докладов:

Статьи в рецензируемых журналах, рекомендованных ВАК:

1. Chertov O.G., Komarov A.S., Nadporozhskaya М.А., Bykhovets S.S., Zudin S.L. ROMUL - a model of forest soil organic matter dynamics as a substantial tool for forest ecosystem modelling // Ecological Modelling, 2001. Vol. 138. No. 1/3. P. 289-308.

2. Быховец C.C., Комаров A.C. Простой статистический имитатор климата почвы с месячным шагом // Почвоведение, 2002. № 4. С. 443-452.

3. Chertov O.G., Komarov A.S., Bykhovets S.S., Kobak K.I. Simulated soil organic matter dynamics in forests of the Leningrad administrative area, northwestern Russia // Forest Ecology and Management, 2002. Vol. 169. No. 1/2.

P.29-44.

4. Komarov A., Chertov O., Zudin S., Nadporozhskaya M., Milchailov A., Bykhovets S., Zudina E., Zoubkova E. EFIMOD 2 - a model of growth and cycling of elements in boreal forest ecosystems // Ecological Modelling, 2003. Vol. 170. No. 2/3. P. 373-392.

5. Shaw C., Chertov O., Komarov A., Bhatti J., Nadporozskaya M., Apps M., Bykhovets S., Mikhailov A. Application of the forest ecosystem model EFIMOD 2 to Jack pine along the Boreal Forest Transect Case Study // Canadian Journal of Soil Science, 2006. Vol. 86. No. 2. P. 171-185.

6. Быховец С.С., Сороковиков В.А., Мартуганов Р.А, Мамыкин В.Г., Гиличинский ДА. История наблюдений за температурой почвы на сети метеорологических станций России // Криосфера Земли, 2007. Т. 11. № 1. С. 7-20.

7. Larionova А.А., Yevdokimov I.V., Bykhovets S.S. Temperature response of soil respiration is dependent on readily decomposable С // Biogeosciences, 2007. Vol. 4. No. 6. C. 1073-1081.

8. Chertov O., Bhatti J.S., Komarov A., Milchailov A., Bykhovets S. Influence of climate change, fire and harvest on the carbon dynamics of black spruce in Central Canada // Forest Ecology and Management, 2009. Vol. 257. No. 3. P. 941-950.

9. Р1адпорожская M.A., Цудлин П., Р1овак Ф., Быховец С.С., Чертов О.Г., Комаров А.С., Михайлов А.В. Применение математической модели ROMUL для анализа устойчивости почв ельников Крконоше в Чехии // Почвоведение, 2009. № 6. С. 708-718.

Монографии и главы в монографиях:

10. Чертов О.Г., Комаров А.С., Надпорожская М.А., Михайлов А.В., Быховец С.С., Зудин C.JI., Зубкова Е.В. Динамическое моделирование процессов трансформации органического вещества почв. Имитационная модель ROMUL: Учебно-методическое пособие / Р1ауч. ред. Б.Ф. Апарин. СПб.: СПбГУ, 2007. 96 с.

11. Чертов О.Г., Быховец С.С., Р1адпорожская М.А., Комаров А.С.,

Ларионова A.A. Оценка скоростей трансформации органического вещества почвы в модели ROMUL // Моделирование динамики органического вещества в лесных экосистемах / Отв. ред. В.Н. Кудеяров. М.: Наука, 2007. С. 83-99.

12. Быховец С.С. Статистическое моделирование почвенного климата // там эюе. С. 105-131.

13. Михайлов A.B., Лукьянов A.M., Быховец С.С., Прйпутина И.В. Влияние уровней выпадения азота и климатических изменений на различные способы ведения лесного хозяйства // там же. С. 297-305.

14. Бхатти Дж., Шоу С., Чертов О.Г., Комаров A.C., Anne М., Михайлов A.B., Быховец С.С., Надпорожская М.А. Калибрация и применение модели EFIMOD для оценки совместного влияния изменеиия климата, рубок, пожаров и инвазий вредителей в бореальных лесах центральной Канады // там же. С. 314-323.

Статьи в сборниках научных трудов;

15. Быховец С.С. Статистический имитатор климата почвы // Новые подходы и методы в изучении природных и природно-хозяйственных систем. Доклады к Международной конференции. Алматы: Казак Университет!, 2000. С. 109-113.

16. Chertov O.G., Komarov A.S., Nadporozhskaya M., Bykhovets S.S., Zudin S.L. Simulation study of nitrogen supply in boreal forests using model of soil organic matter dynamics ROMUL // Plant nutrition - Food security and sustainability of agro-ecosystems through basic and applied research / Ed. by W.J. Horst et al. (Developments in Plant and Soil Sciences. Vol. 92). Dordrecht: Kluwer, 2001. P. 900-901.

17. Комаров A.C., Чертов О.Г., Михайлов A.B., Надпорожская M.A., Быховец С.С., Зудин С.Л., Зудина Е.В., Зубкова E.B. EFIMOD - система имитационных моделей циклов элементов в бореальных лесных экосистемах и ее область применения //' Эмиссия и сток парниковых газов на территории Северной Евразии / Под ред. Н.П. Лаверова. Пущино: ОНТИ ПНЦ РАН,

2004. С. 167-172.

18. Быховец С.С. Климатологическая информация для моделирования баланса углерода в лесных экосистемах // там же. С. 190-195.

19. Bhatti J., Shaw S., Chertov О., Komarov A., Nadporozhskaya M., Apps M., Bykhovets S., Mikhailov A. Influence of climate change, fire and harvest on the soil С dynamics for Jack Pine in Central Canada: simulation approach with the EFIMOD model // The 5th European Conference on Ecological Modelling - ECEM

2005. Proceedings. Pushchino, Russia, 2005. P. 27-28.

20. Bykhovets S. Statistical modeling of soil climate in models of elements cycles in forest ecosystems // там же. P. 31—32.

21. Chertov O., Bhatti J., Komarov A., Apps M., Mikhailov A., Bykhovets S. Difference of ecological strategies of coniferous tree species in Canadian and European boreal forests: simulation modelling analysis // The 6th European Conference on Ecological Modelling (ECEM '07). 2007. Conference Proceedings. Trieste, Italy, 2007. P. 103-104.

Тезисы докладов (основные):

22. Быховец С.С. Статистическое моделирование характеристик климата почвы // Математические методы в экологии. Тезисы докладов Всероссийской научной школы. Петрозаводск, 2001. С. 204—206.

23. Быховец С.С. Климатические сценарии для задач моделирования динамики органического вещества в лесных экосистемах // Конференция «Биосферные функции почвенного покрова» поев. 100-летию со дня рождения чл.-корр. АН СССР В.А. Ковды. Тезисы докладов. Пущино, 2005. С. 15-16.

24. Loukianov A.M., Bykhovets S.S. Carbon exchange under climate changes in the forest soils of the Russian North-West using results of simulation modelling // Climate Changes and their Impact on Boreal and Temperate Forests. Abstracts of International Conference. Ekaterinburg, 2006. P. 59.

Структура и объем диссертации. Диссертация состоит из введения, 6 глав, заключения и библиографии. Она включает 148 страниц текста,

Заключение Диссертация по теме "Метеорология, климатология, агрометеорология", Быховец, Сергей Станиславович

Основные результаты и выводы настоящей работы могут быть сформулированы следующим образом:

1. На основании проведенного анализа моделей круговорота углерода и азота в лесных экосистемах сформулированы требования к разрабатываемой модели климата почвы: быть достаточно простой, иметь месячный шаг по времени, и решать задачу преобразования входных климатических рядов в температуру и влажность лесной подстилки и почвы и, при необходимости, генерировать сами входные ряды с заданными статистическими свойствами; также климатический блок, будучи включенным в систему моделей биологического круговорота углерода и азота в лесных экосистемах EFIMOD, должен позволить сравнение влияния изменений климата и внешних воздействий, таких как сплошные рубки и выпадения соединений азота из атмосферы;.

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

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

4. Разработанная модель ЗСЫББ в составе системы моделей ЕР1МСЮ как источник сценариев климата почвы для моделировании круговорота углерода и азота в лесных экосистемах позволила промоделировать сукцессионные процессы при изменениях климата и антропогенных воздействиях на лесные экосистемы.

5. Вычислительные эксперименты с системой моделей продуктивности ЕР1МСЮ, показали, что:

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

• «сильные воздействия» - рубки и пожары — приводят к значительным потерям органического вещества из экосистемы (как из надземного, так, в конечном счете, и из почвенного пулов), компенсация которых требует длительного времени; при повторяющихся «сильных воздействиях» запасы углерода в экосистеме не восстанавливаются полностью, и остаются существенно ниже, чем в ненарушенных лесах;

• влияние рубок и пожаров более существенно, чем эффект потепления климата; причем, чем интенсивнее рубки, тем меньше влияние как потепления климата, так и увеличения уровней поступления атмосферного азота (но при этом, больше деградация почв и меньше продуктивность древостоев).

ЗАКЛЮЧЕНИЕ

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

1. Аблеева В.А., Терешонок H.A. Метеорологическая характеристика сезонов года и динамика климата в Приокско-Террасном заповеднике // Экосистемы Приокско-Террасного биосферного заповедника: Сб. науч. тр. Пущино, 2005. С. 18-34.

2. Абражко В.И., Говоренков Б.Ф., Пономарева В.В. Водно-физические и химические свойства почв // Структура и продуктивность еловых лесов южной тайги. JL: Наука, 1973. С. 63-70.

3. Абражко В.И. Водный режим сообществ еловых лесов (синэкологические исследования). Дисс. на соискание ученой степени докт. биол. наук. СПб., 1998. 620 с.

4. Агрометеорологические ежегодники по Архангельской области и Коми АССР за 1981-1988 гг. Архангельск, 1982-1989.

5. Адаменко В.Н. Мелиоративная микроклиматология. Л.: Гидрометеоиздат, 1979. 184 с.

6. Александрова Л. Н. Органическое вещество почвы и процессы его трансформации. Л.: Наука, 1980. 287 с.

7. Алексеев В.А., Марков М.В. Статистические данные о лесном фонде и изменение продуктивности лесов России во второй половине XX века. СПб.: Санкт- Петербургский лесной экологический центр, 2003. 273 с.

8. Алисов Б.П. Климат СССР. М.: Высшая школа, 1969. 104 с.

9. Алифанов В.М., Лошакова H.A. Водный режим серых лесных почв // Почвоведение, 1981. № 4. С. 58-70.

10. Анисимов O.A. Оценка влияния ожидаемых изменений климата на режим вечной мерзлоты // Метеорология и Гидрология, 1990. № 3. С. 40-46.

11. Барталев С.А., Ершов Д.В., Исаев A.C., Потапов П.В., Турубанова С.А., Ярошенко А.Ю. Леса России: доминирующие группы древесных пород и сомкнутость древесного полога. Карта. Масштаб. 1:14000000. М., 2004. http://forestforum.ru/info/pictures/rusmap.pdf

12. Бельченко Г.Г. Разработка статистического имитатора погоды // Науч.-тех. бюлл. по агрон. физике, 1989. № 76. С. 51-56.

13. Биопродукционный процесс в лесных экосистемах Севера. СПб.: Наука, 2001. 278 с.

14. Бондарик H.JL, Карпечко Ю.В. Сравнение методов определения среднемноголетней величины испарения с леса // Метеорология и гидрология, 1999. № 9. С. 98-105.

15. Браславская Т.Ю., Тихонова Е.В. Оценка биоразнообразия южнотаежных лесов на северо-востоке Костромской области // Лесоведение, 2006. № 2. С. 34-50.

16. Будыко М.И. Тепловой баланс земной поверхности. JL: Гидрометеоиздат, 1956. 255 с.

17. Будыко М.И. Климат и жизнь. JL: Гидрометеоиздат, 1971. 470 с.

18. Быховец С.С., Комаров A.C. Простой статистический имитатор климата почвы с месячным шагом // Почвоведение, 2002. № 4. С. 443-452.

19. Быховец С.С., Сороковиков В.А., Мартуганов Р.А, Мамыкин В.Г., Гиличинский Д.А. История наблюдений за температурой почвы на сети метеорологических станций России // Криосфера Земли, 2007. Т. XI. № 1. С. 7-20.

20. Васильев И.С. Водный режим подзолистых почв // Материалы по изучению водного режима почв / Труды Почвенного ин-та им. В.В. Докучаева. Т. XXXII. М., 1950. С. 74-296.

21. Васильев И.С. Водный режим дерново-подзолистых почв под лесом и пашней // Современные почвенные процессы в лесной зоне Европейской части СССР. М.: Изд. АН СССР, 1959. С. 3^9.

22. Верхоланцева JI.A. Почвы Усть-Куломского топоэкологического профиля // Вопросы экологии леса (средняя тайга Коми АССР) / Труды Коми филиала АН СССР, № 19. Сыктывкар: Коми кн. изд-во, 1970. С. 95-106.

23. Верхоланцева Jl.А. Водно-физические свойства почв сосняков-зеленомошников // Вопросы экологии сосняков Севера / Труды Коми филиала АН СССР, № 24. Сыктывкар: Коми кн. изд-во, 1972а. С. 42-51.

24. Верхоланцева Л.А. Торфянисто-подзолистые почвы различных типов леса среднетаежной подзоны Коми АССР // Материалы по почвам Коми АССР. Сыктывкар, 19726. С. 28-37.

25. Верхоланцева Л.А. Почвы северотаежных ельников // Экология ельников Севера / Труды Коми филиала АН СССР, № 32. Сыктывкар: Коми кн. изд-во, 1977. С. 52-84.

26. Верхоланцева Л.А., Бобкова К.С. Почвы и растительность сосновых лесов Зеленоборского стационара // Вопросы экологии сосняков Севера / Труды Коми филиала АН СССР, № 24. Сыктывкар: Коми кн. изд-во, 1972. С. 6-20.

27. Винников К.Я., ЛемеШко H.A., Сперанская H.A. Влагосодержание почвы и. сток внетропической части Северного полушария при глобальном потеплении // Метеорология и Гидрология, 1990. № 3. С. 5-10.

28. Волобуев В.Р. Соотношение между тепловым режимом почв и климатом приземного слоя воздуха // Почвоведение, 1983. № 2. С. 52-63.

29. Воронина Т.В. Режим влажности дерново-подзолистых почв Валдайской возвышенности // Почвоведение, 1975. № 2. С. 67-81.

30. Воронков H.A., Кожевникова С.А., Шомполова В.А. Температурный режим почв под лесом и залежью в условиях Подмосковья // Почвоведение, 1979. № 6. С. 90-99.

31. Газизуллин А.Х., Хасаншин Б.Д. Бурые лесные почвы широколиственных и хвойно-широколиственных лесов Среднего Поволжья на двучленных наносах//Почвоведение, 1980. № 11. С. 117-130.

32. Галенко Э.П. Фитоклимат и энергетические факторы продуктивности хвойного леса Европейского Севера. Л.: Наука, 1983. 128 с.

33. Гневко Г.Т., Исаев .A.A. Обобщение рекомендаций по устранению неоднородностей в рядах температуры воздуха и осадков при обработке данных и создании автоматизированных архивов. Обнинск, 1980. 40 с.

34. Голубятников Л.Л., Ямагата Й., Александров Г.А. Стохастический генератор величин месячных осадков и среднемесячных температур // Известия АН. Физика атмосферы и океана, 2002. Т. 38. № 1. С. 47-55.

35. Гусев Е.М., Насонова О.Н. Параметризация тепло- и влагообмена на поверхности суши при сопряжении гидрологических и климатических моделей // Водные ресурсы, 1998. Т. 25. № 4. С. 421-431.

36. Гусев Е.М., Насонова О.Н. Параметризация процессов тепло- и влагообмена в экосистемах бореальных лесов // Известия АН. Физика атмосферы и океана, 2001. Т.'37. № 2. С. 182-200.

37. Гусев Е.М., Насонова О.Н. Методика оценки динамики водного и углеродного балансов -экосистемы хвойного леса // Известия РАН: Физика атмосферы и океана, 2007. Т. 43. № 1. С. 81-92.

38. Данные по климату СССР. Т. 1, 2. Обнинск, 1977.

39. Демченко П.Ф. Динамико-стохастическая модель формирования поверхностного стока // Известия АН. Физика атмосферы и океана, 2003. Т. 39. №2. С. 186-192.

40. Димо В.Н. Расчетный метод определения температуры почвы // Бюлл. Почвенного ин-та им. В.В. Докучаева, 1967. Вып. 1. С. 88-99.

41. Димо В.Н. Тепловой режим почв СССР. М.: Колос, 1972. 360 с.

42. Димо В.Н. Сравнительная характеристика климата дерново-подзолистых, лессивированных и бурых лесных почв Европы // Почвоведение, 1981. № 7. С. 67-74.

43. Дюшофур Ф. Основы почвоведения. М.: Мир, 1970.

44. Елагин И.Н., Изотов В.Ф. Температура почвы в сосняках лесной зоны в различные периоды года//Почвоведение, 1968. № 6. С. 138-142.

45. Еруков Г.В., Власкова Г.В. Гидротермический режим сосновых лесов Карелии. JL: Наука, 1986. 110 с.

46. Жуковский Е.Е., Киселева Т.Л., Мандельштам С.М. Статистический анализ случайных процессов в приложении к агрофизике и агрометеорологии. Л.: Гидрометеоиздат, 1976. 407 с.

47. Забоева И.В. Почвы и земельные ресурсы Коми АССР. Сыктывкар: Коми книжн. изд-во, 1975. 344 с.

48. Заугольнова Л.Б.(ред.) Оценка и сохранение биоразнообразия лесного покрова в заповедниках Европейской России. М.: Научный мир, 2000.

49. Зонн C.B. Водный режим почв дубовых лесов // Труды Ин-та леса. T. VII. М., 1951. С. 9-84.

50. Зубенок Л.И. Испарение на континентах. Л.: Гидрометеоиздат, 1976. 251 с.

51. Иванов Л.А., Силина A.A., Жмур Д.Г., Цельникер Ю.Л. Об определении транспирационного расхода древостоем леса // Ботанический журнал, 1951. Т. 36. № 1. С. 5-20.

52. Изменение климата, 2001 г. Обобщенный доклад. МГЭИК, 2003.

53. Исаев A.A. Статистика в метеорологии и климатологии. М.: Изд-во Моск. ун-та, 1988. 245 с.

54. Исаев A.A. Атмосферные осадки. Часть 1. Изменчивость характеристик осадков на территории России и сопредельных стран. М.: Изд-во Моск. унта, 2002. 192 с.

55. Казакевич Д.И. Основы теории случайных функций в задачах гидрометеорологии. JL: Гидрометеоиздат, 1989.

56. Казимиров Н.И., Морозова P.M., Куликова В.К. Органическая масса и потоки веществ в березняках средней тайги. Л.: Наука, 1978. 216 с.

57. Кайбияйнен Л.К. Сбалансированность системы водного транспорта у сосны обыкновенной. IV. Общие характеристики водного режима в разных экологических условиях// Лесоведение, 1986. № 4. С. 70-75.

58. Карпечко Ю.В. Оценка пространственной и временной неоднородности задержания жидких осадков пологом леса // Лесоведение, 1997. № 4. С. 64-70.

59. Кашанский А.Д. Подзолистые почвы на двучленных наносах Европейского Севера СССР // Современные почвенные процессы. Сб. статей. М., 1974. С. 112-135.

60. Кислов A.B. Стохастическая модель пространственнотвременного распределения влажности почвы на территории СССР // Метеорология и гидрология, 1991. № 3. С. 101-107.

61. Кислов A.B. Учет изменчивости начального состояния в стохастической модели влажности почвы // Метеорология и гидрология, 1991. № 8. С. 109— 111.

62. Кислов A.B. Климат в прошлом, настоящем и будущем. М.: МАИК "Наука/Интерпериодика", 2001. 351 с.

63. Климатология / под ред. O.A. Дроздова и Н.В. Кобышевой. Л.: Гидрометеоиздат, 1989. 568 с.

64. Климат России / под ред. Н.В. Кобышевой. Л.: Гидрометеоиздат, 2001. 655 с.

65. Климатологический справочник СССР: Метеорологические данные за отдельные годы. Вып. 1-27. Части I, II, VII. Л.: Гидрометеоиздат, 19531964.

66. Книзе A.A., Крестовский О.И. Определение влияния структуры лесного фонда на водность рек (на примере бассейнов озер Ладожское и Ильмень). Методические рекомендации. СПб.: СПбНИИЛХ, 1993. 73 с.

67. Книзе A.A., Романюк Б.Д., Федорчук В.Н. Динамика продуктивности таежных и подтаежных лесов европейской части России за длительный период // Труды СПбНИИЛХ. СПб, 2000. Вып. 1 (2). С. 161-169.

68. Кобышева Н.В., Наровлянский Г.Я. Климатологическая обработка метеорологической информации. Л.: Гидрометеоиздат, 1978. 295 с.

69. Кононенко A.B. Гидротермический режим таежных и тундровых почв Европейского Северо-Востока. JL: Наука, 1986. 145 с.

70. Крестовский О.И. Влияние вырубок и восстановления лесов на водность рек. JL: Гидрометеоиздат, 1986. 118 с.

71. Кудрявцев В.А. (ред.) Мерзлотоведение. М.: Изд. Моск. Ун-та, 1981. 239 с.

72. Куртенер Д.А., Чудновский А.Ф. Агрометеорологические основы тепловой мелиорации почв. JL: Гидрометеоиздат, 1979. 231 с.

73. Кучмент JI.C., Демидов В.Н. Моделирование влияния гидрологических процессов на углеродный цикл лесной экосистемы // Метеорология и гидрология, 2004. № 12. С. 71-81.

74. Луговая Д.Л. Разнообразие растительных сообществ после пожаров и рубок в лесах Костромской области // Лесоведение, 2008. № 4. С. 34-43.

75. Лучшев A.A. Данные по влиянию леса на температуру почвы // Труды ВНИИЛХ, 1940, вып. 18. С. 279-287.

76. Мамихин С.В. Воспроизведение температурного и гидрологического режимов почвы в математических моделях сухопутных экосистем // Вестн. Моск. ун-та. Сер. 15. Почвоведение, 1997. № 3. С. 7-10.

77. Материалы наблюдений Валдайского филиала ГГИ. Вып. 9-38 за 1957— 1986 гг. Валдай, 1961-1989.

78. Материалы наблюдений Научно-исследовательской гидрометеорологической обсерватории "Каменная Степь". Вып. 7—25 за 1970-1988 гг. Курск, 1974-1989.

79. Метеорологический ежемесячник. 1965-1990. Вып. 1, 3, 8, 29. Часть 2. №№ 1-13.

80. Методы изучения и расчета водного баланса. Л., Гидрометеоиздат, 1981. 397 с.

81. Методы расчета водных балансов. Международное руководство по исследованиям и практике // Исследования и доклады по гидрологии. Вып. 17 / Под ред. A.A. Соколова и Т.Г. Чапмена. JL: Гидрометеоиздат, 1976. 120 с.

82. Моделирование динамики органического вещества в лесных экосистемах / Отв. ред. В.Н. Кудеяров. М.: Наука, 2007. 380 с.

83. Молчанов A.A. Гидрологическая роль сосновых лесов на песчаных почвах. М.: Изд. АН СССР, 1952. 487 с.

84. Молчанов A.A. Гидрологическая роль леса. М.: Наука, 1960. 487 с.

85. Молчанов A.A. Лес и климат. М.: Наука, 1961. 278 с.

86. Надпорожская М.А., Цудлин П., Новак Ф., Быховец С.С., Чертов О.Г., Комаров A.C., Михайлов A.B. Применение математической модели ROMUL для анализа устойчивости почв ельников Крконоше в Чехии // Почвоведение, 2009. № 6. С. 708-718.

87. Научно-прикладной справочник по климату СССР. Сер. 3. Вып. 1-34. Л.: Гидрометеоиздат, 1988-1992.

88. Павлов A.B. Теплообмен почвы с атмосферой в северных и умеренных широтах территории СССР. Якутск: Якутское книжное издательство, 1975. 302 с.

89. Палагин Э.Г. Математическое моделирование агрометеорологических условий перезимовки озимых культур. Л.: Гидрометеоиздат, 1981. 191 с.

90. Пенман X. Растения и влага. Л.: Гидрометеоиздат, 1968. 162 с.

91. Перунова М.И. Температура почвы // Климат Союза ССР. Ч. VI. Л.: Гидрометеоиздат, 1952. 349 с.

92. Подзолистые почвы Запада Европейской части СССР. М.: Колос, 1977. 288 с.

93. Подзолистые почвы Северо-Запада Европейской части СССР. М.: Колос,1979. 256 с.

94. Подзолистые почвы Центральной и Восточной частей Европейской территории СССР (на суглинистых почвообразующих породах). Л.: Наука,1980. 301 с.

95. Подзолистые почвы Центральной и Восточной частей Европейской территории СССР (на песчаных почвообразующих породах). Л.: Наука,1981. 199 с.

96. Полуэктов P.A. Динамические модели агроэкосистемы. Л.: Гидрометеоиздат, 1991. 312 с.

97. Полуэктов P.A., Терлеев В.В. Моделирование водоудерживающей способности почвы с использованием агрогидрологических характеристик // Метеорология и гидрология, 2005. № 12. С. 98-103.

98. Почвы зоны переброски части стока северных рек. Л.: Наука, 1983. 168 с.

99. Работнов Т.А. Азот в наземных биогеоценозах // Структурно-функциональная организация биогеоценозов. М.: Наука, 1980. С. 69-90.

100. Растворова О.Г. Физика почв (практическое руководство). Л.: Изд. Ленингр. ун-та, 1983. 194 с.

101. Рахманов В.В. Гидроклиматическая роль лесов. М.: Лесная Промышленность, 1984. 241 с.

102. Рудаков В.Е. О независимости соотношения между величинами транспирации и массы хвои у сосны обыкновенной от погодных условий // Изв. ВГО, 1979. Т. 111. Вып. 2. С. 168-170.

103. Руднева E.H. Лесорастительные свойства таежных поверхностно-глееватых почв ельников черничных // Влияние хозяйственныхгмероприятий на лесные почвы Карелии. Петрозаводск, 1983. С. 35-77.

104. Рысин Л.П., Савельева Л.И. Еловые леса России. М.: Наука, 2002.

105. Семко А.П. Гидротермический режим почв лесной зоны Кольского полуострова. Апатиты: Кольский филиал АН СССР, 1982. 142 с.

106. Сиротенко О.Д. Математическое моделирование водно-теплового режима и продуктивности агроэкосистем. Л.: Гидрометеоиздат, 1981, 167 с.

107. Смагин. В., Садовникова Н.Б., Смагина М.В., Глаголев М.В., Шевченко Е.М., Хайдапова Д.Д., Губер А.К. Моделирование динамики органического вещества почв. М.: Изд. МГУ, 2001. 120 с.

108. Смирнова О.В., Восточно-европейские леса, их история в голоцене и современность. Т. 1, 2. М.: Наука, 2004.

109. Справочник по климату СССР. Вып. 1-34. Ч. II, IV. Л.: Гидрометеоиздат, 1964-1969.

110. Справочник по климату СССР: Метеорологические данные за отдельные годы. Вып. 1-34. Ч. I, II, VII, VIII. 1969-1978.

111. Справочник по климату СССР: Устойчивость и точность климатических характеристик. Т. 1,2. Л.: Гидрометеоиздат, 1976.

112. Структура и продуктивность еловых лесов южной тайги. Л.: Наука, 1973. 312 с.

113. Структурно-функциональная организация почв и почвенного покрова Европейского Северо-Востока. СПб.: Наука, 2001. 224 с.

114. Сценарии выбросов. Специальный доклад. МГЭИК. 2000.

115. Топаж А.Г. Моделирование суточных метеоданных как входного сигнала модели продукционного процесса // Почва и растение процессы и модели. Сб. науч. трудов. СПб.: АФИ, 1992. С. 79-86.

116. Федорец Н.Г., Бахмет О.Н. Экологические особенности трансформации соединений углерода и азота в лесных почвах. Петрозаводск: Карельский научный центр РАН, 2003. 240 с.

117. Федоров С.Ф. Исследование элементов водного баланса в лесной зоне Европейской территории СССР. Л.: Гидрометеоиздат, 1977. 264 с.

118. Федоров С.Ф. Определение транспирации еловыми насаждениями по коэффициенту транспирационной активности // Экспериментальные гидрологические исследования на Валдае / Труды ГГИ, вып. 279. Л.: Гидрометеоиздат, 1981. С. 32-34.

119. Фридланд В.М. Почвы Приокско-Террасного государственного заповедника // Труды Почвенного института АН СССР. Т. 46. М.: Изд-во АН СССР, 1955. С. 136-190.

120. Хильми Г.Ф. Энергетика и продуктивность растительного покрова суши. ( Л.: Гидрометеоиздат, 1976. 62 с.

121. Чертов О.Г. Математическая модель экосистемы одного растения // Журнал общей биологии, 1983. Т. 44. С. 406-414.

122. Чертов О.Г. Имитационная модель минерализации и гумификации лесного опада и подстилки. Журнал общей биологии 1985. Т. 46. № 6. С. 794-804.

123. Чертов О.Г., Комаров A.C. Имитационная модель динамики органического вещества почвы // Вестник Санкт-Петербургского университета, 1996. Сер. 3. Вып. 1. С. 104-109.

124. Честных О.В., Замолодчиков Д.Г, Уткин А.И., Коровин Г.Н. Распределение запасов органического углерода в почвах лесов России // Лесоведение, 1999. № 2. С. 13-21.

125. Честных О.В., Замолодчиков Д.Г. Зависимость плотности почвенных горизонтов от глубины их залегания и содержания гумуса // Почвоведение, 2004. № 8. С. 937-944.

126. Чумаченко С.И. Моделирование динамики многовидовых разновозрастных лесных ценозов // Журнал общей биологии, 1998. Т. 59. № 4. С. 363-376.

127. Швер Ц.А. Атмосферные осадки на территории СССР. JL: Гидрометеоиздат, 1976. 302 с.

128. Швер Ц.А. Закономерности распределения количества осадков на континентах. Д.: Гидрометеоиздат, 1984. 285 с.

129. Шеин Е.В. Курс физики почв. М.: Изд-во Моск. ун-та, 2005. 432 с.

130. Шкадова А.К. Температурный режим почв на территории СССР. Д.: Гидрометеоиздат, 1979. 240 с.

131. Шмакин А.В., Попова В.В. Влияние североатлантического колебания на многолетний гидротермический режим Северной Евразии. II. Моделирование внутривековых колебаний теплового и водного балансов // Метеорология и гидрология, 2003. № 6, С. 59-68.

132. Шульгин A.M. Температурный режим почвы. Д.: Гидрометеоиздат, 1957. 242 с.

133. Шульгин A.M. Климат почвы и его регулирование. Д.: Гидрометеоиздат, 1972. 341 с.

134. Aber J.D., Melillo J.M. FORTNITE: a computer model of organic matter and nitrogen dynamics in forest ecosystems // University of Wisconsin Research Bulletin, 1982. No R3130.

135. Alban D.H., Laidly P.R. Generalized biomass equations for jack and red pine in the Lake States Pinus banksiana, USA // Canadian Journal of Forest Research, 1982. Vol. 12. No. 4. P. 913-921.

136. Arp P.A., Yin X.I. Predicting water fluxes through forests from monthly precipitation and mean monthly air temperature records // Canadian Journal of Forest Research, 1992. Vol. 22. No. 6. P. 864-877.

137. Bergstrom S. Development and application of a conceptual runoff model for Scandinavian catchments. SMHI Report No. RH7. Norrkopping, 1976.

138. Blaney H.F., Criddle W.D. Determining water requirements in irrigated areas from climatological and irrigation data / USDA Soil Conservation Service. Techn. Paper No. 96. 1950.

139. Bocock K.L., Jeffers J.N.R., Lindley D.K., Adamson J.K., Gill C.A. Estimating woodland soil temperature from air temperature and other climatic variables // Agricultural Meteorology, 1977. Vol. 18. P. 351-372.

140. Brady N.C. The nature and properties of soils. McMillan: New York, 1990.

141. Burns R.M., Honkala B.H. Silvics of North America: 1. Conifers; 2. Hardwoods. Agriculture Handbook 654. USDA, Forest Service, 1990. Washington, DC. Available at:http://www.na.fs.fed.us/spfo/pubs/silvicsmanual/tableofcontents.htm

142. Chertov O.G., Komarov A.S. SOMM a model of soil organic matter dynamics // Ecological Modelling, 1997. Vol. 94. No. 2/3. P. 177-189.

143. Chertov O.G., Komarov A.S., Karev G.P. Modern Approaches in Forest Ecosystem Modelling / European Forest Institute Research Report No. 8. Leiden: Brill, 1999a. 130 p.

144. Chertov O.G., Komarov A.S., Tsiplianovsky A.M. A combined simulation model of Scots pine, Norway spruce and Silver birch ecosystems in the European boreal zone // Forest Ecology and Management, 1999b. Vol. 116. No. 1/3. P. 189-206.

145. Chertov O.G., Komarov A.S., Nadporozhskaya M.A., Bykhovets S.S., Zudin S.L. ROMUL — a model of forest soil organic matter dynamics as a substantial tool for forest ecosystem modeling // Ecological Modelling, 2001. Vol. 138. № 1/3. P. 289-308.

146. Chertov O.G., Komarov A.S., Bykhovets S.S., Kobak K.I. Simulated soil organic matter dynamics in forests of the Leningrad administrative area, northwestern Russia // Forest Ecology and Management, 2002. Vol. 169. No. 1/2. P. 29-44.

147. Chertov O.G., Komarov A.S., Loukianov A.M., Mikhailov A.V., Nadporozhskaya M.A., Zubkova E.V. The use of forest ecosystem model

148. EFIMOD for research and practical implementation at forest stand, local and regional levels // Ecological Modelling, 2006. Vol. 194. P. 227-232.

149. Chertov O., Bhatti J.S., Komarov A., Mikhailov A., Bykhovets S. Influence of climate change, fire and harvest on the carbon dynamics of black spruce in Central Canada // Forest Ecology and Management, 2009. Vol. 257. No. 3. P. 941-950.

150. Eckersten H., Jansson P.-E., Johnsson H. SOILN Model, Version 9.2, User's Manual / Swedish Agricultural University. Division of Hydrotechnics, Communication 98:6. Uppsala, 1998.

151. Forbes R.D. (Ed.) Forestry handbook. Ronald Press: New York, 1955.

152. Gower S.T., Vogel J., Stow T.K., Norman J.M., Steele S.J., Kucharik C.J.

153. Carbon distribution and above-ground net primary production of upland andlowland boreal forest in Saskatchewan and Manitoba // Journal of Geophysical Research, 1997. Vol. 104. P. 29029-29041.

154. Groisman P.Ya., Koknaeva V.V., Belokrylova T.A., Karl T.R. Overcoming biases of precipitation measurement: A history of the USSR experience // Bull, of American Meteorological Society, 1991. Vol. 72. No. 11. P. 1725-1733.

155. Gundersen P., Callesen I., de Vries W. Nitrate leaching in forest ecosystems is related to forest floor C/N ratios // Environmental Pollution, 1998. 102. No.l. SI. 1. P. 402-407.

156. Halliwell D.H., Apps, M.J. BOReal Ecosystem-Atmosphere Study (BOREAS) biometry and auxiliary sites: locations and descriptions / Natural Resources Canada, Canadian Forest Service, Northern Forestry Centre. Edmonton, 1997a. 120 p.

157. Halliwell D.H., Apps M.J. BOReal ecosystem-atmosphere study (BOREAS) biometry and auxiliary sites: overstory and understory data / Natural Resources Canada, Canadian Forest Service, Northern Forestry Centre. Edmonton, 1997b. 254 p.

158. Halliwell D.Ii., Apps M.J. BOReal ecosystem-atmosphere study (BOREAS) biometry and auxiliary sites: soils and detritus data / Natural Resources Canada, Canadian Forest Service, Northern Forestry Centre. Edmonton, 1997c. 235 p.

159. Halldin S. SOIL water and heat model. I. Synthesis of physical processes // AcTa Universitatis Uppsaliensis, 1980. No. 567. 28 p.

160. Hanson C.L., Cumming K.A., Woolhiser D.A., Richardson C.W. Microcomputer program for daily weather simulation in the contiguous United States / USDA Agriculture Research Service, ARS-114. 1994.34 p.

161. Heino R. Climate in Finland during the period of meteorological observations // FMI Contributions. No. 12. Finnish Meteorological Institute: Helsinki, 1994. 209 p.

162. Holdridge L.R. Simple method for determining potential évapotranspiration from temperature data// Science, 1959. Vol. 130. No. 3375. P. 572.

163. Jansson P.-E. SOIL water and heat model. II. Field studies and applications // AcTa Universitatis Uppsaliensis, 1980. No. 568. 26 p.

164. Jansson P.-E. Simulation model for soil water and heat conditions. Description of the SOIL model. Swedish University of Agricultural Science, Department of Soil Sciences, Report 165. Uppsala, 1991. 73 p.

165. Karjalainen T., Schuck-A. (Eds.) Causes and Consequences of Forest Growth Trends in Europe. Results of the RECOGNTION Project. Leiden: Brill, 2004.

166. Khanina L., Bobrovsky M., Komarov A., Mikhajlov A. Modelling dynamics of forest ground vegetation diversity under different forest management regimes // Forest Ecology and Management, 2007. Vol. 248. No. 1/2. P. 80-94.

167. Klein Tank A.M.G., Wijngaard J.B., Können G.P. et al. Daily dataset of 20* century surface air temperature and precipitation series for the European Climate Assessment // International Journal of Climatology, 2002. Vol. 22. No. 12. P. 1441-1453.

168. Komarov A., Chertov- O., Zudin S., Nadporozhskaya M., Mikhailov A., Bykhovets S., Zudina E., Zoubkova E. EFIMOD 2 a model of growth and cycling of elements in boreal forest ecosystems // Ecological Modelling, 2003. Vol. 170. No. 2/3. P. 373-392.

169. Kreutzer K. Changes in the role of nitrogen in Central European forests // Forest decline in the Atlantic and Pacific Region. Berlin: Springer, 1993. P. 82-95.

170. Kurz W.A., Apps M.J. A 70-year retrospective analysis of carbon fluxes in the Canadian forest sector // Ecological Applications, 1999. Vol. 9. No. 2. P. 526547.

171. Larionova A.A., Yevdokimov I.V., Bykhovets S.S. Temperature response of soil respiration is dependent on readily decomposable C // Biogeosciences, 2009. Vol. 4. No. 6. P. 1073-1081.

172. Li C., Frolking S., Frolking T.A. A model of nitrous oxide evolution from soil driven by rainfall events: 1. Model structure and sensitivity // Journal of Geophysical Research Atmospheres, 1992. Vol. 97. No. D9. P. 9759-9776.

173. Li C., Aber J., Stange F., Butterbach-Bahl K., Papen H. A process-oriented model of N20 and NO emissions from forest soils: 1. Model development // J. of Geophysical Research Atmospheres, 2000. Vol. 105. No. D4. P. 4369-4384.

174. Liski J., Palosuo T., Peltoniemi M., Sievanen R. Carbon and decomposition model Yasso for forest soils // Ecological Modelling, 2005. Vol. 189. No. 1/2. P. 168-182.

175. Luyssaert S., Schulze E.-D., Borner A., Knohl A., Hessenmoller D., Law B.E., Ciais P., Grace J. Old-growth forests as global carbon sinks // Nature, 2008. Vol. 455. No. 7210. P. 213-215.

176. Mekis E., Hogg W.D. Rehabilitation and analysis of Canadian daily precipitation time series // Atmosphere Ocean, 1999. No. 1. P. 53-85.

177. Mitchell T.D. Pattern scaling: an examination of the accuracy of the technique for describing future climates // Climatic Change, 2003. Vol. 60. No. 3. P. 217-242.

178. Monteith J.L. Principles of environmental physics. London: Edward Arnold, 1980. 241 p.

179. Myneni R.B., Keeling C.J., Tucker C. Asrar G., Nemani R.R. Increased plant growth in the northern latitudes from 1981 to 1991 // Nature, 1997. Vol. 386. No. 6626. P. 698-702.

180. New M., Hulme M., Jones P.D. Representing twentieth century space-time climate variability. Part 1: Development of a 1961—1990 mean monthly terrestrial climatology //Journal of Climate, 1999. Vol. 12. No. 5. P. 829-856.

181. New M., Hulme M., Jones P.D. Representing twentieth century space-time climate variability. Part 2: Development of 1901-96 monthly grids of terrestrial surface climate // Journal of Climate, 2000. Vol. 13. No. 13. P. 2217-2238.

182. Parton W.J. Abiotic section of ELM // Grassland simulation model / Ecological Studies, Vol. 26. New York: Springer, 1978. P. 31-53.

183. Parton W.J. Predicting soil temperature in a shortgrass steppe // Soil Science, 1984. Vol. 138. No. 2. P. 93-101.

184. Parton W.J., Shimel D.S., Cole C.V., Ojima D.S. Analysis of factors, controlling soil organic matter levels in Great Planes grasslands // Soil Science Society of America Journal, 1987. Vol. 51. No. 5. P. 1173-1179.

185. Parton W.J., Hartman M., Ojima D., Shimel D. DAYCENT and its land surface submodel: description and testing // Global and Planetary Change, 1998. Vol. 19, No. 1/4. P. 35-48

186. Pastor J., Post W.M. Development of a Linked Forest Productivity Soil Process Model / Oak Ridge National Laboratory. ORNL/TM-9519. 1985. 168 p.

187. Peterson T.C., Vose R.S. An overview of the Global Historical Climatology Network temperature database // Bull, of American Meteorological Society, 1997. Vol. 78. No. 12. P. 2837-2849.

188. Pitman A.J. Review: The evolution of, and revolution in land surface schemes designed for climate models // International Journal of Climatology, 2003. Vol. 23.No. 5. P. 479-510.

189. Price D.T., McKenney D.W., Papadopol P., Logan T., Hutchinson M.F. High resolution future scenario climate data for North America // Proceedings of the

190. American Meteorological Society. 26th Conference on Agricultural and Forest Meteorology. Vancouver, 2004. 13 p. (CD-ROM).

191. Racsko P., Szeidl L., Semenov M. A serial approach to local stochastic weather models // Ecological Modelling, 1991. Vol. 57. No. 1/2. P. 27-41.

192. Rawls W.J., Pachepsky Y.A., Ritchie J.C., Sobecki T.M., Bloodworth H. Effect of soil organic carbon on soil water retention // Geoderma, 2003. Vol. 116. No. 1/2. P. 61-76.

193. Richardson C.W. Stochastic simulation of daily precipitation, temperature, and solar radiation//Water Resources Research, 1981. Vol. 17. No. l.P. 182-190.

194. Saxton K.E., Rawls W.J. Soil water characteristic estimates by texture and organic matter for hydrologic solutions // Soil Science Society of America Journal, 2006. Vol. 70. No. 5. P. 1569-1578.

195. Sevruk B. (Ed.) Correction of precipitation measurements. Zurich, 1986.

196. Shugart H.H. A theory of forest dynamics. New York: Springer, 1984. 278 p.

197. Shvidenko A.Z., Schepaschenko D.G., Nilsson S., Buluy Yu.I. Tables and models of growth and productivity of forest of major forest forming species of Northern Eurasia (standard and reference materials). Moscow, 2006.

198. Spiecker H. Growth trends in European forests Do we have sufficient knowledge? // Karjalainen T., Spiecker H., Laroussinie O. (eds.): Causes and consequences of accelerating tree growth in Europe. EFI Proceedings, 1999. Vol.27. P. 157-169.

199. Stolbovoi V. Carbon in Russian soils // Climatic Change, 2002. Vol. 55. No. 1/2 P. 131-156.

200. Thornthwaite C.W. An approach toward a rational classification of climate // Geographical Review, 1948. Vol. 38. No. 1. P. 84-94.

201. Tiktak A., van Grinsven H.G.M. Review of sixteen forest-soil-atmosphere models // Ecological Modelling, 1995. Vol. 83. No. 1/2. P. 35-53.

202. Toy T.J., Kuhaida A.J. jr., Munson B.E. The prediction of mean monthly soil temperature from mean monthly air temperature // Soil Science, 1977. Vol. 126. No. 3. P. 181-189.

203. Wösten J.H.M., Pachepsky Y.A., Rawls W.J. Pedotransfer functions: bridging the gap between available basic soil data and missing soil hydraulic characteristics // Journal of Plydrology, 2001. Vol. 251. No. 3/4. P. 123-150.

204. Yin X., Arp P.A. Predicting forest soil temperatures from monthly air temperature and precipitation records // Canadian Journal of Forest Research, 1993. Vol. 23. No. 12. P. 2521-2536.