Известия Русского географического общества, 2021, T. 153, № 4, стр. 47-67

Уязвимость полевых культур к переувлажнению в условиях изменений климата в Нечерноземье Европейской России и варианты адаптации

М. В. Николаев *

Агрофизический научно-исследовательский институт
Санкт-Петербург, Россия

* E-mail: clenrusa@mail.ru

Поступила в редакцию 30.11.2020
После доработки 17.05.2021
Принята к публикации 20.05.2021

Полный текст (PDF)

Аннотация

Оценивается воздействие избытка атмосферной влаги на устойчивость растениеводства при ожидаемом в ближайшем будущем изменении климата в Нечерноземье Европейской России. Акцент сделан на полевых культурах, так как они имеют важное продовольственное значение и являются кормовой базой для ведения молочного животноводства. Их уязвимость к переувлажнению охарактеризована повторяемостью лет со значениями гидротермического коэффициента Г.Т. Селянинова (ГТК), превышающими 1.8, которые соответствуют внешним условиям для возникновения полегания. С использованием транзитивных моделей общей циркуляции атмосферы и океана (МОЦАО), обладающих полнотой описания физических механизмов и связей внутри модельных блоков, оценено изменение температуры воздуха, количества осадков и значений ГТК для десятилетия 2021–2030 гг. Показано, что при реализации сценария контролируемого выброса парниковых газов в атмосферу, отвечающего инициативам ООН для обеспечения устойчивого развития до 2030 г., наиболее уязвимой к эффектам переувлажнения оказывается северная часть Нечерноземья Европейской России. Появление резкого переувлажнения посевов также может иметь место и в средней части этой зоны. Обсуждается ряд адаптационных мер и стратегий, направленных на смягчение негативных проявлений изменения климата.

Ключевые слова: Нечерноземье Европейской России, полевые культуры, переувлажнение, транзитивные МОЦАО, оценка уязвимости, варианты адаптации

ВВЕДЕНИЕ

Высокоурожайные полевые культуры – продовольственные (озимая пшеница и озимая рожь), зернофуражные (озимый райграс, некоторые сорта пшеницы и ржи), силосные (кукуруза и подсолнечник на зеленую массу) и стеблевые сеяные травы на сенаж в силу их морфометрических параметров часто подвержены полеганию в избыточно влажных условиях, возникающих в Нечерноземье Европейской территории России (ЕТР) [18]. Переход на короткостебельные виды и сорта лишь отчасти обеспечивает достижение высоких и устойчивых урожаев, поскольку их биологическая продуктивность ниже: например, продуктивность травостоев из клевера и люцерны уступает продуктивности стеблевых злаковых травостоев более чем в два раза [5].

В то же время, многие из перечисленных культур составляют ценную кормовую базу для ведения молочного животноводства, обеспечивающего население высококачественной молочной продукцией. Прослеживающееся повышение температурного фона служит фактором, способствующим продвижению к северу посевов высокоурожайных культур и, следовательно, увеличению потенциала продуктивности сельскохозяйственных земель. Но в целом лимитирующими факторами для выращивания зерновых и кормовых культур в этой зоне, с одной стороны, остаются ограниченные ресурсы тепла, с другой – погодные условия, приводящие к избытку атмосферной влаги. Поэтому по агроклиматической классификации Г.Т.Селянинова, Нечерноземье ЕТР отнесено к зоне избыточного увлажнения со значениями ГТК ≥ 1.3 [21]. Приведенный по данным З.А. Шостак в [23] анализ типов погод в Нечерноземье ЕТР показал, что для озимой ржи в период ее созревания (июль) оптимальное количество осадков составляет 30 мм при средней температуре воздуха около 17.0°С, а для мягкой озимой пшеницы такими показателями являются: сумма осадков 50 мм и средняя температура воздуха немногим выше 16.0°С. Согласно данным из агроклиматических справочников, у силосных культур интенсивное листообразование происходит в июле при средней температуре воздуха от 16.5 до 19.0°С, а интенсивный рост стебля сенажных трав при средней температуре воздуха от 16.5 до 17.5°С, в обоих случаях при осадках от 40 до 50 мм [1]. То есть оптимум условий в период формирования урожая культур лежит в пределах 0.7 < ГТК VII < 1.1.

Однако при нарастающих изменениях климата потери их урожая из-за продолжительных и обильных дождей становятся сопоставимыми с потерями от появления бездождий11, но при этом существенно снижается и качество полученного урожая [3, 14, 23, 34]. Далее приводимые оценки показывают, что в ожидаемых климатических условиях уязвимость полевых культур к эффектам переувлажнения может возрасти. Поэтому понадобится применение комплекса адаптационных мер, чтобы обеспечить устойчивое производство растениеводческой продукции и сведение к минимуму ущерба от негативных проявлений изменения климата.

МЕТОДИКА ИССЛЕДОВАНИЯ

В пределах Нечерноземья ЕТР выбрано десять станций, относящихся к хвойно-широколиственно-лесной (Трубчевск), южно-таежно-лесной (Смоленск, Псков, Кострома и Вологда) и средне-таежно-лесной (Вытегра, Петрозаводск, Шенкурск, Котлас и Сыктывкар) зонам согласно природно-сельскохозяйственному районированию территории России [19]. Длительные метеорологические наблюдения на этих станциях обеспечили возможность составить непрерывные ряды средних месячных значений температуры воздуха и помесячных сумм осадков с 1945 по 2017 гг. На их основе рассчитаны значения гидротермического коэффициента Г.Т. Селянинова [21], рекомендуемого Росгидрометом и ВМО для проведения агроклиматических исследований и представляемого в общем виде:

${\text{ГТК}} = \frac{{\sum P }}{{0.1\sum T }},$
где: $\sum P $ – сумма выпавших осадков (мм) и $\sum T $ – сумма температур воздуха (°С) за период (или отрезок вегетации) со среднесуточной температурой воздуха выше 10°С.

Для оценки отрицательного воздействия избыточного увлажнения на посевы удобно воспользоваться значениями ГТК, превышающими его среднюю многолетнюю величину. По долговременным обследованиям состояния посевов установлено, что надежным индикатором условий возникновения их корнево-стеблевого полегания служит пороговое значение ГТК = 1.8 [18]. Следуя этому критерию, а также исходя из сравнения значений ГТК в июле, когда формируется конечный урожай мелкозернистых озимых злаков и происходит интенсивное накопление зеленой массы силосных и сенажных культур, для характеристики степени полегания нами разработаны градации ГТК: 1.8 < ГТК VII ≤ 2.5 – выраженное полегание; 2.5 < ГТК VII ≤ 3.5 – сильное полегание; 3.5 < ГТК VII ≤ 4.5 – очень сильное полегание [15]. Градации отражают усиливающееся воздействие обильных осадков на степень полегания культур из-за возникновения условий, которые одновременно приводят к разжиженному состоянию поверхности почвы и увеличению механической нагрузки на стебель. Такому состоянию почвы способствует преобладание суглинистых почв со свойственной им слабой водопроницаемостью и часто близкое к поверхности расположение зеркала грунтовых вод [20].

Эмпирический анализ данных с применением математических методов их обработки позволил выявить наблюдаемые тенденции в изменении температуры воздуха и количества осадков в теплое время года в пределах Нечерноземья ЕТР. Так, в изменяющемся климате наблюдается увеличение средней месячной температуры воздуха в июле. Для временного отрезка 1982–2017 гг. (по сравнению с временным отрезком 1945–1981 гг.) такое увеличение для выбранных станций составляет от 0.6 до 1.3°C, тогда как для других летних месяцев отличия в средних значениях температур менее выражены. Для станций также проведено сравнение хода 7-летних скользящих средних значений температуры воздуха в июле за весь рассматриваемый период. Такое сравнение показало, что условия тепло-обеспеченности на северных станциях в последние два десятилетия приближаются к условиям теплообеспеченности на более южных станциях в ранние десятилетия. Например, по величине сумм температур условия в Петрозаводске приближаются к условиям в Костроме, в Вытегре к условиям в Смоленске, а в Пскове к условиям в Трубчевске [17].

Вместе с тем, в Нечерноземье ЕТР отмечается учащение экстремально высокого количества выпадающих осадков, включая июльские их суммы, из-за усиления циклонической активности в северной Атлантике [710, 29]. Это находит преломление в наблюдаемом расширении территорий, подверженных атмосферному переувлажнению и изменении их конфигурации [11, 13]. Сопоставление гистограмм отклонений сумм сезонных осадков от тренда показало, что такие гистограммы для летних осадков в последующие после 1979 г. десятилетия отличаются положительной асимметрией, а рост величины стандартного отклонения обусловлен возрастающей частотой экстремально высоких их сумм [15]. С помощью функционального анализа, при котором ГТК представляется в виде функции от двух переменных, выявлен усиливающийся вклад обильных осадков в высокие значения ГТК в июле по сравнению с вкладом сумм температур. В свою очередь, это проявляется в увеличении тесноты регрессионных связей между возрастающими суммами июльских осадков и высокими значениями ГТК вследствие снижения остаточной дисперсии [17]. Однако в северной части Нечерноземья ЕТР очень высокие значения ГТК за июль фиксируются при суммах температур заметно выше их средних многолетних показателей [16]. Из сводных таблиц, помещенных в [15, 17], также следует, что процент лет с эффектами сильного и очень сильного полегания, отвечающими пороговым значениям ГТК выше 2.5 и 3.5, значительно возрос в среднетаежной зоне, для которой характерны подзолистые суглинистые почвы с признаками оглеения.

Оценка ожидаемой повторяемости переувлажнения посевов в следующем десятилетии основывается на результатах численного моделирования климата по транзитивным моделям общей циркуляции атмосферы и океана (МОЦАО). При выборе МО-ЦАО предпочтение отдано таковым, которые отличаются полнотой описания физических механизмов и связей внутри модельных блоков. Отметим, что такими качествами обладают семейства моделей, созданных в известных центрах моделирования глобального климата – Метеорологическом отделении Института им. М. Планка (Германия), Метеорологическом офисе Центра им. Дж. Хэдли (Великобритания) и Канадском центре климатического моделирования и анализа, а также модель Европейского Сообщества Earth3, широко применяемая в странах Европы в последние годы. Размеры и форма ячеек регулярных сеток моделей несколько отличаются. Последние версии моделей Института им. М. Планка имеют размер ячейки регулярных сеток по градусам широты и долготы 0.935° × 0.9375°; Канадского центра – 2.79° × 2.8125°; Центра им. Дж. Хэдли – 1.25° × 1.875° (вытянутая ячейка). У европейской модели размер ячейки наименьший и составляет 0.702° × 0.703°. Между тем, использование моделей, созданных в других странах мира, оказывается затруднительным, так как их регулярные сетки лишь частично покрывают исследуемую нами территорию: например, модель CESM2 – WACCM Национального центра атмосферных исследований (NCAR, США), модель AСCESS – ESM1.5 Организации Содружества Научных и Индустриальных Исследований (CSIRO, Австралия) и модель MRI – ESM-2-0 Метеорологического научно-исследовательского института (Meteorological Agency, Япония).

Подробное описание выбранных моделей приводится в [3032, 35]. Интерполяция модельных оценок из узлов регулярных сеток моделей в координаты станций осуществлена с помощью билинейной интерполяции. Модельные оценки отнесены для сценария контролируемого выброса парниковых газов в атмосферу RCP4.522 [25].

РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ

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

• MPI – ESM1-2 – HR (модель Института им. М. Планка, модельные эксперименты 1 и 2);

• Нad GEM3- GC31 – LL (модель Центра им. Дж. Хэдли, модельный эксперимент 1);

• EC – Earth3 ( модель Европейского Сообщества, модельные эксперименты 9 и 15);

• Can ESM5 (модель Канадского центра климатического моделирования и анализа, модельные эксперименты 14, 23).

Предпочтение в выборе конкретных модельных экспериментов обусловлено полнотой представления оценок в моделируемых рядах текущего климата (1945–2014 гг.), достаточной длительностью моделируемых рядов на перспективу и очень незначительными, но все же существующими различиями в вариантах описания физических связей в модельных блоках. По британской модели пока выполнен один эксперимент, который позволяет получить оценки на перспективу после 2020 г. Но он дает возможность оценить значения температуры воздуха и ГТК в Смоленске и Трубчевске только для 6 лет внутри десятилетия 2021–2030 гг. по причине пропусков данных о температуре в узлах регулярной сетки. Так как по модели MPI – ESM1-2 – HR представлены результаты только двух экспериментов, для сбалансированности получаемых оценок также выбрано по два эксперимента по моделям EC – Earth3 и Can ESM5, но исходя из упомянутых выше критериев их отбора. Иначе говоря, какой-либо из моделей изначально не отдается приоритета.

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

Рис. 1.

Тенденции в изменении и вариабельности средней месячной температуры воздуха, количества осадков и значений ГТК в июле в Петрозаводске, оцененные по модели MPI – ESM1-2-HR (Институт им. М. Планка), и их проекции исходя из сценария контролируемого выброса парниковых газов в атмосферу. а – смоделированные ряды с 1945 по 2014 гг.; б – проекция на 2015–2036 гг. по модельному эксперименту 1; в – проекция на 2015–2036 гг. по модельному эксперименту 2; точки – ежегодные оценки; линии – 3‑летние скользящие средние оценки. Fig. 1. Trends in change and variability of mean monthly air temperatures, precipitation totals and HTC values in July at Petrozavodsk simulated by MPI – ESM1-2-HR model (Max Planck Institute) and their projections according to the scenario of controlled greenhouse gas emission into the atmosphere. a – modeled series since 1945 till 2014; б – projection from 2015 to 2036 by model run 1; в – projection from 2015 to 2036 by model run 2; dotes – year to year estimates; lines – 3-yr moving average estimates.

Результаты обоих экспериментов показывают, что условия увлажнения в июле в Петрозаводске на десятилетии 2021–2030 гг. становятся более экстремальными по сравнению с 70-летним периодом с 1945 по 2014 гг., но с преобладанием “выбросов” в сторону очень высоких сумм осадков и значений ГТК при повышающемся температурном фоне.

Поскольку модельные эксперименты основаны на единых физических принципах описания связей в модельных блоках, а в рядах климатических переменных присутствует стохастичность, целесообразно объединить результаты экспериментов для каждой модели, тем самым увеличив количество случаев, и представить оценки в вероятностной интерпретации. На рис. 2–4 изображены гистограммы (с их сглаженным методом сверстки видом) смоделированных оценок для двух периодов, отличающихся по степени антропогенного влияния на климат [2, 4, 6, 24] и для десятилетия 2021–2030 гг. с учетом выбранного эмиссионного сценария. Примеры расчетов взяты произвольно, но таким образом, чтобы наглядно продемонстрировать оценки по разным моделям, охватив всю природно-сельскохозяйственную зональность Нечерноземья ЕТР.

Рис. 2.

Гистограммы распределений средней месячной температуры воздуха, количества осадков и ГТК в июле для станции Котлас, полученные по модели Can ESM5 (Канадский центр климатического моделирования и анализа): а – период 1945–1979 гг.; б – период 1980–2014 гг.; в – десятилетие 2021–2030 гг. (объединение экспериментов 14 и 23). Fig. 2. Histograms of mean monthly air temperatures, precipitation totals and HTC values in July at Kotlas based on the simulation by Can ESM5 model (Canadian Center for Climate Modelling and Analysis): а – interval 1945–1979; б – interval 1980–2014; в – slice 2021–2030 (conjunction of run 14 and run 23).

Рис. 3.

Гистограммы распределений средней месячной температуры воздуха, количества осадков и ГТК в июле для станции Вологда, полученные по модели MPI – ESM1-2-HR (Институт им. М. Планка): а – период 1945–1979 гг.; б – период 1980–2014 гг.; в – десятилетие 2021–2030 гг. (объединение экспериментов 1 и 2). Fig. 3. Histograms of mean monthly air temperatures, precipitation totals and HTC values in July at Vologda based on the simulation by MPI – ESM1-2-HR model (Max Planck Institute): а – interval 1945–1979; б – interval 1980–2014; в – slice 2021–2030 (conjunction of run 1 and run 2).

Рис. 4.

Гистограммы распределений средней месячной температуры воздуха, количества осадков и ГТК в июле для станции Трубчевск, полученные по модели EC – Earth3 (Европейское Сообщество): а – период 1945–1979 гг.; б – период 1980–2014 гг.; в – десятилетие 2021–2030 гг. (объединение экспериментов 9 и 15). Fig. 4. Histograms of mean monthly air temperatures, precipitation totals and HTC values in July at Trubchevsk based on the simulation by EC–Earth3 model (European Community): а – interval 1945–1979; б – interval 1980–2014; в – slice 2021–2030 (conjunction of run 9 and run 15).

На рис. 2 приведены гистограммы смоделированных для периодов 1945–1979 гг., 1980–2014 гг. и десятилетия 2021–2030 гг. распределений средней месячной температуры воздуха, количества осадков и ГТК в июле по модели Канадского центра климатического моделирования и анализа Can ESM5 для станции Котлас (среднетаежная зона). Распределения для десятилетия 2021–2030 гг. построены путем объединения результатов двух модельных экспериментов – 14 и 23.

Из графиков видно, что распределения средней месячной температуры воздуха в июле для всех трех временных отрезков имеют почти симметричный вид, но при смещении моды в сторону более высокой температуры с появлением диапазона температурного максимума для десятилетия 2021–2030 гг. В распределении количества осадков для периода 1980–2014 гг. прослеживается асимметрия в сторону их увеличения. Такая асимметрия становится выраженной для десятилетия 2021–2030 гг. с одновременным появлением диапазона очень высоких сумм июльских осадков. Для распределения ГТК на этом же временном отрезке характерно заметное увеличение частоты экстремально высоких значений ГТК (до 0.2). Для периодов 1945–1979 и 1980–2014 гг. их частота соответственно составляет 0.07 и 0.1 .

Нужно отметить, что характер распределений отражает изменяющийся во времени вклад влажностного и термического факторов в величину экстремально высоких значений ГТК, и это подтверждается из данных наблюдений. Так, на станции Котлас в течение периода 1945–1979 гг. экстремально высокое значение ГТК = 4.18 в июле 1951 г. фиксируется при сумме температур, равной 450°C (85% средней многолетней величины) и сумме осадков, равной 188 мм (261% нормы). В течение периода 1980–2014 гг. экстремально высокое значение ГТК = 4.14 в июле 2000 г. достигается при сумме температур, равной 603°C (113% средней многолетней величины), и сумме осадков, равной 250 мм (347% нормы). Увеличение количества осадков на повышенном термическом фоне по сравнению с их количеством на пониженном термическом фоне, по-видимому, объяснимо дополнительным вкладом конвекции. Но все же роль адвективных осадков, обусловленных циклоничностью, остается определяющей, поскольку высокие значения ГТК в июле чаще отмечаются на пониженном термическом фоне.

Аналогичным образом, на рис. 3 показаны гистограммы смоделированных распределений средней месячной температуры воздуха, количества осадков и ГТК в июле по модели Института им. М. Планка MPI – ESM1-2-HR для станции Вологда (южно-таежно-лесная зона) с объединением для десятилетия 2021–2030 гг. результатов двух модельных экспериментов – 1 и 2.

В распределениях средней месячной температуры воздуха в июле также выявляется смещение моды в сторону более высокой температуры. Все распределения количества осадков показывают увеличение частоты высоких июльских осадков, а для десятилетия 2021–2030 гг. вероятность появления очень высоких сумм осадков. Это также отражается и в асимметрии распределения значений ГТК в июле для десятилетия 2021–2030 гг. с появлением диапазона значений ГТК > 3.5, характеризующегося вероятностью 0.1.

По данным наблюдений в течение периода 1945–1979 гг. появление экстремально высокого значения ГТК = 3.06 в июле 1954 г. соответствовало сумме температур 588°C (111% средней многолетней величины) и сумме осадков 180 мм (247% нормы). В течение периода 1980–2014 гг. экстремально высокое значение ГТК = 3.18 в июле 1998 г. отвечало сумме температур 516°C (99% средней многолетней величины) и сумме осадков 165 мм (226% нормы).

Рис. 4 иллюстрирует гистограммы смоделированных распределений средней месячной температуры воздуха, количества осадков и ГТК в июле по модели Европейского Сообщества EC – Earth3 для станции Трубчевск (зона хвойно-широколиственных лесов) с объединением для десятилетия 2021–2030 гг. результатов двух модельных экспериментов – 9 и 15.

В распределениях средней месячной температуры воздуха в июле в целом присутствует смещение моды в сторону более высокой температуры. Однако вид гистограммы для периода 1980–2014 гг. существенно отличается от ее вида для десятилетия 2021–2030 гг. В первом случае наблюдается преобладание диапазонов с температурой воздуха ниже 20°C, тогда как во втором случае, наоборот – с температурой воздуха выше 20°C при отсутствии диапазонов с температурой воздуха ниже 18°C. Распределения количества июльских осадков показывают смещение моды в сторону высоких их сумм при прослеживающейся бимодальности для десятилетия 2021–2030 гг. за счет заметного увеличения частоты сумм осадков в диапазоне от 20 до 40 мм и в диапазоне от 80 до 100 мм. Распределение ГТК для десятилетия 2021–2030 гг. имеет более плоскую вершину вследствие равенства частот для диапазонов ГТК < 0.5 и 1.0 < ГТК < 2.0.

В течение периода 1945–1979 гг. экстремально высокое значение ГТК = 3.31 в июле 1979 г. отмечалось при сумме температур, равной 477°C (83% средней многолетней величины), и сумме осадков, равной 158 мм (193% нормы). В течение периода 1980–2014 гг. экстремально высокое значение ГТК = 3.55 в июле 1984 г. соответствовало сумме температур, равной 489°C (85% средней многолетней величины), и сумме осадков, равной 174 мм (212% нормы). Близкое количество осадков, равное 170 мм, также наблюдалось в июле 1999 г., но при сумме температур, равной 642°C (111% средней многолетней величины).

На повышающемся температурном фоне тенденции к росту частоты высоких осадков и значений ГТК для десятилетия 2021–2030 гг. (по сравнению с предшествующими двумя 35-летними периодами) прослеживаются и для других станций. Так, по модели Can ESM5 для Шенкурска обнаруживается появление дополнительного диапазона осадков от 200 до 220 мм и значений ГТК от 3.5 до 4.0 (при вероятности 0.1). По модели MPI – ESM1-2-HR для Петрозаводска предсказывается появление двух дополнительных диапазонов осадков от 180 до 200 мм и от 200 до 220 мм и значений ГТК от 3.0 до 3.5 (при вероятности 0.3) и от 3.5 до 4.0 (при вероятности 0.1). По модели EC– Earth3: для Сыктывкара прогнозируется появление дополнительного диапазона осадков от 180 до 200 мм и увеличение частоты значений ГТК от 3.0 до 3.5 (при вероятности 0.3); для Вытегры – увеличение частоты осадков в диапазоне от 140 до 160 мм и значений ГТК от 2.0 до 2.5 (при вероятности 0.5); для Костромы – увеличение частоты осадков в диапазоне от 120 до 140 мм и значений ГТК от 2.0 до 2.5 (при вероятности 0.4).

Но ряд модельных экспериментов не позволяет выявить тенденции к росту частоты высоких сумм осадков и значений ГТК для десятилетия 2021–2030 гг. Например, по модели Нad GEM3-GC3.1 – LL в Вологде предсказывается учащение сумм осадков в диапазоне от 60 до 80 мм и значений ГТК от 1.0 до 1.5 (при вероятности 0.7). Модель Can ESM5 для Пскова выявляет рост повторяемости еще более низких сумм осадков в диапазоне от 40 до 60 мм и значений ГТК от 0.5 до 1.0 (при вероятности 0.8). Модель MPI – ESM1-2-HR для Трубчевска также показывает учащение повторяемости сумм осадков в диапазоне от 40 до 60 мм и значений ГТК от 0.5 до 1.0 (при вероятности 0.5).

Для того, чтобы судить о достоинствах или недостатках каждой модели, которые состоят в степени отличия описываемого текущего климата от действительного, осуществлена верификация выбранных моделей. Для оценки точности воспроизведения моделями характеристик реального климата в табл. 1–3 приведены разности между смоделированными и фактическими значениями осредненных климатических переменных за период с 1945 по 2014 гг. (за 70 лет) для станций Нечерноземья ЕТР.

Таблица 1.

Точность воспроизведения нормы июльских осадков моделями (разность между смоделированными и фактическими величинами; минимальные отклонения по станциям выделены жирным шрифтом) Table 1. The reproduction accuracy of the July precipitation norm by models (the difference between simulated and real values; the minimum deviations by station are shown in bold)

Станция Широта Долгота Реальная норма (P), мм Разность
P(P(MPI-ESM1-2 – HR) – P, мм
Разность
P(EC- Earth3) – P, мм
Разность
P(Нad GEM3- GC3.1- LL) – P, мм
Разность
P(Can ESM5) – P,
мм
Петрозаводск 61.8° с.ш. 34.3° в.д. 73.65 13.73 –2.16 33.77 21.95
Вытегра 61.0° с.ш. 36.4° в.д. 76.63 23.33 1.06 36.40 22.09
Шенкурск 62.1° с.ш. 42.9° в.д. 68.21 26.49 19.00 30.33 28.15
Котлас 61.2° с.ш. 46.7° в.д. 72.15 23.71 13.32 25.89 31.78
Сыктывкар 61.7° с.ш. 50.8° в.д. 71.33 26.59 11.79 24.90 24.41
Вологда 59.3° с.ш. 39.9° в.д. 72.67 27.77 16.17 27.07 37.41
Кострома 57.7° с.ш. 40.8° в.д. 73.32 38.15 3.96 20.74 41.44
Псков 57.8° с.ш. 28.4° в.д. 75.71 26.04 0.60 14.34 18.74
Смоленск 54.5° с.ш. 32.9° в.д. 90.99 10.01 –14.73 –12.82 22.80
Трубчевск 52.6° с.ш. 33.8° в.д. 82.01 2.59 –20.70 –16.78 13.85
Таблица 2.  

Точность воспроизведения средней многолетней температуры воздуха в июле моделями (разность между смоделированными и фактическими величинами; минимальные отклонения по станциям выделены жирным шрифтом) Table 2. The reproduction accuracy of the mean annual air temperature in July by models (the difference between simulated and real values; the minimum deviations by station are shown in bold)

Станция Широта Долгота Реальная средняя (T), °C Разность
T(MPI M1-2 – HR) – – T, °C
Разность
T(EC- Earth3) – T, °C
Разность
T(Нad GEM3- GC3.1- LL) – T, °C
Разность
T(Can ESM5) – T, °C
Петрозаводск 61.8° с.ш. 34.3° в.д. 16.38 0.44 –1.65 –1.54 0.68
Вытегра 61.0° с.ш. 36.4° в.д. 16.93 –0.33 –1.85 –1.83 0.12
Шенкурск 62.1° с.ш. 42.9° в.д. 17.37 –0.90 –2.60 –3.07 0.39
Котлас 61.2° с.ш. 46.7° в.д. 17.15 –0.52 –1.94 –1.88 0.39
Сыктывкар 61.7° с.ш. 50.8° в.д. 17.08 0.32 –2.15 –1.48 0.56
Вологда 59.3° с.ш. 39.9° в.д. 17.14 0.39 –0.89 –1.42 0.67
Кострома 57.7° с.ш. 40.8° в.д. 18.07 –0.79 –0.96 –1.64 0.32
Псков 57.8° с.ш. 28.4° в.д. 17.78 –1.29 –1.36 –2.04 0.32
Смоленск 54.5° с.ш. 32.9° в.д. 17.38 0.48 0.93 0.43 1.37
Трубчевск 52.6° с.ш. 33.8° в.д. 18.62 0.10 0.96 –0.72 1.03
Таблица 3.  

Точность воспроизведения среднего многолетнего значения ГТК в июле моделями (разность между смоделированными и фактическими значениями; минимальные отклонения по станциям выделены жирным шрифтом) Table 3. The reproduction accuracy of the mean annual HTC value in July by models (the difference between modeled and real values; the minimum deviations by station are shown in bold)

Станция Широта Долгота Реальное среднее ГТК(H) Разность
H(MPI M1-2 – HR) – – H
Разность
H(EC- Earth3) –H
Разность
H(Нad GEM3- GC3.1- LL) – H
Разность
H(Can ESM5) – H
Петрозаводск 61.8° с.ш. 34.3° в.д. 1.48 0.22 0.13 0.92 0.34
Вытегра 61.0° с.ш. 36.4° в.д. 1.48 0.50 0.23 0.99 0.40
Шенкурск 62.1° с.ш. 42.9° в.д. 1.30 0.61 0.66 1.04 0.55
Котлас 61.2° с.ш. 46.7° в.д. 1.38 0.54 0.48 0.74 0.54
Сыктывкар 61.7° с.ш. 50.8° в.д. 1.39 0.54 0.46 0.59 0.37
Вологда 59.3° с.ш. 39.9° в.д. 1.40 0.58 0.42 0.68 0.61
Кострома 57.7° с.ш. 40.8° в.д. 1.36 0.76 0.16 0.53 0.66
Псков 57.8° с.ш. 28.4° в.д. 1.40 0.62 0.15 0.47 0.36
Смоленск 54.5° с.ш. 32.9° в.д. 1.72 0.13 –0.32 –0.21 0.25
Трубчевск 52.6° с.ш. 33.8° в.д. 1.46 0.02 –0.40 –0.31 0.14

Как видно из данных табл. 1, наименьшие расхождения между оцененной и фактической нормой июльских осадков достигаются по модели EC – Earth3 (для некоторых станций погрешность расчета составляет менее 5%). Модель MPI – ESM1-2-HR обеспечивает наибольшую близость смоделированных норм июльских осадков их фактическим нормам только для двух южных станций.

Если сравнить смоделированное количество осадков с наблюдаемым их количеством в последующие годы: 2015, 2016 и 2017 гг., то также можно обнаружить, что наиболее часто очень близкое (менее ±10 мм) или полное совпадение смоделированных и наблюденных сумм осадков достигается моделью EC–Earth3 – 14 случаев. Для модели MPI – ESM1-2-HR число таких случаев – 9; для модели Can ESM5 – 8 и для модели Нad GEM3-GC3.1 – LL.

Табл. 2 показывает, насколько отличаются смоделированные и реальные величины средней многолетней температуры воздуха в июле. Из данных этой таблицы следует, что наименьшие расхождения в величинах температуры воздуха обеспечиваются моделями Can ESM5 и MPI – ESM1-2-HR (в среднем погрешность находится в пределах ±0.5–0.7°C). Но если модель Can ESM5 показывает преобладание положительных отклонений, то модель MPI – ESM1-2-HR – наоборот, отрицательных отклонений. Модели Нad GEM3-GC3.1 – LL и EC – Earth3 температуру воздуха несколько занижают, в особенности модель Нad GEM3-GC3.1 – LL.

Сравнение по трем последующим годам после 2014 г. также показало, что очень близкое (менее 0.5°C) или полное совпадение между смоделированными и наблюдаемыми величинами июльской температуры воздуха достигаeтся моделью Can ESM5 в 7 случаях, а моделью MPI – ESM1-2-HR в трех случаях. По модели EC – Earth3 в пяти случаях расхождения составляют менее 1°C , в то время как модель Нad GEM3-GC3.1 – LL дает более выраженную погрешность.

В табл. 3 сведены невязки значений ГТК – показателя влаго- (тепло) обеспеченности культур, вклад количества осадков в величину которого является определяющим. Видно, что наименьшие расхождения в значениях ГТК выявляются по моделям EC–Earth3 и MPI – ESM1-2-HR.

За 2015, 2016 и 2017 гг. несущественные отличия (менее 0.1) или совпадения смоделированных и рассчитанным по данным наблюдений значений ГТК за июль обеспечиваются моделью EC–Earth3 в 13 случаях, моделью Can ESM5 в 9 случаях, моделью MPI – ESM1-2-HR в 6 случаях и моделью Нad GEM3-GC3.1 – LL в трех случаях.

В целом из приведенных в табл. 1–3 отклонений можно заключить, что моделями преимущественно воспроизводятся более влажные климатические условия, чем в реальности. В этом отношении представляет интерес проверить, в какой мере описываемое моделями экстремально высокое количество осадков в июле соответствует наблюдаемому максимуму июльских осадков за весь рассматриваемый период с 1945 по 2017 гг. Проверка такого соответствия показала, что для Котласа, Вытегры и Смоленска максимумы смоделированных осадков всегда оказываются меньше наблюденных максимумов. Для Петрозаводска в четырех случаях, а для Костромы и Трубчевска в одном случае смоделированные максимумы превышают действительные максимумы только по модели Нad GEM3-GC31 – LL. Завышения фактических максимумов сумм осадков для остальных станций распределяются следующим образом: моделью EC–Earth3 – 2 раза (по одному для Шенкурска и Вологды); моделью MPI – ESM1-2-HR – 5 раз (для Шенкурска и Сыктывкара – дважды, для Пскова – один раз); моделью Can ESM5 – 5 раз (для Пскова – дважды, для Шенкурска, Сыктывкара и Вологды – по одному разу) и моделью Нad GEM3-GC31 – LL – 4 раза (для Шенкурска – дважды, для Вологды и Пскова – по одному разу).

Для оценок уязвимости полевых культур к переувлажнению в ближайшем будущем выбраны пороговые значения ГТК за июль, которые отвечают условиям возникновения разной степени полегания – соответственно: ГТК VII > 1.8, ГТК VII > 2.5 и ГТК VII > 3.5 . Затем рассчитана доля лет с такими условиями, характеризующая повторяемость эффектов переувлажнения посевов. В табл. 4 представлено сравнение повторяемости лет с эффектами переувлажнения по данным, основанным на наблюдениях, с повторяемостью, оцененной по МОЦАО и спроецированной для десятилетия 2021–2030 гг.

Таблица 4.  

Повторяемость эффектов переувлажнения посевов в периоды 1945–1979, 1980–2014 гг. и для десятилетия 2021–2030 гг. в разных частях Нечерноземья ЕТР (% лет) Table 4. Recurrence of crop over-wetting effects within intervals of 1945–1981, 1982–2017 and for slice 2021–2030 in different parts of the Non-chernozem region of the European Russia (% of years)

Станция;
координаты; субъект РФ
Значения
ГТКVII
Реальные значения Значения, рассчитанные по моделям
MPI-ESM1-2-HR EC-Earth3 HadGEM3-GC3.1-LL CanESM5
1945–1979 1980–2014 1945–1979 1980–2014 2021–2030 1945–1979 1980–2014 2021–2030 1945–1979 1980–2014 2021–2030 1945–1979 1980–2014 2021–2030
Петрозаводск >1.8 29 26 37 54 50 39 35 50 77 59 50 47 51 40
61.8° с.ш.; 34.3° в.д. >2.5 11 9 9 17 20 10 10 20 49 35 20 17 11 8
Республика Карелия >3.5 3 0 0 0 5 2 1 0 9 9 10 1 0 0
Вытегра >1.8 26 40 63 51 55* 41 42 45 89 71 44 49 56 35
61.0° с.ш.; 36.4° в.д. >2.5 11 14 20 17 35 10 22 10 63 29 22 13 16 5
Вологодская обл., север >3.5 0 0 6 3 0 1 2 0 11 13 0 0 0 0
Шенкурск >1.8 20 26 51 46 35 55 54 40 69 59 67* 43 47 40
62.1° с.ш.; 42.9° в.д. >2.5 3 6 26 31 10 19 23 20 46 38 22 20 13 8
Архангельская обл. >3.5 0 0 0 6 0 2 3 0 9 6 0 1 0 2
Котлас >1.8 20 23 54 43 35 54 48 40 61 50 40 54 51 55
61.2° с.ш.; 46.7° в.д. >2.5 14 9 31 23 10 16 19 10 45 25 20 13 16 18
Архангельская обл., юг >3.5 3 3 6 9 0 2 5 0 6 7 20 0 0 2
Сыктывкар >1.8 23 37 54 60 40 56 46 50* 72 40 33 47 41 48
61.7° с.ш.; 50.8° в.д. >2.5 6 9 31 23 10 18 19 25 28 16 0 14 4 8
Республика Коми >3.5 3 0 6 3 0 2 3 0 7 4 0 1 0 0
Вологда >1.8 23 31 60 57 40 49 48 30 65 52 29 69 60 40
59.3° с.ш.; 39.9° в.д. >2.5 9 20 26 26 20 17 19 10 32 21 0 14 20 8
Вологодская обл., юг >3.5 0 0 0 0 5 1 6 0 6 0 0 3 1 0
Кострома >1.8 29 31 66 60 55 34 32 40 54 47 33 73 61 50
57.7° с.ш.; 40.8° в.д. >2.5 14 6 40 26 20 11 12 5 20 26 0 23 13 5
Костромская обл. >3.5 6 3 3 3 0 0 1 0 0 9 0 0 0 0
Псков >1.8 37 20 54 63 45 36 31 35* 49 46 20 51 39 32
57.8° с.ш.; 28.4° в.д. >2.5 9 9 29 20 20 10 12 5 23 17 10 17 11 10
Псковская обл. >3.5 0 0 3 0 5 0 4 0 0 3 0 4 1 0
Смоленск >1.8 43 43 54 46 40 24 34 30 33 29 0 70 51 55*
54.5° с.ш.; 32.9° в.д. >2.5 23 17 31 11 10 9 10 5 8 7 0 20 14 12
Смоленская обл. >3.5 6 9 3 0 0 1 1 0 8 0 0 0 0 2
Трубчевск >1.8 34 31 37 23 30* 10 15 15 17 13 0 47 24 22
52.6° с.ш.; 33.8° в.д. >2.5 20 11 6 6 5 1 5 5 6 0 0 7 4 2
Брянская обл. >3.5 0 0 0 0 0 0 0 0 0 0 0 0 0 0

Примечание: жирным шрифтом выделены проценты лет с эффектами переувлажнения, превышающие для десятилетия 2021–2030 гг. таковые для двух прошлых периодов; знаком * выделены проценты лет с эффектами переувлажнения, превышающие для десятилетия 2021–2030 гг. таковые для периода 1980–2014 гг.

Результаты сравнения показывают, что по разным моделям достигается разное расхождение между полученными оценками повторяемости лет с эффектами переувлажнения и реальной их повторяемостью. Например, модель Центра им. Хэдли Нad-GEM3-GC3.1 – LL дает значительное завышение доли таких лет для северных станций и значительное занижение для южных станций. В то же время, модель Европейского Сообщества EC – Earth3 обеспечивает достаточно высокую согласованность оцененного и реального процента лет для ряда станций в разных частях Нечерноземья ЕТР.

В табл. 4 выделены случаи, когда для десятилетия 2021–2030 гг. повторяемость лет с эффектами переувлажнения превышает их повторяемость для периодов 1945–1979 гг. и 1980–2014 гг. при одновременной близости оцененного по моделям и реального процентов лет для этих периодов. Для порового значения ГТК > 1.8 такие случаи выявляются по модели Европейского Сообщества EC – Earth3 для двух северных станций (Петрозаводск – 50% лет; Вытегра – 45% лет) и Костромы (40% лет). Кроме того, по модели Канадского центра климатического моделирования и анализа Can ESM5 предсказывается увеличение процента таких лет в Котласе (до 55% лет) и Сыктывкаре (до 48% лет).

Возрастающую повторяемость лет с эффектами резкого переувлажнения (ГТК > 2.5) на десятилетии 2021–2030 гг. демонстрируют три модели. Для Петрозаводска модель Европейского Сообщества EC – Earth3 и модель Института им. М. Планка MPI – ESM1-2 – HR дают одинаковую оценку: 20% лет. Для других станций оценки соответственно составляют: для Вытегры по модели Института им. М. Планка MPI – ESM1-2 – HR – 35% лет; для Сыктывкара по модели Европейского Сообщества EC – Earth3 – 25% лет; для Котласа по модели Канадского центра климатического моделирования и анализа Can ESM5 – 18% лет.

По модели Института им. М. Планка MPI – ESM1-2 – HR для десятилетия 2021–2030 гг. также индицируется пятипроцентная повторяемость лет с очень резким переувлажнением (ГТК > 3.5) в Петрозаводске, Вологде и Пскове, а по модели Канадского центра климатического моделирования и анализа Can ESM5 двухпроцентная повторяемость таких же лет в Шенкурске. Модель Центра им. Хэдли Нad GEM3-GC3.1 – LL демонстрирует возрастание повторяемости сверхвлажных лет в Петрозаводске и Котласе, но нужно напомнить, что оценки по этой модели основываются на результате лишь одного модельного эксперимента.

Для двух южных станций увеличение повторяемости лет с эффектами переувлажнения для десятилетия 2021–2030 гг. (по сравнению с двумя прошлыми периодами) не прослеживается, хотя модель Европейского Сообщества EC – Earth3 и показывает равенство процента лет в Трубчевске для периода 1980–2014 гг. и десятилетия 2021–2030 гг.

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

ВЫВОДЫ

Результаты выполненного исследования свидетельствуют о возможностях использования МОЦАО в качестве инструмента для оценок устойчивости полеводства в условиях будущего климата в целях упреждающей адаптации. Различия в полученных по моделям оценках объясняются тем, что выбранные модели отличаются между собой как пространственной разрешающей способностью, так и детальностью описания механизмов обратных связей, которые уточняются в процессе проведения модельных экспериментов. Моделируемые ими условия влагообеспеченности культур также имеют отличие от действительных – вероятно, по причине неточностей в описании процессов осадкообразования и переноса влаги. Это в итоге приводит к тому, что оцененная по моделям повторяемость лет с возникновением переувлажнения имеет расхождения с их повторяемостью по данным наблюдений. Но для высоких пороговых значений ГТК, которые индицируют резкое переувлажнение, достигается значительно лучшая сопоставимость.

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

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

В этой связи, ориентируясь на отечественный и зарубежный опыт ведения рациональной практики возделывания полевых культур в гумидных регионах [5, 12, 22, 2628, 33] и принимая во внимание полученные на перспективу оценки, к мерам по смягчению негативных проявлений изменения климата следует отнести:

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

• эффективное применение дренажных систем – увеличение пропускной способности дренажа закрытого типа и широкое применение систем открытого дренажа вследствие более частых и сильных дождей;

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

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

Одновременно в качестве адаптационных стратегий могут выступать:

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

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

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

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

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

  1. Агроклиматические ресурсы Архангельской, Брянской, Вологодской, Костромской, Псковской, Смоленской областей; Карельской АССР и Коми АССР (справочники) Л.: Гидрометеоиздат, 1971–1974.

  2. Второй оценочный доклад Росгидромета об изменениях климата и их последствиях на территории Российской Федерации. Общее резюме. М.: Росгидромет, 2014. 58 с.

  3. Грингоф И.Г., Клещенко А.Д. Основы сельскохозяйственной метеорологии. Обнинск: ФГБУ “ВНИИГМИ-МЦД”, 2011. Т. I. С. 601–609.

  4. Доклад Росгидромета об особенностях климата на территории Российской Федерации за 2018 г. М.: Росгидромет, 2019. С. 10–14, 18–20.

  5. Зотов А.А., Косолапов В.М., Кобзин А.Г., Трофимов И.А., Уланов А.Н., Шевцов А.В. Шельменкина Х.Х., Щукин Н.Н. Сенокосы и пастбища на осушенных землях Нечерноземья. Москва–Астана, 2012. 1198 с.

  6. Катцов В.М., Школьник И.М., Ефимов С.В. Перспективные оценки изменений климата в российских регионах: детализация в физических и вероятностных пространствах // Метеорология и гидрология. 2017. № 6. С. 68–81.

  7. Киктев Д.Б., Круглова Е.Н., Куликова И.А. Влияние крупномасштабных мод атмосферной циркуляции на режим температуры и осадков в Арктике // Метеорология и гидрология. 2018. № 1. С. 5–21.

  8. Куликова И.А., Круглова Е.Н, Киктев Д.Б. Крупномасштабные моды атмосферной изменчивости. Ч. II. Их влияние на пространственное распределение температуры и осадков на территории Северной Евразии // Метеорология и гидрология. 2015. № 4. С. 5–16.

  9. Муравьев А.В., Куликова И.А Взаимосвязь суммарных осадков над Евразией с центрами действия атмосферы Северного полушария и главными модами изменчивости температуры поверхности Северной Атлантики // Метеорология и гидрология. 2011. № 5. С. 5–16.

  10. Нестеров Е.С. О формировании взрывных циклонов в северо-восточной части Атлантического океана // Метеорология и гидрология. 2010. № 10. С. 44–53.

  11. Николаев М.В. Оценка смещения границ зон рискованного земледелия в условиях изменения климата // Изв. РГО. 2015. Т. 147. Вып. 1. С. 54–65.

  12. Николаев М.В. Использование метода пространственно-временных аналогов для оценок уязвимости земледелия к воздействию климатических изменений // Изв. РГО. 2015. Т. 147. Вып. 2. С. 1–12.

  13. Николаев М.В. Климатический мониторинг для оценок уязвимости сельскохозяйственных территорий к эффектам переувлажнения в Нечерноземной зоне Европейской России // Изв. РГО. 2017. Т. 149. Вып. 5. С. 4–16.

  14. Николаев М.В. Уязвимость и адаптация полеводства в регионах России к факторам атмосферной засухи и избыточного увлажнения в условиях изменяющегося климата // Материалы Международной научной конференции, посвященной 85-летию Агрофизического НИИ “Тенденции развития агрофизики: от актуальных проблем земледелия и растениеводства к технологиям будущего”, Санкт-Петербург, 27–29 сентября 2017 г. СПб: ФГБНУ АФИ, 2017. С. 513–520.

  15. Николаев М.В. Оценка изменяющегося вклада обильных осадков в рискованность земледелия в Нечерноземье Европейской России // Изв. РГО. 2018. Т. 150. Вып. 6. С. 1–14. https://doi.org/10.1134/S0869607118060010

  16. Николаев М.В. Геостатистическая оценка рискованности земледелия России в условиях изменения климата// Материалы 2-й Международной научной конференции “Тенденции развития агрофизики: от актуальных проблем земледелия и растениеводства к технологиям будущего”, посвященной памяти акад. Е.И. Ермакова, Санкт-Петербург, 2–4 октября 2019 г. СПб: ФГБНУ АФИ, 2019. С. 631–639.

  17. Николаев М.В. Влияние климатических изменений на продуктивность мелиорируемых земель // Материалы Всероссийской научной конференции с международным участием “Вклад агрофизики в решение фундаментальных задач сельскохозяйственной науки”, Санкт-Петербург, 1–2 октября 2020 г. СПб: ФГБНУ АФИ, 2020. С. 478–483.

  18. Пасечнюк А.Д. Погода и полегание зерновых культур. Л.:Гидрометеоиздат, 1990. 212 с.

  19. Природно-сельскохозяйственное районирование земельного фонда СССР. М.: Агропромиздат, 1986. 71 с.

  20. Роде А.А. Водный режим почв и его регулирование. М.: Изд-во АН СССР, 1963. 117 с.

  21. Селянинов Г.Т. Принципы агроклиматического районирования СССР // Вопросы агроклиматического районирования СССР. М.: Изд-во МСХ СССР, 1958. С. 7–13.

  22. Система ведения сельского хозяйства северо-западной зоны РСФСР. Л.: Колос, 1968. 592 с.

  23. Справочник агронома по сельскохозяйственной метеорологии. Нечерноземная зона Европейской части РСФСР / Под ред. И.Г. Грингофа. Л.: Гидрометеоиздат, 1986. 527 с.

  24. Climate change 2014 Synthesis Report / Ed. by The Core Writing Team, R. K. Pachauri, L. Meyer. WMO: Geneva, 2014. 132 p.

  25. Climate Extremes and Their Implications for Impact and Risk Assessment / Ed. by J. Sillman, S. Sippel, S. Russo. Elsevier Inc., 2020. 355 p.

  26. Eitzinger J., Utset A., Trnka M., Zalud Z., Nikolaev M., Uskov I. Weather and climate and optimization of farm technologies at different input levels // Managing Weather and Climate Risks in Agriculture / M.V.K. Sivakumar, R. Motha (eds.). Springer, 2007. P. 141–170.

  27. Eitzinger J., Thaler S., Kubu G., Alexandrov V., Utset A., Mihailovic D.T., Lalic B., Trnka M., Zalud Z., Semiradova D., Ventrella D., Anastasiou D.P., Medany M., Altaher S., Olejnik J., Lesny J., Nemeshko N., Nikolaev M., Simota C., Cojocaru G. Vulnerability and Adaptation Options of European Agriculture // Global Environmental Change: Challenges to Science and Society in Southeastern Europe / V. Alexandrov, M.F. Gajdusek, C.G. Knight, A. Yotova (eds.). Springer, 2010. P. 139–161.

  28. Eitzinger J., Kubu G., Thaler S., Glauninger J., Alexandrov V., Utset A., Mihailovic D.T., Lalic B., Trnka M., Zalud Z., Semiradova D., Ventrella D., Anastasiou D.P., Medany M., Altaher S., Olejnik J., Lesny J., Nemeshko N., Nikolaev M.V., Simota C., Cojocaru G. Adaptation Options to Climate Change Impacts in European Agriculture // Climate Change Adaptation: Ecology, Mitigation and Management / Adam L. Jenkins (ed.). Nova Science Publishers, 2011. P. 151–162.

  29. Franzke C., Feldstein S.B. The continuum and dynamics of Northern Hemisphere teleconnection patterns // J. Atmospheric Science. 2005. V. 62. Iss. 9. P. 3250–3267.

  30. Kuhibrodt T. et al. The Low-Resolution Version of HadGEM3 GC 3.1 model: Development and Evaluation for Global Climate // J. Advances in Modeling Earth Systems. 2018. V. 10. Iss. 11. P. 2865–2888. doi.org/https://doi.org/10.1029/2018MS001370

  31. Massonnet F. et al. Replicability of the EC-Earth3 Earth System Model under a Change in Computing Environment // J. Geoscientific Model Development. 2020. V. 13. P. 1165–1178. https://doi.org/10.5194/gmd-13-1165-2020

  32. Muller W.A. et al. A Higher-resolution Version of the Max Planck Institute Earth System Model (MPI-ESM 1.2-HR) // J. Advances in Modeling Earth Systems. 2018. V. 10. Iss. 7. P. 1383–1413. doi.org/https://doi.org/10.1029/2017MS001217

  33. Nikolaev M.V. Impact of Climate Change on Agriculture in North-West Russia and Adaptation Options // Advances in Environmental Modeling and Measurements / D.T. Mihailovic, B. Lalic (eds.). Nova Science Publishers, 2010. Ch. 20. P. 223–231.

  34. Nikolaev M.V. Integrated Assessment of Change in Contribution of Excessive Moisture to Farming Risks in the Humid Zone of Western Russia // J. Meteorology, Hydrology and Water Management. 2020. V. 8. Iss. 1. P. 46–53. https://doi.org/10.26491/mhwm/111543

  35. Swart N.C. et al. The Canadian Earth System Model version 5 (Can ESM5.03) // J. Geoscientific Model Development. 2019. V. 12. Iss. 11. P. 4823–4873. https://doi.org/10.5194/gmd-12-4823-2019

Дополнительные материалы отсутствуют.