Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Лушев, Ю. Г. Физика верхней атмосферы Земли учебник

.pdf
Скачиваний:
10
Добавлен:
22.10.2023
Размер:
12.92 Mб
Скачать

ГЛАВА 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

Соседние файлы в папке книги из ГПНТБ