|
Характеристики слоев ионосферы |
Таблица 15.1 |
|
|
|
|
|
|
|
|
|
Название |
Высота, |
Ne, |
Число соуда- |
Условия суще- |
слоя |
км |
эл/см3 |
рений ν, 1/с |
ствования |
|
D |
60−80 |
103–104 |
107 |
Днем |
|
E |
100 |
1,5 105–днем |
105 |
Всегда |
|
|
|
5 103–ночью |
|
|
|
F (F1) |
200−230 |
(2–4) 105 |
104 |
Днем, летом |
|
F (F2) |
300−400 |
(4–9) 105 |
103–104 |
Днем и ночью, |
|
|
|
|
|
летом |
|
F2 |
220− 240 |
(0,6–2) 106 |
103–104 |
Днем, зимой |
|
|
300−330 |
2,5 105 |
|
Ночью, зимой |
|
Из данных таблицы следует, что слой D — нижний слой — существует лишь днем и зимой и летом. Ночью свободные электроны в нем полностью рекомбинируют с ионами. Слой Е существует всегда, хотя ночью, с заходом Солнца, электронная концентрация в нем уменьшается. Слой Е является самым устойчивым ионосферным слоем, и процессы в нем хорошо описываются современной теорией. Слой F зимой ведет себя подобно слою Е, но электронная концентрация в нем гораздо больше. Летом слой F распадается на два слоя — F1 и F2, причем слой F1 существует только днем. Необходимо иметь в виду, что число электронов даже в максимумах электронной концентрации гораздо меньше числа нейтральных молекул, концентрация которых в области слоя D составляет 1015 см–3 и уменьшается до 1010 см–3 в слое F2. Этим объясняется такое большое число соударений электронов с нейтральными молекулами за одну секунду (ν).
В заключение описания основных физических свойств ионосферы следует отметить, что гипотеза о существовании ионизированной области в верхних слоях атмосферы высказывалась задолго до изобретения радио. Свое косвенное подтверждение она получила в 20-х годах XX века, когда была установлена радиосвязь на коротких волнах между Америкой и Европой. Теория дифракции не могла это объяснить из-за сравнительно малой длины волны. Однако прямое доказательство существования ионосферы было получено только в 30-х годах после вертикального зондирования ионосферы радиоимпульсом и приема на Земле отраженного сигнала. Для
проведения таких исследований в разных районах Земли было построено много ионосферных станций. Одна из первых в СССР такая станция была построена в г. Томске при Томском государственном университете. Поскольку сначала было неизвестно общее число ионосферных слоев, то им присваивали названия не с первой буквы, а с некоторым запасом. Первым был открыт слой Е, а ниже его оказался лишь один слой D. По данным ионосферных станций можно определить значения максимальной электронной концентрации в каждом слое. Конечно, эти данные были уточнены с началом полетов искусственных спутников Земли.
Существование ионосферы оказывает определяющее влияние на распространение радиоволн всех диапазонов, кроме УКВ, но и в нем ионосфера влияет на космическую связь, связь со спутниками, радиоастрономию и т.п.
15.2. Электрические параметры ионосферы
Наличие в верхних слоях атмосферы кроме нейтральных молекул свободных электронов и ионов определяет электрические параметры ионосферы. В целом ионосфера является электрически нейтральной средой, поскольку число положительных и отрицательных частиц в ней одинаково. В физике такой ионизированный газ называется плазмой. Рассмотрим распространение плоской монохроматической волны в плазме, считая ее однородной.
При прохождении радиоволны через плазму, в данном случае ионосферу, положительные ионы и свободные электроны начинают двигаться упорядоченно в такт с изменением поля волны. Движущиеся в ионосфере заряды можно рассматривать как некоторый конвекционный ток. Плотность конвекционного тока электронов определяется концентрацией электронов в газе Ne и скоростью их
движения V :
где е — заряд электрона.
В этом уравнении неизвестна скорость движения электронов V . Определим ее, используя уравнение движения электрона
m dV |
+ mνV = eE, |
(15.3) |
dt |
|
|
где m — масса электрона; ν — число столкновений электронов с ионами и молекулами в единицу времени; еЕ — сила, действующая на электрон со стороны электрического поля волны.
Поскольку векторы E и V параллельны, знаки векторов в формуле (15.3) опускаем.
Величина mνV характеризует изменение количества движения электронов при соударениях за единицу времени и имеет смысл силы трения.
Решение уравнения (15.3) ищем в виде
где ω — круговая частота электромагнитной волны. В результате подстановки (15.4) в (15.3) получим
V = eE , m(ν + jω)
и выражение для конвекционного тока электронов преобразуется к виду
|
Nee2E |
|
Nee2E |
|
δэл = |
|
= |
|
(ν − jω). |
(15.5) |
m( jω+ ν) |
m(ω2 + ν2 ) |
Очевидно, что конвекционный ток положительных ионов будет определяться таким же выражением, только параметры m, e и ν будут относиться к иону. Поскольку масса иона во много раз больше массы электрона, то величина тока ионов будет во столько же раз меньше тока электронов и им можно пренебречь.
Представим плазму как вакуум, в котором находятся только заряженные частицы. Присутствием нейтральных молекул можно пренебречь, так как даже в тропосфере, где плотность воздуха гораздо больше, показатель преломления отличается от единицы на
миллионные доли. Таким образом, кроме конвекционного тока в ионосфере будет существовать также ток смещения в вакууме
|
|
|
|
|
|
δсм = jωε0E |
|
|
|
|
|
|
(15.6) |
и полный ток составит |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Nee2 |
|
|
|
Nee2ν |
|
|
δ=δ |
см |
+δ |
эл |
= jωε |
0 |
E 1 |
− |
|
|
+ |
|
E. (15.7) |
|
|
|
|
|
|
|
|
|
ε0m(ν |
2 2 |
|
|
|
m(ν |
2 2 |
) |
|
|
|
|
|
|
|
|
|
+ω |
) |
|
+ω |
|
В формуле (15.7) разделены вещественная и мнимая части. Известно [7, 10, 11], что в любом веществе полный ток может
быть представлен суммой токов проводимости и смещения:
δ = δпр + δсм = σE + jωε0εr E. |
(15.8) |
Сравнивая выражения (15.7) и (15.8), получим формулы для проводимости и относительной диэлектрической проницаемости ионосферы:
|
|
|
|
|
|
|
|
|
|
σ = |
|
Nee2ν |
|
; |
|
(15.9) |
|
m( |
ν2 + ω2 ) |
|
|
εr =1− |
|
Nee2 |
. |
(15.10) |
|
ε0m(ν2 |
+ ω2 ) |
|
|
|
|
|
Формулам (15.9) и (15.10) можно придать более удобный вид, если подставить численные значения входящих в них констант (е, m, ε0 ):
σ = 2,82 10−8 |
|
|
N |
ν |
|
|
|
|
|
|
|
e |
|
|
; |
(15.11) |
|
ν |
2 |
|
2 |
|
|
|
|
+ ω |
|
|
|
|
εr =1−3190 |
|
|
Ne |
|
. |
|
(5.12) |
ν2 |
+ ω2 |
|
|
|
|
|
|
Если частота ω много больше числа соударений (ω >> ν), то формулы (15.11) и (15.12) упрощаются:
|
|
|
|
|
|
σ = 2,82 10−2 |
|
Neν |
; |
(15.13) |
|
2 |
|
|
ω |
|
|
εr =1−80,8 |
Ne |
. |
|
(15.14) |
|
|
|
f 2 |
|
|
В формулах (15.13) и (15.14) электронная концентрация выражена в электронах на кубический сантиметр, а частота — в килогерцах.
Из последних соотношений видно, что εr и σ зависят от часто-
ты. Следовательно, ионосфера является диспергирующей средой. Частота, при которой показатель преломления плазмы обращается в нуль, называется плазменной частотой или частотой Ленгмюра:
Соударения электрона с нейтральными молекулами при этом не учитываются.
Используя соотношение для плазменной частоты, запишем выражение для показателя преломления плазмы:
n = εr |
= 1− |
80,8Ne |
. |
(15.16) |
|
|
|
f 2 |
|
Формула (15.16) имеет основополагающее значение при рассмотрении распространения радиоволн в ионосфере, хотя в ней не учтены соударения электрона с нейтральными молекулами. Однако для большинства диапазонов частот, как следует из табл. 15.1, выполняется условие ω >> γ и поэтому учет столкновений вносит несущественные поправки.
Основные особенности распространения радиоволн в реальной ионосфере можно объяснить, полагая в первом приближении, что концентрация электронов в ионосфере, а следовательно, и ее плазменная частота зависят только от высоты над поверхностью Земли и не зависят от широты и долготы места.
Изменение диэлектрической проницаемости в ионосфере существенно отличается от изменения этого параметра в тропосфере. Показатель преломления в тропосфере близок к единице и меняется
с высотой незначительно. В ионосфере, как видно из выражения (15.16), при частотах радиоволны, меньших плазменной частоты (f < f0), относительная диэлектрическая проницаемость εr оказыва-
ется меньше нуля, а показатель преломления — мнимой величиной:
n = − j εr .
Отсюда следует, что при f < f0 электромагнитные колебания в ионосфере не распространяются и затухают при увеличении расстояния r по экспоненциальному закону. Нетрудно заметить, что выражение (15.16) совпадает с формулой для показателя преломления в волноводах без учета в нем потерь, так что плазменной частоте соответствует критическая частота в волноводе. Очевидно, что все особенности полей в волноводах, изучаемые в курсе электродинамики, переносятся на ионосферу. Отличие между ними состоит в том, что показатель преломления волновода характеризует его в целом, а в ионосфере он изменяется по высоте из-за изменения электронной концентрации.
15.3. Отражение и преломление радиоволн в ионосфере
В отличие от однородной плазмы, рассмотренной в предыдущем подразделе, распределение электронной концентрации в реальной ионосфере имеет сложный характер, зависящий от высоты, времени суток и года, земных координат. Таким образом, ионосфера является неоднородной средой, точнее, неоднородной плазмой. Анализ основных особенностей распространения радиоволн в ионосфере проводят для идеализированной модели в виде простого слоя, распределение электронной концентрации в котором приведено на рис. 15.2. В пределах полутолщины слоя В оно может быть представлено параболическим законом [3]
|
|
h − h |
2 |
|
|
Ne (h) = Nem 1 |
− |
m |
|
. |
(15.17) |
|
B |
|
|
|
|
|
|
|
|
|
|
|
|
|
Как и в тропосфере, анализ распространения радиоволн в ионосфере проводится в приближении геометрической оптики. Будем пренебрегать потерями в ионосфере. Будем также считать Землю и ионосферу плоскими, делая в окончательных выводах поправки на их сферичность. Изменение электронной концентрации с высотой приводит согласно соотношению (15.16) к изменению показателя преломления и искривлению траектории, т.е. к рефракции волны. Основные формулы, описывающие рефракцию радиоволн, остаются такими же, как и в тропосфере, но зависимость n(h) выражена гораздо сильнее.
h
hm B
Nem Ne
Рис. 15.2. Распределение электронной концентрации по высоте в простом слое
В качестве исходных соотношений будем использовать (15.16), (14.10) и (14.13), подчеркивая в них зависимость параметров ионосферы от высоты. Тогда показатель преломления на разных высотах h в простом слое будет определяться формулой
n(h) = 1− |
80,8Ne (h) |
. |
(15.18) |
|
|
f 2 |
|
Уравнение траектории волны имеет вид |
|
n(h)sin θ(h) = n0 sin θ0 , |
(15.19) |
где n0 и θ0 — значения показателя преломления и угла падения на нижней границе простого слоя.
Можно принять, что на нижней границе n0 = 1. Радиус кривизны луча на разных высотах h
|
ρ(h) = |
|
|
n(h) |
|
. |
|
(15.20) |
|
− |
∂n sin θ(h) |
|
|
|
|
|
|
∂h |
|
|
|
|
Вычисляя производные |
|
|
|
|
|
|
|
dn = − 40,4 |
dNe |
, |
|
dNe |
= |
2Nem |
(h |
− h) |
|
|
|
dh |
nf 2 dh |
|
dh |
B2 |
m |
|
|
|
|
и подставляя их в равенство (15.20), окончательно для радиуса кривизны луча получим выражение
ρ(h) = |
n(h)3 f 2B2 |
|
. |
(15.21) |
80,8N |
em |
(h |
− h)sin θ |
0 |
|
|
m |
|
|
|
Остановимся ниже на следующих вопросах для модели простого слоя:
–условие отражения вертикально направленной волны;
–условие отражения наклонно направленной волны;
–траектория радиоволн в простом слое;
–время распространения радиоволн в простом слое.
Условие отражения вертикально направленной волны.
В простом слое плазменная частота (15.15) изменяется с высотой вместе с изменением электронной концентрации. Максимальное значение плазменной частоты в ионосферном слое называется его критической частотой. Очевидно, что критическую частоту слоя можно определить как плазменную частоту максимума ионизации
Рассмотрим распространение вертикально направленной волны в простом слое в зависимости от соотношения частоты волны и критической частоты слоя. При этом возможны три случая
(рис. 15.3).
|
|
|
n |
|
|
n |
|
|
|
|
|
0 |
1 |
0 |
1 |
0 |
1 |
Re(n) |
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
б |
|
|
|
в |
|
Рис. 15.3. Зависимость n(h) для вертикально направленной волны:
а – f1 > fкр, б – f2 = fкр, в – f3 < fкр
На частоте f1 > fкр (рис. 15.3,а) показатель преломления нигде
не обращается в нуль (n(h) > 0) и волна проходит через слой, лишь испытывая изменение скорости. На частоте f2 = fкр (рис. 15.3,б)
показатель преломления в точке максимума электронной концентрации равен нулю (n(hm) = 0). При приближении к этой точке скорость движения энергии волны (групповая скорость) замедляется до нуля, распространение волны вверх прекращается и она поворачивает вниз. Будем называть этот процесс отражением волны, хотя, конечно, он отличается от классического отражения волны от резкой границы раздела двух сред.
На частоте f3 < fкр (рис. 15.3,в) показатель преломления обра-
щается в нуль в двух точках по высоте и волна отражается в нижней из них. Между этими точками показатель преломления — мнимая величина. Чем ниже частота волны f, тем ниже оказывается высота первой точки отражения.
Формула (15.22) для критической частоты простого слоя позволяет определить критические частоты реальных ионосферных слоев. Так, согласно табл. 15.1:
для слоя D — Nem =104 эл/см3 , |
fкрD = 80,8 104 = 900 кГц; |
для слоя Е — Nem =104 эл/см3 , |
fкрE = 3,5 ГГц; |
для слоя F1 — Nem =104 |
эл/см3 |
, |
fкрF |
= 5,5 ГГц; |
для слоя F2 — Nem =104 |
эл/см3 |
|
1 |
|
, |
fкрF |
=12 ГГц. |
|
|
|
2 |
|
Условие отражения наклонно направленной волны. Из уравнения луча (15.19) следует, что с продвижением волны вглубь простого слоя при h < hm угол падения θ увеличивается (поскольку n уменьшается), а при h > hm — уменьшается. Об этом же свидетельствует анализ радиуса кривизны луча (см. формулу (15.20)). Из формулы (15.20) следует, что при h < hm радиус кривизны ρ>0, при h > hm ρ<0. Таким образом, ниже hm траектория волны обращена выпуклостью от Земли, а выше hm — выпуклостью к Земле. Если при этом в какой-либо точке угол падения станет равным 90° (sinθ = 1), то в ней произойдет поворот волны к Земле, то есть отражение волны. Очевидно, что отражение волны может произойти там, где ρ>0, т.е. ниже максимума ионизации. Из уравнения траектории волны (15.20) следует, что условием отражения наклонного луча является выполнение равенства
n(hотр) = sin θ0. |
(15.23) |
Таким образом, отражение наклонного луча в ионосфере происходит на высоте, при которой показатель преломления n(h) равен синусу начального угла падения на ионосферу θ0 . Возведем выра-
жение (15.23) в квадрат и, используя равенство (15.18), придадим ему вид
|
|
|
|
|
|
|
|
|
1− |
80,8Ne (hотр) |
= sin2 |
θ0 или |
80,8Ne (hотр) |
= cos2 θ0. |
|
|
|
f 2 |
|
|
|
f 2 |
Величина |
80,8Ne (hотр) = f0 |
(hотр) |
определяет частоту волны, |
отражающуюся |
на высоте |
hотр |
при |
вертикальном падении, а |