Бесплатный автореферат и диссертация по сельскому хозяйству на тему

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

ВАК РФ 06.01.03, Агропочвоведение и агрофизика

Автореферат диссертации по теме "Многолетняя изменчивость климата и прогнозирование системы погода-урожай"

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

СМИРНОВА Инна Иннокентьевна

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

Специальность: 06.01.03 - агропочвоведение, агрофизика

Автореферат

на соискание ученой степени ии^4ВВО 14 кандидата сельскохозяйственных наук

Улан-Удэ, 2009

003488014

Работа выполнена на кафедре мелиорации и охраны земель Института землеустройства, кадастров и мелиорации ФГОУ ВПО «Бурятская государственная сельскохозяйственная академия им. В.Р.Филиппова»

Научный руководитель: доктор биологических наук, профессор,

Куликов Анатолий Иннокентьевич

Официальные оппоненты: доктор сельскохозяйственных наук, профессор

Чимитдоржиева Галина Доржиевна кандидат сельскохозяйственных наук, доцент Сордонова Маргарита Николаевна

Ведущая организация: Институт почвоведения и агрохимии СО РАН

Защита диссертации состоится " 25 " декабря 2009 г. в 9 час. на заседании диссертационного совета К 220.006.02 при Бурятской государственной сельскохозяйственной академии им. В.Р. Филиппова по адресу: 670034, г. Улан-Удэ, ул. Пушкина, 8. www.bgsha.ru факс: (3012) 44-21-33

С диссертацией можно ознакомиться в библиотеке Бурятской государственной сельскохозяйственной академии им. В.Р. Филиппова и на сайте www.bgsha.ru

Автореферат разослан 25 ноября 2009 г.

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

Л-/

р Корсунова Т.М.

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

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

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

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

Задачи исследования

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

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

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

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

Защищаемые положения:

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

менения и позволяют установить климатические индикаторы опустынивания;

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

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

Апробация работы. Основные результаты работы доложены на научных, научно-практических конференциях и симпозиумах: международных (Улан-Удэ, 2003; Улан-Батор, 2008), региональных (Улан-Удэ, 2004; Кемерово, 2005) на заседаниях кафедры мелиорации и охраны земель Бурятской ГСХА им. В.Р. Филиппова в 2002-2006 гг. и нашли отражение в 8 опубликованных работах, в том числе 2 в рецензируемых изданиях.

Структура и объем диссертации. Диссертация изложена на 120 страницах компьютерного текста и состоит из введения, 5 глав и выводов. Содержит 9 таблиц, 29 рисунков и 2 приложений. Список литературы включает 180 наименований, в том числе 7 иностранных авторов.

Глава 1. Природные условия степных ландшафтов Западного Забайкалья

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

Глава 2. Методика исследования временных рядов

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

Методика определения динамики обеспеченности метеоэлементами в зависимости от циклов солнечной активности основана на обнаружении и статистическом подтверждении связей между элементами климата и солнечной активностью.

Глава 3. Цикличность и трендовые изменения климата

Цикличность изменения климата в Забайкалье изучалась рядом исследователей (Жуков, 1965; Киссис, 1969; Куликов и др., 1997; Андреев, 2004; Шестср-нев, Шейнкман, 2008; Волошин и др., 2007).

Цикличность и тренды температуры воздуха. По данным метеостанции Улан-Удэ за 103 года (1900-2003 гг.), потепление климата оценивается ростом температуры воздуха на 2,5°С (рис. 1). При этом в Новоселенгинске среднегодовая температура воздуха повысилась на 1,8, в Кяхте — на 1,6°С. Получены уравнения линейного тренда среднегодовой температуры воздуха (У), которое для Улан-Удэ имеет вид: У = 0,0247 Х-49,383. Наиболее резкие изменения среднегодовой температуры воздуха происходят с 70-х гг. XX в., что характерно для северного полушария и согласуется с данными Всемирной метеорологической организации (Изменение климата..., 2003), а также данными для ряда пунктов России (Перльштейн и др., 2006; Мельников и др., 2007; Куликов, Куликов, Смирнова, 2008 и др.).

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

Рис. 1. Многолетняя динамика и тренды среднегодовой температуры воздуха в Западном Забайкалье. 1 - Улан-Удэ, 2 - Новоселенгинск, 3 - Кяхта

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

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

1900 1920 1940 1960 1980 2000 Годы

Рис. 2. Синхронно-асинхронная изменчивость многолетней динамики среднегодовой температуры воздуха (1) и солнечной активности (2) (сглаженные ряды с шагом 5 лет)

Посезонный анализ позволяет выяснить структуру явления потепления климата, установить количественный вклад каждого сезона и обнаружить критический период, ответственный за интегральный эффект. Исследования временного ряда с 1970 по 2003 гг. показывают, что рост средних температур зимнего периода составил 1,1°С, примерно в такой же степени происходит потепление лета. При этом для весеннего периода отмечается похолодание примерно на 0,5°С. Осень за последние 33 года потеплела незначительно, приблизительно на 0,3°С. Это указывает на то, что повышение среднегодовых температур происходит в большей степени за счет потепления зимнего и летнего сезонов.

В многолетней динамике суммарных значений температуры воздуха временной темп роста суммы температур воздуха выше 0 и выше 10°С больше, чем темп нарастания суммы температуры воздуха выше 5°С (рис. 3). Так, из угловых коэффициентов уравнений следует, что £р-0°С и 10°С имеют темп нарастания соответственно 10,3 и 11,4°С/год, тогда как 2л> 5°С во времени увеличивается с темпом 9,4°С/год.

Рис. 3. Многолетняя изменчивость сумм температур выше 0°С (/), 5°С (2), 10°С (5) с линиями и уравнениями тренда

Сумма активных температур является важным показателем теплообеспе-ченности, которая в свою очередь определяет потенциальные природные ресурсы сельского хозяйства, обусловливающие набор сельскохозяйственных культур по их требованию к теплу, а также формирование их продуктивности. В Улан-Удэ 2>10°С = 1897°С, при экстремумах 1548 (1995 г.) и 2371°С (2002 г.).

Тренд-анализ показывает (Убугунов, Куликов, Смирнова, Мангатаев, 2008), что показатель ^О10°С на территории Западного Забайкалья возрастает в разных пунктах с длинными (более 100 лет) рядами наблюдений на 370-450°С. Если воспользоваться эмпирической формулой Г.Т. Селянинова:

Ис= X? > 10 " С =1 ММ;

10

где Ис - испаряемость воды, то вытекает, что потенциал испарения возрос на 37-45 мм. Из этого следует, что по дефициту испарения климатическую норму орошения необходимо принимать на 370-450 м3/га больше существующих. Этот вывод имеет существенное значение и вкладывает новое агромелиоративное содержание в глобальные изменения. Тем более это важно, что исследования многолетнего ряда (с 1900-2003 гг.) не показывают какого-либо статистически достоверного прироста количества атмосферных осадков. Величина роста климатической нормы орошения весьма существенна, если учесть, что в настоящее время в Бурятии орошается 151257 га земель, а потенциальный мелиоративный фонд составляет 2306 тыс. га.

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

Глобальные изменения отражаются на продолжительности сезонов (рис. 4). Установлено увеличение весеннего, летнего и осеннего сезонов, соответственно зимний период сокращается. Если в начале 1970-х годов продолжительность периодов с положительными и отрицательными температурами была примерно равна (180-185 дней), то в начале 2000-х годов продолжительность периода с температурами выше 0°С составляет более 200 дней.

Годы

Рис. 4. Динамика продолжительности периода между датами перехода через 0°С и линия тренда.

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

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

Средний уровень количества выпадающих атмосферных осадков в Западном Забайкалье незначительный, колеблется от 220 (Новоселенгинск) до 346 мм/год (Кяхта), причем на протяжении почти 80-ти лет он остается практически неизменным по средней величине. Однако весьма значительна неустойчивость количества осадков по годам, что является одним из важных условий рискованности земледелия. Так, абсолютный максимум осадков, приходящийся на 1959 г., составляет 413,3 мм, а абсолютный минимум, наблюдавшийся в 1989 г., равняется 106,5 мм.

Сглаженный ряд по полиномиальной модели 6-й степени позволяет выделить 3 цикла большого временного диапазона. Длительность больших циклов составляет примерно 35 лет, что укладывается в брикнеровскую цикличность, период которой равняется 30-50 лет, а в среднем 35 лет (Климатология, 1989). На брикнеровские циклы наложены более мелкие циклы. При фильтрации временного ряда с пятилетним шагом выявляется до 9 циклов с меньшим временным диапазоном.

С 1935 г. выявляется достаточно длительный (8 лет) цикл повышенного атмосферного увлажнения, сопровождающийся пессимумом до 1947 г. Довольно низко увлажненность падает в обертоне 1951-1957 гг. Последующий цикл 19581972 гг. характеризуется достаточно высоким увлажнением, но с мелкими раз-нознаковыми изменениями на всем его протяжении. Такие же частотные, но уже с гораздо большей амплитудой колебания наблюдаются в последние годы, начиная с 1991 г. Этот период включает в себя одно-, двухгодичные обертоны, которые невозможно считать полноценно значимыми циклами ввиду их непродолжительности. В силу этого период с 1991 по 2003 гг. принимается как единый цикл.

Из выделенных 9 циклов, 5 - повышенного увлажнения и соответственно 4 - пониженного увлажнения. По длительности циклы различны: от 3 до 15 лет (табл. 1) и совпадают с известными 11 -летними циклами солнечной активности.

Стандартное отклонение, как мера того, насколько широк разброс эмпирических данных относительно их среднего, показывает, что во времени атмосферное увлажнение становится все менее стабильным. Так, если в первом периоде (1935-1942) вариабельность количества годовых осадков измерялась стандартным отклонением, равным ±32,2 мм (13%), то в настоящее время отклонс-

ния от климатической нормы возросли и стандартное отклонение стало составлять уже ± 65-83 мм, т.е. отклонения от средней достигли 25-32%. '

Таблица 1. Характеристика циклов повышенного (+) и пониженного (-) атмосферного увлажнения

№ цикла Климациклы по годам ± Т Хтах у Z

I 1935-1942 + 8 283,0 336,6 230,0 32,2

II 1943-1947 - 5 237,8 289,0 160,0 48,4

III 1948-1950 + 3 265,8 305,8 232,3 -

IV 1951-1957 - 7 212,9 271,1 158,2 40,3

V 1958-1972 + 15 268,9 413,3 152,7 60,3

VI 1973-1981 - 9 206,7 365,9 138,9 64,7

VII 1982-1985 + 4 312,1 327,9 298,2 -

VIII 1986-1990 - 5 235,1 312,2 106,5 82,8

IX 1991-2003 + 13 277,1 407,1 186,4 65,0

Среднее 1935-2003 69 256,0 413,3 106,5 -

Т- длительность цикла, год; ± циклы повышенного или пониженного атмосферного увлажнения; ~х- среднее значение суммы осадков, мм; Хтш, Хт.п -экстремумы количества осадков, мм; - стандартное отклонение количества осадков, мм

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

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

Видно, что чем позже происходит переход среднесуточных температур через определенный предел, тем меньше сумма температур за данный период. Аналогичные данные получены и по другим пунктам Западного Забайкалья. Тем самым, полученные данные по Западному Забайкалью согласуются с концепцией Ф.Ф. Давитая (1964) и с результатами расчетов А.И.Куликова и др. (1997).

Прогноз теплообеспеченности предстоящего сельскохозяйственного года в виде суммы активных температур воздуха (2/>Ю°С) проводят, зная дату (Д) перехода средней температуры в мае через 10°С (можно узнать на ближайшей ме-

теостанции), по регрессионной модели У>> 10°С = -64,72Д+ 904,27. Тем самым уже в мае можно располагать информацией об ожидаемых тепловых ресурсах вегетационного сезона.

о 1200 Т

о

-1—1-1-Г—Г-1-Г—Г-Г—1-1—Т-Г~1-1 1 1 1,1 I-1-1-1-г—1-1-1-1-1-1-1-1-1

£ 0 1 2 3 4 5 5 7В 9 10 11 12 13 1 4 1516 17 18 19 20 21 22 23 24 25 2Б272В 29 30 31

Дни месяца

Рис. 5. Связь между датой весеннего перехода температуры через 10 °С (май) и суммой температур выше 10°С

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

450 т-

1935 1940 1945 1950 1955 1960 1965 1970 1975 1930 1985 1990 1995 2000

Голы

Рис. 6. Сопряженная изменчивость годовой суммы осадков (1) и солнечной активности (2)

При этом на каждый пик возрастания активности солнца приходится снижение количества осадков ниже среднемноголетнего уровня, вероятность совпадения составляет 100%, причем в 4-х случаях из 6-ти (60,6%), количество осадков снижается до значений менее 200 мм (острозасушливая ситуация), в 2-

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

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

Ё 400

§ 350

к к 300

250

О И.-! 200

12 Я 150

а а) 100

■во 50

у 0

4!

О I III | III

ГТТ г Г 1

^Г МИШИ"

Рис. 7. Динамика сумм годовых осадков (Г) и шкала динамики тенденций (1Г)

При допустимой ошибке 5% построен график (рис. 7) динамики атмосферного увлажнения по методике И.П. Дружинина (1977), где годы с тенденцией повышения увлажнения отмечены как столбики вверх от средней линии, а годы с тенденцией понижения увлажнения, соответственно, как столбики вниз. Общее число тенденций-39, из них 1-годичных-23 (59%), 2-годичных- 14 (36%), 3-годичных - 1 (2,5%), 4-годичных -1 (2,5%). Залитыми столбиками отмечены годы совпадения солнечных реперов с годами начала смены тенденции.

Большая вероятность таких совпадений (16 из 19 реперов - 84 %) свидетельствует об их неслучайности, что позволяет использовать данные о солнечных циклах для прогноза количества атмосферных осадков, особенно экстремальных по величине.

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

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

№ подхода Прогнозные подходы Вероятность, %

1 Вероятность уменьшения количества осадков в годы максимума солнечной активности 100

2 Вероятность а) уменьшения количества осадков на второй год после максимума солнечной активности 60

б) снижения количества осадков до 200 мм и менее 60

в) снижения количества осадков до 150 мат и менее 30

3 Вероятность ¿7) повышения количества осадков после окончания пика солнечной активности 100

б) наступления двухгодичной тенденции повышения увлажнения 83

4 Вероятность смены текущей тенденции в предстоящем году солнечного репера 84

5 Вероятность смены текущей тенденции на следующий год 59

Прогноз № 1 имеет абсолютную надежность (100%) и вытекает из факта, что в год пика солнечной активности наблюдается спад количества атмосферных осадков. Пик солнечной активности продолжается 1-3 года. Подход № 2а 60%-ной вероятности основан на том, что наибольший спад количества атмосферных осадков приходится именно на 2-й год после пика солнечной активности. Подходы № 26 и 2в имеют уточняюще-дополняющую функцию к подходу 2а. Прогноз по модели № 3 основан на том, что после окончания периода с высокой солнечной активностью наблюдается повышение количества осадков, особенно на второй год спада солнечной активности. Эта модель свидетельствует о том, что увеличение количества осадков в период спада солнечной активности происходит в каждом из годослучаев, т.е. вероятность равна 100%. С вероятностью 83% (подход № 36) можно утверждать, что тенденция увеличения количества осадков продлится два года. Кроме этого на характер увлажнения влияют годы солнечных реперов (повышенные значения чисел Вольфа между смежными годами) - вероятность 84% (прогнозный подход № 4).

Подход № 5 основан на том, что 1-летние годослучаи достаточно представлены в статистическом ряду. Вероятность того, что следующий год будет противоположен по знаку тенденции настоящего года равен 59%.

Тем самым, модели, основанные на качественном подходе нередко более надежны по сравненшо с аналитическими.

Количественные показатели климатической составляющей опустынивания. Наиболее благоприятные условия для опустынивания складываются в днищах котловин Западного Забайкалья при радиационном балансе, превышающем 1600 МДж/м2 и значениях RJLr > 2,0 (R, - радиационный баланс, £ - теплота парообразования, г - количество атмосферных осадков). В осях Н-j (Я -абсолютная высота, j - широта) с координатами < 900 м - < 52° формируется региональное «ядро» опустынивания, где R/Lr > 2,5, а по внешнему кольцу -вторичный ареал. Сдвиг вторичного ареала в сторону высоких широт связан с наличием в регионе высокоширотных степных изолятов в таких котловинах, как Баргузинская и Муйско-Куандинская.

Климатическую параметризацию опустынивания также проведена, сравни-

вая наши данные за 1970-2003 гг. и данные, помещенных в «Справочниках по климату СССР» (1966,1968), где обработан временные ряды с 1900 или 1935 гг. по 1960 или 1965 гг. (для разных метеопараметров временные отрезки различаются). Исходили из допущения, что временной период, охваченный «Справочниками. ..», не характеризовался столь сильно выраженными процессами опустынивания, как последние три десятилетия.

Сравнение параметров преддезертификационного (1900-1965 гг.) и современного (1970-2003 гг.) периодов однозначно указывает на климатическую предрасположенность усиления процессов аридизации и дезертификации (табл. 3). Климатическим условием аридизации и опустынивания является рост суммы температур выше 10°С, продолжительности периода с этими температурами.

Таблица 3. Разностные показатели метеорологических параметров, полученные при сравнении авторских данных и данных из «Справочников по климату СССР» (1966, 1968), метеостанция Улан-Удэ

Метеопараметр Разность

Средняя годовая температура воздуха, °С +1.2

Средняя температура воздуха января, °С +6

Сумма температур воздуха > 10°С +400

Продолжительность периода с 1в > 10°С, дни +13

Сумма атмосферных осадков, мм/год +6

Относительная влажность воздуха, % -8

ГТК -0,15

Количество дней с ветром более 5 м/с +12

Количество дней с суховеями +9

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

Сравнение данных, полученных по временным рядам с 1900 по 2003 гг., с данными из «Справочников по климату СССР» (1966,1968) показывает их значительное, а по некоторым параметрам очень значительное расхождение. Всем организациям, пользующимся метеорологической информацией, в первую очередь, сельскохозяйственным органам, а также энергетикам, строителям и др. необходимо учитывать это обстоятельство.

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

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

(2003) и продиктован решениями Всемирной конференции ООН по природным опасностям (Yokohama..., 1995).

По значениям удельного физического риска видно, что из всех типов опустынивания наибольшую опасность для сельскохозяйственных земель Бурятии представляет дефляция. Площадная дефляция земель происходите интенсивностью 2,2 год"1. Водная эрозия имеет интенсивность 1,5 год-1. Полный физический риск потери сельскохозяйственных земель в результате дефляции для всей Бурятии составляет254,5 га/год. Дефляция приносит материальный ущерб, равный 89 га/год. Экономический риск в зависимости от типа опустынивания колеблется в пределах 430-700 руб./год. Экономическая уязвимость земель наибольшая при дефляционном разрушении.

Уязвимость населения от дефляции оценивается величиной 15 чел./год, а неблагоприятное действие водной эрозии и водной эрозии совместно с дефляций ежегодно испытывают соответственно 12 и 11 человек соответственно, а уязвимость от засоления составляет 7,2 чел/год.

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

Глава 4. Статистические связи в системе погода-урожай Статистический анализ многолетней динамики урожайности пшеницы. Урожай является функцией многих факторов, среди которых погодные условия занимают одно из основных мест. Среднее значение урожайности пшеницы за период 1968-2004 гг., по данным ОПХ «Иволгинское», составляет 17,5 ц/ га. Линейный тренд указывает на заметное (около 7 ц/га) снижение урожайности. Тренд описывается регрессионным уравнением вида У = - 0,13Х+ 269,44.

В многолетней изменчивости количества урожая выделяются определенные циклы. В 1968-1973 гг. урожайность приближается к среднему значению. В 1974-1976 гг. урожайность стабильно снижается, что связано, по-видимому, с периодом пониженной влагообеспеченности. Наибольшие урожаи получены в начале 80-х годов XX в. В дальнейшем наблюдается устойчивое снижение урожайности, которое описывается регрессионной моделью У--0,9711'-»- 1953,6. Тренд за период 1981 -2004 гг. характеризуется высоким значением углового регрессионного коэффициента, который показывает, что урожайность пшеницы падает с темпом 1 ц/га в год.

Такая отрицательная динамика урожайности характерна для всего Западного Забайкалья, Бурятии и всей России. Она вызвана теми социально-экономическими потрясениями, которые постигли страну в годы перестройки и постперестроечное время. Тем самым в динамике урожайности с учетом социально-экономического фактора выделяется два этапа: предстагнационный и стагнационный. Логично также выделить пока еще неизученный третий этап—современный, связанный с новыми рыночно-экономическими условиями и земельными реформами. Удостовериться в этом можно на объективной основе. Проведенные иссле-

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

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

35 30

.£ 25

и £ 20

5 15

10

/ Л \ Г г А >

/1 \ п Л. Ч

"к \ ¡\ к / г * X

V / ( С V у V Г ■

V. > г. - V У -

140

га

100 ¡3

о «з

-ВО

20

1968

1972 1976 1980 1984 1988 1992 1996 2000 Годы

Рис. 8. Динамика урожайности и солнечной активности в числах Вольфа (сглаженные ряды с шагом в 5 лет). 1 - фактическая урожайность, 2 - расчетная урожайность, 3 - числа Вольфа

Зависимость урожайности от чисел Вольфа в предстагнационный период, когда она формировалась в функции от погодных условий и биологических возможностей культуры описывается моделью У = 26,784 - 0,125а". Используя эту зависимость, проведено математическое моделирование урожайности при условии, что если хозяйственно-экономический комплекс сохранился бы на уровне предстагнационного этапа. По расчетной кривой 2 (см. рис. 8), видно, что урожайность пшеницы была бы выше, особенно в последние годы.

Корреляционно-регрессионные функции урожая от погодных условий.

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

Таблица 4. Корреляционно-регрессионные связи урожайности (у) и количества осадков (х)

Метеорологический показатель корреляция г уравнение регрессии

сумма осадков за год г = 0,6225 у= -1,3452+ 0,0765х

сумма осадков за лето г = 0,5904 у = 2,8443 +0,0906х

сумма осадков за июнь г = 0,5742 у = 8,9063 + 0,2447х

сумма осадков за период от уборки урожая предыдущего года до уборки урожая текущего года г = 0,6527 у = - 2,9627 + 0,0840х

сумма осадков за 3 декаду июня г = 0,4237 у = 13,846+ 0,2187х

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

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

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

Расчет ГТК проводился по формуле:

ГТК= 5>10°С/0,1Д,

где 2>10°С - сумма активных температур выше 10°С, /?- количество осадков за тот же период.

Установлена достаточно плотная связь урожайности с ГТК для периода от посевов до всходов (табл. 5).

Таблица 5. Корреляционно-регрессионные связи урожайности (У) с ГТК{Х)

Метеорологический показатель г уравнение регрессии

ГТК периода от посева до уборки урожая г = 0,6538 у = 0,3838+ 16,492х

ГТК за июнь г = 0.6394 \= 8,4889+ 12.141х

ГТК периода между датами перехода через 10°С г = 0,5583 у= -0,5839 + 19.786х

ГТК 3-й декады июня г = 0,4401 у = 16,888 + 3,3641.\

Связь между урожайностью и ГГК в третьей декаде июня, т.е. в наиболее важный период для развития яровой пшеницы, ниже, чем связь урожайности и ГТК за весь июнь. Это может свидетельствовать о том, что осадки конца июня не могут восполнить потери урожая при неблагоприятных условиях 1-3 июня, когда растения находятся в тех фенологических фазах, которые наиболее требовательны к воде.

Разработана трехфакторная модель урожайности: у = -0,621 + О,24\х1 + 0,129x2, где х] - сумма осадков в период от уборки предыдущего года до посева в текущем году, х2 - ГТК за июнь. Коэффициент множественной корреляции Я = 0,535, который указывает на то, что можйо надежно прогнозировать до 28% го-дослучаев (Я2 = 0,285). Проверочные расчеты показывают, что уравнение не отличаясь большой точностью, может очень надежно предсказывать, какова будет урожайность в терминах «низкая», «средняя», «выше средней».

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

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

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

ПЭи = 12,5 (2-ЮП Дс Е'-с^ло'С(КУ-Р) {

Л'ЛГ + 100

где ПЭи - почвенно-экологический индекс, V- плотность (объемная масса) почвы (в среднем для метрового слоя); 2 - максимально возможная плотность

почв при их предельном уплотнении, г/см3; 77- «полезный» объем почвы (в метровом слое); Дс - дополнительно учитываемые свойства почв; КУ - коэффициент увлажнения; Р - поправка к этому коэффициенту; КК-коэффициент конти-нентальности; А - итоговый агрохимический показатель.

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

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

ПЭи=12,5(2- У)П Y,t^™0C(Ky-P), як+ 10

Дк R

где КУ~ ^ t> j 0 о с + 500 ; (Дк - дополнительный коэффициент; R - средне-

360(t -t . )

годовая сумма осадков; КК=-шах mm

ср + ш

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

Самый низкий индекс принадлежит Баргузину и Улан-Удэ. В районах, характеризуемых этими пунктами, свойственны низко плодородные песчаные почвы. Здесь климат наиболее засушливый и континентальный. Тот факт, что Кяхта имеет довольно высокий почвенно-экологический индекс объясняется достаточно высоким атмосферным увлажнением (345 мм/год), а значит, здесь относительно высок коэффициент увлажнения, а пониженный коэффициент континентальное™ свидетельствует об относительной мягкости климата. Противоречивое (см. формулу ПЭи) взаимодействие этих двух параметров приводит к повышению величины ПЭи.

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

Как видно из представленных данных, количественные изменения во времени почвенно-экологического индекса значительны, от 33 до 57. Из этого факта вытекает, что при бонитировочных работах с использованием почвенно-экологического индекса необходимо учитывать его изменчивость во времени. Пользование только средними климатологическими значениями при расчетах

60

г 36

а

Г)

о

о

И 30

20 И

1970-1974 1975-1979 1930-1984 1985-1989 1990-1994 1995-1999 2000-2004 Пятилетия

Рис. 9. Временное изменение почвенно-экологического индекса (/), коэффициента увлажнение (2), коэффициента континентальности (3). Метеостанция Улан-Удэ

коэффициентов увлажнение и континентальное™ может привести к неправильным выводам.

Глава 5. Микроклиматическая неоднородность земель в расчлененном рельефе. Наряду с вскрытыми выше закономерностями временной динамики климата, его элементы имеют ярко выраженную динамику по еще одной координате - пространственной. Среди известных разномасштабных в пространстве изменений климата особое значение имеет микроклиматическая изменчивость. Микроклиматические ресурсы территории до сих пор используются неполно, и в этом отношении микроклимат служит резервом для роста продуктивности сельского хозяйства. Так, только стандартная агрометеорологическая служба путем информационного обеспечения способствует снижению потерь урожая от неблагоприятных явлений на 20-30%.

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

Ландшафтная структура Тапхарской депрессии. Ландшафтные исследования в Западном Забайкалье проводились рядом исследователей (Типы местности, 1959; Рещиков, 1961; Фадеева, 1963; Предбайкалье и Забайкалье, 1965; Ландшафты юга Восточной Сибири, 1977; Михеев, 1988, 1992; Дамбиев, 1991; Иметхенов, 1997; Мельник, 1999; Сымпилова и др., 2003 и др.). Трудами уче-ных-ландшафтоведов приведены в известность структурно-функциональные особенности геосистем Западного Забайкалья.

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

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

Составлена ландшафтная карта Тапхарской депрессии М. 1 : 25000. Легенда к ней содержит 12 ландшафтных урочищ. Наибольшее распространение имеют склоновые каменистые сухостепные мелкодерновинные на горных маломощных каштановых почвах и пологонаклонные равнинные сухостепные сухостепные твердоватоосоковые полынно-разнотравные на среднемощных каштановых почвах мучнистокарбонатных почвах урочшца. В подножии Ганзуринского хребта также довольно широко распространены наклонные предгорные сухостепные полынно-злаково-разнотравные на маломощных каштановых почвах урочища. Наличие мелкодолинных овражно-балочных и ложбинных древнего и современного неустойчивого стока мелкодерновинных разнотравно-злаковых урочищ свидетельствует о развитии неблагоприятных эрозионных явлений.

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

Микроклиматические условия Тапхарской депрессии. Микроклиматические исследования в Западном Забайкалье, равно как в Бурятии и Забайкалье, проведены только фрагментарно. По региону известны работы А.П.Волошиной (1962), Н.А.Мячковой (1962). Н.Б.Бадмаев и др. (1996) изучили тепловлагоо-беспеченность склоновых земель одной из сопок Еравнинской котловины. Исследованиями А.И.Куликова и др. (1997) установлены основные закономерности радиационного режима, перераспределения физических, физико-химических свойств почв и формирования ресурсов влаги и тепла на склонах разной экспозиции, при этом установлены топоэкологические коэффициенты контрастности склоновых ландшафтов.

В дневные часы температура почвы на глубине 0,2 м стабильно выше на 2-3°С на склоне южной экспозиции по сравнению со склонами северной экспозиции. В ночные часы нередко происходит пространственная инверсия и северные склоны оказываются более теплыми на 0,5-0,7°С. На фоне общей очень низкой влажности почвы (2-9% от массы) инсолируемые склоны больше иссушены в среднем на 1-1,5%. Наименьшая влажность характерна для маломощной сильно хрящеватой почвы вершины увала.

Составлена карта микроклимата Тапхарской депрессии, где выделено 9 микроклиматических участков, различающиеся по параметрам (табл. 6). Самые теплые сильно аридизированные склоны в виду инсолируемых позиций обеспечены энергией фотосинтетически активной радиации (ФАР) на 7 млн. МДж/га больше, чем фоновый участок. При этом они на территории депрессии наиболее континентальны, т.к. в ночные часы сильно выхолаживаются в виду выпуклости

теплоотдающей поверхности, а днем сильно нагреваются, о чем свидетельствуют знаки разностей экстремальных температур воздуха (?тт и ? ). Почвы инсо-лируемыхучастков лучше прогреваются, на глубине 0,2 м температура (70,)выше, чем на фоновом на 4,5°С. Влагозапас в слое 0-0,5 м меньше, чем в почвах фонового участка на 40 мм, что объясняется тепловым режимом и повышенными возможностями испарения влаги. В Иволгинской котловине потенциально может испариться более 600 мм воды, а на склонах южной экспозиции испаряемость (Ео) еще выше на 200 мм. Вследствие указанных причин период с температурой выше 5°С (г (г>5°С) на этих склонах увеличивается на 30 дней.

Таблица 6. Микроклиматическая характеристика Тапхарской депрессии (май-сентябрь)

Микроклиматические участки ДФДР •105, МДяс/га йСш ,„ °с ЛС„ л, °С "С А1Г„. 0.5. мм ДЕо, мм Лт (/>5» С), дни

1. Самые теплые очень засушливые 70 -2 4,5 3,5 -40 200 30

2. Самые теплые сухие 40 -1 3 2.5 -30 100 20

3. Теплые засушливые* 140 5,6 21,2 14,7 60 620 148

4. Теплые умеренно засушливые 0... -5 -I 0 -1,5 25 -50 -10

5. Умеренно теплые умеренно увлажненные -25 2,5 -3 -2 25 -50 -5

6. Прохладные умеренно увлажненные -30 0 -5 -2 70 -100 -15

7. Прохладные влажные -40 0,5 -10 -7 >10 0 -150 -15

8. Заморозкоопасные замкнутые -50 -5 -20 -10 100 -250 -20

9. «Фитоклимат стволов» <8-15 раз -5 -10 -7 10 -100 -20

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

ВЫВОДЫ

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

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

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

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

5.Установлены почвенно-экологические индексы и показана их пространственно-временная изменчивость.

6. Выявлены микроклиматические условия Тапхарской депрессии и проведено ландшафтно-микроклиматическое картографирование.

Предложения производству

1. Для заблаговременного определения обеспеченности земледелия суммами активных температур (2л> 10°С) рекомендуется простая и легко реализуемая на обычном калькуляторе регрессионная модель: £г>10°С = -64,72Д - 904,27, где Д-дата устойчивого перехода температуры воздуха через 10~С в мае.

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

3. Для того, чтобы получить достоверное представление о видах на предстоящий урожай в категориях «низкая», «средняя», «выше средней» предлагается трехфакторная модель урожайности: у = -0,621 + 0,241х/ + 0,129x2, где х! -сумма осадков в период от уборки предыдущего года до посева в текущем году, х2 — ГТК за июнь.

Список опубликованных работ.

1. Амшеев Р.М., Куликов А.И., Смирнова И.И. О тенденциях изменения климата г. Улан-Удэ за последние 3 0 лет // Материалы научно-практической конференции, посвященной 50-летию агрономического факультета БГСХА, 2002.-е. 13-16.

2. Куликов А.И., Куликов М.А., Смирнова И.И. О глубине протаива-иия почв при изменениях климата // Вестник БГСХА - Выпуск 14. - Улан-Удэ: Изд-во БГСХА, 2009.-С.121-126.

3. Куликов А.И., Куликов М.А., Смирнова И.И. О проблеме водного дефицита// Вестник БГСХА - Выпуск 13. - Улан-Удэ: Изд-во БГСХА, 2008.-С.112-115.

4. Смирнова И.И. Об изменениях климата степной зоны Республики Бурятия за последние 50 лет // Научное обеспечение устойчивого развития АПК в Сибири / Материалы конференции молодых ученых Сибирского федерального округа (7-11 июля 2004 г.). - Улан-Удэ, 2004. - С. 68-69.

5. Смирнова И.И. Прогноз обеспеченности осадками степных ландшафтов Республики Бурятия // Инновационное развитие аграрного производства в Сибири / Сборник материалов Третьей конференции молодых ученых вузов «Аг-рообразования» СФО. - Кемерово, 2005. - С. 80-83.

6. Куликов А.И., Смирнова И.И.. О тенденциях увлажнения степной зоны Байкальского региона // Устойчивое землепользование в экстремальных условиях / Материалы международной научно-практической конференции (26-28 ноября 2003). - Улан-Удэ, 2003. - С. 42-43

7.Убугунов Л. Л., Куликов А.И., Смирнова И.И., Мангатаев А.Ц. Региональные особенности опустынивания экосистем Забайкалья и подходы к оценке их устойчивости. // Глобальные и региональные особенности трансформации экосистем Байкальского региона./ Материалы симпозиума. - Улан-Батор, 2008. - С. 198-203.

Куликов А.И., Куликов М.А., Смирнова И.И. Термическое состояние деятельного слоя в криолитозоне Байкальского региона в контексте глобального потепления // Изменение климата Центральной Азии: социально-экономические и экологические последствия / Материалы международного симпозиума. -Чита, 2008. С. 171-178.

Подписано в печать 23.11.2009. Бумага офс. №1. Формат 60x84 1/16. Усл.печ.л. 1,5. Тираж 100. Заказ Цена договорная.

Издательство ФГОУ ВПО «Бурятская государственная сельскохозяйственная академия им. В. Р. Филиппова» 670034, г. Улан-Удэ, ул. Пушкина, 8 e-mail: rio-bqsha@maiI. ru