Региональные особенности трансформации структуры водного баланса в ходе восстановительных сукцессий на вырубках темнохвойных лесов Енисейского кряжа

Полный текст

Введение

Подход к лесопользованию с позиций поддержания и восполнения водных ресурсов на реках Средней Сибири ориентируется на изучение последствий рубок на гидрологический режим территории. Удаление лесного полога в результате рубки сказывается на трансформации гидрологического режима, и, прежде всего, на соотношении основных элементов водного баланса – суммарного испарения и стока. Если данные по стоку для некоторых рек доступны благодаря регулярным наблюдениям на гидропостах гидрометеорологической сети, то данные по суммарному испарению представлены фрагментарными натурными исследованиями на лесных водосборах (Антипов, Антипова, 1984; Буренина, 1982; Лебедев, 1982; Молчанов,1960; Раунер, 1972; Федоров, 1979, 1981; Bosch and Hewlett, 1982; Turner, 1991) или основаны на расчетных методах и эмпирических зависимостях (Будаговский, 1964; Будыко, 1956; Карпечко, Мясникова, 2014). В Советском Союзе большое внимание уделялось исследованию данных натурных наблюдений, но при этом многие исследователи отмечали большую трудоемкость прямого измерения суммарного испарения и невозможность экстраполяции полученных данных на большие территории.

Из расчетных методов наиболее часто используемый – метод водного баланса, который применялся при наличии данных сетевых наблюдений за атмосферными осадками и речным стоком. Широкое распространение получили расчетные схемы М.И. Будыко (1956), А.И. Будаговского (1964), Х.Л. Пенмана (1968). Они основаны на теории турбулентной диффузии и теплового баланса подстилающей территории. Одной из наиболее распространенных моделей для изучения суммарного испарения благодаря удобству расчета и легкости получения необходимых метеорологических данных является Формула Пенмана-Монтейта, которую многие исследователи модифицируют с учетом факторов, влияющих на процесс эвапотранспирации (Yang, 2007; Zhang et al., 2016).

Некоторые модели, такие как SEBAL (Surface Energy Balance Algorithm for Land) (Bastiaanssen et al., 1998) и METRIC (Mapping Evapotranspiration with Internalized Calibration) (Allen et al., 2007), определяют энергетический баланс на земной поверхности с использованием спутниковых снимков. Это позволяет рассчитывать как фактическую, так и потенциальную эвапотранспирацию на пиксельной основе. SEBAL и METRIC могут отображать эти ключевые индикаторы во времени и пространстве, по дням, неделям или годам. С появлением спектрометра среднего пространственного разрешения MODIS космическое агентство США NASA для оценки эвапотранспирации использует метод Пенмана–Монтейта (Monteith, 1965), учитывающий климатические характеристики исследуемой территории, типы подстилающей поверхности, а также разнообразие растительного покрова (Расулова, 2021; Mu, 2007, 2011).

Как показали исследования (Лебедев, 1982; Побединский, 1976; Федоров, 1979, 1981; Bates and Henry, 1928; Hibbert, 1967 и др.), в результате масштабных рубок на водосборах резко меняется структура водного баланса, но по мере восстановления лесной растительности на вырубках соотношение эвапотранспирации и стока постоянно меняется в зависимости от направленности лесовосстановительного процесса. На основе комплексного анализа изменения структуры лесов после рубки и динамики элементов водного баланса за несколько десятилетий, выполненного для лесной зоны Европейской территории СССР, О.И. Крестовский (1984) установил, что восстановление водного режима территории до рубки леса и, соответственно, водности рек произойдет не раньше чем через 50 лет.

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

Объекты и методы

Исследования проводились в бассейне р. Сухой Пит (рис. 1), которая является левым притоком р. Большой Пит, принадлежащей к бассейну р. Енисей. Сухой Пит протяженностью 156 км и площадью бассейна 1262 км² формирует свой сток на западном склоне Енисейского кряжа, который, несмотря на его незначительную высоту, является орографическим барьером на пути господствующего переноса влагонесущих воздушных масс и существенно влияет на распределение осадков в регионе. С повышением местности над уровнем моря количество осадков увеличивается. По данным метеорологических станций Усть-Пит и Северо-Енисейск, на наветренных склонах кряжа выпадает от 540 мм до 1002 мм осадков в год. На восточных склонах (метеостанции Южно-Енисейск и Раздолинск) количество осадков уменьшается и составляет 406–490 мм (http://meteo.ru/it/178-aisori). По лесорастительному районированию изучаемая территория относится к Приенисейской провинции, округу темнохвойных и темнохвойно-лиственничных лесов (Кутафьев, 1970). Преобладают темнохвойные леса, представленные пихтовыми насаждениями, спелые и перестойные VI–VIII классов возраста и низкой продуктивности (III–IV классов бонитета). Сомкнутость полога 0.7–0.9. Запас древесины не превышает 200–250 м³/га. Подлесок представлен рябиной, черемухой, красной смородиной, подрост куртинный – из пихты, осины и березы. В живом напочвенном покрове преобладает кислица, майник, звездчатка Бунге. Зеленые мхи образуют сплошной покров (90%). Проективное покрытие травяно-кустарничкового яруса – 40–60%, средняя высота 10 см. Подстилка рыхлая, толщиной до 6 см. Почва сильно оподзоленная, влажная, суглинистая, средней мощности. За последние десятилетия в результате хозяйственного использования лесов лесной покров исследуемого района претерпел значительную антропогенную трансформацию (Лапшина и др., 1971).

 

Рис. 1. Участки, нарушенные вырубками и пожарами на водосборе р. Сухой Пит за период 2001–2020 гг.

 

Территория водосбора р. Сухой Пит была определена при помощи функций гидрологического анализа (ESRI ArcGIS 10 Spatial Analysis) на основе цифровой модели рельефа (SRTM 1 Arc-Second Global). Динамика площадей вырубок для водосбора определялась по данным спутниковой съемки Landsat TM/ETM/OLI из открытых архивов USGS (https://earthexplorer.usgs.gov) за период 2001–2020 гг. При этом для погодичной агрегации вырубок использовался продукт Global Forest Change (https://www.science.org/doi/10.1126/science.1244693). Изображения Landsat, отобранные для территории водосбора, дешифрировались последовательным применением методов сегментации Multiresolution Segmentation (Trimble Cognition Developer 8) и неуправляемой классификации Isoclaster (ESRI ArcGIS 10) с последующим уточнением классов по методу Maximum Likelihood Classification (ESRI ArcGIS 10). Выделенные классы вырубок для каждой спутниковой сцены идентифицировались, датировались и дополнялись на основе данных Global Forest Change.

Для выделения классов вырубок были привлечены данные лесоводственных исследований 2011 г., полученные в результате полевых исследований. Пробные площади, где проводились обследования состояния лесовозобновления на вырубках, находятся на арендных территориях Лесосибирского ЛДК № 1, где Институтом леса ведется мониторинг за процессом лесовосстановления на площадях, нарушенных рубкой леса. Бассейн Сухого Пита вследствие небольшой площади и незначительных перепадов высоты характеризуется относительно однородными лесорастительными условиями и, соответственно, растительным покровом.

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

В работе были использованы данные многолетних наблюдений за стоком и осадками по гидрологическом посту Сухой Пит (http://www.textual.ru/gvr/). Данные по испарению по месяцам были получены по спутниковым данным MODIS – эвапотранспирация с шагом агрегации 8 дней (продукт MOD16A2 – пространственное разрешение 500 м/пикс., https://modis.gsfc.nasa.gov/data/dataprod/mod16.php). Также были использованы данные по индексу листовой поверхности на основе Terra MODIS (MOD15A2H, 500 м/пикс., https://lpdaac.usgs.gov/products/mod15a2hv006/). Значения среднемесячной температуры воздуха (TMP) и суммарных месячных осадков (PRE) были взяты нами из пространственной растровой базы данных CRU TS 4.01 (разрешение 0.5 ~ 50 км/пикс., https://crudata.uea.ac.uk/cru/data/hrg/cru_ts_4.01/).

Нами использовался продукт MODIS Global Terrestrial Evapotranspiration (ET) Product (MOD16A2/A3) рассчитанный на основе спутниковых данных MODIS по алгоритму Year-end Gap-filled (MOD16A2GF/A3GF) NASA Earth Observing System MODIS Land Algorithm (https://lpdaac.usgs.gov/documents/931/MOD16_User_Guide_V61.pdf). Приведенный метод расчета данных ET основан на уравнении Пенмана–Монтейта (Monteith, 1965):

ET = (Δ(R_n – G) + ρ_ac_p(δe)g_a) / ((Δ + γ(1 + g_a/g_s))L_v),

где ET – объем эвапотранспирированной воды (мм с^–1); Δ – скорость изменения удельной влажности насыщения с температурой воздуха (Па К^–1); R_n – чистая освещенность (Вт м^–2), внешний источник потока энергии; G – тепловой поток почвы (Вт·м^–2); c_p – удельная теплоемкость воздуха (Дж кг^–1 К^–1); ρ_a – плотность сухого воздуха (кг м^–3); δe – дефицит давления пара (Па); g_a – проводимость воздуха, проводимость атмосферы (м с^–1); g_s – проводимость стомы, поверхностная или устьичная проводимость (м с^–1); γ – психрометрическая константа (γ ≈ 66 Па К^–1).

Алгоритм расчета параметров уравнения подробно описан в документации продуктов MOD16A2GF/A3GF и включает параметры, определяемые дистанционно (по данным съемки MODIS), такие как индекс листовой поверхности (LAI), альбедо и тип земной поверхности, а также данные о температуре и влажности воздуха. В расчетах мы использовали набор изображений. MOD16A2 – (суммарную) за 8 дней эвапотранспирацию (ET, мм) с пространственным разрешением 500 × 500 м в пикселе.

Средневзвешенные по суммарной площади пикселей тестовых участков значения эвапотранспирации, индекса листовой поверхности, осадков и температура воздуха были рассчитаны для трех участков: объединенных полигонов вырубок за период 2001–2020 гг.; фонового полигона внутри водосбора (без нарушений) и для всего водосбора целиком. При этом 8-дневные значения эвапотранспирации суммировались для каждого года по трем диапазонам дат: 17 мая – 14 сентября (летний период), 8 февраля – 30 марта (зимний период) и для всего года в целом. Интервалы дат были выбраны с учетом наличия данных в архивах продуктов MOD16A2 и MOD15A2H. Все указанные выше растровые данные приводились к единой картографической проекции (WGS 1984 UTM Zone 46N) с разрешением растровой сетки 500 м/пикс. Внутри трех тестовых полигонов (1 – вырубки, 2 – ненарушенный фон, 3 – весь водосбор) рассчитывалось среднее значение величин ET и LAI (8-дневный шаг дат), а также TMP и PRE (месячный шаг дат) из диапазона 2001–2020 гг. Сформированные временные ряды значений ET, LAI, TMP и PRE (для полигонов 1, 2 и 3) использовались далее для определения среднегодовых, сезонных значений (летний и зимний периоды) и сравнительного статистического анализа. Для оценки уровня связи между показателями использовался коэффициент корреляции Пирсона при уровне значимости p < 0.05 (5%).

Обработка и анализ данных выполнялись с использованием процедур открытого сервиса Google Earth Engine и функций зональной статистики пакета ESRI ArcGIS 10.

Результаты и обсуждение

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

 

Таблица 1. Характеристика мониторинговых объектов ИЛ СО РАН, использованных в работе (данные на 2011 г.)

Тип растительных сообществ	Лесоводственно-таксационные характеристики древостоя					Подрост		Травяно-кустарниковый ярус	Координаты
	Состав	Возраст древостоя, лет	Полнота	Высота, м	Диаметр, см	Состав	Густота, шт./га
Вырубки
Трехлетняя вырубка вейниково-разнотравная	–	–	–	–	–	–	–	светолюбивое опушечное разнотравье, луговые и рудеральные виды	58°27´49 93°16´59
Трехлетняя вырубка –кипрейно-вейниковый малинник	–	–	–	–	–	–	–	лугово-лесное разнотравье, злаки и рудеральные виды	58°27´51 93°17´61
Семилетняя гарь кипрейно-малинниковая (по вырубке)	–	–	–	–	–	8Ос2Б, ед. К, П	9750	пирофильные рудеральные виды, луговые и боровые виды, кустарники	58°27´22 93°15´31
Молодняки
Березняк малинниково-кипрейно-вейниковый	5Б5Ос	8	0.3	0.5	1.2	9Б1П, ед. К, Е	35500	кустарники, лесное разнотравье, таежные рудеральные и виды	58°26´52 93°13´10
Пихтарник с березой чернично-разнотравный	5П5Б	10	0.3	0.5	1.5	–	–	кустарники, бореальные кустарнички, таежное мелкотравье, лесное разнотравье	58°26´58 93°13´01
Березняк вейниковый с красной смородиной	10Б+Ос, Е, П	13	0.7	1.8	1.7	4П3Б2Ос1Е	7800	кустарники, лугово-лесное разнотравье, таежные виды	58°26´33 93°10´42
Жердняки
Березняк с осиной разнотравный	8Б2П	23	0.3	2.0	1.6	8П2Б	6300	лугово-лесное разнотравье, таежное мелкотравье, лесные папоротники	58°26´21 93°16´19
Пихтарник мелкотравно- зеленомошный	7П3Б	30	0.8	6.4	5.6	–	–	таежное мелкотравье, зеленые мхи, лесное разнотравье, бореальные кустарнички	58°26´23 93°06´18
Контроль
Пихтарник чернично-зеленомошный	10П	120	0.8	19.5	17.8	–	–	таежное мелкотравье, зеленые мхи, бореальные кустарнички	58°27´26 93°15´28
Пихтарник вейниково-мелкотравный	10П+Ос+Б	100	0.9	18.0	17.0	–	–	лесное разнотравье, таежные виды, зеленые мхи	58°27´54 93°16´57

 

Проведение снегомерных съемок на объектах мониторинга в течение нескольких лет позволили проследить динамику формирования снежного покрова на вырубках в пространственном и временном аспектах, а также рассчитать испарение со снежного покрова на вырубках с учетом стадии лесовосстановления (Burenina et al., 2017). Испарение со снега на вырубках текущего года варьирует от 12 до 34 мм в зависимости от экспозиции склона, по мере зарастания вырубок оно уменьшается, но в пихтовых жердняках отмечается некоторое повышение испарения со снега за счет перехваченных осадков сомкнувшимся пологом древостоя.

Из наших работ и публикаций других авторов (Лебедев, 1982; Fu, 1981; Morton, 1983, 1984) известно, что основной расход влаги на испарение растительным покровом происходит за счет транспирации и испарения осадков, задержанных растительным пологом. С использованием зависимости расхода влаги на транспирацию и испарение задержанных осадков от лесоводственно-биометрических параметров (Грибов, 1997; Лебедев, 1982; Онучин, 2001; Onuchin and Burenina, 2008) были рассчитаны составляющие суммарного испарения на исследуемых вырубках. Для расчетов использовались литературные данные по интенсивности транспирации древесных пород и травяно-кустарникового яруса (Бейдеман, 1983), а также фондовые материалы лаборатории техногенных лесных экосистем Института леса СО РАН по запасам фитомассы травяно-кустарничковой растительности и листвы и хвои древесного яруса.

Располагая расчетными данными по эвапотранспирации, мы проанализировали структуру суммарного испарения в зависимости от сукцессий растительного покрова на вырубках темнохвойных лесов Енисейского кряжа (рис. 2). Снижение суммарного испарения от 3-летних вырубок до 13–14-летних лесных молодняков связано с уменьшением транспирации и незначительным задержанием атмосферных осадков древостоем. Переход от травянисто-кустарниковой стадии на вырубке в стадию молодняков ведет к уменьшению запасов транспирирующей фитомассы (Буренина, 1982; Burenina et al., 2021). Значительные различия суммарного испарения между 23-летним березняком и 25-летним пихтарником – это результат того, что хвойные древостои задерживают не только жидкие, но и твердые атмосферные осадки. Полученные данные характеризуют локальные условия лесорастительных условий и не всегда могут быть использованы для оценки суммарного испарения на бассейновом уровне. Для анализа динамики эвапотранспирации по стадиям лесовосстановления на бассейновом уровне использованы спутниковые данные Terra (MOD16A2).

 

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

 

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

 

Таблица 2. Эвапотранспирация на вырубках и фоновой территории за вегетационный сезон (15.05–15.09)

Год	Эвапотранспирация, мм		Средняя температура, °C	Сумма осадков, мм
	условный фон	вырубки
2001			12.9	302
2002	164.4	163.1	13.6	311
2003	187.2	188.4	13.2	288
2004	191.3	189.9	13.3	216
2005	187.1	185.1	12.5	302
2006	195.2	194.4	13.2	317
2007	186.1	189.0	12.7	221
2008	185.6	189.8	12.3	323
2009	182.2	187.5	12.3	313
2010	180.3	187.6	12.0	321
2011	177.4	185.4	11.5	332
2012	186.3	194.5	12.8	322
2013	180.8	193.2	13.5	228
2014	172.7	177.5	11.9	364
2015	178.6	182.6	11.9	352
2016	196.3	208.9	13.2	359
2017	194.9	197.7	13.2	265
2018	196.6	203.6	12.6	350
2019	195.6	201.3	12.6	256
2020	188.8	203.8	13.4	276

 

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

Для того чтобы понять, насколько общие модели применимы на региональном и локальном уровнях, мы оценили связь между годовой эвапотранспирацией и среднегодовой температурой воздуха, годовой суммой осадков (несмещенная и с опережением эвапотранспирации на 1 и 2 года) и накопленной площадью вырубок на территории водосбора р. Сухой Пит за период 2001–2020 гг. (длина выборки n = 20 лет). В табл. 3 приведены значения коэффициентов корреляции Пирсона (R) для ET и указанных характеристик, отобранных по уровню значимости p < 0.05 (5%).

 

Таблица 3. Оценка уровня корреляции годовой эвапотранспирации (ET) с основными метеорологическими характеристиками и площадью вырубок на водосборе р. Сухой Пит за период 2001–2020 гг.

Характеристика водосбора	Коэффициент корреляции (R)	Уровень значимости (p)
Среднегодовая температура воздуха	0.62	0.003
Сумма годовых осадков (несмещенная)	0.36	0.05
Сумма годовых осадков (с опережением ET на 1 год)	0.82	0.00004
Сумма годовых осадков (с опережением ET на 2 года)	0.44	0.04
Накопленная с 2000 г. площадь вырубок	0.71	0.0002

 

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

Сравнение межгодовой вариабельности эвапотранспирации с климатическими характеристиками по вегетационным сезонам показало, что максимальный уровень эвапотранспирации наблюдается в годы с оптимальным сочетанием температуры воздуха и осадков. В годы с осадками менее 300 мм и температурой воздуха выше 12.5°С наблюдается максимальное испарение как на фоновых участках водосбора р. Сухой Пит, так и на вырубках. При величине осадков более 300 мм даже при температуре воздуха 13.0°–13.9°С наблюдается снижение суммарного испарения. Минимальные значения эвапотранспирации характерны в годы с максимальными значениями осадков (более 350 мм) и температурой воздуха менее 12°С (см. табл. 2).

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

Как следует из материалов диаграммы (рис. 3), эвапотранспирация на вырубках и в целом по бассейну, несмотря на различия в абсолютных величинах, изменяется синхронно, что соответствует вариабельности метеорологических условий по годам. Но при этом до 2005 г. средневзвешенные значения эвапотранспирации на вырубках и в целом по бассейну практически не различаются. С 2009 г. на вырубленных участках эвапотранспирация увеличилась на 10–20 мм по сравнению с ее средневзвешенными значениями в бассейне и на фоновых участках. Эти данные значительно расходятся с результатами расчетов, выполненных для наших мониторинговых объектов. Объяснить это можно как использованием при расчетах средних многолетних показателей, так и тем, что контрольные пробные площади характеризуют коренные древостои региона, тогда как фон и водосбор в целом – это совокупность различных категорий насаждений и других природных комплексов. Эвапотранспирация этих участков – это интегральная оценка (средневзвешенная) для данной территории в отличие от точечных данных, каковыми являются конкретные пробные площади.

 

Рис. 3. Динамика летней эвапотранспирации (ET) для вырубок и водосбора р. Сухой Пит.

 

Согласно данным рис. 3 на временном отрезке 2001–2020 гг. эвапотранспирация на вырубленных участках и на водосборе в целом имеет тенденцию к увеличению. На водосборе суммарное испарение увеличилось на 22 мм, что можно объяснить устойчивым годовым приростом среднегодовой температуры воздуха для этого региона (Корец и др., 2020). Но более интенсивный прирост эвапотранспирации на вырубленных участках (40 мм за 20 лет) при тех же климатических условиях, скорее всего, связан с особенностями лесовосстановления, т.е. с ростом суммарного испарения на сукцессионных стадиях.

В бассейне р. Сухой Пит вырубки проводились неравномерно от 0.4 до 3% от площади водосбора в год. Максимальные площади были пройдены рубками в 2011–2012 и 2018–2020 гг. На рис. 4 показано, как в бассейне р. Сухой Пит менялось соотношение различных возрастных послерубочных категорий. Даже в годы, когда отмечались максимальные рубки на водосборе, преобладали площади, вырубленные пять лет назад, что могло внести определенный вклад в формирование суммарного испарения.

Изменение летней эвапотранспирации по годам с учетом дифференциации вырубок по возрастным категориям показало, что максимальные значения эвапотранспирации приурочены к максимальным площадям вырубок 5-летнего возраста, но в 2018–2020 гг. отмечается снижение суммарного испарения, возможно связанное с увеличением площадей молодняков и жердняков на месте вырубок предыдущих лет. Более высокий расход влаги в летний период на 5-летних вырубках (травяно-кустарниковая стадия) по сравнению с другими производными фитоценозами вполне обоснован их высокой продуктивностью и высокими показателями транспирации произрастающих здесь видов растений (Буренина и др., 2021).

 

Рис. 4. Динамика площадей вырубок различной давности и летней эвапотранспирации (ET) для территории водосбора р. Сухой Пит.

 

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

Констатируя, что вырубка леса и последующие стадии лесовосстановления вносят определенный вклад в суммарное испарение в бассейне р. Сухой Пит, определенный интерес представляет вопрос, в какой степени это отражается на речном стоке.

Анализ динамики стока р. Сухой Пит с учетом изменения площадей свежих вырубок показал, что наиболее значительные пики подъема воды в реке приурочены не к годам проведения наиболее масштабных рубок, а являются результатом, когда 2–3 года подряд вырубаются значительные площади (2004–2006 гг.), т.е. на увеличении стока отражается аккумулирующий эффект нарушенности земель на водосборе (рис. 5). Резкий скачок модуля стока с 2012–2014 гг. связан как с аномальными погодными условиями (в 2010–2011 и 2013 гг. выпало на 50% больше нормы осадков, а в 2014 г. почти в два раза больше среднемноголетней величины), так и с тем, что в 2012–2013 гг. ежегодно вырубалось более 2% от площади водосбора.

 

Рис. 5. Динамика площади вырубок на водосборе и модуля стока р. Сухой Пит.

 

Из вышесказанного следует, что при ежегодно вырубаемой площади лесов менее 2% от площади водосбора лесные рубки вызывают незначительные изменения в величине суммарного испарения и, соответственно, речного стока. Следует отметить, что большее влияние на гидрологический режим промышленных лесозаготовок проявляется в аномальные годы: если годовая сумма осадков превышает среднюю многолетнюю норму на 50–60%, уже при вырубке 2–2.5% от площади водосбора отмечается увеличение годового стока до 30% с высоким уровнем воды в реке в половодье и летние паводки. При этом следует принимать во внимание факт, что при накоплении вырубленных площадей в сочетании с экстремальными осадками на небольших водосборах можно ожидать катастрофических наводнений. В ряде публикаций (Sun et al., 2006; Wang, 2008; Wei et al., 2003) утверждается, что переход лесных земель в результате рубок или пожаров в травяные сообщества приводит к возникновению пиков паводков на локальном уровне, но не являются главной причиной высоких половодий. Из выводов, сделанных Центром Международных лесных исследований (FAO, 2006) на основе обобщения данных многолетних наблюдений, следует, что леса не могут предотвращать катастрофические наводнения, вызванные метеорологическими явлениями. Отдельные гидрологи (Vertessy and Bessard, 1999; Zhang et al., 2001, 2012) придерживаются концепции, что лесистость водосборов является более значимым фактором, чем метеорологические параметры.

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

Заключение

Таким образом, оценка испарения с поверхности суши по данным дистанционного зондирования Земли (продукт MOD16A2) может быть использована при воднобалансовых расчетах в различных по масштабам водосборных бассейнах, особенно при отсутствии базы данных по стоку, а также при экологической экспертизе последствий нарушенности лесов на водосборе в результате воздействия природных и антропогенных факторов. Анализ динамики показателей эвапотранспирации, полученных с использованием спутниковых данных MODIS для бассейна р. Сухой Пит, позволил констатировать, что на бассейновом уровне вариация эвапотранспирации по годам в первую очередь связана с изменчивостью температуры воздуха и осадков в сочетании с динамикой свежих вырубок и уже облесившихся лесосек. Ранжирование площадей вырубок с учетом восстановительно-возрастной динамики растительного покрова позволило оценить вклад различных стадий восстановительных сукцессий в суммарное испарение в бассейне р. Сухой Пит и показать, как соотношение площадей свежих вырубок с площадями уже облесившихся лесосек отражается на величине суммарного испарения. Сравнение трендов эвапотранспирации на фоновых участках бассейна и площадях, нарушенных вырубками, показывает, что на изменение испарения на водосборе влияют не только климатические характеристики, но и устойчивый тренд увеличения площадей молодняков и жердняков на месте вырубок конца XX и начала XXI в. Результаты исследований показали, что при поддержании оптимальной лесистости и выполнении технологии рубки с соблюдением сроков примыкания и площадью ежегодной рубки не более 2% от площади водосбора возможно поддержание стабильности гидрологического режима на малых и средних реках Сибири с учетом климатических трендов, зафиксированных на данный момент.

Финансирование

Работа выполнена при поддержке базового проекта ИЛ СО РАН “Теоретические основы сохранения экологического и ресурсного потенциала лесов Сибири в условиях возрастающего антропогенного пресса и климатических аномалий”, № АААА-А17-117101940014-9 (0356-2019-0027) и гранта РФФИ “Исследование особенностей гидрологического цикла в речных бассейнах таежной зоны Средней Сибири при усилении антропогенного пресса на лесные экосистемы в условиях меняющегося климата” № 20-05-00095.

Funding

This work was supported by the basic project of the Institute of forest Siberian Branch of RAS “Theoretical Foundations of Preserving the Ecological and Resource Potential of Siberian Forests under the Conditions of Increasing Anthropogenic Press and Climate Anomalies,” no. AAAA-A17-117101940014-9 (0356-2019-0027) and RFBR grant “Study of the features of the hydrological cycle in the river basins of the taiga zone of Central Siberia with increased anthropogenic pressure on forest ecosystems in a changing climate” no. 20-05-00095.

Об авторах

Т. А. Буренина

Институт леса им. В.Н. Сукачева Сибирского отделения Российской академии наук – обособленное подразделение ФИЦ КНЦ СО РАН (ИЛ СО РАН)

Автор, ответственный за переписку.
Email: burenina@ksc.krasn.ru
Россия, Красноярск

М. А. Корец

Институт леса им. В.Н. Сукачева Сибирского отделения Российской академии наук – обособленное подразделение ФИЦ КНЦ СО РАН (ИЛ СО РАН)

Email: mik@ksc.krasn.ru
Россия, Красноярск

Ж. Р. Сулейманова

Институт леса им. В.Н. Сукачева Сибирского отделения Российской академии наук – обособленное подразделение ФИЦ КНЦ СО РАН (ИЛ СО РАН)

Email: janetta_syleiman@mail.ru
Россия, Красноярск

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

Дополнительные файлы