CC BY
15
НАУКИ О ЗЕМЛЕ
«НАУКА. ИННОВАЦИИ. ТЕХНОЛОГИИ», № 4, 2018
25.00.29
ФИЗИКА АТМОСФЕРЫ И ГИДРОСФЕРЫ
ДАНИЛОВА Н.Е., СЕМЕНОВА Ю.А., СМЕРЕК Ю.Л., ЗАКИНЯН Р.Г.
УДК 551.513.22
Северо-Кавказский федеральный университет, г. Ставрополь, Россия
zakinyan@mail.ru
ВЛИЯНИЕ ПОДОБЛАЧНОЙ КОНВЕКЦИИ НА РАЗВИТИЕ ОБЛАЧНОЙ КОНВЕКЦИИ
Введение: Облачная конвекция является причиной опасных явлений погоды [23]. Так
сильные дожди и последующие наводнения являются прямым следствием конвективных штормов. Явленияшторма часто сильно локализованы, поэтому прогнозирование того, почему, где и когда шторм развивается, является одной из самых сложных проблем в прогнозе погоды [10]. На сегодняшний день оперативные численные модели часто не могут предсказать местоположение и время начала облачной конвекции при инициировании ее процесса-мив подоблачном слое, поскольку эти процессы недостаточно хорошо представлены в моделях. Однако влияние этих процессов на инициирование облачной конвекции часто отмечается в исследованиях [25, 39], так что для прогнозирования начала облачной конвекции требуется знание термодинамики и динамикиподоблачного слоя.
Материалы и методы: Для прогноза параметров облачной конвекции используется данные радиозондирования, которые наносятся на аэрологическую диаграмму. Как правило, при этом делается допущение, что скорость восходящих потоков и перегрев на уровне конденсации равны нулю. Но это не всегда так. Поэтому представляется важным определить, как параметры приземного слоя влияют на значения параметров облачной конвекции на уровне конденсации.
Результаты
исследования: в статье разработана математическая модель конвекции влажного подоблач-
ного слоя атмосферы. Математическая модель представляет собой систему уравнений движений, теплопроводности и неразрывности, записанные в двумерной вертикальной плоскости. Получено аналитическое решение указанной системы. Установлены выражения для уровня конвекции, уровня максимальной скорости, для профиля вертикальной и горизонтальной составляющей скорости, вертикальное распределение температуры. Из указанных выражений выведен критерий развития подоблачной конвекции. Полученные выражения позволяют в свою очередь вывести выражения для параметров конвекции на уровне конденсации. А это в свою очередь позволяет вывести критерий развития облачной конвекции. Установлено, что развитие облачной конвекции зависит, как от абсолютных значений дефицита точки росы в приземном слое, так и в большей степени от вертикального градиента массовой доли водяного пара. В частности при критическом значении вертикального градиента массовой доли водяного пара наблюдается «взрывной рост» конвекции.
Обсуждение и
заключение: Получены аналитические выражения для дефицитов точки росы в приземном
слое, определяющие значения функции перегрева и скорости восходящих потоков на уровне конденсации. Показано, что на развитие подоблачной конвекции существенное влияние оказывает значение массовой доли водяного пара. Установлено, что при значениях градиента массовой доли водяного пара, приближающихся к критическому значению, наблюдается «взрывной» рост конвективных ячеек в подоблачном слое. Также показано, что существует критическое значение градиента массовой доли водяного, способствующее развитию облачной конвекции.
Ключевые слова: приземный слой атмосферы, подоблачная конвекция, уровень конденсации, облачная конвекция, двумерная модель конвекции, аналитическое решение/
Danilova N.Ye., Semyonova Yu.A., SmerekYu. L, Zakinyan R.G.
North-Caucasian Federal University, Stavropol, Russia,
zakinyan@mail.ru
Influence under-cloudy convections on development of a cloudy convection
Introduction:
Materials and methods:
Results:
Discussion and conclusion:
Keywords:
Cloud convection causes hazardous weather [23]. So heavy rains and subsequent floods are a direct consequence of convective storms. Storm phenomena are often highly localized, so predicting why, where and when a storm develops is one of the most difficult problems in weather forecasting [10]. Today, operational numerical models often cannot predict the location and time of the onset of cloud convection when it is initiated in the skid layer, since these processes are not well represented in the models. However, the influence of these processes on the initiation of cloud convection is often noted in studies [25, 39], so that to predict the beginning of cloud convection requires knowledge of the thermodynamics and dynamics of the sub-cloud layer. To predict the parameters of cloud convection, radiosonde data is used, which are plotted on the upper-air diagram. As a rule, it is assumed that the rate of ascending flows and overheating at the level of condensation are zero. But it's not always the case. Therefore, it is important to determine how the parameters of the surface layer affect the values of the parameters of cloud convection at the level of condensation. In the article the mathematical model of a convection moist under-cloud layer atmosphere is developed. The mathematical model represents system of the equations of movements, heat conductivity and the continuity, recorded in a two-dimensional vertical plane. The analytical solution of the indicated system is received. Expressions for level of convection, level of the maximum velocity, for vertical and horizontal components ofvelocity, vertical distribution of temperature are established. From the indicated expressions the criterion of development under cloudy convections is derived. The received expressions allow derived in turn expressions for convection parameters on a condensation level. And it in turn allows derived criterion of development of a cloudy convection. It is established that development of a cloudy convection depends, as from absolute values of deficit of dew point in a ground layer, and in a greater degree from a vertical gradient of a mass fraction of a water vapor. In particular at critical value of a vertical gradient of a mass fraction of water vapor «explosive growth» convections is observed.
Analytical expressions are obtained for the dew point deficiencies in the surface layer, which determine the values of the overheating function and the ascending flow rate at the condensation level. It has been shown that the development of subliminal convection is significantly affected by the value of the mass fraction of water vapor. It has been established that with the values of the gradient of the mass fraction of water vapor approaching the critical value, there is an "explosive" growth of convective cells in the sublayer. It is also shown that there is a critical value of the gradient of the mass fraction of water, contributing to the development of cloud convection. surface layer of atmosphere, under-cloudy convection, a condensation level, a cloudy convection, two-dimensional model of convection, the analytical solution
Введение
Очевидно, что характеристики приземного слоя влияют на развитие конвекции в подоблачном слое, что в свою очередь сказывается на развитии облачной конвекции. Вопрос заключается в том, какие характеристики приземного слоя в большей или меньшей степени влияют на развитие, в конечном счете, облачной конвекции. Так, в статье [11] приводятся экспериментальные подтверждения влияния вертикального профиля влажности на развитие конвекции. Исследованию влияния вертикального распределения
НАУКИ О ЗЕМЛЕ
Влияние подоблачной конвекции на развитие облачной конвекции Данилова Н.Е., Семенова Ю.А., Смерек Ю.Л., Закинян Р.Г.
температуры и влажности приземного слоя на состояние свободной тропосферы посвящена статья [33].
Надо заметить, что, хотя конвективное движение является одним из самых распространенных видов движений в атмосфере, как физика этого процесса, так и его математическое описание остаются открытыми. Считается, что конвекция начинается с формированием термиков [14]. Физика образования термиков аналогична фазовому переходу. Здесь аналогично при определенных условиях возникает термик критического размера. Если размер тер-мика меньше критического, то он неустойчив и со временем исчезает. Если же размер термика больше критического размера, то он устойчив. Далее он отрывается и растет. Вопрос о размере критического термика с точки зрения теории остается открытым. Как правило, его определяют из баланса силы плавучести и силы сопротивления Стокса.
Есть еще один аспект, который делает представленное в настоящей работе исследование актуальным. В настоящее время эффекты, связанные с мелкомасштабными процессами, к которым относится мелкая конвекция параметризуются в общей циркуляционной, региональной или мезомасштаб-ной моделях. От того, насколько хорошо мы понимаем с точки зрения физики и умеем описывать эти процессы, будет зависеть качество прогнозов в глобальных и региональных моделях [17, 28]. Подоблачная конвекция определяет вертикальное распределение влажности, теплоты и импульсаи играет клю-чевуюроль в глобальном гидрологическом и энергетическом цикле [2, 7, 16, 24, 34, 41]. Для углубления понимания роли процессов подоблачного слоя в климатических моделях использовались простые одномерные модели [3, 3, 8, 12, 15, 26, 33] и их совершенствование является актуальной задачей.
Облачная конвекция является причиной опасных явлений погоды [23]. Так сильные дожди и последующие наводнения являются прямым следствием конвективных штормов. Явленияшторма часто сильно локализованы, поэтому прогнозирование того, почему, где и когда шторм развивается, является одной из самых сложных проблем в прогнозе погоды [10]. На сегодняшний день оперативные численные модели часто не могут предсказать местоположение и время начала облачной конвекции при инициировании ее процесса-мив подоблачном слое, поскольку эти процессы недостаточно хорошо представлены в моделях. Однако влияние этих процессов на инициирование облачной конвекции часто отмечается в исследованиях [25, 39], так что для прогнозирования начала облачной конвекции требуется знание термодинамики и динамикиподоблачного слоя.
Проведенные ранее полевые исследования дают представление о процессах, вызывающих инициирование облачной конвекции [18, 35, 37, 39]. Хотя имеется хорошее понимание общей динамики рассматриваемых процессов, начальная стадия развития конвективных ячеек остается недостаточно изученной [4].
В ряде исследований предполагалось [9, 11], что инициирование облачной конвекции очень чувствительно к вариациям температуры и влажности в подоблачном слое.В частности, поле влажности подоблачного слоя является ключевым фактором для понимания инициирования осадкообразующей конвекции [20]. Небольшие изменения в массовой доле водяного пара могут управлять процессом и определять, происходят ли штормы или нет. И что особенно важно, изменения приблизительно в 1 g/kg в массовой доле водяного пара могут дать диаметрально противоположные прогнозы:отсутствие конвекции или интенсивная конвекция [11]. Много исследований проведено по изучению влияния влажности поверхности на развитие конвекции [31, 44]. Главный механизм положительной обратной связи заключается в том, что влажный грунт обогащает воздух пограничного слоя влажностью [5, 13].
Таким образом, понимание и успешное предсказание штормового инициирования зависят от точногоопределенияначальныхтермодинамических и кинематических полей подоблачного слоя, с особым вниманием к полям влажности водяного пара. Вообще, типичное разрешение обычных систем наблюдения, таких как радиозонды, недостаточно для того, чтобы наблюдать изменчивость параметров в подоблачном слое атмосферы. В ряде исследований [37] предполагалось, что критическим фактором, ограничивающим прогноз конвективных осадков, является неопределенность измерения вертикального распространения водяного пара. Радиозонды, являющиеся традиционными средствами измерения содержания водяного пара, недостаточны для этих целей. Это связано с тем, что они предоставляют информацию о вертикальном профиле в широко распределенных пунктах зондирования и обычно доступны только два раза в день, иногда содержа существенные погрешности [32, 36].
Таким образом, влиянию характеристик приземного слоя атмосферы на развитие, как крупномасштабной циркуляции [1, 27], так и развитие облачной конвекции посвящено много работ [6].
В настоящей работе и подоблачная и облачная конвекция описываются с помощью двухмерной модели. Находятся соответствующие аналитические решения, а также выявляются критерии влияния параметров приземного слоя атмосферы на развитие облачной конвекции.
Материалы и методы исследования
Для прогноза параметров облачной конвекции используется данные радиозондирования, которые наносятся на аэрологическую диаграмму. Как правило, при этом делается допущение, что скорость восходящих потоков и перегрев на уровне конденсации равны нулю. А это не всегда так. Поэтому представляется важным определить, как параметры приземного слоя влияют на значения параметров облачной конвекции на уровне конденсации.
НАУКИ О ЗЕМЛЕ
Влияние подоблачной конвекции на развитие облачной конвекции Данилова Н.Е., Семенова Ю.А., Смерек Ю.Л., Закинян Р.Г.
В работе [42] было получено аналитическое решение двумерной задачи тепловой конвекции сухого воздуха в атмосфере. Причем воздух считался абсолютно сухим. А в работе [43] получено аналитическое решение двумерной модели конвекции влажногонасыщенного воздуха. Последний случай для краткости изложения будем называть «облачной» конвекцией. Таким образом, у нас есть аналитическое решение двумерной модели конвекции сухого воздуха, которую будем называть «подоблачной» конвекцией, и аналитическое решение «облачной» конвекции. Целью настоящей работы является применить полученные в работах [42, 43] результаты для анализа влияния характеристик приземного слоя атмосферы на развитие облачной конвекции.
Для описания облачной конвекции нужно знать параметры конвекции, такие как функция перегрева, скорость восходящих потоков и массовая доля водяного пара на уровне конденсации. Эти параметры могут быть найдены из рассмотрения подоблачной конвекции. Она в свою очередь определяется характеристиками приземного слоя. Таким образом, характеристики приземного слоя определяют характер облачной конвекции. Задача настоящей работы состоит в том, чтобы эту связь установить. Ранее, как было отмечено во введении, такие исследования выполнялись численно и экспериментально или же использовались одномерные модели конвекции. Здесь же мы исследуем данную задачу аналитически в рамках двумерной нелинейной модели конвекции. Данный анализ, как уже было сказано, основан на полученных нами ранее аналитических решениях конвекции сухого и влажного насыщенного воздуха.
ходимо, предварительно, описать конвекцию влажного ненасыщенного воздуха. По аналогии с работой [42], где мы получили аналитическое решение двумерной модели конвекциисухого воздуха,получим теперь решение для конвекции влажного ненасыщенного воздуха и рассмотрим влияние параметров приземного слоя атмосферы на развитие облачной конвекции, т.е. на значения параметров конвекции на уровне конденсации.
В настоящем разделе мы шаг за шагом повторим вывод, проделанный нами в работе [42], но учтем, что воздух теперь влажный. Уравнение движения рассматривается в плоскости (х, z). Давление представляется в виде суммы давлений в состоянии статики и его возмущения: р = р + р'. Плотность влажного воздуха определим выражением:
уравнения тепловой конвекции влажного воздуха
Подоблачный воздух не является сухим. Поэтому нам необ-
(1)
где а=ж - коэффициент теплового расширения воздуха (Т0 = 273 К);
АГ(г) = Т- (г)-Ге (г) - функция перегрева;
ДТ , Те - температура воздуха, соответственно, возмущенной и невозмущенной атмосферы; я - массовая доля водяного пара;
Ля = ^ — яе- функция пересыщения;
p=§í-l=a608 ;
,М(1
:Му
Md = 29 §/шо1- молярная масса сухого воздуха;
Му = 18 ^шо1 - молярная масса водяного пара.
Здесь р; - плотность воздуха возмущенной атмосферы;
ре - плотность воздуха в состоянии статикиатмосферы. Пара-
метры окружающей атмосферы мы рассматриваем как невозмущенное состояние. Для атмосферы аДТ < 1. Считается, что подъем воздушной частицы происходит адиабатически. Другими словами, при приближении (1) предполагается, что плотность воздуха зависит от температуры и массовой доли водяного пара и не зависит от давления [14, 21].
Функция перегрева равна:
AT (z ) = A0T -Ду z, (2)
где AoT,=7jo-2^o - значение функции перегрева вблизи поверхности земли;
Ду = Ya - Y, где у - градиент температуры воздуха невозмущенной атмосферы в состоянии статики; уа - сухоадиабатический градиент температуры;
Tío - температура воздуха вблизи поверхности земли возмущен-
ной атмосферы;
Te0 - температура воздуха вблизи поверхности земли невозму-
щенной атмосферы в состоянии статики.
Для адиабатически поднимающейся воздушной частицы массовая доля водяного пара остается постоянной sí(z) = si0 , где si0 - массовая доля водяного пара вблизи поверхности земли. Пусть в невозмущенной атмосфере в состоянии статики массовая доля водяного пара убывает по линейному закону
¿e (z) = %) -bz , (3)
где b - градиент массовой доли водяного пара. Тогда для функции
пересыщения получим выражение
Д (z) = A 0 ¿ + bz, (4)
где Д o s = si0 - se0
НАУКИ О ЗЕМЛЕ
Влияние подоблачной конвекции на развитие облачной конвекции Данилова Н.Е., Семенова Ю.А., Смерек Ю.Л., Закинян Р.Г.
- пересыщение у поверхности земли. Выражение для плотности воздушной частицы запишется в виде
Pi = Ре [1 -(аДоГ + РД о 5 ) + (аДу-рЬ ) ^ ]. (5)
Отсюда найдем уровень выравнивания плотностей воздуха в воздушной частице и невозмущенной атмосфере
2 = аАоТ_+РДо£^. (6)
Р аДу - РЬ
Из формулы (6), формально принимая в = 0, получим выражение для уровня выравнивания температур которое следует непосредственно из формулы (2). Отсюда видно, что уровень выравнивания плотностей воздуха для влажной частицы выше, чем уровень выравнивания температур.
Принимая во внимание уравнение статики, которое имеет место для невозмущенной атмосферы, система уравнений, описывающих тепловую конвекцию влажного ненасыщенного воздуха в стационарном состоянии, представится в виде:
ди ди 1 др'
и— + т— =---—, (/)
дх дх ре дх
и дТ + М' дТ = ^ (аДоГ + в Д о 5 - (аДу - РЬ ) • х), (8)
^= о. (9)
дх Эх
Здесь мы при записи уравнения (8) сделали допущение, такое же, как в работе [22]:
-—1Г<< 8 (аДТ + РД?) .
Ре дх
Из уравнения неразрывности (9) следует, что можно ввести функцию тока у следующим образом:
и = ——, т = —. (Ю)
дх дх
Сравнивая эту систему с системой уравнений тепловой конвекции сухого воздуха [42], замечаем, что они отличаются только лишь правой частью в уравнении (8). Поэтому решение системы уравнений (7) - (9) находится аналогично системе уравнений для сухого воздуха.
Из формулы (6) следует, что, если градиент массовой доли водяного пара будет приближаться к критическому значению, равному Ьсг = аДу/р , уровень выравнивания плотностей будет стремиться к бесконечности, произой-
дет так называемое «взрывное» усиление конвекции. Данный факт является отличительной особенностью конвекции во влажной атмосфере по отношению к конвекции сухого воздуха.
Таким образом, мы видим, что критическим параметром для конвекции является не абсолютное значение массовой доли водяного пара вблизи поверхности земли, а его вертикальный градиент.
Найдем значение критического градиента массовой доли водяного пара при условии, что у = 0:
(ьсг) =^! = ш .ю-5 т-1. (10)
^ " 'тах р
Для случая конвекции влажного ненасыщенного воздуха функция тока имеет вид
^ (аАу-рь) г-г-)
у = ----^ (2zр - z)cos кх. (11)
Отсюда для проекций скорости получим выражения:
„ >(aAy-Pb) ■ ^Z cos k*, (12)
k ф ( 2Zp- z) w = g (aAy-pb) z(2zp - z)sin kx. (13)
Отсюда из условия w = 0 для уровня конвекции влажного ненасыщенного воздуха получим выражение
„aA0T + RA 0 5
zw = 2 zp = 2—0-(14)
w p aAy-pb • (14)
Из (14) видно, что уровень конвекции зависит, как от перегрева и пересыщения в приземном слое, так и от вертикальных градиентов температуры и массовой доли водяного пара. Оценки параметров, входящих в формулу (14) показывают, что фактор перегрева Д0 T сильнее влияет на развитие конвекции, чем фактор пересыщения Д0 5. Как видно из формулы (14) наиболее существенным фактором является градиент массовой доли водяного пара b. Если градиент массовой доли водяного пара приближается к критическому значению, то происходит «взрывное» усиление конвекции, т. е. уровень конвекции растет и стремится к бесконечности. По поводу «взрывного» роста следует заметить, что в нашей модели конвекции не учитывается вязкость воздуха. Поэтому, учет вязкости ограничит рост уровня конвекции.
Уровень максимальных скоростей найдем из условия: dw/dz = 0. Взяв производную по переменной z в выражении (13), получим
НАукИ о ЗЕмлЕ
Влияние подоблачной конвекции на развитие облачной конвекции Данилова Н.Е., Семенова Ю.А., Смерек Ю.Л., Закинян Р.Г.
— = NBV Zp Z sin kx = 0, (15) f ( 2 zp- z )
где величину Nbv = ^Jg (aAy - P») (16)
назовем частотой Брента - Вяйсяля колебаний влажного ненасыщенного воздуха. Из этой формулы видно, что при критическом градиенте массовой доли водяного пара частота Брента - Вяйсяля равна нулю, а значит, период колебаний стремится к бесконечности. Из формулы (15) следует, что уровень максимальных скоростей совпадает с уровнем выравнивания плотностей: max = Zp.
Подставляя выражение для величины zp в формулу (13), для максимальной скорости восходящих потоков влажного ненасыщенного воздуха получим равенство
wmax = ^BVzpsin kc = ^aAyg_p¿ (aA0T + PA0S)sin kx. (17)
Запишем выражение для квадрата скорости восходящих потоков влажного ненасыщенного воздуха:
w2 =[g (aAy —pb) z (2zp — z)] sin2kx. (18)
Отсюда следует, что, если zp = 0, т. е. плотности воздуха возмущенной и невозмущенной атмосферы у поверхности земли совпадают, то условием развития конвекции во влажной атмосфере является условие b > bcr. В частности, для сухой атмосферы отсюда следует известное условие развития конвекции: у > уа. Если же атмосфера влажная, то конвекция развивается даже при условии у < уа.
Тогда для квадрата скорости восходящих потоков получим выражение
w = -^Pg (b — bcr) z • sin kx, (19)
т. е. скорость восходящих потоков будет расти с высотой по линейному закону и уровень конвекции не будет достигнут, даже если у < уа. Ясно, что полученный результат является идеализацией. В реальной атмосфере инверсионные слои и вязкость воздуха (напомним, что в нашей модели рассматривается уравнение движение идеальной жидкости) не дадут реализоваться указанному сценарию.
Выражение для амплитуды максимальной скорости восходящего потока во влажной атмосфере имеет вид
"0 р» (aA0" + PA0(20)
Таким образом, из полученного выражения (20) для амплитуды скорости восходящего потока видно, что вертикальная скорость зависит от неустойчивости атмосферы, градиента массовой доли водяного пара в окружающей атмосфере и перегрева у поверхности земли. Из формулы (20) также видно, что при критическом значении массовой доли водяного пара максимальная скорость стремится к бесконечности («взрывной рост»). Значения параметров конвекции на уровне конденсации Определим уровень конденсации по формуле [21]:
„ _ Т0 - т0 _ ¿0
с Та-Тт ДТт' где т0 - температура точки росы у поверхности земли;
¿0 - дефицит точки росы вблизи поверхности земли ^ = 0);
градиент температуры точкиросы.
(21)
У
т
Как известно, формула (21) дает заниженные значения уровня конденсации, так как при ее выводе не учитывается вовлечение поднимающегося воздуха с окружающей средой. Для дальнейшего анализа имеет значение только лишь зависимость уровня конденсации от дефицита точки росы вблизи поверхности земли.
На уровне конденсации zc выражение для скорости восходящих потоков будет иметь вид:
м>с2=2£(аЛ0Г+рД05)|й-
^ аДу~рЬЛ
вт2^. (22)
2(аА0Г+рД05) Дух
Подставляя выражение (21) в формулу (2), получим формулу для функции перегрева на уровне конденсации
ЛСГ = Л0Т - ¿0 . (23)
АУт
Для функции пересыщения на уровне конденсации получим выражение
ЛС5 = Л05 + ь-^. (24)
ЛУт
Из выражений (22) и (23) видно, что существуют критические значения дефицита точки росы вблизи поверхности земли, при которых соответственно равны нулю функция перегрева на уровне конденсации ДсТ и скорость восходящих потоков воздуха (на уровне конденсации):
ЛсТ = 0: (¿0 )сг1 =Л0Т А^, (25)
^с = 0: (¿0 )сг2 = 2 (аЛ0Т + РЛ 0 ^^Л^ • (26)
НАУКИ О ЗЕМЛЕ
Влияние подоблачной конвекции на развитие облачной конвекции Данилова Н.Е., Семенова Ю.А., Смерек Ю.Л., Закинян Р.Г.
Из формулы (26) следует, что влажность подоблачного слоя будет способствовать развитию облачной конвекции. Более того, из этой формулы следует, что при значении градиента массовой доли водяного пара, равного критическому значению Ьа= ^р , критическое значение дефицита точки росы вблизи поверхности земли (^0)сг стремится к бесконечности, а это значит, что подоблачная конвекция будет всегда «пробивать» уровень конденсации, т. е. ис > 0.
Полученные критерии определяют начальные условия облачной конвекции.
1. Если дефицит точки росы вблизи поверхности земли равен второму критическому значению й0 = (¿0)с12 , то на уровне конденсации перегрев будет отрицательным:
а скорость восходящего потока на этом же уровне будет равна нулю (ис = 0). Следовательно, в этом случае подоблачный слой не будет способствовать развитию облачной конвекции. То есть атмосфера при таком дефиците точки росы вблизи поверхности земли недостаточно влажная, и подоблачный слой не будет пробиваться восходящим воздухом (рис. 1, график 1).
Формально принимая в = 0, получим, что для «сухой подоблачной атмосферы» функция перегрева на уровне конденсации будет равна АсТ = -А0Т < 0. Данный результат известен из адиабатической модели конвекции сухого воздуха [21].
2. Если дефицит точки росы вблизи поверхности земли
(^0)сг1 < ^о < (^0)сг2, то перегрев на уровне конденсации будет отрицательным, но по модулю меньше максимального значения -(АсТ)тах <АсТ < 0, а скорость восходящего потока на этом же уровне положительна. В этом случае уровень конденсации будет пробиваться восходящей подоблачной струёй, что способствует развитию облачной конвекции (рис. 1, график 2).
3. Если дефицит точки росы вблизи поверхности земли равен первому критическому значению й0 < (^0)сг1 , то на уровне конденсации перегрев равен нулю ДсТ = 0, а скорость восходящего потока будет равна ис ^ (меньше максимального значения), т. е. равна значению на уровне zt выравнивания температур:
(аАу + РЬ ) А0Т + 2рАуА 0 5
(АсТ) < 0,
V ^ / тах
(27)
аАу-рЬ
или
(28)
(29)
Заметим, что в этом случае уровень конденсации совпадает с уровнем выравнивания температур: zc = zt. Аналогично, функция пересыщения примет значение, равное значению на этом же уровне zt:
Лс5 = Л 05 + Ь = Л 05 + bzt. (30)
Лу
Следовательно, в этом случае будут существовать условия для развития облачной конвекции (рис. 1, график 3).
4. Если дефицит точки росы вблизи поверхности земли меньше первого критического значения дефицита < (¿0)СгЬ то перегрев на уровне конденсации положителен Д0Т > 0, а скорость восходящего потока будет положительной, но меньше максимального значения, что также может способствовать развитию облачной конвекции (рис. 1, график 4).
МЧ^ ьст=-Мтах, ¿0<К)сг2, -(Д0г)тах<Лсг<0,
^тах = 0 ^с>0, 5и>/&<0
Из выражения (22) видно, что, если градиент массовой доли водяного пара в окружающей атмосфере принимает критическое значение Ьсг , то скорость восходящих потоков на уровне конденсации при любом значении дефицита точки росы вблизи поверхности земли положительна и равна:
= 2 Я (аЛ0Т + РЛ 0 5 (31)
ЛУт
Итак, учет влажности подоблачного слоя приводит не только к количественным изменениям критериев, определяющих параметры облачной конвекции на уровне конденсации, но и качественным изменениям. Мы видим, что существенное влияние на параметры облачной конвекции на уровне конденсации имеет не столько значение абсолютной влажности подоблачного слоя, сколько градиент массовой водяного пара в подоблачном слое.
Таким образом, полученное в работеаналитическое выражение для двумерной модели конвекции позволило провести анализ влияния характеристик приземного слоя атмосферына развитие облачной конвекции.
Результаты исследования и их обсуждение
В работе [43] получено аналитическое решение двумерной модели конвекции влажного насыщенного воздуха. Приведем основные положения и допущения, на которых было основано аналитическое решение задачи в работе [43].
НАУКИ О ЗЕМЛЕ
Влияние подоблачной конвекции на развитие облачной конвекции Данилова Н.Е., Семенова Ю.А., Смерек Ю.Л., Закинян Р.Г.
Zc _ Zw I z
*4 V 4 4 N 4 4 4 N 4 \ \ \ V l Wmax
ДТ, w -^-►
До Т
График 1
d0 < (d0)cr1 , ДСТ = 0, dw/dz = 0,
A Z
zp zc zw max
1 Wm
ДТ, w
ДоТ
1 z
Zw
Zc 4 V 4 4 N
\ 4 4 N 4 \ \ \ V l Wmax
ДТ, w
График 2 ДоТ
d0 < (d0)cr1 , ДСТ > 0, dw/dz > 0,
v z
Zw 4 4 4 4
zp _ zw max \ 4 V \ \ \ * Wmax
Zc
ДТ, w
ДоТ
График 3
График 4
Рис. 1.
Распределение скорости восходящих потоков и функции перегрева с высотой в зависимости от различных значений дефицита точки росы вблизи поверхности земли.
Fig. 1. The distribution of the ascending flow velocity and the overheating function with height depending on various values of the dew point deficit near the earth's surface.
z
w
Градиент температуры влажного насыщенного воздуха представлен нами в виде [43]:
__d7[_ L dsm
Yma — , — У а ,
uz с„ dz
: Ymac +e(z_zc)'
(32)
где
Ymac У a.
L ds„
cp dz
- влажноадиабатический градиент температуры на уровне конденсации - известная функция температуры и давления
^ < 0 I [21];
L -
&
удельная теплота конденсации;
удельная теплоемкость при постоянном давлении; массовая доля насыщенного водяного пара;
L I I -7о 2
e=cdA'|«3-10 С/т - константа;
К=
dX
dz2
_ - некоторая функция от величины 5тс - массовой доли водяного пара на уровне конденсации, также принимается постоянной.
Функция перегрева запишется в виде:
AT (z ) = ACT + Ay mac (z - zc)--( z - zc )2, значение функции перегрева на уровне конденсации;
(33)
где ЛсТ -
АУшас=У_Утас - разность градиентов температуры воздуха невозмущенной атмосферы и поднимающейся частицы на уровне конденсации. На аэрологической диаграмме Дутас определяет угол между кривой состояния и кривой стратификации на уровне конденсации.
Для функции пересыщения получено выражение:
As ( z) = Ac 5 +
f Л
ds b | u m
dz
V z=zc /
K.
( z - zc ) + ^(z - zc )
с 2 v -c , (34)
где
Лс5 = 5тс - 5ес - пересыщение на уровне конденсации
Формула для плотности воздушной частицы запишется в виде
\2"
Pi = Pe
1 -а0 -ai (z - zc ) + a2 (z - zc )
(35)
s
m
cp -
НАукИ о ЗЕмлЕ
Влияние подоблачной конвекции на развитие облачной конвекции Данилова Н.Е., Семенова Ю.А., Смерек Ю.Л., Закинян Р.Г.
где
а0 = аДсТ + рДе5,
а1 = аДТ mac + Р
ь +
а2 = 2 (ае-РК ).
Уровень выравнивания плотностей воздуха в облачной воздушной частице и окружающей атмосфере
а1 + ^а2 + 2 (ае - РК) (аДеГ + РДС 5)
ае-рК
(36)
Уровень выравнивания температур определяется выражением: ЛУтас^(ЛУтас)2+2еДс7
(37)
Тепловая конвекция влажного насыщенного воздуха описывается системой уравнений
ди ди 1 др'
и--— =---,
дх дг ре дх
Ау(г) = Дс5 +
/ л
. (15™
ь+—^
Аг
'су
ди дм . — + — = 0. дх дг
Решая эту систему, для функции тока получим выражение
(38)
(39)
(40)
(41)
(42)
¥
(г-гс) к1( \2 2 31 с)
сое кх,
(43)
г
с у
где
ё
аАУтас+Р
сЬ
с у
(44)
- частота Брента - Вяйсяля для влажного насыщенного воздуха;
ае-р/С
аАУтас+Р аДсГ + р Дс5
г \
Ъ + **
аАУтас+Р
<ь
(45)
су
Для вертикальной скорости получим выражения:
2-2,
2 3 К с>
•зт2Ьс
(46)
Толщина конвективной ячейки («облака») находится из условия м> = 0 и для случая, когда скорость на уровне конденсации равна нулю (^с = 0), равна
16,
(47)
сктт
Так как согласно то из формул (47)
видно, что и в случае облачной конвекции существует критическое значение градиента массовой доли водяного пара окружающей среды, которое способствует развитию облачной конвекции:
Кг =
(к,
т 62
'с У
(48)
Более того, из выражений (45) и (47) видно, что когда Дутас < 0, т. е. атмосфера влажно устойчива и согласно традиционным представлениям конвекция не возможна, то при значениях градиента массовой доли водяного пара, больших критического значения, и в этом случае конвекция возможна.
На рисунке 2 функция тока для подоблачной конвекции рассчитывалась по формуле (11), а для облачной конвекции по формуле (43).
НАукИ о ЗЕмлЕ
Влияние подоблачной конвекции на развитие облачной конвекции Данилова Н.Е., Семенова Ю.А., Смерек Ю.Л., Закинян Р.Г.
Рис. 2.
Взаимодействие подоблачной и облачной конвекции. На рисунке по осям x и z отложены метры.
Fig. 2. Interaction of under-cloud and cloud convection. In the figure along the axes x and z laid off meters.
Из рисунка видно, что увеличение влажности в приземном слое приводит к развитию облачной конвекции.
Таким образом, из полученных выражений для параметров облачной конвекции видно, что она зависит от значений функции перегрева ДсТ и скорости восходящих потоков wc на уровне конденсации, которые сами зависят от характеристик приземного слоя атмосферы.
Выводы
На основании полученных в работах [42, 43] аналитических решений уравнений конвекции сухого и влажного насыщенного воздуха, а также аналитического решения конвекции влажного ненасыщенного воздуха, показано влияние на развитие облачной конвекции параметров приземного слоя атмосферы. Получены аналитические выражения для дефицитов точки росы в приземном слое, определяющие значения функции перегрева и скорости восходящих потоков на уровне конденсации. Показано, что на развитие подоблачной конвекции существенное влияние оказывает значение массовой доли водяного пара. Установлено, что при значениях градиента массовой доли водяного пара, приближающихся к критическому значению, наблюдается «взрывной» рост конвективных ячеек в подоблачном слое. Также показа-
но, что существует критическое значение градиента массовой доли водяного, способствующее развитию облачной конвекции.
Таким образом, основной результат, следуемый из предложенной теоретической модели, о влиянии влажности приземного слоя атмосферы на развитие конвекции, находится в согласии с многочисленными экспериментальными и теоретическими исследованиями, обзор которых был проведен во введении. Отличительной особенностью нашей модели является то, что, во-первых, мы говорим не просто об абсолютном значении влажности приземного слоя, а о критическом значении градиента массовой доли водяного пара. Кроме того, для анализа мы находили решение двумерной нелинейной модели конвекции. В качестве недостатка нашей модели отметим ее стационарность и отсутствие вязкости.
Библиографический список
1. Тихонов А.Н., Самарский А.А. Уравнения математической физики. М.: Изд-во МГУ. 6-е издание. 1999. 799 с.
2. Alekseev V.V., Gusev A.M. Free convection in geophysical processes. Sov. Phys. Uspekhi 1983, 26, 906-922.
3. Arakawa A. The cumulus parameterization problem: Past, present, and future, J. Climate, 17, 2493-2525, 2004.
4. Bellon G., and Stevens B., 2005: On bulk models of shallow cumulusconvection. J. Atmos. Sci., 62, 3286-3302.
5. Bennett L.J., Browing K.A., Blyth A.M., Parker D.J., and Clark P.A., 2006: A reviewof the initiation of precipitating convection in the United Kingdom. Q.J.R. Meteorol. Soc., 132, 1001-1020.
6. Betts A.K., Ball J.H., Beljaars A. C.M., Miller M.J., and Viterbo P.A.: The land surface-atmosphere interaction: A review based on observational and global modeling perspectives, J. Geophys. Res.-Atmos., 101, 7209-7225, 1996.
7. Bluestein H.B. Severe Convective Storms and Tornadoes. -Chichester, UK: Springer, 2013. 481 p.
8. Bony S., and Dufresne J.L., 2005: Marine boundary layer clouds at the heart of tropical cloud feedback uncertainties in climate models. Geophys. Res. Let., 32, L20806.
9. Bretherton C.S., and Park S., 2008: A new bulk shallow-cumulusmodel and implications for penetrative entrainment feedbackon updraft buoyancy. J. Atmos. Sci., 65, 2174-2193.
10. Brooks H.E., Doswell C.A., and Wicker L.J., 1993: STORM-TIPE: A forecastingexperiment using a three-dimensional cloud model. Wea. Forecasting, 8, 352-362.
11. Browning K., Blyth A., Clark P., Corsmeier U., Morcrette C, Agner J., Bamber D., Barthlott C., Bennett L., Beswick K., Bitter M., Bozier K.,. Brooks B,. Collier C, Cook C., Davies F., Deny B., Engelhardt M., Feuerle T., Forbes R., Gaffard C., Gray M., Hanken R., Hewison T., Huckle R., Kalthoff N., Khodayar S., Kholer M., Kraut S., Kunz
НАУКИ О ЗЕМЛЕ
Влияние подоблачной конвекции на развитие облачной конвекции Данилова Н.Е., Семенова Ю.А., Смерек Ю.Л., Закинян Р.Г.
M., Ladd D., Lenfant J., Marsham J., McGregor J., Nicol J., Norton E., Parker D., Perry F., Ramatschi M., Roberts H., Russel A., Schulz H., Slack E., Vauhan G., Waight J., Watson R., Webb A., Wieser A., and Zinz K., 2007: The Convective Storm Initiation Project. Bull. Am. Meteorol. Soc., 1-17.
12. Crook N. A. Sensitivity of Moist Convection Forced by Boundary Layer Processes to Low-Level Thermodynamic Fields. Monthly Weather Review. Volume 124, 1996, p. 1767-1785.
13. Deardorff J.W., 1979: Prediction of convective mixed-layer entrainment forrealistic capping inversion structure. J. Atmos. Sci., 36, 424-436.
References
1. Tikhonov A.N., Samara A.A. Equations of mathematical physics. M.: MGU Publishing House. 6th edition. 1999. 799 p.
2. Alekseev V.V., Gusev A.M. Free convection in geophysical processes. Sov. Phys. Uspekhi 1983, 26, 906-922.
3. Arakawa A. The cumulus parameterization problem: Past, present, and future, J. Climate, 17, 2493-2525, 2004.
4. Bellon G., and Stevens B., 2005: On bulk models of shallow cumu-lusconvection. J. Atmos. Sci., 62, 3286-3302.
5. Bennett L.J., Browing K.A., Blyth A. M., Parker D.J., and Clark P.A., 2006: A reviewof the initiation of precipitating convection in the United Kingdom. Q.J.R. Meteorol. Soc., 132, 1001-1020.
6. Betts A.K., Ball J.H., Beljaars A.C.M., Miller M.J., and Viterbo P.A.: The land surface-atmosphere interaction: A review based on observational and global modeling perspectives, J. Geophys. Res.-Atmos., 101, 7209-7225, 1996.
7. Bluestein H.B. Severe Convective Storms and Tornadoes. Chichester, UK: Springer, 2013. 481 p.
8. Bony S., and Dufresne J.L., 2005: Marine boundary layer clouds at the heart of tropical cloud feedback uncertainties in climate models. Geophys. Res. Let., 32, L20806.
9. Bretherton C.S., and Park S., 2008: A new bulk shallow-cumulus-model and implications for penetrative entrainment feedbackon up-draft buoyancy. J. Atmos. Sci., 65, 2174-2193.
10. Brooks H.E., Doswell C.A., and Wicker L.J., 1993: STORM-TIPE: A forecastingexperiment using a three-dimensional cloud model. Wea. Forecasting, 8, 352-362.
11. Browning K., Blyth A., Clark P., Corsmeier U., Morcrette C, Agner J., Bamber D., Barthlott C., Bennett L., Beswick K., Bitter M., Bozier K.,. Brooks B,. Collier C, Cook C., Davies F., Deny B., Engelhardt M., Feuerle T., Forbes R., Gaffard C., Gray M., Hanken R., Hewison T., Huckle R., Kalthoff N., Khodayar S., Kholer M., Kraut S., Kunz M., Ladd D., Lenfant J., Marsham J., McGregor J., Nicol J., Norton E., Parker D., Perry F., Ramatschi M., Roberts H., Russel A., Schulz H., Slack E., Vauhan G., Waight J., Watson R., Webb A., Wieser A., and Zinz K., 2007: The Convective Storm Initiation Project. Bull. Am. Meteorol. Soc., 1-17.
12. Crook N. A. Sensitivity of Moist Convection Forced by Boundary Layer Processes to Low-Level Thermodynamic Fields. Monthly Weather Review. Volume 124, 1996, p. 1767-1785.
13. Deardorff J.W., 1979: Prediction of convective mixed-layer entrapment forrealistic capping inversion structure. J. Atmos. Sci., 36, 424-436.
об авторах
Данилова Нина Евгеньевна, аспиранта кафедры общей и теоретической физики Северо-Кавказского федерального университета. Тел. 8-918-769-22-07 E-mail: ninochka841@mail.ru Семенова Юлия Александровна, аспиранта кафедры общей и теоретической физики Северо-Кавказского федерального университета. Тел. 8-918-753-65-56 E-mail: brilliance_wave@mail.ru Смерек Юлия Леонтьевна, кандидат физико-математических наук, доцент кафедры общей и теоретической физики Северо-Кавказского федерального университета. Тел. 8(918)758-71-29 E-mail: smerek@mail.ru, Scopus ID: 6503952686 ResearcherlD is: Q-5159-2018 Закинян Роберт Гургенович, доктор физико-математических наук, профессор кафедры общей и теоретической физики Северо-Кавказского федерального университета. Тел. 8(918)77-88-675 E-mail: zakinyan@mail.ru
About the authors
Danilova Nina E., Postgraduate Student, Department of General and Theoretical Physics, North-Caucasian Federal University, E-mail: ninochka841@mail.ru Semenova Yulia Alexandrovna, Postgraduate Student, Department of General and Theoretical Physics, North-Caucasian Federal University, Email: brilliance_wave@mail.ru Smerek Yulia Leontevna, Candidate of Physical and Mathematical Sciences, Associate Professor, Chair of General and Theoretical Physics, North Caucasus Federal University Phone 8 (918) 758-71-29, E-mail: smerek@mail.ru. Scopus ID:6503952686 ResearcherlD is: Q-5159-2018 zakinyan Robert Gurgenovich, The doctor of physical and mathematical sciences, Professor of the Department of General and Theoretical Physics IMNS NCFU
Telephone 8(918) 77-88-675 E-mail: zakinyan@mail.ru Scopus ID: 7801621547
