Научная статья на тему 'Численное моделирование параметров полярной ионосферы'

Численное моделирование параметров полярной ионосферы Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
196
52
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ИОНОСФЕРА / ЭЛЕКТРОДИНАМИКА / ЭЛЕКТРОННАЯ ПЛОТНОСТЬ / ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ / ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ / ИТЕРАЦИОННЫЕ МЕТОДЫ

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Лукьянова Рената Юрьевна

В статье дается обзор разработанных трех численных моделей для расчета электродинамических параметров высокоширотной ионосферы Земли. Модель глобального распределения ионосферного электрического потенциала, построенная на основе решения краевой задачи о растекании ионосферных токов, позволяет рассчитать траектории конвекции ионосферной плазмы в северном и южном полушариях. Модель высокоширотной ионосферы позволяет оперативно рассчитывать трехмерную структуру электронной плотности в диапазоне высот 120-500 км при различных гелио-геофизических условиях. Учитывается определяющая роль электрических полей магнитосферного происхождения. Концентрация основных ионосферных ионов определяется решением уравнения фотохимического баланса и конвективно-диффузионного уравнения вдоль траектории конвекции плазменных трубок с учётом параметров термосферы. Разработана методика и алгоритмы численного расчета распределения магнитного поля над ионосферой, которое создается электрическими токами магнитосферного происхождения. Модель базируется на решении уравнения для векторного магнитного потенциала и позволяет рассчитать двумерную картину магнитных вариаций.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Лукьянова Рената Юрьевна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Numerical modeling of the polar ionosphere parameters

An overview of the three numerical models for calculating the electrodynamic parameters of the Earth’s high-latitude ionosphere is presented. The model of the global distribution of the ionospheric electric potential, constructed on the basis of the solution of the boundary value problem on the spreading of ionospheric currents, makes it possible to calculate the trajectories of convection of the ionospheric plasma in the northern and southern hemispheres. The model of the high-latitude ionosphere allows calculating the three-dimensional structure of the electron density in the altitude range of 120-500 km under various helio-geophysical conditions. Importance of the electric fields of magnetospheric origin is stressed. Concentration of the main ionospheric ions is determined by the solution of the photochemical balance equation and the convective-diffusion equation along the trajectory of the plasma tube convection, taking into account the parameters of the thermosphere. A methodology and algorithms for calculating the distribution of magnetic field over the ionosphere, which is created by electric currents of magnetospheric origin, are developed. The model is based on the solution of the equation for the vector magnetic potential and makes it possible to calculate a two-dimensional picture of the magnetic variations.

Текст научной работы на тему «Численное моделирование параметров полярной ионосферы»

ЧЕБЫШЕВСКИИ СБОРНИК

Том 19. Выпуск 4

УДК 51-73:537, 519.63

DOI 10.22405/2226-8383-2018-19-4-91-102

Численное моделирование параметров полярной ионосферы 1

Лукьянова Рената Юрьевна — доктор физико-математических наук, главный научный сотрудник, Геофизический центр РАН. e-mail: [email protected]

В статье дается обзор разработанных трех численных моделей для расчета электродинамических параметров высокоширотной ионосферы Земли. Модель глобального распределения ионосферного электрического потенциала, построенная на основе решения краевой задачи о растекании ионосферных токов, позволяет рассчитать траектории конвекции ионосферной плазмы в северном и южном полушариях. Модель высокоширотной ионосферы позволяет оперативно рассчитывать трехмерную структуру электронной плотности в диапазоне высот 120-500 км при различных гелио-геофизических условиях. Учитывается определяющая роль электрических полей магнитосферного происхождения. Концентрация основных ионосферных ионов определяется решением уравнения фотохимического баланса и конвективно-диффузионного уравнения вдоль траектории конвекции плазменных трубок с учётом параметров термосферы. Разработана методика и алгоритмы численного расчета распределения магнитного поля над ионосферой, которое создается электрическими токами магнитосферного происхождения. Модель базируется на решении уравнения для векторного магнитного потенциала и позволяет рассчитать двумерную картину магнитных вариаций.

Ключевые слова: Ионосфера, Электродинамика, Электронная плотность, Численное моделирование, Дифференциальные уравнения, Итерационные методы

Библиография: 15 названий. Для цитирования:

Р. Ю. Лукьянова. Численное моделирование параметров полярной ионосферы // Чебы-шевский сборник, 2018, т. 19, вып. 4, с. 91-102.

Аннотация

1 Исследование выполнено при поддержке Российского научного фонда (проект 16-17-00121)).

CHEBYSHEVSKII SBORNIK Vol. 19. No. 4

UDC 51-73:537, 519.63 DOI 10.22405/2226-8383-2018-19-4-91-102

Numerical modeling of the polar ionosphere parameters

Lukyanova Renata Yurievna — D.Sc., Principal research scientist, Geophysical Center RAS. e-mail: [email protected]

Abstract

An overview of the three numerical models for calculating the electrodynamic parameters of the Earth's high-latitude ionosphere is presented. The model of the global distribution of the ionospheric electric potential, constructed on the basis of the solution of the boundary value problem on the spreading of ionospheric currents, makes it possible to calculate the trajectories of convection of the ionospheric plasma in the northern and southern hemispheres. The model of the high-latitude ionosphere allows calculating the three-dimensional structure of the electron density in the altitude range of 120-500 km under various helio-geophysical conditions. Importance of the electric fields of magnetospheric origin is stressed. Concentration of the main ionospheric ions is determined by the solution of the photochemical balance equation and the convective-diffusion equation along the trajectory of the plasma tube convection, taking into account the parameters of the thermosphere. A methodology and algorithms for calculating the distribution of magnetic field over the ionosphere, which is created by electric currents of magnetospheric origin, are developed. The model is based on the solution of the equation for the vector magnetic potential and makes it possible to calculate a two-dimensional picture of the magnetic variations.

Keywords:

Bibliography: 15 titles.

For citation:

R. Yu. Lukianova, 2018, "Numerical modeling of the polar ionosphere parameters" , Chebyshevskii sbornik, vol. 19, no. 4, pp. 91-102.

1. Введение

При взаимодействии плазмы солнечного ветра и межпланетного магнитного поля с земным диполем в ионосфере и магнитосфере Земли генерируются мощные электрические поля и токи. Электродинамические процессы, происходящие в ионосферно-магнитосферной системе, оказывают значительное влияние на функционирование современных технологических систем на Земле и в космосе. Они являются источниками спорадических возмущений геомагнитного поля, особенно интенсивных в полярных регионах, в Арктике и Антарктике. Изучение электродинамики верхней атмосферы Земли важно для понимания солнечно-земных связей, а разработка инструментов прогнозирования имеет прикладное значение для определения, например, условий распространения радиоволн в высоких широтах и предупреждения сбоев в системах электроснабжения вследствие индукционных токов.

Последние десятилетия характеризуются интенсивным развитием космических и наземных систем наблюдения за электродинамическими параметрами ионосферы и магнитосферы. Статистический анализ большого количества полученных данных позволяет строить эмпирические модели, описывающие усредненные состояния окружающей среды [1, 2, 3, 4]. С другой стороны, эффективным способом количественного описания глобального распределения параметров является математическое моделирование, основанное на решении уравнений, описывающих физические процессы, происходящие в ионосферно-магнитосферной системе. В отличие от эмпирических математические модели позволяют не только воспроизводить фоновые, достаточно спокойные условия, но и дают возможность проводить численные эксперименты, в которых входные параметры задаются в широких пределах, что важно для прогнозирования экстремальных ситуаций.

В течение нескольких лет в ГЦ РАН проводились работы по развитию семейства численных физико-математических моделей для исследования электродинамики высокоширотной ионосферы. Основными компонентами этого семейства являются (*) модель распределения магнитного поля на верхней границе ионосферы, которое создается электрическими токами магнитосферного происхождения; (*) модель распределения электрического потенциала в ионосфере Земли в глобальной постановке, т.е. при учете электродинамической связи между ионосферными оболочками северного и южного полушарий; (*) модель распределения концентрации электронов в диапазоне высот 120-500 км. В статье дается обзор полученных результатов.

2. Общая схема расчетной области

Общая схема расчетной области, аппроксимирующей ионосферу Земли, представлена на Рис. 1. Ионосфера представляет собой сферическую оболочку конечной толщины, нижняя поверхность которой находится на высоте 120 км от поверхности Земли, а верхняя - на высоте 500 км. Продольные электрические токи, которые генерируются в пограничных слоях магнитосферы (вне расчетной области), втекают и вытекают из ионосферы в высокоширотных областях и растекаются по всей нижней поверхности ионосферной оболочки, проводимость которой неоднородна. С высотой проводимость ионосферы не меняется, на нижней поверхности она задается в виде интегральной по высоте величины. Продольные токи, текущие вдоль силовых линий геомагнитного поля, являются источником электрического потенциала на нижней границе ионосферы. Эти же токи являются источником магнитного потенциала на верхней границе. На рис. 1 нижняя ионосферная оболочка разделена на три области: северная и южная полярные шапки, где силовые линии геомагнитного поля разомкнуты и вытянуты далеко в хвост магнитосферы, и остальная часть сферы, где силовые линии опираются на сопряженные точки в северном и южном полушариях.

:

■с

Рис. 1: Схема расчетной области для ионосферы. Двумя концентрическими окружностями показаны верхняя и нижняя границы ионосферы. Пунктирные .линии обозначают замкнутые в среднеширотной области и разомкнутые в полярных шайках силовые линии магнитного диполя. Толстые серые стрелки в северной и южной полярных шайках обозначают продольные токи, втекающие и вытекающие из ионосферы. Сплошными горизонтальными линиями показаны границы полярных шапок, приэкваториальные границы ереднеширотных областей и магнитный экватор. Буквы, обозначающие границы и заштрихованные области относятся к модели электрического потенциала (см. раздел 4)

3. Моделирование магнитного поля, создаваемого продольными токами над ионосферой

Магнитное ноле, которое создается продольными токами на верхней границе ионосферы, измеряется низколетящими спутниками с полярной орбитой. Обычно это три компоненты ноля: зональная, меридиональная и вертикальная. Спутник может измерять ноле только в точках своей траектории [5]. Для получения общей картины распределения магнитного ноля и верификации спутниковых измерений была разработан алгоритм численного расчета распределения магнитного потенциала. Модель базируется на решении уравнения

где В вектор индукции магнитного ноля, которое создается электрическим током, А векторный магнитный потенциал.

Уравнение 1 записывается в сферических координатах в, <f, г, где в - геомагнитная широта, р - долгота, г - радиус Земли + 120 км, и дискретизируется на сетке с шагом ДО, Др. Продольный ток j, создающий магнитное ноле В, в высоких широтах земного шара направлен практически перпендикулярно к ионосферной оболочке. Тогда из выражения В = rot А((0, (р)ег широтная и долготная компоненты индукции определяются как

Для прямоугольного контура с центром в точке сетки, периметром L, площадью s и плот-

В = rotA при дополнительном условии divA = 0

(1)

ностью тока внутри контура j можно записать:

1дА в _ 1 дА г д9 ' г sin 9 д<р

Bv = -1 дА, Во = дА (2)

1 [

-В1И--Пйз = 0 (3)

Ц- с .)

Ь 8

Компоненты В для узла (¿,.1), находящегося в центре единичного контура с угловыми точками (1234), дискретизируются через компоненты А:

(В„)1 = -

Aíj Aí— 1 j

гАв

А

В)2 = -

(Во )з =

Аг+1 з — Aíj гАв 1 Aíj — Aíj—i г sin вг Аip 1 Aíj+i — Aíj

(Во )з = -i-r íJ+" íJ (4)

sin í А p

Итерационное уравнение для А записывается как

А

(п+1)

= - АА^

í3

(п,п+1)

А

(п,п+1)

í3

=

fíj (КА1 + 1,

А(п+1) А(п)

í3

í3

дА

(п)

fíj (КА1+1,

А

( п+1)

Aí-13 Ап)

(5)

где / - левая часть уравнения , w - параметр релаксации, п - номер итерации. Двумерное распределение А однозначно определяет значения широтной и долготной компонент В в каждой точке расчетной сетки. Задавая произвольное распределение продольных токов j, можно получить картину создаваемого этими токами магнитного поля на высотах верхней границы ионосферы. Для решения задачи был разработан алгоритм, реализованный в программном модуле на языке Fortran. Задача решалась на сетке с шагом АО, Ар =0.5°, покрывающей полярную область от полюса до магнитной широты 60° при различных распределениях входных параметров.

д

з

4. Численная модель глобального распределения электрического потенциала в ионосфере Земли

Трехмерное уравнение непрерывности электрических токов, обеспечивающих магнито-сферно-ионосферное взаимодействие, записывается как

¿тЗ= ] ■ втх (6)

где 3 - интегральные горизонтальные токи, текущие в ионосфере, у - продольные токи, текущие вдоль силовых линий геомагнитного поля и в высоких широтах практически перпен-

х

Согласно постановке задачи, ионосфера разделена на три области (см. рис. 1): северная и южная полярные шапки и средние широты, где силовые линии геомагнитного поля замкнуты

[6]. На рис. 1 границы этих областей обозначены как $2 и Оз, а сами области заштрихованы (из среднеширотных - только северная) и обозначены как а = 1, 2, 3. Индекс 1 соответствует северной полярной шапке с экваториальной границей на в = вц индекс 2 - южной полярной шапке с границей на в = 0% индекс 3 - среднеширотной области с границами на в = в\т в = 9$. Постановка задачи учитывает электродинамическое сопряжение ионосфер противоположных полушарий. Внутри северной и южной полярных шапок процесс растекания горизонтальных ионосферных токов описывается своим уравнением непрерывности, с присущим данной шапке распределением проводимости и источников (продольных токов), в то время как вне области полярных шапок замкнутые магнитные силовые линии земного диполя эффективно выравнивают электрический потенциал в сопряженных точках противоположных полушарий. Выполнение граничных условий, отражающих неразрывность общей токовой цепи и выравнивание потенциала на границах полярных шапок, приводит к взаимозависимости распределения потенциала внутри шапок и влиянию обеих шапок на распределение потенциала в среднеширотной области.

Уравнение (6) в координатах {в, р) решается для северной, южной и среднеширотной областей отдельно, а соответствующие граничные условия связывают эти три области в единую систему. В средних широтах уравнение решается на половине (для определенности, северной) области с низкоширотной границей, несколько отстоящей от экватора, а в качестве интегральной проводимости и источника используются сумма проводимостей и сумма источников в сопряженных точках обоих полушарий.

Краевая задача имеет вид

йтЗ 1 = ^для^ < в1 (7)

йтЗ 2 = --в1 < в <ж (8)

йтЗз = ]здля^ < в < вз (9)

Граничные условия:

и 1(61, ф) = из(в1,ф) = и2(62, ф) (10)

31(в1,р) - -3 з(в1,р) = 3 2(в2, <р) (11)

3 з(вз,<р) = 0, (12)

где иа, Заи ионосферный электрический потенциал, интегральный ионосферный ток и источник возбуждения потенциала в виде продольного тока в соответствующих областях (а = 1, 2, 3). Граничное условие (10) означает отсутствие скачка потенциала поперек границы заданной шапки и между границами противоположных шапок. Условие (11) означает, что возможные разрывы нормальной компоненты тока в плоскости {в, ф) па границах северной и южной шапок компенсируют друг друга за счет токов перетекания из одного полушария в другое. Условие (12) соответствует отсутствию тока через экваториальную границу.

Параметры 1!аи /«связаны законом Ома:

За = Яа(-уиа), (13)

где - тензор интегральной по высоте ионосферной проводимости [7], т.е.

Е = ( ^ А (14)

= ХР/ sin2x, Е^ = Хн/ sin2x, = Sp sinx = 2 ■ cos в / (1 + 3 ■ cos2 0)1/2,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

где Хри Ея™- ионосферные проводимости Педерсена и Холла.

Уравнения (7) - (12) представляют собой краевую задачу эффективным численным методом решения которой является метод конечных разностей. Для перехода к итерационной схеме вводится переменная Za:

Zin+1) = (U(n+1) -uin))!т{п+1)

(15)

где т - релаксационный параметр [8], п - номер итерации, и(и и(п) - решения в соответствующих областях (а = 1, 2, 3), соответственно на итерациях п и (п+1). Из (15) видно, что если итерационный процесс сходится (т.е. п+1"> —> 0), то иа(иа(п) иа, а {/„является решением исходной задачи (7) - (12).

Уравнения (7) - (9) для переменной ^^в итерационном виде имеют вид

di v(-VZ^+1)) + divJin) = j0

(16)

Граничные условия 10 - 12 удовлетворяется на каждой итерации. Условие 10 в терминах (вертикальная линия обозначает границу на соответствующей границе в) записывается как

(17)

(18)

Z 1а+1) в1 Z (п+1) Z2 в2 = (u(n) в2 u(n) /Т(п+1)

Z (n+1) в1 Z (п+1) Z3 в2 = (u(n) в1 u(n) J/ /т(п+1)

Граничные условия (11) и (12) в итерационном виде выражаются как

dZ1

( П+1)

дв

dz:

( п+1)

в1

д

д Z

(п+1)

02

д

J

(п)

в1

J

(п)

91

J

(п)

02

91

д Z

(п+1)

д

J

(п)

в3

в3

(19)

(20)

Далее применяется метод разделения переменных. В левую часть уравнения (16) вводится регуляризатор К, который представляет собой угловую (в, ф) часть оператора Лапласа [9]:

R = div(-V) = -Дв1р

Для переменной Дополучаем

RZ(an+1) = di v(-VZ(an+1)) = ^Zfr+V Д Z(^+1) = divJ(an) - ja

(21)

(22) (23)

Лапласиан допускает разделение переменных, и решение (16) можно представить в виде ряда Фурье

т=+ж

Z^+1)(9,<р) = £ Сат(9) е

тир

(24)

т=—оо

где Сат{&) - комплексные коэффициенты Фурье, зависящие только от 0,т номер узла сетки по долготе.

Подстановка (24) в (16) - (20) приводит к одномерной краевой задаче для Сат(в):

1 д ... дС, (sin0

sin д

д

) - ■ Саш = Fm{ -ja + divJin)}

sin2

Clm\gi С2т\ Q2 — ^

а

(п)

02

— U

(п)

01

-(п+1)

(26)

С\т\д\ Сзт\д2 — Fi

U

(п)

U

(п)

01

01

-(п+1)

(27)

dCim

дв

дС2т

ei

дв

— Fn

J

(п)

02

ei

- J

(п)

02

— J

(п)

01.

дСзт

д

— п

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

03

ни

(28)

(29)

где Fm{-} - зависящий от в коэффициент Фурье от выражения в скобках для т-го узла сетки по долготе, а Сатесть Сат(в).

В конечно-разностном виде уравнение (25) имеет вид

1

А в ■ sin вг

sin0 1

Ci+1 — Сi А

- sin 1 i-/2

Ci — Ci-1 А

1

А в ■ sin в,

■Сг — Ft[ja + divjW

(30)

где А в - шаг по кошироте, г - помер узла сетки по кошироте. Значки а и т далее опущены, так что Сат(в) в г-м узле сетки обозначено как Ci.

Три системы разностных уравнений вида (30) для областей, дополненные разностными аналогами граничных условий (26) - (29), приводят к связности этих систем по значениям в

1 2 3

лить отдельную замкнутую подсистему. После ее решения каждая из трех основных систем разностных уравнений решается методом прогонки.

Изложенные выше алгоритмы были реализованы в виде набора программ на языке Фортран. Расчеты проводились с двойной точностью. Программное средство было протестировано на ряде модельных задач [10, 11]. Сходимость оценивалась по нормам невязок

divjW - ja I divJ^) - j.

(31)

a

a

5. Численная модель распределения электронной концентрации в Е области полярной ионосферы

Модель рассчитывает распределение электронной концентрации (Же) в полярной ионосфере в области выше 50° магнитной широты [12, 13]. Основными факторами, определяющими состояние высокоширотной ионосферы, являются динамика плазмы и процессы образования и химических потерь ионов. Соответственно с этим модель состоит из двух блоков. В первом блоке с целью определить траектории, вдоль которых дрейфует плазменная трубка, рассчитываются геометрические и электродинамические параметры картин конвекции ионосферной плазмы. Скорость движения заряженной частицы определяется суммой скорости коротации и скорости дрейфа в скрещенных электрическом и магнитном полях. Так как па высотах >120 км плазма замагничена, линии конвекции совпадают с изолиниями ионосферного электрического потенциала. В рамках модели распределение электрического потенциала параметризовано и описывается соответствующими аналитическими формулами согласно [14]. Во втором блоке решается задача интегрирования по высоте конвективно-диффузионного уравнения и уравнения фотохимического равновесия в вертикальной плазменной трубке, дрейфующей под действием электрических полей магнитосферной конвекции и коротации.

Уравнения траектории конвекции записываются в геомагнитных сферических координатах в, р. г. Система обыкновенных дифференциальных уравнений, описывающая конвективные траектории, имеет вид:

йв/<И = -2 ■ Е1р(в,ф). (32)

йр/<И = Ев (в, р). (33)

где Ед и Е,р - компоненты ионосферного электрического поля, которое возбуждается продольным током магнитосферного происхождения, М - магнитный момент Земли, ш - частота углового вращения Земли вокруг геомагнитной оси, К = Д#+200 км, Ке — радиус Земли, £ -время.

Конфигурация и плотность продольного тока j аппроксимировались аналитическими функциями в зависимости от параметров солнечного ветра.

з(= 12 [з'т ■ й1п(т (34)

т

3 (р) = 30 + ^ [Зш ■ сон(т<р)]т = 1,.... (35)

т

Базовая картина конвекции задавалась распределением электрического потенциала Ф, симметричным относительно меридиана полдень-полночь:

Ф(в ,р) = -Т (0)Е (<р), (36)

где Е(р)=вт(р-ро) — определяет поворот картины конвекции относительно меридиана

ции распределения электрического поля (Еу__в) при условии равенства падения потенциала в

полярной шапке Фрс и в авроральпой зоне. Аналитическое описание базируется на непрерывном изменении профиля Еу--в, в зависимости от параметров солнечного ветра, геомагнитной активности, сезона.

Ионосферный блок описывает трёхмерное распределение Ие. В точках, расположенных вдоль траекторий конвекции, решается задача интегрирования нестационарного конвективно-диффузионного уравнения для иона кислорода 0+:

В В ( В \

-п{о+) = - (л■ -п(о+) + я■ п(о+)) + д0+ -Ь (37)

где п{0+) - ионная концентрация, £ - время, 2 - высота, Лий - коэффициенты, соответствующие вертикальной скорости 0+, зависящей от скорости новообразования, параметров термосферы, частот соударений, географических координат, значений электрического и магнитного поля, (¡}о+ и I - скорость ионизации и потерь ионов, соответственно. Также интегрируется уравнение фотохимического равновесия для концентрации обобщенного молекулярного иона М+:

Ям + + Ь = а ■ п(М+) ■ (п(М+) + п(0+)), (38)

где п(М)+ - ионная концентрация, а — эффективный коэффициент рекомбинации.

На нижней границе ионосферы задаётся концентрация п{0+) из условия фотохимического равновесия, а на верхней - величина уходящего потока. Потери ионов определяются скоростью рекомбинации в химических реакциях, часть из этих реакций зависит от температуры. Для иона 0+ учитываются следующие реакции:

о+ + о2 ^ о2+ + о

О + 2 ^ Nо+ + N

(39)

В уравнении (38) учтено образование иона М+ под действием излучения (корпускулярного, связанного с высыпаниями энергичных частиц из магнитосферы, а также прямого и рассеянного солнечного УФ) и в реакции потерь 0+.

Скорость ионообразования за счет высыпающихся частиц рассчитывается по эмпирическим статистическим моделям для потоков электронов. Скорость образования ионов за счет солнечной радиации на различных длинах волн определяется по эмпирическим формулам. Для сокращения времени расчета параметры термосферы не моделируются внутри ионосферного блока, а задаются на основе эмпирической модели [15]. Таким образом, в ионосферном блоке при заданных гелиогеофизических условиях, дне года и моменте мирового времени рассчитывается высотное распределение Ие в любом заданном наборе точек на сетке, покрывающей высокоширотную область.

В результате расчетов для заданных входных параметров рассчитывается высотные профили Иев точках расчетной сетки, покрывающей полярную область, а также строятся карты изолиний максимума Ие.

6. Заключение

Разработаны численные модели для расчета распределения электродинамических параметров высокоширотной ионосферы Земли. Модель глобального распределения электрического потенциала, построенная на основе решения краевой задачи о растекании ионосферных токов, позволяет рассчитать траектории конвекции ионосферной плазмы в северном и южном полушариях и выявляет возможность эффективного проникновения электрических полей, генерируемых в полярной шапке продольными токами магнитосферного происхождения, в средние широты и в противоположную полярную шапку за счет электродинамической связи полушарий.

Расчет магнитного потенциала, создаваемого продольными токами над ионосферой, позволяет моделировать общую картину распределения горизонтальных компонент магнитной индукции и может служить для верификации измерений спутниковых магнитометров, которые регистрируют компоненты поля в точках вдоль своей орбиты.

Численная модель высокоширотной ионосферы позволяет оперативно рассчитывать трехмерную структуру Р-области при различных гелио-геофизических условиях. Траектории конвекции моделируются в блоке расчёта электрических полей с непрерывной зависимостью их геометрии и интенсивности от входных параметров в виде значений межпланетного магнитного поля, скорости солнечного ветра, дня года, мирового времени и уровня геомагнитной и солнечной активности. В ионосферном блоке рассчитывается высотное распределение Ие. Учитываются обобщенный молекулярный ион М+ и тон 0+, концентрация которых определяется соответственно решением уравнения баланса и конвективно-диффузионного уравнения вдоль траектории конвекции с учётом поля конвекции и коротации и параметров термосферы. Все модели реализованы в виде программных средств, удобных для практических приложений.

п(М+) = п^О+) + п(О2) + п^ 2+)

(40)

СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ

1. Heppner J.P. Maynard N.C. Empirical high-latitude electric field model // J.Geophvs.Res. 1987. Vol. A2, № 5, P.4467-4489.

2. Papitashvili, V.O., Christiansen F., Neubert T. A new model of field-aligned currents derived from high-precision satellite magnetic field data // Geophvs. Res. Lett. 2002. Vol. 29, № 14, 1683. http://dx.doi.org/10.1029/2001GL014207.

3. Ruohoniemi J.M., Greenwald R.A. Dependencies of high-latitude plasma convection: consideration of IMF, seasonal, and UT factors in statistical patterns //J. Geophvs. Res. 2005. Vol. 110, A09204. http://dx.doi.org/10.1029/2004JA010815.

4. Bilitza D., McKinnell L.-A. Reinisch B.W., Fuller-Rowell T. The international reference ionosphere today and in the future // J. Geodesy. 2010, Vol. 85, P. 909-920, doi: 10.1007/s00190-010-0427-x.

5. Lukianova R., Troshichev OA., Galperin Y.I., Jorjio N.V. Reconstruction of the FAC structure along spacecraft trajectory with regard for edge effects of current sheets // Int. J.Geomag. Aeronomv. 2001, Vol. 2, № 3. P.159-172. .

6. Лукьянова Р.Ю. Распределение электрического потенциала в ионосфере Земли с учетом электродинамической сопряженности полушарий: численная модель // Матем. моделирование, 2017. Т. 29, № 5. С. 122-132.

7. Брюнелли Б.Е., Намгаладзе А.А. Физика ионосферы // Л.: Наука, 1988. 526с.

8. Владимиров B.C. Уравнения математической физики // М.: Наука, 1981. 512 с.

9. Самарский А.А., Николаев Е.С. Методы решения сеточных уравнений // М.: Наука, 1978, 592 с.

10. Lukianova R., Christiansen F. Modeling of the global distribution of ionospheric electric field based on realistic maps of field-aligned currents //J. Geophvs. Res. 2006. Vol. Ill, A03213, doi: 10.1029/2005J A011465.

11. Lukianova R., Hanuise C., Christiansen F. Asymmetric distribution of the ionospheric electric potential in the opposite hemispheres as inferred from the SuperDARN observations and FAC-based convection model // J. Atmos. Solar-Terr. Phvs. 2008. Vol. 70. P. 2324-2335, doi:10.1016/j.jastp.2008.05.015.

12. Uvarov V.M., R.Yu. Lukianova, Numerical modeling of the polar F region ionosphere taking into account the solar wind conditions // Adv. Space Research, 56, 2563-2574 10.1016/j.asr.2015.10.004, 2015.

13. Lukianova R., Uvarov V.M., Coisson P. Evolution of the high-latitude F region large-scale ionospheric irregularities under different solar wind and zenith angle conditions // Adv. Space Res. 2017. Vol. 59, № 2. P. 557-570, doi:10.1016/j.asr.2016.10.010.

14. Уваров B.M., Барашков П.Д. Типы распределения электрических полей и соответствующие им типы конвекции в полярной ионосфере. Модель // Геомагнетизм и аэрономия. 1989. Т. 29. № 4. С. 621-628.

15. Hedin, А. Е. MSIS-86 Thermospheric Model // J. Geophvs. Res. 1987. Vol. 92, № 5. P. 46494662.

REFERENCES

1. Heppner J.P. Mavnard N.C. 1987, "Empirical high-latitude electric field model" , J. Geophys. Res. Vol. A2, № 5, P. 4467-4489.

2. Papitashvili, V.O., Christiansen F., Neubert T. 2002, "A new model of field-aligned currents derived from high-precision satellite magnetic field data" , Geophys. Res. Lett. Vol. 29, № 14, 1683. http://dx.doi.org/10.1029/2001GL014207.

3. Ruohoniemi J.M., Greenwald R.A. 2005, "Dependencies of high-latitude plasma convection: consideration of IMF, seasonal, and UT factors in statistical patterns", J. Geophys. Res. Vol. 110, A09204. http://dx.doi.org/10.1029/2004JA010815.

4. Bilitza D., McKinnell L.-A. Reinisch B.W., Fuller-Rowell T. 2010, "The international reference ionosphere today and in the future", J. Geodesy. Vol. 85, P. 909-920, doi: 10.1007/s00190-010-0427-x.

5. Lukianova R., Troshichev O.A., Galperin Y.I., Jorjio N.V. 2001, "Reconstruction of the FAC structure along spacecraft trajectory with regard for edge effects of current sheets", Int. J. Geomag. Aeronomy. Vol. 2, № 3. P. 159-172.

6. Lukianova R.Yu. 2017, "Distribution of the electric potential in the Earth's ionosphere taking into account the electrodvnamic coupling of the hemispheres: numerical moel" , Math. Modeling. Vol. 29, № 5. P. 122-132*.

7. Brunelli B.E., Namgaladze A.A. 1988, Physics of the ionosphere L.: Nauka, 526 p.

8. Vladimirov V.S. 1981 Equations of mathematical physics M.: Nauka, 512 p.

9. Samarskv A.A., Nikolaev E.S. 1978, Methods of solution of network equations M.:Nauka, 592p.

10. Lukianova R., Christiansen F. 2006, "Modeling of the global distribution of ionospheric electric field based on realistic maps of field-aligned currents", J. Geophys. Res. Vol. Ill, A03213, doi: 10.1029/2005J A011465.

11. Lukianova R., Hanuise C., 2008, "Christiansen F. Asymmetric distribution of the ionospheric electric potential in the opposite hemispheres as inferred from the SuperDARN observations and FAC-based convection model" , J. Atmos. Solar-Terr. Phys. Vol. 70. P. 2324-2335, doi:10.1016/j.jastp.2008.05.015.

12. Uvarov V.M., R.Yu. Lukianova, 2015, "Numerical modeling of the polar F region ionosphere taking into account the solar wind conditions" , Adv. Space Research, 56, 2563-2574 10.1016/j.asr.2015.10.004.

13. Lukianova R., Uvarov V.M., Coisson P. 2017, "Evolution of the high-latitude F region large-scale ionospheric irregularities under different solar wind and zenith angle conditions" , Adv. Space Res. Vol. 59, № 2. P. 557-570, doi:10.1016/j.asr.2016.10.010.

14. Uvarov V.M., Barashkov P.D. 1989 "Types of electric field distribution and corresponding types of convection in the polar ionosphere. A model" , Geomagnetizm i aeronomy. Vol. 29. № 4. P. 621-628.

15. Hedin, A. E. 1987, "MSIS-86 Thermospheric Model", J. Geophys. Res. Vol. 92, № 5. P. 46494662.

Получено 27.07.2018

Принято в печать 22.10.2018

i Надоели баннеры? Вы всегда можете отключить рекламу.