Известия РАН. Серия физическая, 2023, T. 87, № 2, стр. 254-261

Связь высоты стабильных авроральных красных дуг с солнечной активностью. Результаты наблюдений на меридиане Якутска

И. Б. Иевенко 1*, С. Г. Парников 1, А. В. Сокольников 1

1 Федеральное государственное бюджетное учреждение науки Институт космофизических исследований и аэрономии имени Ю.Г. Шафера Сибирского отделения Российской академии наук
Якутск, Россия

* E-mail: ievenko@ikfia.ysn.ru

Поступила в редакцию 17.09.2022
После доработки 05.10.2022
Принята к публикации 26.10.2022

Полный текст (PDF)

Аннотация

Представлены первые результаты триангуляционных измерений высоты красных дуг на меридиане Якутска в 2015–2020 гг. Получены значения высоты нижней границы свечения в дугах 350–370 км вблизи максимума и 250–290 км во время минимума солнечной активности в 24 цикле. Уменьшение высоты красных дуг в годы минимума обусловлено значительным падением плотности верхней нейтральной атмосферы.

ВВЕДЕНИЕ

Стабильные авроральные красные (SAR) дуги возникают вследствие перекрытия кольцевого тока с внешней плазмосферой (плазмопаузой) во время магнитных бурь, которые определяются по вариации геомагнитного индекса Dst [13]. Потоки энергичных ионов нагревают плазмосферные электроны согласно механизмам, предложенным Коул [4, 5] а также Корнуолл и др. [6]. Возникающий нисходящий поток сверхтепловых электронов вдоль линий магнитного поля увеличивает окружающую электронную температуру на высотах области F2 ионосферы в виде субаврорального пика электронной температуры (Te пика). В результате происходит усиление интенсивности красной линии атомарного кислорода в SAR-дуге, отображающей плазмопаузу [3].

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

Однако спутниковые исследования, проведенные в период 1980−1990 гг., иногда указывали на то, что энергичные частицы кольцевого тока проникают во внешнюю плазмосферу во время главной фазы магнитной бури и/или во время отдельных суббурь [3, 7]. Результаты наблюдений SAR-дуг с помощью цифровых камер всего неба (All-sky imagers, ASI) также указывают на проникновение энергичных ионов кольцевого тока во внешнюю плазмосферу в начале главной фазы магнитной бури во время суббурь [810].

Наши многолетние наблюдения субаврорального свечения на меридиане Якутска (62° с.ш., 130° в.д.) на большом количестве событий детально показали, что SAR дуги появляются и/или уярчаются во время фазы расширения суббурь. Формирование SAR-дуги начинается в окрестности экваториальной границы диффузного сияния (ДС). В случае длительной суббуревой активности SAR-дуга отделяется от ДС и движется в экваториальном направлении [1115]. В последних статьях мы впервые сделали сопоставление наземной регистрации SAR-дуги ASI “Keo Sentry” (см. раздел 2) с данными одновременных измерений параметров плазмы и потоков частиц на борту спутника Van Allen Probe (B) в начале большой магнитной бури [16, 17]. В нашей работе [18] впервые представлены результаты одновременных наблюдений ASI динамики сияний и формирования SAR-дуги на авроральных и субавроральных широтах, соответственно.

В настоящее время принято считать, что высота излучения красной линии атомарного кислорода в SAR дугах имеет среднее значение 400–450 км [39]. В то же время высота максимальной скорости возбуждения эмиссии 630.0 нм в дуге потоком сверхтепловых электронов в область F2 ионосферы должна зависеть от высотного распределения плотности атомарного кислорода и коэффициента дезактивации возбужденных атомов O(1D). Известно значительное изменение высотного распределения плотности верхней нейтральной атмосферы в течение солнечного цикла. Мы сделали оценку вариаций плотности атомарного кислорода (O) по модели MSISE-90 для магнитно-спокойных дней в феврале в 2000 и 2008 гг. (для максимума и минимума 23 солнечного цикла) в локальную полночь на широте 64° с.ш. на меридиане Якутска. Было получено, что на высотах 350 и 450 км концентрация О уменьшалась в ~4 и ~9 раз, соответственно, в минимуме 23 солнечного цикла. Это дает основание полагать, что высота SAR-дуг вероятнее всего падает во время низкой солнечной активности.

По данным многолетних наблюдений на меридиане Якутска мы знаем, что в условиях низкой и средней магнитной активности (Dst ≥ –50 нТл; Kp ≤ 3–4) SAR-дуги наблюдаются на геомагнитных широтах 55°–60° [14]. В 2018–2020 гг. мы выполнили триангуляционные измерения высоты красных дуг в минимуме 24 солнечного цикла по данным одновременных наблюдений субаврорального свечения на двух станциях. Также была определена высота дуг методом триангуляции в начале спада этого цикла в 2015 г. по данным наблюдений на одной станции и спутниковой регистрации Те пиков. В этой статье представлены первые результаты этих исследований.

МЕТОДЫ НАБЛЮДЕНИЙ И ОПРЕДЕЛЕНИЯ ВЫСОТЫ SAR-ДУГ

Мониторинговые наблюдения субаврорального свечения и полярных сияний проводятся на оптической станции Маймага (MMG: 63° с.ш., 130° в.д.) с исправленными геомагнитными координатами 58° CGMLat и 202° CGMLon для эпохи 2015 г. Станция располагается севернее на 1 градус широты г. Якутска (62° с.ш., 130° в.д.) на одном географическом меридиане. Полночь по магнитному локальному времени (MLT) на меридиане Якутска соответствует ~15:50 UT. Наблюдение сияний и SAR-дуг с 2013 г. выполняются цифровой камерой всего неба (All-sky imager (ASI)) “Keo Sentry” производства фирмы “Keo Scientific Ltd”. ASI оснащен CCD камерой PI Acton ProEM 1024B и пяти полостными интерференционными фильтрами с полушириной полосы пропускания 2.0 нм. Фильтры центрированы на длины волн 557.7 и 630.0 [OI], 470.9 $\left( {{\text{N}}_{2}^{ + }} \right),$ 486.1 (Hβ), 620.0 (континуум), 480.0 нм (континуум). Данные ASI представлены в этой работе в виде кеограмм, а также отдельных изображений в проекции на поверхность Земли в исправленных геомагнитных координатах на основе моделей геомагнитного поля DGRF/IGRF.

В 2018 г. на авроральной станции Жиганск (ZGN: 66.8° с.ш., 123.4° в.д.) с геомагнитными координатами 62° CGMLat и 196° CGMLon был установлен сканирующий вдоль меридиана цифровой фотометр для регистрации субаврорального свечения и сияний в эмиссиях 557.7 и 630.0 нм [15]. Для событий одновременных наблюдений SAR-дуг на двух станциях были определены высоты их нижней границы методом триангуляции. Для выполнения измерений строились меридиональные сканограммы станций Жиганск и Маймага по исходным данным с угловым шагом 0.5 и 0.17 градуса, соответственно.

Красные дуги обычно имеют протяженность 150–350 км вдоль меридиана и до ~400 км по высоте. Поэтому определение высоты свечения в SAR-дуге по зенитным углам максимумов ее интенсивности на двух станциях дает не однозначные результаты, так как метод триангуляции применим к более локализованным источникам излучения. Например, этим методом проводились измерения высоты дискретных форм полярных сияний с линейными размерами от сотен метров до нескольких километров. Учитывая большие пространственные размеры красной дуги, мы выполнили триангуляционные измерения высоты ее нижней границы. Схема таких измерений на экваториальном крае дуги приведена на рис. 1а. Вычисления искомой высоты выполнялись по формулам (1) и (2):

(1)
$h{\text{\;}}: = \left( {\sqrt {{{c}^{2}} + {{R}^{2}} + 2cR \cdot {\text{cos}}\left( {z1} \right)} } \right) - R{\text{,}}$
(2)
$с{\text{\;}}: = \frac{{2R \cdot \sin \left( {\frac{{{\alpha }}}{2}} \right) \cdot \cos \left( {z2 - \frac{{{\alpha }}}{2}} \right)}}{{\sin \left( {z1 + z2 - {{\alpha }}} \right)}},$
где: h – искомая высота, R – радиус Земли, α – угол между станциями в градусах широты, z1 и z2 – зенитные углы регистрации нижней границы свечения на станциях Маймага и Жиганск, соответственно.

Рис. 1.

Схемы определения высоты нижней границы свечения в SAR-дуге. Измерение высоты по данным одновременных наблюдений на двух станциях (а); по данным наблюдений на одной станции и одновременной регистрации Те пика на спутнике (б).

Формулы были выведены для схемы измерений с положением SAR-дуги в области широт между двумя станциями. Поэтому, в случае регистрации красной дуги экваториальнее этой области зенитный угол z1 определяется со знаком минус. В ситуации наблюдения дуги полярнее этого интервала широт будет отрицательным угол z2. Между станциями Жиганск и Маймага угол α = = 3.8°.

Схема показывает, что для зенитных углов больше z1 и z2 на двух станциях регистрируется фоновое свечение верхней атмосферы вне области SAR-дуги. При значениях углов меньше z1 и z2 начинается регистрация усиления свечения выше фонового уровня в области нижней границы красной дуги. На схеме в градациях серого условно показано увеличение интенсивности эмиссии 630 нм в дуге на высотах F2 слоя ионосферы. Такие же измерения можно выполнить на северном крае SAR-дуги, когда она наблюдается в окрестности зенита полярной станции.

На рис. 1б приведена схема измерения высоты нижней границы свечения в дуге по одной наземной станции с использованием данных спутниковой регистрации Те пика вблизи меридиана наблюдений. В этом методе угол α определяется как разность между географическими широтами наземной станции и экваториальной (полярной) границы Те пика. Также, зенитный угол z1 = 0 в выражениях (1) и (2). Измерения электронной температуры на спутниках DMSP с полярной орбитой выполняются значительно выше (~850 км) нижней границы свечения в SAR-дугах. Поэтому географическую широту экваториального (полярного) края Te пика нужно предварительно определить на подобранной высоте нижней границы свечения в дуге посредством сопряжения вдоль силовой линии геомагнитного поля (см. рис. 1б). На этой высоте проекции на поверхность Земли границы свечения и Te пика должны примерно совпадать (см. рис. 3б и 4б). В принципе, такое сопряжение двух границ позволяет определить высоту нижней границы свечения в красной дуге. Но более точно искомая высота вычисляется по формулам (1) и (2) с использованием параметров α и z2.

Рис. 2.

Определение высоты SAR-дуги по наблюдениям на двух станциях 9 февраля 2019 г. Кеограммы в эмиссии 630.0 нм станций Жиганск (ZGN) и Маймага (MMG) в проекции на поверхность Земли для высоты свечения 270 км (а). Z – зенит станций; одновременные сканограммы SAR-дуги на двух станциях (б). Стрелками указаны отсчеты зенитных углов z1 и z2 на экваториальном крае дуги; оригинальное изображение (в автоконтрасте) всего неба в эмиссии 630.0 нм во время измерения высоты красной дуги (в); обработанное изображение в проекции на поверхность Земли для высоты свечения 270 км в исправленных геомагнитных координатах (г). Шкала интенсивности эмиссии в градациях серого приведена до 0.1 килорелея (кРл) для более детального отображения границы ДС и SAR-дуги.

Рис. 3.

Измерение высоты свечения в SAR-дуге по данным наблюдений на одной станции и спутниковой регистрации Те пика 15 марта 2018 г. Кеограммы в эмиссиях 557.7 и 630.0 нм ст. Маймага в проекции на поверхность Земли для высот свечения 110 и 270 км, соответственно (а). Шкалы интенсивности эмиссий приведены в килорелеях. Z – зенит станции; одновременные сканограмма SAR-дуги и измерения Te пика в проекции на поверхность Земли для высоты 270 км (б). Стрелкой указана широта экваториального края дуги и Te пика. Ниже, сканограмма дуги с двумя отсчетами угла z2 экваториального края дуги; как на рис. 2 (в); стрелкой показан проход спутника DMSP F15 через область красной дуги (г). Остальное как на рис. 2.

РЕЗУЛЬТАТЫ НАБЛЮДЕНИЙ С ОПРЕДЕЛЕНИЕМ ВЫСОТЫ SAR-ДУГ

На рис. 2 представлен пример определения высоты SAR-дуги по данным одновременных наблюдений на станциях Жиганск и Маймага 9 февраля 2019 г. в минимуме 24 солнечного цикла. Кеограммы на рис. 2а в эмиссии 630.0 нм [OI] двух станций показывают экваториальное движение SAR-дуги в интервале геомагнитных широт 62°–58° от момента ее возникновения в окрестности границы ДС во время суббури. На кеограмме ст. Маймага в ~1930 UT (стрелка 2 на кеограмме) центр дуги регистрировался в окрестности зенита на широте 58°. Положение проекции дуги на поверхность Земли в зените станции наблюдений не зависит от ее высоты. Это дает возможность примерно определить высоту свечения в дуге как 270 км, при которой центр ее проекции на кеограмме ст. Жиганск также находится на этой широте.

Далее, определим высоту нижней границы свечения в SAR-дуге более точно по сканограммам в эмиссии 630.0 нм [OI] двух станций в ~1858 UT на рис. 2б (стрелка 1 рис. 2а). Шкалы углов на рисунке приведены с началом отсчета от южного горизонта станций. На сканограмме ст. Жиганск стрелкой указан зенитный угол z2 экваториальной граница дуги, на которой начался резкий рост интенсивности красной линии атомарного кислорода. На ст. Маймага интенсивность красной дуги была ~35 Релей. Экваториальный край дуги отмечен на двух зенитных углах z1 в начале плавного и резкого роста интенсивности. В данном случае два отсчета z1 определяют возможную ошибку измерения высоты на границе свечения. Таким образом, по формулам (1) и (2) при α = = 3.8°, z2 = 63.0° были получены высоты нижней границы свечения в SAR-дуге 257 и 292 км для отсчетов z1 = –7.1° и –17.1°, соответственно. Среднее значение высоты составляет 274.5 км, что хорошо согласуется с заданной высотой красной дуги на кеограммах двух станций (см. рис. 2а). Рассмотренный пример указывает на достоверность триангуляционных измерений согласно схеме на рис. 1а. Изображения всего неба в эмиссии 630.0 нм [OI] на рис. 2в и 2г показывают положение красной дуги в окрестности зенита ст. Маймага во время определения ее высоты. В западном и восточном направлениях интенсивность дуги увеличивается в следствие эффекта ван Райна. Красная дуга располагалась вблизи экваториальной границы ДС.

Рисунок 3 представляет событие с измерением высоты SAR-дуги по данным наблюдений на одной станции и спутниковой регистрации Те пика 15 марта 2018 г. также во время минимума 24 цикла. Кеограммы ст. Маймага на рис. 3а показывают экваториальное расширение ДС в эмиссиях 557.7 и 630.0 нм [OI] до широты ~59° во время усиления магнитосферной конвекции после поворота Bz ММП к югу (на рис. 3 не показано). С ~16 UT началось формирование SAR-дуги в эмиссии 630.0 нм [OI] в окрестности границы ДС на геомагнитной широте ~60°. Далее, красная дуга отделилась от ДС и двигалась к экватору через зенит станции наблюдений. В ~1905 UT спутник DMSP F15 пересек SAR-дугу на высоте 853 км вблизи меридиана наблюдений (см. рис. 3г).

На рис. 3б экваториальные края SAR дуги и Te пика в проекции на поверхность Земли примерно совпали на подобранной высоте 270 км (см. предыдущий раздел). Экваториальная граница Te пика была определена с исходным широтным разрешением ~0.22° по резкому подъему электронной температуры на ~1000 К после отсчета, отмеченного стрелкой. На этой широте (60.72°) начинается и более быстрый рост интенсивности эмиссии 630.0 нм в SAR-дуге. На сканограмме (нижняя панель рис. 3б) стрелками показаны два отсчета зенитного угла z2 на экваториальной границе красной дуги согласно схеме на рис. 1б. Два значения z2 обуславливают возможную ошибку определения высоты на границе свечения в этом случае. Подстановкой значений параметров α = = 2.3°, z1 = 0 в формулы (1) и (2) для нижней границы свечения в дуге были получены высоты ~244 и ~262 км при значениях z2 = 48.0° и 46.0°, соответственно. Среднее значение высоты составляет 253 км, что на 17 км (~6%) ниже полученной высоты при сопряжении границ Te пика и дуги. Интенсивность красной линии атомарного кислорода в дуге во время определения высоты была ~30 Релей. Оригинальное и обработанное изображения всего неба в эмиссии 630.0 нм [OI] на рис. 3в и 3г показывают положение красной дуги относительно ДС во время определения ее высоты. Интенсивность дуги в окрестности меридиана станции наблюдений падает в следствие эффекта ван Райна.

Событие наблюдений SAR-дуги на ст. Маймага 10 декабря 2015 г. вблизи максимума 24 солнечного цикла приведено на рис. 4. Определение высоты свечения в дуге было выполнено также с использованием спутниковой регистрации Те пика. Кеограммы ст. Маймага на рис. 4а показывают экваториальное расширение ДС в эмиссиях 557.7 и 630.0 нм [OI] во время суббури. С ~1615 UT началось формирование SAR-дуги в эмиссии 630.0 нм в окрестности границы ДС. Далее, красная дуга отделилась от ДС и двигалась к экватору через зенит станции наблюдений.

Рис. 4.

Определение высоты свечения в SAR-дуге по данным наблюдений на одной станции и спутниковой регистрации Те пика 10 декабря 2015 г. Кеограммы как на рис. 3 (а). Высота красной линии – 350 км; как на рис. 3 для принятой высоты 350 км (б); как на рис. 3 (в); стрелкой показан проход спутника DMSP F16 через область красной дуги (г). Высота свечения 350 км. Остальное как на рис. 2 и 3.

В ~2026 UT спутник DMSP F16 пересекает SAR-дугу на высоте 857 км вблизи меридиана оптических наблюдений (см. рис. 4г). В этом случае экваториальные края SAR-дуги и Te пика в проекции на поверхность Земли совпали на подобранной высоте 350 км, как показано на рис. 4б (верхняя панель). Рост интенсивности дуги на экваториальном крае неравномерен, поэтому для определения широты общей границы было взято два значения по отсчетам электронной температуры с интервалом 0.22° (показаны стрелками). На сканограмме (нижняя панель) указан зенитный угол z2 резкой экваториальной границы дуги согласно схеме на рис. 1б. По формулам (1) и (2) с установленными параметрами z1 = 0, z2 = 59.85° были получены высоты нижней границы свечения в дуге ~367 и ~349 км при значениях α = 5.0° и 4.78°, соответственно. Среднее значение высоты равно ~358 км, что незначительно выше подобранной высоты при сопряжении границ Te пика и дуги. Интенсивность свечения в дуге во время определения высоты была ~20 Релей. Изображения всего неба в эмиссии 630.0 нм на рис. 4в свидетельствуют о регистрации слабой SAR-дуги. Обработанное изображение на рис. 4г показывает положение красной дуги на геомагнитной широте 53°–55° во время определения ее высоты. Распределение интенсивности эмиссии 630.0 нм в дуге вдоль параллели также определяется эффектом ван Райна.

ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ

Мы рассмотрели данные для трех событий наблюдения SAR-дуг и получили высоту нижней границы свечения на их экваториальном крае. Кроме этого, вблизи максимума 24-го солнечного цикла была определена высота свечения в событии 16 ноября 2015 г. по данным наблюдений на станции Маймага и регистрации Te пика на спутнике SWARM-B. Во время минимума солнечной активности было также выполнено измерение высоты нижней границы свечения в дуге по данным наблюдений на станциях Маймага и Жиганск 21 декабря 2020 г. В табл. 1 приведены все пять случаев определения высоты нижней границы излучения красной линии атомарного кислорода в SAR-дуге. Значения индекса SYM-H в таблице свидетельствуют о низком уровне буревой магнитной активности в рассмотренных событиях. Измеренные высоты показывают, что во время минимума 24 солнечного цикла нижняя граница свечения в дугах была ниже на ~100 км чем в максимуме. Средние значения этого параметра равны ~260 и ~360 км для минимума и максимума, соответственно. Уменьшение высоты красных дуг в годы минимума, наиболее вероятно, обусловлено значительным падением плотности верхней нейтральной атмосферы согласно модели MSISE-90. Для подтверждения этого вывода мы планируем проведение одновременных наблюдений SAR-дуг на двух станциях в годы максимума 25 солнечного цикла.

Таблица 1.

Результаты определения высоты нижней границы свечения в SAR-дугах на меридиане Якутска

Дата наблюдений дд.мм.гг. Метод определения высоты Высота нижней границы, км Минимум SYM-H, нТл Фаза солнечного цикла
16.11.2015
10.12.2015
15.03.2018
09.02.2019
21.12.2020
1 ст. + Те пик
1 ст. + Те пик
1 ст. + Те пик
2 станции
2 станции
352–360
349–367
244–262
257–292
242–273
–45
–10
–30
–15
–10
Максимум
Максимум
Минимум
Минимум
Минимум

(2 станции) – триангуляционные измерения высоты свечения в дуге по данным одновременных наблюдений на двух станциях; (1 ст. + Те пик) – определение высоты свечения методом триангуляции по данным наблюдений на одной стации и спутниковой регистрации Те пика.

ЗАКЛЮЧЕНИЕ

1) Рассмотрены схемы измерения высоты нижней границы свечения на южном (северном) крае SAR-дуги методом триангуляции, в том числе с использованием спутниковой регистрации Те пика. 2) Выполнены триангуляционные измерения высоты нижней границы свечения по наблюдениям на двух станциях в 2018–2020 гг. В минимуме 24 солнечного цикла получены значения высоты в интервале 240–290 км. 3) Измерены высоты нижней границы свечения в красной дуге методом триангуляции по наблюдениям на одной станции и данных спутниковой регистрации Те пика. Вблизи максимума 24 цикла в 2015 г. были получены высоты в интервале 350–370 км. 4) Высота нижней границы свечения в SAR-дугах вероятнее всего изменяется в течение цикла солнечной активности как следствие вариаций плотности верхней атмосферы согласно модели MSISE-90.

Вычисления по модели MSISE-90 были выполнены на сайте UK Solar System Data Centre (http://www.wdc.rl.ac.uk/wdcc1/msis.html). Данные измерений электронной температуры и ионосферного дрейфа на спутниках DMSP F15,16 из Space Sciences at the University of Texas at Dallas были получены в Madrigal database (http://cedar.openmadrigal.org/single/). Исправленные геомагнитные координаты рассчитывались в NASA service (https:// omniweb.sci.gsfc.nasa.gov/vitmo/cgm.html). Геомагнитные индексы были получены в WDC C2 for Geomagnetism, Kyoto (http://swdcwww.kugi.kyoto-u.ac.jp/index.html).

Работа выполнена в рамках темы государственного задания (номер гос. регистрации № 122011700172-2) при частичной поддержке Российского фонда фундаментальных исследований (проект № 21-55-50013ЯФ_а).

Список литературы

  1. Rees M.H., Akasofu S. // Planet. Space Sci. 1963. V. 11. No. 1. P. 105.

  2. Зайцева С.А., Пудовкин М.И., Дряхлов В.В., Дьяченко В.Н. // Геомагн. и аэроном. 1971. Т. 11. № 5. С. 853.

  3. Kozyra J.U., Nagy A.F., Slater D.W. // Rev. Geophys. 1997. V. 35. No. 2. P. 155.

  4. Cole K.D. // J. Geophys. Res. 1965. V. 70. No. 7. P. 1689.

  5. Cole K.D. // Ann. Geophys. 1970. V. 26. No. 1. P. 187.

  6. Cornwall J.M., Coroniti F.V., Thorne R.M. // J. Geophys. Res. 1971. V. 76. No. 19. P. 4428.

  7. Horwitz J.L., Brace L.H., Comfort R.H., Chappell C.R. // J. Geophys. Res. 1986. V. 91. No. A10. Art. No. 11203.

  8. Shiokawa K., Miyoshi Y., Brandt P.C. et al. // J. Geophys. Res. 2013. V. 118. No. 1. P. 256.

  9. Mendillo M., Baumgardner J., Wroten J. // J. Geophys. Res. Space Phys. 2016. V. 121. No. 1. P. 245.

  10. Martinis C., Baumgardner J., Mendillo M. et al. // J. Geophys. Res. Space Phys. 2019. V. 124. No. 6. P. 4658.

  11. Иевенко И.Б. // Геомагн. и аэроном. 1993. Т. 33. № 5. С. 42; Ievenko I.B. // Geomagn. Aeron. 1994. V. 33. No. 5. P. 599.

  12. Иевенко И.Б. // Геомагн. и аэроном. 1995. Т. 35. № 3. С. 37; Ievenko I.B. // Geomagn. Aeron. 1995. V. 35. No. 3. P. 331.

  13. Иевенко И.Б. // Геомагн. и аэроном. 1999. Т. 39. № 6. С. 26; Ievenko I.B. // Geomagn. Aeron. 1999. V. 39. No. 6. P. 697.

  14. Иевенко И.Б., Алексеев В.Н. // Геомагн. и аэроном. 2004. Т. 44. № 5. С. 643; Ievenko I.B., Alekseev V.N. // Geomagn. Aeron. 2004. V. 44. No. 5. P. 592.

  15. Ievenko I.B., Parnikov S.G., Alexeyev V.N. // Adv. Space Res. 2008. V. 41. No. 8. P. 1252.

  16. Иевенко И.Б., Парников С.Г. // Геомагн. и аэроном. 2020. Т. 60. № 6. С. 751; Ievenko I.B., Parnikov S.G. // Geomagn. Aeron. 2020. V. 60. No. 6. P. 737.

  17. Ievenko I.B. // J. Atmos. Sol. Terr. Phys. 2020. V. 209. Art. No. 105386.

  18. Иевенко И.Б., Парников С.Г. // Геомагн. и аэроном. 2022. Т. 62. № 2. С. 171; Ievenko I.B., Parnikov S.G. // Geomagn. Aeron. 2022. V. 62. No. 1–2. P. 32.

Дополнительные материалы отсутствуют.