Пространственно-временная изменчивость углеродного баланса наземных экосистем ЕТР по данным измерений и математического моделирования (географический факультет)НИР

Spatial and temporal variability of carbon balance of terrestrial ecosystems of the European part of Russia based on field measurements and mathematical modeling (faculty of geography)

Соисполнители НИР

МГУ имени M.В.Ломоносова, Московский университет или МГУ Соисполнитель

Источник финансирования НИР

Правительство Российской Федерации, Программа развития Московского университета

Этапы НИР

# Сроки Название
1 1 августа 2023 г.-31 декабря 2023 г. Пространственно-временная изменчивость углеродного баланса наземных экосистем ЕТР по данным измерений и математического моделирования (географический факультет)
Результаты этапа: Измерения потоков СО2 с поверхности почвы на территории полигона проводились, начиная с февраля 2023 года на шести мониторинговых площадках в разных типах ландшафтов (хвойные и хвойно-широколиственные леса, агроландшафты с многолетними травами и черным паром, луга), с использованием портативного газоанализатора PS-9000 (LICA, Китай). Измерительная система работает по закрытой схеме, при которой скорость почвенной эмиссии определяется по скорости изменения содержания СО2 в закрытой почвенной камере. Определение скорости эмиссии с поверхности почвы на каждой выбранной площадке проводилось в 8-кратной повторности 2 раза в месяц, а также в ходе интенсивных полевых компаний. Параллельно с измерениями потоков определялись температура и влажность почвы. Результаты измерений эмиссии СО2 с поверхности почвы показали значительную временную и пространственную изменчивость потоков. Минимум почвенной эмиссии СО2 наблюдался в зимний период (февраль) около 1.24+-1.52 мкмоль СО2/м2с, и увеличивался в летние месяцы (август) в среднем до 7.69+-2.48 мкмоль СО2/м2с. Выделение СО2 из почвы в атмосферу в осенние месяцы было в среднем в 3 раза ниже, а зимой  почти в 10 раз ниже, чем в летние месяцы. Поступление СО2 из почвы в атмосферу статистически значимо зависело от изменения температуры и объемной влажности почвы. Для масштабирования точечных измерений потоков парниковых газов методами турбулентных пульсаций и экспозиционных камер используется трехмерная математическая модель переноса парниковых газов между неоднородной подстилающей поверхностью и атмосферой. Модель основана на поиске стационарного решения системы уравнений Навье-Стокса с использованием осреднения Рейнольдса и в предположении о постоянстве плотности воздуха. Для определения входных параметров модели, описывающих структуру растительности и рельеф, использовались данные дистанционного зондирования, включающие спутниковые снимки (Landsat 8), и данные лидарной и мультиспектральной съемки с БПЛА. Для задания параметров ветра на верхней границе модельной области использовались данные реанализа ERA5. Результаты модельных экспериментов для участка карбонового полигона Чашниково (Рисунок 2) показали значительную пространственную неоднородность вертикальных потоков CO2 в приземном слое атмосферы для обоих участков, определяемых расчлененностью рельефа, мозаичной растительностью и хаотично расположенными поселковыми зданиями и сооружениями. Модельные расчеты показали, что в дневное время при малооблачной погоде в летние месяцы средняя скорость поглощения СО2 для всей территории полигона Чашниково составляла 4.4 мкмоль/м2с на высоте 5 м (внутри полога древесной растительности и над травянистой растительностью), тогда как на высоте 50 метров средняя величина поглощения СО2 экосистемами полигона за счет вклада древесного яруса составляла уже около 11 мкмоль/м2с. Разработанная модель может являться связующим звеном между измерениями потоков в локальном и экосистемном масштабе и потоками, моделируемыми моделью деятельного слоя суши (ДСС) ИВМ РАН-МГУ. Кроме того, модель служит инструментом для оценки чувствительности потоков СО2 и других парниковых газов к изменению условий внешней среды и точности задания ключевых входных биофизических и биохимических параметров, описывающих свойства растительности и почвы в региональных и глобальных моделях. На основании разработанной модели был предложен новый подход к оценке потоков парниковых газов над горизонтально-неоднородной поверхностью, основанный на измерении концентрации газов на нескольких уровнях над поверхностью (например, измеренных с помощью портативных газоанализаторов, закрепленных на БПЛА) и решении обратной задачи по восстановлению соответствующего потока вблизи верхней границы растительности. Для концентрации исследуемого парникового газа формулируется уравнение "диффузии-адвекции", в котором используются предварительно вычисленные с помощью трехмерной модели поля скорости ветра и коэффициента турбулентности над неоднородной поверхностью. Поток в задаче для концентрации выступает в роли входного параметра в граничном условии. Уровень, на котором он задается, выбирается вблизи верхней границы растительного покрова для того, чтобы избежать детального описания источников и стоков парниковых газов внутри растительности. После вычисления пространственного распределения концентрации газа при заданном потоке на входе осуществляется сравнение модельных данных и результатов измерения концентрации парникового газа на нескольких уровнях, расположенных над растительностью в пределах атмосферного пограничного слоя. Затем производится направленное изменение входного потока в узлах сетки таким образом, чтобы минимизировать норму разности измеренной и вычисленной концентрации. Для оценки начального приближения потока определяющим фактором являются значения измеренных значений концентрации парниковых газов, по крайней мере, на двух уровнях в пределах атмосферного пограничного слоя. Они используются для аппроксимации производных по высоте в выражениях для адвективного и турбулентного вертикального потока исследуемого парникового газа. Проведенное исследование обратной задачи в модельном случае позволило получить предварительные оценки точности восстановления потоков в зависимости от выбора высот для измерения содержания парниковых газов, а также от точности измерения их концентраций.
2 1 января 2024 г.-31 декабря 2024 г. Пространственно-временная изменчивость углеродного баланса наземных экосистем ЕТР по данным измерений и математического моделирования (географический факультет)
Результаты этапа: Для оценки пространственной и временной изменчивости потоков диоксида углерода с поверхности почвы, а также параметризации их зависимости от факторов внешней среды в четвертом квартале 2024 года были продолжены регулярные наблюдения за потоками диоксида углерода с поверхности почвы камерным методом. Измерения проводились на 6 заранее выбранных площадках в разных типах ландшафтов на карбоновом полигоне в Чашниково (Московская область): Во вторичном елово-осиновом щитовниково-медунично-звездчатковом лесу. Почва среднедерново-неглубокоподзолистая поверхностно-слабоглееватая языковатая среднесуглинистая (Epygleyic Glossic Retisol (Loamic)); На многолетних злаковых травах 2-го года использования. Почва дерново-мелкоподзолистая освоенная языковатая мелкопахотная среднесуглинистая слабосмытая (Glossic Antric Retisol (Loamic)); Во вторичном папоротниково-снытиево-кисличном ельнике. Почва глубокодерново-глубокоподзолистая легкосуглинистая (Retisol (Loamic)); На суходольном тимофеевко-таволговом разнотравном лугу. Почва глубокодерново-глубокоподзолистая поверхностно-оглеенная среднесуглинистая (Epygleyic Retisol (Loamic)); В переувлажненном ельнике. Почва торфяная болотно-низинная типичная торфяно-глеевая (Eutric Histosol); На черном паре. Почва дерново-неглубокоподзолистая освоенная языковатая среднепахотная легкосуглинистая (Glossic Antric Retisol (Loamic)). Измерения проводились, начиная с февраля 2023 года, с периодичностью каждые 2-3 недели. Для измерений использовался портативный газоанализатор PS-9000 (LICA, Китай), соединенный с автоматической камерой. На каждой площадке измерения проводились в 8-ми кратной повторности. Анализ временной изменчивости почвенного дыхания на разных площадках показал значительную временную и пространственную изменчивость потоков. Минимум почвенной эмиссии СО2 на всех площадках наблюдался в зимний период 2023 и 2024 года (декабрь-февраль), и увеличивался в летние месяцы. Максимальные значения почвенной эмиссии за исключением площадки в переувлажненном ельнике отмечались в летние месяцы 2023 достигая на суходольном разнотравном лугу в июне до 11.22+-3.19 мкмоль СО2/м2с. Минимальные скорости эмиссии СО2 из почвы в атмосферу в зимние месяцы отмечались в 2024 году, варьируя от 0.32+-3.19 мкмоль СО2/м2с в январе в переувлажненном ельнике до 0.92+-0.27 мкмоль СО2/м2с на суходольном разнотравном лугу. Сравнение скорость эмиссии СО2 с поверхности почвы на травянистых участках (точки 2 и 4) с участком без растительности (пар, точка 6) позволяет оценить средний вклад корневого дыхания в интегральный поток СО2 с поверхности почвы. Результаты сравнения показывают, что, что скорость эмиссии СО2 с поверхности черного пара в летние месяцы 2023 года была в 1.5 – 2 раза ниже скорости эмиссии полевых и луговых участков, и не превышала в июне 2023 года - 4.6 мкмоль/м2с. В летние месяцы 2024 года на участке с паром отмечалась эмиссия СО2 на 20 - 40% ниже, чем на суходольном разнотравном лугу. Скорость почвенного дыхания на участке с многолетними злаковыми травами отличалась значительной изменчивостью в летние месяцы 2024 года: если в июне скорость почвенной эмиссии СО2 была в среднем на 80% ниже почвенного дыхания с поверхности пара, то в августе их соотношение изменилось на противоположное (скорость эмиссии СО2 с поверхности почвы превышала эмиссию СО2 с поверхности пара более чем в 1.5 раза). Данная изменчивость может быть обусловлена различиями температурного и влажностного режима выбранных участков, а также разной чувствительностью автотрофного и гетеротрофного почвенного дыхания к изменению факторов внешней среды. Для параметризации зависимости скорости почвенного дыхания (R_s) от внешних факторов нами использовалась модифицированный алгоритм Hashimoto et al. (2015), учитывающий зависимость скорости почвенного дыхания от влажности почвы (Balogh et al., 2011). Модель была параметризована с использованием данных измерений на черном паре (точка 6) . Для аппроксимации модели использовался алгоритм Левенберга — Марквардта — метод оптимизации, направленный на решение задач о наименьших квадратах. Результаты аппроксимации показали достаточно хорошее соответствие выбранной модели и экспериментальных данных (R2=0.65, p<0.05). Данный алгоритм для параметризации зависимости скорости почвенного дыхания от внешних параметров планируется интегрировать в модель деятельного слоя суши (ДСС) ИВМ РАН-МГУ. Трехмерная локальная математическая модель переноса парниковых газов между неоднородной подстилающей поверхностью и атмосферой была разработана для описания переноса парниковых газов над неоднородной подстилающей поверхностью. Последняя версия модели была также адаптирована для масштабирования точечных измерений потоков парниковых газов методами турбулентных пульсаций и экспозиционных камер, и определения региональных значений потоков. Модель состоит из взаимосвязанных блоков, описывающих атмосферный перенос и турбулентность, перенос прямой и рассеянной солнечной радиации в растительном покрове, перенос парниковых газов (СО2, СН4) в атмосферном пограничном слое. Модель рассматривается как связующее звено между измерениями потоков в локальном и экосистемном масштабе и потоками, моделируемыми моделью деятельного слоя суши (ДСС) ИВМ РАН-МГУ. Она, с одной стороны, позволит описать пространственную (подсеточную) изменчивость потоков парниковых газов для конкретных атмосферных условий, а также оценить эффект неоднородной структуры растительности и рельефа на вертикальные потоки. В ходе реализации проекта была выполнено: • Развитие и оптимизация численной схемы трехмерной гидродинамической модели для описания поля ветра и потоков парниковых газов в атмосферном пограничном слое. • Калибровку входных параметров трехмерной модели, включая структуры растительного покрова, и данные по потокам парниковых газов у поверхности почвы. • Проведение численных экспериментов с помощью трехмерной математической модели по расчету полей вера и вертикальных потоков парниковых газов на карбоновом полигоне Чашниково для разных направлений ветра и стратификации атмосферы.
3 1 января 2025 г.-31 декабря 2025 г. Пространственно-временная изменчивость углеродного баланса наземных экосистем ЕТР по данным измерений и математического моделирования (географический факультет)
Результаты этапа:

Прикрепленные к НИР результаты

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