Файл: Лушев, Ю. Г. Физика верхней атмосферы Земли учебник.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 111

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Приводим (по данным статьи1) на рис. 5.10 и 5.11 результа­ ты измерения плотности и температуры на о. Хейса с помощью

О

200

400

600

вО

Т°К

Рис. 5.11.

Вертикальные

профили

температуры над

о. Хейса:

1 — измерения с помощью ионизационных маномет­ ров в сентябре 1966 г.; 2 — то же в январе 196? г.; 3— измерения с помощью тепловых манометров в де­ кабре 1967 г.; 4 — измерения с помощью масс-спек­ трометра; 5 — данные малой метеоракеты в сентябре

1966 г.; 6 — по Лори-Чанину (1965)

; 7 — по Хедину

и Ниру

(1966); 8 — по Ллойду и Шеппарду (1965);

9 — по

Стандартной

атмосфере

(1964); 10 — по

CIRA-65 в сентябре;

И — по CIRA-65 в январе

манометров и масс-спектрометров. Для сравнения на том же гра­ фике нанесены данные других авторов. Как видно, температура

') И. Н. И в а н о в а, Г. А. К о к и н, А. Ф. Ч и ж о в. Температура и плот­ ность мезосферы и термосферы Центральной Арктики. — «Метеорология и гид­ рология», 1968, № 5.

9 Зак . 5025

129

мезосферы в Арктике довольно низкая (около 200 °К), минимум ее (мезопауза) наблюдается на высоте около 115 км.

Наблюдения над искусственными светящимися облаками, вы­ полненные в рамках программы франко-советских исследований термосферы высоких широт, позволили определить температуру на высоте 170 км в осенние (пуски ракет 9 и 10 октября 1967 г.) и весенние (пуски 28 февраля, 16, 19 и 20 марта 1968 г.) месяцы. Согласно этим данным ■), температура воздуха на высоте 170 км составляет 900— 1000 °К в осенние месяцы и 400—500 °К в начале весны (при сходных условиях геомагнитной и солнечной активно­ сти). Таков эффект охлаждения термосферы Арктики в течение зимы. Однако если геомагнитная активность повышена (как было 16 марта 1968 г.), то температура и весной может приближаться к 1000 °К- Подмечено также, что по мере погружения Солнца под горизонт температура может довольно значительно понизиться (так, 28 февраля 1968 г. температура упала от 500 до 400 °К, а 16 марта 1968 г. — от 1100 до 950°К при изменении угла погру­ жения Солнца о т —9 д о —10°).

') Ж. Е. Б л а м о н, М. Л. Шанэн , М. М а й я р, С. М. П о л о с к о в, Г. Ф. Т у л и н о в. Измерение температуры термосферы Центральной Арктики с помощью искусственных светящихся облаков. — «Метеорология и гидроло­ гия», 1970, № 11.

130



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