Численное исследование влияния уменьшения удельной влажности в вертикальном столбе атмосферы на параметры развивающегося конвективного облака
Численное исследование влияния уменьшения удельной влажности в вертикальном столбе атмосферы на параметры развивающегося конвективного облака
Аннотация
В работе представлены результаты численного моделирования и количественной оценки чувствительности динамических, микрофизических и макроструктурных параметров изолированного внутримассового конвективного облака к вариациям начального влагосодержания атмосферы. На основе трехмерной негидростатической модели конвективного облака с детальным учетом термодинамических, микрофизических и электрических процессов проведена серия из шести численных экспериментов, в которых начальный профиль удельной влажности во всем вертикальном столбе атмосферы равномерно уменьшался на определенный процент: от 0% (контрольный эксперимент) до 1%, 2%, 3%, 5% и 10%. Количественный анализ временных рядов (продолжительностью 90 минут) показал наличие выраженной нелинейности и порогового характера деградации облачной системы. Если слабые возмущения влагосодержания до 2% компенсируются внутренней термодинамической перестройкой облака, то при достижении критического порога около 3% запускается механизм каскадного ослабления конвекции. При дефиците влажности в 10% максимальная восходящая скорость падает на 51,2%, а масса наиболее чувствительной ледяной фракции (града и крупы) сокращается на 78%. Снижение интенсивности конвекции обусловлено падением скрытой теплоты конденсации на 21,6%, что подавляет развитие восходящего потока и снижает плавучесть, а также активирует негативное влияние вовлечения сухого воздуха. Полученные зависимости важны для оптимизации параметризации глубокой конвекции в климатических моделях и прогнозирования опасных погодных явлений.
1. Введение
Конвективные облака являются одним из наиболее значимых объектов исследования в динамической метеорологии и физике атмосферы. Они оказывают определяющее влияние на глобальный и региональный водный цикл, перенос тепла, влаги и импульса, а также на формирование экстремальных погодных явлений (шквалистые ветра, интенсивные ливни, крупный град)
, . Среди факторов окружающей среды, определяющих интенсивность и глубину конвекции, особую роль играет влагосодержание тропосферы.Влияние распределения водяного пара в тропосфере на конвективную активность подробно исследовалось теоретически и в наблюдениях
, . Ключевым механизмом этого влияния является втягивание сухого окружающего воздуха в конвективный восходящий поток, что снижает плавучесть поднимающегося воздуха и ослабляет восходящие движения , , для которого предлагается термин — энтрейнмент, от англ. entrainment (вовлечение, увлечение; размыв). Чем суше свободная тропосфера, тем интенсивнее испарение облачных капель при их смешении с окружающим воздухом и тем быстрее теряется плавучесть восходящего термика. Вообще вариации параметров конвективного облака тесно связаны с процессами градообразования , , например, современные нейросетевые модели позволяют оценивать максимальный размер града именно по термодинамическим параметрам, среди которых ключевую роль играет удельная влажность .Проведённые ранее численные исследования, как правило, рассматривали влияние влажности свободной тропосферы в контексте идеализированных экспериментов с постоянными значениями относительной влажности
либо анализировали зависимость от влажности нижнего слоя атмосферы и слоя смешения . Анализ источников показывает, что количественная оценка чувствительности конвективных облаков к изменениям вертикального профиля удельной влажности имеет важное прикладное значение. Во-первых, иссушение атмосферы в отдельных регионах, связанное с изменением климата или адвекцией сухих воздушных масс, существенно модифицирует характеристики конвекции . Во-вторых, понимание этой зависимости критически важно для развития технологий активного воздействия на облака и осадки. В-третьих, корректная параметризация чувствительности конвекции к влажности остаётся одной из основных задач при верификации моделей общей циркуляции атмосферы .В контексте климатических изменений и антропогенного воздействия количественная оценка чувствительности развивающихся облачных систем к вариациям удельной влажности приобретает стратегическое значение для понимания темпов изменений. Например, это необходимо как для повышения точности гидрологических прогнозов, так и для верификации конвективных схем в моделях прогноза погоды и климата
. А значит актуальным становится изучение отклика облачных параметров на однородное пропорциональное иссушение всего вертикального столба атмосферы, которое моделирует масштабные адвективные процессы или региональные климатические сдвиги.Таким образом, целью настоящей работы является комплексное численное исследование и количественный анализ чувствительности динамических, микрофизических и макроструктурных параметров развивающегося конвективного облака к равномерному уменьшению удельной влажности в вертикальном столбе атмосферы в диапазоне от 1% до 10% относительно естественного хода.
2. Методы и принципы исследования
2.1. Описание численной модели
Для проведения численных экспериментов использовалась трехмерная негидростатическая модель конвективного облака с детальным учетом термодинамических, микрофизических и электрических процессов, базирующаяся на уравнениях гидротермодинамики в приближении Буссинеска
. Модель воспроизводит пространственно-временную эволюцию векторов скорости, возмущений давления, потенциальной температуры и массовых долей влаги.Уравнение движения записывается в векторной форме:
где — вектор скорости, — базовая плотность воздуха, p' — возмущение давления, FD— сила турбулентной вязкости и аэродинамического сопротивления гидрометеоров, k — единичный вектор в вертикальном направлении, B — сила плавучести, определяемая выражением:
Здесь θ′ — возмущение потенциальной температуры, θ0 — базовая потенциальная температура стратифицированной атмосферы, а q′v,qc,qr,qi,qs,qg,qh — удельные массовые содержания водяного пара, облачной воды, дождевых капель, первичных ледяных кристаллов, снега, крупы и града соответственно.
Микрофизический блок модели реализует детальное описание фазовых переходов капель воды и ледяных кристаллов в процессах нуклеации, диффузионного роста частиц (включая эффект Бержерона — Финдейзена), фазовых переходов (замерзание, плавление, сублимацию и десублимацию) и контактных взаимодействий частиц, таких как коагуляция, аккреция капель и ледяных кристаллов и агрегация кристаллов
. Микроструктура облака описывается с помощью кинетических уравнений для функций распределения частиц по массам (или радиусам). Модель детально учитывает спектры размеров основных влагосодержащих компонентов, разделяя их на следующие категории: жидкая фаза — облачные капли (qc) и частицы осадков - дождь (qr); твердая (ледяная) фаза — ледяные кристаллы (qi), снег (qs), крупа (qg) и град (qh).В общем математическом виде для каждого типа облачных частиц i уравнение микрофизического блока записывается как:
где: fi(x,y,z,m,t) — функция распределения частиц i-го типа по массам m;
Динамика изменений данных функций рассчитывается с учетом их пространственно-временного переноса воздушными потоками, турбулентного перемешивания, а также скорости седиментации (гравитационного осаждения) частиц каждой фракции под действием силы тяжести. Электрические силы значительно изменяют коэффициенты столкновения и слияния частиц (особенно малых фракций), что существенно ускоряет укрупнение гидрометеоров и стимулирует раннее градообразование на стадии активного роста конвективного облака.
Расчетный домен имеет геометрию 60 Х 60 Х 18 км. Пространственная дискретизация выполнена на регулярной сетке с шагом по вертикали Δz = 200 м и по горизонтали Δx и Δy = 500 м. Шаг интегрирования по времени принят равным Δt = 5 с. На боковых границах заданы условия непротекания, на верхней границе реализован демпфирующий поглощающий слой для предотвращения отражения акустико-гравитационных волн.
2.2. Начальные условия и параметры экспериментов
Начальное термодинамическое состояние атмосферы задавалось на основе реального профиля радиозондирования из архива ФГБУ «ВГИ» (архивный файл input_sounding_2025_08_13_1200Z.txt), характеризующегося высоким запасом конвективной нестабильности. Поверхностная потенциальная температура составляла T0 = 300 К, а базовая удельная влажность у земли достигала q0 = 12.52 г/кг, почти экспоненциально убывая с высотой. Конвекция инициировалась искусственным введением термического возмущения (теплового пузыря) в нижнем 2-километровом слое с радиальным размером 10 км и температурным перегревом Δθ = 3 К в геометрическом центре области моделирования.
Набор численных экспериментов включает один контрольный и пять модифицированных запусков, в которых профиль удельной влажности qv(z) масштабировался коэффициентом (1 - δ), где δ — относительный дефицит влаги:
z01 Исходный базовый профиль, δ = 0%.
z02-z06 Равномерное иссушение на δ = 1%, 2%, 3%, 5%, 10%, соответственно.
Начальное температурное поле, профиль давления и вектор фонового сдвига ветра оставались строго одинаковыми во всех итерациях численного моделирования, что позволило выделить вклад интегрального влагосодержания. Полная продолжительность симуляций составила 90 минут, фиксация диагностических полей производилась с дискретностью 1 минута.
3. Основные результаты
Анализ динамических полей продемонстрировал жесткую, но нелинейную зависимость интенсивности вертикальных движений от дефицита влаги. В контрольном эксперименте z01 развитие кучево-дождевой системы характеризуется классической двухпиковой структурой временного хода максимальной вертикальной скорости (wmax). Первый интенсивный восходящий импульс фиксируется на 38-й минуте, достигая 45.6 м/с. Этот импульс обусловлен взрывным ростом плавучести при начальной конденсации водяного пара приземного влажного объема воздуха. Вторичный пик формируется на этапе зрелости облачной системы (56-я минута), генерируя максимальные значения восходящих потоков, сопряженные с тепловыделением при кристаллизации и аккреции.
При переходе к экспериментам z02 и z03 (иссушение на 1% и 2%) пиковые значения wmax снижаются слабо — до 45.2 м/с и 43.8 м/с соответственно, а время их наступления совпадает с контролем (38-я минута). Однако дальнейшее уменьшение влажности приводит к резкому перелому тренда. В эксперименте z04 (δ = 3%) пиковое значение скорости падает до 41.4 м/с (относительное снижение на 9,3%), а сам динамический максимум смещается во времени на 41-ю минуту, что указывает на замедление темпов облачной эволюции.
В экстремальных сценариях z05 (5%) и z06 (10%) наблюдается глубокая деградация динамического потенциала. Максимум восходящей скорости в эксперименте z06 падает более чем в два раза — до 22.29 м/с, а временной сдвиг пика достигает 24 минут.
Интенсивность нисходящих потоков, оцениваемая в нашем случае по модулю минимальной вертикальной скорости (|wmin|), демонстрирует аналогичную тенденцию к затуханию. Пиковое значение нисходящего потока уменьшается с -14.31 м/с (z01) до -8.74 м/с (z06). Ослабление нисходящих движений физически предопределено снижением массы выпадающих осадков и замедлением процессов испарительного охлаждения в подоблачном слое.
Максимальная радиолокационная отражаемость (Zmax) выступает в качестве важнейшего индикатора зрелости облачной системы и наличия крупнодисперсной ледяной фазы. В экспериментах z01, z02 и z03 пиковая отражаемость стабилизирована на отметке 66.7 дБZ, что однозначно свидетельствует о формировании мощного градового ядра (таблица 1). Начиная с эксперимента z04 (δ = 3%), значение Zmax начинает монотонно снижаться, достигая 60.9 дБZ в сценарии z06. Снижение отражаемости на ~9 дБZ физически означает резкое падение концентрации и эффективного диаметра градин и крупных капель дождя.
Сводные количественные характеристики экспериментов приведены в таблице 1.
Таблица 1 - Динамические параметры, максимальная радиолокационная отражаемость, макроструктурные характеристики облака, интегральные массы гидрометеоров, общее выделение скрытой теплоты конденсации (le) и суммарный объем осадков на поверхность
Эксп. | Дефицит qv, % | Пик wmax (wmin), м/с | Время пика, мин | Относ. изм. пика wmax, % | Макс. высота Hвг, км Пик (Zmax), дБZ | Время пика, мин | Относ. изм. Zmax, % | Интегральная масса qc+qr, (qi+qs+qg+qh), 106 кг | Изменение интегральной массы qg+qh, % | Интеграл le, 1015 Дж Суммарные осадки (90 мин), мм | Относ. изм. осадков, % |
z01 | 0 | 45,64(14,31) | 38 (38) | 0,0 | 13,8 (66,68) | 40 (56) | 0,0 | 11,274 (14,5) | 0,0 | -5,19 (13,1) | 0,0 |
z02 | 1 | 45,21 (14,11) | 38 (39) | -0,9 | 13,8 (66,7) | 40 (56) | 0,0 | 11,272 (14,5) | 0,0 | -5,18 (13,1) | 0,0 |
z03 | 2 | 43,81 (13,77) | 38 (39) | -4,0 | 13,6 (66,7) | 41 (56) | 0,0 | 11,273 (14,3) | -1,4 | -5,18 (13,1) | 0,0 |
z04 | 3 | 41,40 (13,74) | 41 (41) | -9,3 | 13,6 (66,4) | 41 (56) | -0,5 | 10,568 (12,0) | -17,2 | -5,14 (13,0) | -0,8 |
z05 | 5 | 35,39 (11,93) | 44 (45) | -22,5 | 13,4 (64,6) | 42 (56) | -3,1 | 10,106 (9,7) | -33,1 | -4,96 (12,6) | -3,8 |
z06 | 10 | 22,29 (8,74) | 52 (52) | -51,2 | 12,6 (60,9) | 45 (57) | -8,7 | 6,770 (3,2) | -78,0 | -4,07 (10,7) | -18,3 |
Макроструктурная эволюция облачной геометрии также выявила деградацию кучево-дождевого облака. Максимальная высота вертикального развития (верхняя граница облака Hвг) в контрольном пуске составляет 13.8 км, что указывает на пробой тропопаузы и формирование развитой наковальни. При иссушении атмосферного столба на 10% (z06) верхняя граница падает до 12.6 км (снижение на 1,2 км). Одновременно с этим нижняя граница облачности испытывает подъем с 1.5 км до 2.3 км вследствие дефицита влаги в подоблачном слое и соответствующего падения уровня конденсации. В результате суммарная геометрическая мощность облака уменьшается более чем на 2 км.
Пороговый характер отклика облачной системы прослеживается и при анализе интегральных по области моделирования масс отдельных микрофизических фракций влаги (таблица 1). В то время как пиковая интегральная масса водяного пара (qv) в силу геометрии начального возмущения остается квазипостоянной (0,4805 х 109 кг), массы конденсированных фаз демонстрируют стремительное нелинейное сокращение.
Масса мелкодисперсной облачной воды (qc) снижается на относительно умеренную величину: с 10,93 х 106 кг в контрольном эксперименте до 6,62 х 106 кг в эксперименте z06 (снижение на 39,4%). В то же время масса укрупненной фракции жидких осадков (qr — дождь) испытывает выраженный коллапс после преодоления 3%-го барьера, сокращаясь в эксперименте z06 более чем в 2,3 раза (до 1,5 х 107 кг).
Наиболее уязвимым звеном облачного микрофизического «конвейера» оказалась плотная ледяная фаза — град и крупа (qh и qg) соответственно. До уровня иссушения в 2% масса града стабильна (14,3 х 106 кг). Однако в эксперименте z04 (δ = 3%) она скачкообразно падает до 12,0 х 106 кг (уменьшение на 17,2%), в эксперименте z05 теряет треть массы, а при 10%-м иссушении (z06) сокращается на 78%, опускаясь до 3,2 х 106 кг. Столь катастрофический спад указывает на то, что механизмы градообразования требуют строго определенного уровня влагосодержания тропосферы, ниже которого быстрый рост твердой фазы блокируется.
Суммарный приток энергии за счет фазового перехода первого рода — выделения тепла при конденсации пара (le) — монотонно и нелинейно убывает. В контроле суммарный интеграл выделенного скрытого тепла составляет 5,19 х 1015 Дж, тогда как в эксперименте z06 он падает до 4,07 х 1015 Дж (падение на 21,6%, таблица 1).
Количество суммарных накопленных осадков на поверхности к 90-й минуте моделирования демонстрирует устойчивость к слабым флуктуациям влажности (вплоть до δ = 2% объем осадков неизменно составляет 13.1 мм). Перелом наступает при преодолении 3-5%-го барьера иссушения. В финальном эксперименте z06 количество выпавших осадков снижается до 10.7 мм, что соответствует сокращению осадков на 18,3%. То есть, интегральные параметры энергетического обмена и влагоотдачи полностью подтверждают системный характер подавления конвекции.
4. Обсуждение
Фундаментальный характер нелинейного отклика динамики и микрофизики облака на уменьшение влажности находит объяснение в рамках классической теории вовлечения и концепции «буферного сортировщика плавучести» (buoyancy sorting)
. Согласно этим представлениям, учтенным в термодинамическом блоке модели, восходящий конвективный поток представляет собой гетерогенный ансамбль элементарных воздушных объемов («пузырей»), которые в процессе подъема непрерывно взаимодействуют со средой, испытывая различную степень турбулентного перемешивания и вовлечения сухого окружающего воздуха.При уменьшении начальной удельной влажности во всем вертикальном столбе атмосферы парциальное давление водяного пара падает, увеличивая дефицит насыщения тропосферы. В процессе вовлечения — втягивание более сухого воздуха вызывает мгновенное испарение части облачной капельной фракции для восстановления насыщения. Процесс испарения сопровождается поглощением скрытой теплоты, что резко охлаждает вовлеченный объем воздуха. В результате виртуальная потенциальная температура элементарного объема падает ниже температуры окружающей среды, плавучесть становится отрицательной, и данный объем выбывает из восходящего ядра, формируя облачные нисходящие потоки. При иссушении атмосферы доля конвективных объемов, способных сохранить положительную плавучесть и достичь уровня свободной конвекции, падает нелинейно, что объясняет резкое падение интегральных параметров восходящего потока.
Ключевым результатом данной работы является четкая фиксация бифуркационного порога в районе Δqv ≈ 3%. В диапазоне дефицита влаги 122% конвективная система успешно компенсирует нехватку влаги. Механизм компенсации реализуется за счет локального сужения восходящего ядра, что временно снижает эффективный периметр вовлечения и защищает центральную часть восходящего потока от разрушительного испарения капель.
Блок-схема каскадного разрушения обратных связей глубокой конвекции при преодолении порогового дефицита влажности
Микрофизические процессы спада градообразования при Δqv = 10% показывают, что интегральная масса qg+qh падает на 78%. Это объясняется тем, что рост зародышей града требует строго определенной критической водности облачной среды. При дефиците пара скорость конденсации снижается, и массовая доля переохлажденных капель падает ниже критического предела. Ледяные ядра, не встречая достаточного количества жидкой влаги для аккреции, прекращают рост и сублимируют или выпадают в виде мелкой крупы и снега, что подтверждается резким падением радиолокационной отражаемости на 9 дБZ. То есть, снижение уровня нейтральной плавучести ниже изотермы 0 °C заблокировало развитие глубокой конвекции, ограничив эволюцию облака стадией мощных кучевых (cumulus congestus).
Обнаруженный пороговый переход согласуется с нелинейной зависимостью интенсивности осадков от влагосодержания столба атмосферы, описанной в наблюдениях тропической конвекции
, . Физически этот порог связан с тем, что при достаточно низкой влажности уровень нейтральной плавучести опускается ниже уровня замерзания, и эффективность микрофизических механизмов роста ледяных частиц резко снижается. Тем самым прерывается цепочка обратных связей, поддерживающая глубокую конвекцию (конденсация → скрытая теплота → увеличение плавучести → рост облака → дополнительная конденсация).С точки зрения активного воздействия на облака и осадки результаты работы показывают, что при иссушении атмосферы более чем на 10-20% относительно нормы эффективность воздействия реагентами (например, йодистым серебром) существенно снижается, поскольку уже сама реализующаяся высота облачной вершины не превышает уровня замерзания. Полученные данные согласуются и с результатами наблюдений
, указывавших на подавление конвективной активности в периоды мезомасштабных адвекций сухого воздуха с востока.
5. Заключение
На основе результатов серии численных экспериментов, выполненных с помощью трехмерной негидростатической модели конвективного облака с детальным учетом термодинамических, микрофизических и электрических процессов, исследована чувствительность развивающегося конвективного облака к равномерному дефициту удельной влажности в вертикальном столбе атмосферы. Сформулированы следующие основные выводы:
- имеет место выраженный нелинейный характер отклика динамических и микрофизических параметров облака на иссушение атмосферы с фиксацией критического порога перелома зависимостей в районе Δqv ≈ 3%. При дефиците влаги до 2% облачная система сохраняет пиковые характеристики за счет внутренней структурной перестройки.
- преодоление 3%-го порога влажностного дефицита запускает каскадный механизм подавления конвекции. В финальной точке расчетной матрицы (Δqv = 10%) максимальная скорость восходящего потока падает на 51,2% (до 22.29 м/с), а интегральная вертикальная активность сокращается в 2,7 раза. Эволюция облака существенно замедляется, сдвигая динамический пик на 14 минут позже относительно контроля.
- селективная деградация плотной ледяной фракции гидрометеоров: при 10%-м иссушении интегральная масса града и крупной крупы по области моделирования падает на 78%, тогда как масса облачных капель сокращается лишь на 39,4%. Это указывает на высокую уязвимость фазовых переходов второго рода (вода-лед) к влажностному режиму тропосферы.
- макроструктурные параметры облака демонстрируют монотонный спад: под влиянием дефицита влаги верхняя граница облачной вершины снижается на 1,2 км, уменьшая общую вертикальную мощность кучево-дождевой системы. Суммарные приземные накопленные осадки уменьшаются на 18,3%.
Полученные количественные закономерности и выявленные аналитические зависимости имеют практическое значение для прогнозирования конвективных осадков при мезомасштабном иссушении атмосферы, для планирования операций по активному воздействию на облака, а также для верификации параметризаций конвекции в глобальных климатических моделях в части описания чувствительности к профилям влажности в тропосфере.