книги из ГПНТБ / Лушев, Ю. Г. Физика верхней атмосферы Земли учебник
.pdfГЛАВА VI
МОДЕЛИ РАСПРЕДЕЛЕНИЯ С ВЫСОТОЙ ПАРАМЕТРОВ СОСТОЯНИЯ ВЕРХНЕЙ АТМОСФЕРЫ
§ 1. ЗА М ЕЧ А Н И Е О БАРОМ ЕТРИ ЧЕСКО Й Ф О РМ УЛ Е ВЫ СО К И Х СЛ О ЕВ АТМ ОСФ ЕРЫ
Неоднократно ставился вопрос о пределах применимости ос новного уравнения статики и барометрической формулы в высо ких слоях атмосферы. Совершенно очевидно, что это уравнение— приближенное. Однако хорошо известно, что отбрасываемые в третьем уравнении движения члены в большинстве случаев (ис ключение составляют области с большими вертикальными ускоре ниями) ничтожно малы по сравнению с силой тяжести и верти кальным градиентом давления. В высоких слоях атмосферы (на высотах более 100—120 км) этот вопрос не столь очевиден, как в более низких слоях. Все дело в том, что состояние высоких слоев атмосферы в значительной степени определяется процесса ми, обусловленными солнечной и космической радиацией (хими ческие превращения молекул, диссоциации, рекомбинации, воз буждения и т. п.). На больших высотах атмосфера испытывает значительные суточные колебания. Если время между соударе ниями частиц соизмеримо с сутками, то атмосфера может стать неравновесной (оценки показывают, что такая неравновесность может возникнуть лишь на высотах, значительно превышающих 800 км). Более существенным источником неравновесности слу жат поглощаемая на больших высотах радиация и обусловлен ные ею фотохимические реакции. Попытка оценить влияние этих процессов на вид барометрической формулы предпринята А. И. Ивановским [10, с. 103—109].
Под влиянием радиации в газе возникают добавочные напря жения Рц, так что общее изотропное напряжение Р равно сумме:
Р = р + Ч '
9* |
131 |
где р — давление газа, связанное с температурой Т уравнением состояния
Р = |
k |
m Рт. |
Это уравнение определяет температуру Т, за которой сохраняет ся свойство характеризовать лишь состояние газа вне зависимо сти от характеристик химических реакций.
Уравнение статики при наличии химических реакций (когда Рп Ф 0) принимает вид
- ^ (P + P z z ) = gP • |
|
|
|
|
Введем эффективную температуру Т * = Т |
+ |
PPzz = |
дг*. |
|
Тогда уравнение статики перепишется в виде |
|
R*P |
|
|
|
|
|
||
д (R* |
_ |
|
|
|
|
р* |
|
|
|
Интеграл этого уравнения |
|
|
|
|
Т* |
Pgdz |
\ |
|
|
р(2) =Ро т^ехр |
R*T* (г) I ’ |
|
|
о
где ро, Т* — значения плотности и эффективной температуры на
исходном уровне.
Последняя формула совпадает с обычной формулой для плот ности воздуха, но с заменой температуры Т на эффективную тем пературу Т*. Приращение ДТ* пропорционально (Д Г *~Q) коли
честву солнечной радиации, поглощаемой |
за |
единицу |
времени |
|||
* |
|
„ |
ДГ* |
3 |
Q |
|
в единичном объеме воздуха. |
Отношение |
~~ф~ = |
^—- гДе |
|||
Qcm— энергия, |
размешиваемая |
в том же |
|
* |
~ |
Wcm |
единичном |
объеме |
|||||
упругими столкновениями молекул или атомов |
( Q cm — kTN cm, |
|||||
N cm— число столкновений, k —постоянная |
Больцмана). |
|
||||
Таким образом, барометрическая формула нарушается тогда, |
||||||
когда Q становится соизмеримой с Q cm. Оценка, |
выполненная по |
|||||
ДТ * |
10-2, показала, |
что обычная |
барометрическая |
|||
критерию - = - = |
формула обеспечивает необходимую точность расчета до высоты около 450 км. Правда, при этом не принято во внимание погло щение корпускулярного излучения Солнца, которое может играть заметную роль в полярных районах.
132
§ 2. ТЕО РЕТ И ЧЕСКИ Е ПРЕД ПОСЫ ЛКИ , |
П О Л О Ж ЕН Н Ы Е В О СН О ВУ |
М ЕЖ Д УН АРО Д Н О Й СПРАВО ЧН ОЙ |
АТМ ОСФ ЕРЫ (C IR A -65) |
Анализ первых наблюдений за торможением ИСЗ, позволив ших обнаружить колебания плотности на высотах более 200 км, привлек исследователей к объяснению этого явления. Наиболее ранние работы по этой проблеме выполнил М. Николе1). Ему принадлежит заслуга выяснения физического механизма влия ния различных факторов на плотность воздуха и построения про стейших решений уравнений. Более общие уравнения проанали зировали, привлекая численные методы и ЭВМ, И. Харрис и В. Пристер23*). Работа этих авторов положена в основу расчета Международной справочной атмосферы (CIRA-65).
Исходным уравнением служит уравнение притока тепла
dT |
р |
dp |
|
(2.1) |
|
'v dt ~ |
7 |
I t |
~ 6 • |
||
|
Привлекая уравнение неразрывности
1 |
dp _ (ди |
dv |
dw |
р |
dt ~ |
ду |
дг |
и уравнение состояния р = RpT, приведем уравнение (2.1) к виду
дТ |
, |
дТ , |
дТ |
дТ\ |
RT (ди |
dv |
d w \ |
е |
|
Tt |
+ |
“ dx + |
V ду |
W дг ) |
(dx |
ду |
дг ) — cv ' |
^ |
Поскольку рассматривается движение верхней атмосферы, осредненное по горизонтали, то все элементы зависят только от времени и высоты. Представим температуру в виде суммы
T(z, t) = T(z) + Ti(z, t), где T(z) — средняя (на высоте z) тем пература, Г1 — отклонение температуры от средней, которое ма
ло по сравнению с T(t). С учетом этих двух предположений урав нение (2.2) принимает вид
дТг , дТ |
%— 1 dw |
е |
~W+ wTz |
а дг |
(2.3) |
С„ |
где х = —с р И а —
1
Cm
Вертикальная скорость w определена из следующего усло вия. В процессе сжатия (или растяжения) атмосферы воздушная
1) |
М . |
N i с*о l e t . |
S tructure of |
th e |
T h erm o sp h ere . P lan . |
a. S p a c e |
Sc., |
5’ |
|||
1— 32, |
1961. |
|
|
|
|
|
|
|
|
|
|
3) |
1. |
H a r r i s |
and |
W. |
P r i e s t e r . T im e -d e p e n d e n t |
stru ctu re |
of |
the |
|||
u p p er |
atm o sp h ere . |
J . |
A tm . |
Sc., |
19, |
Nu 4, |
286— 301, 1962. |
|
|
|
133
частица, находившаяся на высоте z в момент t, оказывается на высоте z + w Ы в момент t + Ы. Поскольку масса столба возду ха, расположенного выше этой границы, при этом не изменяется, то справедливо равенство
p ( z + w M , t + M ) = p { z , t ) . |
(2.4) |
Воспользовавшись барометрической формулой, последнее ра венство перепишем в виде
|
z + w M |
gdz' |
' |
Ро ехр н |
gdz' |
||
RT (/ -f- At) |
A) exp hi RT(t) |
(2-5) |
где ро — давление на уровне моря.
Представляя |
Т (z, t+ M ) |
|
|
дТ |
|
и ^ |
|
1 |
|||
= Т (z, t) + - ^ М |
|
: |
|||||||||
1 |
|
1 |
дТ .Л |
и ограничиваясь |
в |
соотношении |
|||||
~ T { z , *)(* |
T ( z , t ) d t |
Дг I |
|||||||||
|
|
|
|
(относительно М), |
|||||||
(2.5) малыми |
величинами первого |
порядка |
|||||||||
приведем его |
к |
виду |
|
|
|
|
|
|
|
|
|
“"г<г'' > Jг ^ |
|
|
|
|
z |
1 |
дТ |
|
|
||
|
|
|
|
|
dz'. |
( 2.6) |
|||||
Т |
) |
iz' ~ т (г) T 2(z) dt |
|||||||||
|
|
о |
|
|
|
|
|
|
|
|
|
Отсюда |
|
|
|
|
|
|
|
|
|
|
|
dw |
дТ |
f |
1 |
дТ |
|
1 |
дТ |
|
(2.7) |
||
дг |
dz |
J f |
2(z) |
dt aZ + |
f ( z ) |
dt |
|
|
|||
|
|
|
|||||||||
Вставляя выражения |
(2.6) и |
(2.7) |
в уравнение |
(2.3), |
получаем |
||||||
|
дТ |
|
|
|
1 |
dT dz1 |
|
|
|
( 2.8) |
|
|
dt |
|
|
о Тг (г) dt |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|||
|
|
|
dTx |
dT |
|
|
|
|
|
|
|
при этом учтено, что -щ- = |
. |
|
|
|
|
|
|
||||
Вошедший в это уравнение приток |
тепла |
е представляется |
|||||||||
в виде суммы: |
|
|
|
|
|
|
|
|
|
|
|
|
|
е = |
+ |
е уф + |
е О + |
> |
|
|
|
(2.9) |
|
|
|
|
|
|
|
где еЛ, еуф, е0, е' — притоки тепла, обусловленные соответствен
но молекулярной теплопроводностью, поглощением ультрафиоле товой радиации, излучением атомного кислорода в инфракрас
134
ной области и корпускулярным излучением. Для определения последних привлекаются соотношения, вытекающие из кинетиче ской теории газов и теории переноса излучения в газообразной среде. В частности, формула для притока тепла под влиянием молекулярной теплопроводности имеет вид
|
|
|
|
^ = T zk { T ) T z ’ |
|
|
|
<2Л0> |
||||
|
|
Е |
А(П( (Z) ± |
|
|
|
|
|
|
|
||
где |
к (Т ) = |
— —— Т 2 (г) — коэффициент |
температуропро- |
|||||||||
|
|
2j яд *) |
|
|
|
|
|
|
|
|
||
водности |
(О |
|
|
п{ — число |
молекул или атомов »-го |
|||||||
воздуха. Здесь |
||||||||||||
газа в единице объема; |
А {— постоянные, для которых |
приняты |
||||||||||
значения: |
4(H ) = |
2,1 • 103, |
А (Не) = |
9,0-102. |
Л (0 ) = |
3,6-Юг и |
||||||
A (0 2N2) = 1,8.10». |
|
|
(2.8) |
принимает окончательный |
||||||||
|
С учетом (2.10) |
уравнение |
||||||||||
вид: |
|
|
|
|
|
|
|
|
|
|
|
|
|
„ |
дТ _ |
д , 1ТлдТ |
|
я = . . дТ |
|
дТ dz' + |
|
||||
|
Cpf> dt ~ |
d z k ^ |
*•>дг |
|
Ср*Т ^ |
дг |
Jрт ц (z) dt |
|
|
|||
|
|
|
|
+ Р ( е уф + е о + О • |
|
|
( 2. 11) |
|||||
|
Решение этого уравнения найдено с помощью численных ме |
|||||||||||
тодов, а |
расчет |
выполнен |
на ЭВ$Л |
(при |
временном |
шаге |
||||||
Дt |
= 0,25 ч и пространственном Az= 1 км). |
задавались |
на вы |
|||||||||
|
Граничные условия для уравнения (2.11) |
соте 120 км, где температура и плотность, по предположению, не изменяются во времени, и на высоте 1000 км, где обращается в нуль вертикальный градиент температуры. В качестве профи
лей средней температуры Т (г) испытано несколько моделей. Од нако расчет показал, что по истечении 4—5 суток влияние про
филя Т (г) практически перестает сказываться. После того, как на основе решения уравнения (2.11) определена температура Т (г, t), расчет других параметров выполняется с помощью сле дующих соотношений:
а) давление p t (z, t) i-го газа — по барометрической формуле
Pi (z, о = |
Z |
|
|
Pi (Z0) exp | — j* |
dz' ’ |
(2.12) |
|
|
Zo |
|
|
где zQ= 120 k m , f t |
— молекулярный вес f-го |
газа; |
|
135
б) |
плотность ?i (г, |
t) и число |
молекул |
или |
атомов tit (г, t) |
||||
в единице объема — по уравнению |
состояния |
|
|
|
|||||
|
Pi (г, 0 = |
V-iPi (г, |
t) . |
nt {z,.t) |
P i ( 2 , t ) |
’ |
(2.13) |
||
|
R*T(z,t) ’ |
kT(z, T) |
|||||||
где R* — универсальная |
Тазовая |
постоянная, k — постоянная |
|||||||
Больцмана; |
|
|
и плотность |
p(z, t) |
воздуха — по |
||||
в) |
общее давление p(z, t) |
||||||||
закону Дальтона |
|
|
|
|
|
|
|
|
|
|
|
Р = Е Р / 5 |
P = 2 pi- |
|
|
(2Л4) |
|||
|
|
|
( О |
|
( О |
|
|
|
|
Формула (2.12) предполагает, что в пределах всего слоя, для которого строится модель (от 120 до 1000 км), справедливо ос новное уравнение статики, а состав воздуха изменяется с высо той в согласии с гипотезой диффузного разделения газов.
Первые расчеты Харрис и Пристер выполнили в предположе нии, что колебания температуры и плотности обусловлены лишь поглощением ультрафиолетовой радиации Солнца. При этом ока залось, что на высоте 600 км максимум температуры (1959°К)
иплотности наступает в 17 ч, а минимум ее (849 °К) — в 6 ч местного времени. По данным же торможения ИСЗ максимум температуры и плотности наблюдается в 14 ч.
Для согласования результатов расчета с опытными данными
ибыл введен поток тепла е', обусловленный поглощением кор пускулярного излучения Солнца.
Природа этого потока не вполне ясна. Однако его удается связать с излучением Солнца в диапазоне ультракоротких радио
волн (длина волны 10,7 и 20 см). Приток е' по величине сравним
с еуФ‘
Результаты расчета основных физических параметров атмо сферы по изложенной методике сведены в специальные таблицы, которые носят название «Международной справочной атмосфе ры (CIRA-65)». Последующие исследования показали, что мо дель CIRA-65 находится в согласии с опытными данными в широт ной зоне от 40° с. ш. до 40° ю. ш. В более высоких широтах поток корпускулярного излучения превышает те величины, которые бы ли приняты при разработке CIRA-65. Некоторые улучшения в модель CIRA-65 внесены исследованиями Смола, Джекобса, Вриеса и др .!).
Недостаток модели CIRA-65 состоит в том, что высотно-вре менные колебания параметров представлены в ней в виде таб
*) R. Мое. A review of atmospheric models in the altitude range 100 to
1000 km. A1AA Paper, No 69—50, 1969.
136
лиц и графиков. Использование их при расчете и прогнозе орбит ИСЗ на ЭВМ затруднено. В этом смысле более предпочтительны формулы.
Такие формулы, описывающие зависимость плотности и дру гих параметров от высоты и времени, за последние 10 лет пред ложены Л. Яккия, В. Пристером и X. Мартином, И. М. Шваревым, Б. В. Кугаенко, А. С. Суворовым и др. Модели строения атмосферы, в которых учтена зависимость параметров состояния от времени, широты и др., носят название динамических (в от личие от стандартной — статической — модели).
ГЛАВА VII
СЕРЕБРИСТЫЕ ОБЛАКА И ВОДЯНОЙ ПАР В ВЕРХНИХ СЛОЯХ АТМОСФЕРЫ
§ 1. ОСНОВНЫЕ РЕЗУЛЬТАТЫ НАБЛЮДЕНИИ ЗА СЕРЕБРИСТЫМИ ОБЛАКАМИ
Со времени первого наблюдения серебристых облаков1), ко торое провел в ночь с 12 на 13 июня 1885 г. астроном Московской обсерватории В. К- Цераский, многие исследователи обращались к измерению высот этих облаков. Уже 26 июня 1885 г., спустя всего две недели после открытия серебристых облаков, тот же Цераский и известный русский астроном А. А. Белопольский вы полнили первые базисные определения высот этих облаков и на шли, что средняя высота их (из четырех измерений) составляет
79 |
км. |
Все последующие |
измерения (О. Иессе, К. Штермер, |
М. |
И. |
Буров, Г. Витт и др.) |
позволили сделать заключение, что |
средняя высота серебристых облаков достаточно постоянна — по всем данным (включающим 1578 измерений) она составляет 82,4 км с отклонениями (у разных исследователей) не более 1 км в ту и другую сторону. Однако индивидуальные высоты (Нс) се ребристых облаков могут существенно отличаться от среднего значения — диапазон колебаний охватывает высоты от 73 до 95 км (и тоже практически одинаков по данным многих авто ров). В некоторые годы диапазон колебаний Нс может сильно сужаться. Так, в 1958 г. все наблюденные значения Нс были за ключены в интервале 81—85 км.
Систематические наблюдения за серебристыми облаками в Советском Союзе проводятся с середины 30-х годов, когда уси лиями ряда опытных наблюдателей (О. В. Деминева, И. С. Аста
повича, |
Е. Л. |
Кринова, Г. О. Затейщикова, В. А. |
Бронштэна |
и др.) |
была |
создана сеть постоянно действующих |
пунктов для |
регулярного патрулирования этих облаков. Новый этап в изуче нии серебристых облаков наметился в 1950 г., когда под руковод-)*
*) В последние годы все более широкое распространение получает новое название — мезосферные облака.
138
ством И. А. Хвостикова была разработана программа комплекс ных исследований их, включавшая проведение полевых работ с широким использованием измерительной и регистрирующей аппаратуры, лабораторные эксперименты и обработку наблюде
ний. Большой |
вклад в осуществление этой программы внесли |
Н. И. Гришин, |
В. В. Шаронов, М. А. Дирикис, Ч. И. Виллман |
и др. Много новых пунктов по наблюдению серебристых облаков было организовано в период Международного геофизического года (1957—1958 гг.) и последующие годы.
Серебристые облака до настоящего времени удавалось на блюдать в широтных зонах от 50 до 80° и только в теплую поло вину года (с апреля по октябрь). Максимум частоты появле ний этих облаков в среднем наступает через 10—20 дней после летнего солнцестояния. Наблюдения за 1885—1955 гг. показы вают, что 51% облаков появляется в июле и 94% — в летние (июнь—август) месяцы. С увеличением широты максимум числа появлений облаков смещается на более позднее время — в сред нем он приходится на 8 июля в широтной зоне 50—55°, на 18 ию ля в зоне 55—60° и на 28 июля в зоне 60—65°. Наблюдается суточный ход появления серебристых облаков—они бывают ярче и появляются чаще (в 1,5—2 раза) в утренние сумерки, чем в ве черние.
Серебристые облака можно наблюдать только в период суме рек. Однако в период так называемых гражданских сумерек (когда Солнце находится под горизонтом наблюдателя на глу бине, не превышающей 6°, и солнечные лучи освещают всю тро посферу и нижнюю поверхность тропосферных облаков) интен сивность рассеянного света столь значительна и богата разнооб разием ярких красок, что никаких облаков в стратосфере и мезо сфере обнаружить не удается. Наиболее благоприятные условия для обнаружения серебристых облаков возникают в период аст рономических сумерек (когда Солнце находится на глубине от 6 до 12° под горизонтом наблюдателя). В это время земная тень закрывает нижние наиболее запыленные слои атмосферы, а сол нечные лучи освещают только верхние, начиная с мезосферы, слои, проходя под углом около 10° к нижней поверхности сереб ристых облаков. Рассеянный в мезосфере солнечный свет, посту пающий к наблюдателю, и позволяет обнаружить свечение обла ков, цвет которых различные наблюдатели определяют как жемчужно-серебристый с голубоватым отливом или бело-голу
бой.
Внешний вид серебристых облаков имеет некоторое сходство с перистыми, перисто-кучевыми и высококучевыми облаками. В них также видны узкие длинные полосы, волокна и более ши рокие полосы, волны (сходные внешне с рябью на поверхности воды при слабом ветре), округлые темные «кратеры», «флер» в виде туманной клочковатой пелены, местами прикрывающий
139
вышележащие волнистообразные облака. В период наибольшего развития серебристые облака кажутся достаточно плотными. Вместе с тем даже в самых плотных и ярких местах эти облака довольно прозрачны — сквозь них хорошо видны звезды и Луна.
Можно указать следующие признаки, позволяющие отличать серебристые облака от тропосферных облаков верхнего и сред него ярусов: 1) серебристые облака с переходом к гражданским сумеркам и особенно на восходе Солнца исчезают, в то время как облака верхнего яруса в этих условиях становятся видимы луч ше; 2) серебристые облака всегда бывают светлее неба даже тогда, когда они видны на фоне наиболее ярких, окрашенных в желто-оранжевые тона частей зоревого сегмента. Тропосферные же облака выглядят светлее неба только на фоне менее ярких, серовато-голубых зон сумеречного неба, при этом их освещен ность создается светом ярких зон зари; 3) если часть облачной системы, которая проектируется на фоне ярких зон зари, выгля дит темной или, если облака имеют одинаковый вид в области зоревого сегмента и в других частях неба, то это — тропосферные облака.
Т а б л и ц а 7.1
|
|
Частота появления серебристых облаков |
|
|||
Год, месяц |
Тс н |
tc ч |
ч °/о |
Тс —тс ч |
м % |
|
1957 |
г. |
|
|
|
|
|
Май |
|
465,2 |
4,7 |
1,0 |
290,45 |
1,6 |
Июнь |
|
572,6 |
10,4 |
1,8 |
433,31 |
2,4 |
Июль |
|
1579,4 |
41,2 |
2,6 |
1217,54 |
3,4 |
Август |
|
1828,5 |
23,4 |
1,3 |
1201,42 |
1,9 |
1958 |
г. |
|
|
|
|
|
Апрель |
|
1061,42 |
0,58 |
0,05 |
602,09 |
0,1 |
Май |
|
1464,81 |
3,72 |
0,25 |
836,22 |
0,4 |
Июнь |
|
1294,06 |
32,25 |
2,5 |
830,27 |
3,9 |
Июль |
|
1728,17 |
29,07 |
1,7 |
1197,87 |
2,4 |
Август |
|
2470,14 |
28,79 |
1,2 |
1451,27 |
2,0 |
Как же часто появляются серебристые облака? В табл. 7.1 приведены сведения о частоте появления этих облаков по данным наблюдений на всех станциях АН СССР за период Международ ного геофизического года. Здесь Тс — суммарная длительность сумеречных периодов суток за данный месяц, tc— суммарная
140