Усовершенствованная схема параметризации гидрологического цикла суши для моделей Земной системы и систем прогноза погодыНИР

An advanced scheme for parameterization of land surface hydrology in Earth System models and numerical weather forecast systems

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

грант РФФИ

Этапы НИР

# Сроки Название
1 1 января 2020 г.-31 декабря 2020 г. Усовершенствованная схема параметризации гидрологического цикла суши для моделей Земной системы и систем прогноза погоды
Результаты этапа: На территории Метеорологической обсерватории (МО) МГУ продолжен автоматизированный мониторинг термовлажностного режима деятельного слоя почвы, начатый 1 мая 2019 г. На лабораторных установках факультета почвоведения МГУ получены данные по агрегатному и гранулометрическому составу почвы МО МГУ, плотности почвы, плотности твёрдой фазы, пористости, содержанию органического материала, кислотности, проведена аппроксимация гидравлической проводимости и потенциала влажности почвы по параметризации Муалема - ван Генухтена. Полученные физические параметры и переменные состояния почвы используются как входные данные и данные для валидации модели почвы ИВМ РАН-МГУ. Создана новая версия модели тепловлагопереноса в почве ИВМ РАН-МГУ. Параметризации теплопроводности Йохансена, Коте и Конрада, и фрактальная параметризация внедрены в программный код модели. Для проверки чувствительности модели к выбору параметризации привлечены данные измерений в семи различных типах почвы в Тункинской котловине (республика Бурятия). Показано, что параметризация Пильке, используемая в базовой версии модели ИВМ РАН-МГУ, уступает в точности другим параметризациям теплопереноса. В частности, выбор параметризации Коте и Конрада позволяет получить наиболее точные результаты для суглинка и песков. На глобальной сетке модели Земной системы ИВМ РАН-МГУ предлагается переход к параметризациям Йохансена или Коте и Конрада. Параметризация переноса жидкой влаги в почве Муалема-ван Генухтена включена в модель деятельного слоя суши ИВМ РАН-МГУ в качестве альтернативы традиционной параметризации Клаппа-Хорнбергера. Проведены численные эксперименты, в которых в качестве метеорологических данных использованы измерения МО МГУ. Наибольшее различие в результатах расчётов при использовании указанных параметризаций отмечается в слое суточных колебаний температуры. Существенно, что изменение режима влажности существенно сказывается на температурном режиме почвы. Усовершенствовано представление снежного покрова в модели деятельного слоя суши ИВМ РАН-МГУ. Внедрены параметризации задержания снежной массы кронами деревьев и баланса массы снега внутри лесного полога (включая испарение и таяние с поверхности лесного снежного покрова). Валидационные численные эксперименты с обновлённой моделью деятельного слоя ИВМ РАН - МГУ проведены на разрешении 0.5º×0.5º для области, охватывающей речной бассейн р. Северной Двины. В качестве граничных условий использовался реанализ ERA-Interim. Модель удовлетворительно воспроизводит сезонный ход снегозапасов; при этом, имеет место занижение моделью водного эквивалентна снега в среднем по ячейке на ~2 см на протяжении всего холодного периода (исключая период снеготаяния). Межгодовая изменчивость максимальных (мартовских) снегозапасов также удачно воспроизводится моделью; однако, согласие с данными наблюдений при этом лучше на севере, чем на юге бассейна Северной Двины. Модель воспроизводит наблюдаемую в данном речном бассейне закономерность: водный эквивалент снега в лесу выше по сравнению с открытыми участками. Для обеспечения модели ИВМ РАН-МГУ данными о распределении типов подстилающей поверхности высокого качества, проведён анализ существующих баз данных. В качестве оптимального решения выбрана и подключена к модели карта GLCCv2. Карта GLCCv2 представляется не только более точной альтернативой используемой в базовой версии модели карте Вилсона и Хендерсон-Селлерс, но и единственной из широко известных карт, представляющей данные согласно применяемой в модели классификации типов поверхности SiB. Проведён детальный обзор литературы в части методов расчёта параметра термической шероховатости – ключевого параметра, описывающего физические механизмы распространения водяного пара в приповерхностном слое воздуха.
2 1 января 2021 г.-31 декабря 2021 г. Усовершенствованная схема параметризации гидрологического цикла суши для моделей Земной системы и систем прогноза погоды
Результаты этапа: Гидрологический цикл суши играет важнейшую роль в формировании погоды и климата. Его представление в моделях Земной системы и системах прогноза погоды требует непрерывного совершенствования в процессе их развития, связанного с постоянным ростом доступных вычислительных ресурсов. Настоящий проект направлен на создание новой версии модели деятельного слоя суши ИВМ РАН-МГУ, соответствующей уровню ведущих моделей в части описания физических механизмов формирования водного баланса речных бассейнов в умеренных и высоких широтах. Основные результаты, полученные в 2021 г., заключаются в следующем: - продолжены экспериментальные исследования физических параметров и мониторинг температуры и влажности в почве на территории Метеорологической обсерватории (МО) МГУ; экспериментально получено вертикальное распределение основной гидрофизической характеристики, теплофизические характеристики почвы; - проведен расчёт термовлажностного режима почвы моделями почвы SWAP и Hydrus-1D; существенная ошибка моделей в воспроизведении влажности почвы связана со значительной неоднородностью почвы по вертикали; - в почвенный блок модели суши ИВМ РАН-МГУ включен следующие процессы: эффект криовсасывания жидкой влаги, задержка просачивания воды при наличии льда, перенос жидкой влаги за счёт градиента температуры, расчёт диффузии водяного пара в поровом пространстве согласно обобщённой степенной модели; показано, что наибольшее влияние на термовлажностный режим модельной почвы оказывает выбор ОГХ и коэффициента гидравлической проводимости, в то время как неизотермическим переносом влаги можно пренебречь; - в модель деятельного слоя суши ИВМ РАН-МГУ внесены усовершенствования для уточнения расчёта эвапотранспирации: для распределения типов поверхности внедрена карта Барталева (2016), включён сезонный пространственно распределённый годовой ход листового индекса по данным MODIS, в параметризацию устьичного сопротивления добавлена связь с интенсивностью фотосинтеза; в результате существенно улучшена годовая норма испарения в модели относительно наблюдений для трёх выбранных небольших водосборах Европейской территории России; - на основе одномерной RANS-модели пограничного слоя атмосферы проведено исследование на чувствительность потоков тепла, импульса в приземном слое атмосферы к параметрам метелевого переноса; - проведён обзор существующих методов построения поля направлений речного стока на сетках грубого разрешения по данным высокого разрешения для моделей деятельного слоя суши.
3 1 января 2022 г.-31 декабря 2022 г. Усовершенствованная схема параметризации гидрологического цикла суши для моделей Земной системы и систем прогноза погоды
Результаты этапа: Гидрологический цикл суши играет важнейшую роль в формировании погоды и климата. Его представление в моделях Земной системы и системах прогноза погоды требует непрерывного совершенствования, связанного с постоянным ростом доступных вычислительных ресурсов. Настоящий проект направлен на создание новой версии модели деятельного слоя суши ИВМ РАН-МГУ, соответствующей уровню ведущих моделей в части описания физических механизмов формирования водного баланса речных бассейнов в умеренных и высоких широтах. Модель в настоящее время используется для описания взаимодействия атмосферы с подстилающей поверхностью в рамках модели Земной системы ИВМ РАН и системы прогноза погоды ПЛАВ. В ходе проекта в модели были внедрены следующие блоки параметризаций: – усовершенствованные параметризации вертикального теплопереноса (4 параметризации) и влагопереноса в почве, включающие эффекты криовсасывания, задержки гравитационного переноса воды при наличии льда, переноса жидкой влаги за счёт градиента температуры, диффузии водяного пара в поровом пространстве согласно обобщённой степенной зависимости коэффициента диффузии от полной пористости и объёмного содержания газа; – блок описания процессов снегонакопления в высокой растительности, учитывающий разделение твёрдых осадков между кронами и пологом, специфику радиационного режима лесного полога; – новая схема подготовки данных о маршрутизации речных потоков на пространственных разрешениях, характерных для моделей погоды и климата (с шагом сетки >0.1 град.). Кроме того, в модели усовершенствованы схемы расчёта устьичного и аэродинамического сопротивления, применяемые для параметризации потоков явного и скрытого тепла. Уточнены глобальные карты по входным параметрам модели: листовому индексу (в новой версии модели используются данные MODIS) и площади типов растительности (по карте ИКИ РАН). Проверка и калибровка модели осуществлена как с привлечением данных специально организованных in-situ измерений термического и влажностного режима почвы в Метеорологической обсерватории МГУ, так и по осреднённым на различных пространственных масштабах составляющим гидрологического цикла. В качестве последних использован сток равнинных рек России (Северная Двина, Печора, Онега, Мезень), наблюдаемый на гидрометрических постах, а также спутниковые данные о тенденции влагосодержания деятельного слоя проекта GRACE. Продемонстрировано качественное улучшение метрик ошибки модели по составляющим водного баланса бассейнов этих рек и аномалиям влагосодержания деятельного слоя. Подготовлена конфигурация модели для массовых расчётов по протоколу международного проекта сравнения моделей речного стока ISIMIP-Water (https://www.isimip.org).

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

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