Новый подход в прогнозе полярной и глобальной геомагнитной активности: теоретическое и экспериментальное обоснование Тамара Кузнецова, Алексей Лаптухов и Валерий Петров ИЗМИРАН 2-я научная конференция ПОЛАР-2012, мая 2012 года, ИЗМИРАН

Тезисы доклада (что представляем) Необходимо направить усилия на понимание физики взаимодействия солнечного ветра с магнитосферой и методы мониторинга ключевых параметров космической погоды. Оспорим распространенное мнение, что Bz ММП 0) в GSM, (или другой с.к.) - ключевой параметр пересоединения и контролирует процесс взаимодействия. На основе результатов (теории и эксперимента) покажем, что это не так. Предложим геоэффективные независимые инвариантные параметры, учитывающие изменение геометрии взаимодействия при орбитальном и суточном движениях Земли. Применим их для решения проблемы годовой и UT вариации ГА, стоящей с 1856 г (Sabine), и представим объяснение этой вариации. Представим результаты исследования распределения в году значений Кр=8-9, имеющего наряду с классическими пиками вблизи равноденствий максимумы в июле и ноябре, что важно для оценки вероятности возникновения больших возмущений.

ПРЕДЛАГАЕМЫЙ ПОДХОД: КЛАССИФИКАЦИЯ ПРИЧИН ГЕОМАГНИТНОЙ АКТИВНОСТИ Результаты, которые использованы в докладе были опубликованы в работах Кузнецова и др. 2002, 2006a, 2006b, 2009a, 2009b, 2010, Соавторы: А. Лаптухов, В. Кузнецов, В. Петров. Зависящие от времени вариации параметров солнечного ветра и ММП. Источник этих изменений – изменения на солнце. Временные вариации взаимной ориентации скорости солнечного ветра и ММП относительно геомагнитного момента, возникающие при орбитальном и суточном движении Земли. Изменение этой геометрии может привести к изменению эффективности взаимодействия с магнитным полем Земли (и ГА) без каких- либо изменений величин параметров солнечного ветра и ММП. Учет геометрии взаимодействия и ее проявление в геомагнитной активности – основная цель наших исследований

ОСНОВНЫЕ МЕХАНИЗМЫ, УЧИТЫВАЮЩИЕ ВЛИЯНИЕ ГЕОМЕТРИИ ВЗАИМОДЕЙСТВИЯ Russell-McPherron механизм, 1973 (RM). магнитные поля, лежащие в плоскости солнечного экватора имеют максимальную компоненту, южного направления в GSM вблизи max ГА. параметр Bz ММП в GSM Ключевой угол : между ММП и Ме Теоретическая основа: механизм пересоединения Dangey (1961) Основные трудности: неспособность объяснить количественно наблюдения, неверная фаза UT вариации Equinoctical механизм (EM) Bartels (1932). основан на изменениях эффективности взаимодействия от ориентации диполя по отношению к V солн.ветра. Ключевой угол : между V и Ме Теория: механизм квазивязкого вязкого взаимодействия Boller and Stolov (1970) Основные трудности: те же самые, что и для RM

КРИТИКА Bz ММП КАК ГЕОЭФЕКТИВНОГО ПАРАМЕТРА Bm = (B,M)=BxMx+ByMy+BzMz= B M cos (1) Bm is a component of the IMF B directed along the dipole magnetic moment vector M, is an angle between the vectors of M and B; B is module of B. В GSM My=0 (колебания вектора М в YZ-плоскости отсутствуют), значит ByMy=0 in (1), т.е. влияние By ММП не учитывается в GSM (By-дополнение к правой с.к.). Изменения геоэффективности By ММП при орбитальном и суточном движениях магнитного момента Земли Ме должны учитываться!

ЭЛЕКТРИЧЕСКОЕ ПОЛЕ КАК ОСНОВНОЙ ПАРАМЕТР СКОРОСТИ ПЕРЕСОЕДИНЕНИЯ Bz не является основным параметром в физических моделях пересоединения. Основным параметром, который определяет скорость пересоединения, является электрическое поле Е Е не является результатом процесса, а внешним независимым параметром, который практически полностью определяет скорость пересоединения независимо от условий задачи. (Vasyliunas, 1975 Pudovkin, 1976) Ориентация электрического поля солнечного ветра Е относительно геомагнитного момента - основной эффективный параметр в нашем исследовании

Инвариантные независимые от системы координат параметры пересоединения We assume that reconnection in general case of arbitrary oriented fields is determined by three parameters (Kuznetsova and Laptukhov, 2002) Bm = (B,M)=BxMx+ByMy+BzMz= B M cos, (1) Bm is a component of the IMF B directed along the dipole magnetic moment vector M, is an angle between the vectors of M and B; B is module of B. Here and below M =1 for simplicity. Em = (E,M) = V(MzBy - MyBz)= E [(My) 2 +(Mz) 2 ] 1/2 cos(Uem), (2) Em is a component of the solar wind electric field E along the geomagnetic moment vector M.V is directed along X (GSE), see (2a); Ume is angle between vectors E and M Emv= (E,[M,V])/( M V )=V(MyBy + BzMz)= E [(My) 2 +(Mz) 2 ] 1/2 cos(Uemv)(3) E=[B,V] is vector of the electric field of the solar wind, V is vector of velocity of the solar wind and V= V is its value, Emv is a parameter that is proportional to a component of the solar wind electric field E along the [M,V] vector. cos, cosUem, cosUemv можно рассчитать на основе измерений, а их влияние на ГА можно оценить количественно

РЕЗУЛЬТАТЫ: ЗАВИСИМОСТЬ Кр ОТ УГЛА МЕЖДУ ММП И ГЕОМАГНИТНЫМ МОМЕНТОМ При повороте ММП от севера (cos = -1) к югу (cos = 1) Кр изменяется от Кр=15 до Кр=40; Кр=25при изменении cos Kp~25, Kp объясняет 40% изменений. Кривые Кр=F(cos ) для различных лет коррелируют друг с другом с Сс=0.99, отличаясь уровнем Кр: cos - геометрический фактор, повторящийся год от года и объясняющий ~ % изменений Кр.

ЗАВИСИМОСТЬ Кp ОТ Emv (КОМПОНЕНТЫ Е, ПОПЕРЕЧНОЙ Ме В ПЛОСКОСТИ YZ GSE с.к.) Emv объясняет ~95% Kp изменений Kp несимметрично относительно знака Emv: Быстрое увеличение Kp с увеличением Emv при Emv>0 (поле Е с утра на вечер) и медленное при Emv

Зависимость Кр от угла между электрическим полем и геомагнитным моментом Фиг. Зависимость Кр от угла Uem между векторами M и E для разных диапазонов величины E. При Uem=0 0 (E M) и Uem=180 0 (E M) Кр принимает средние значения При Uem=90 0 (E M и направлено с вечера на утро) Кр имеет минимум При Uem=270 0 (E M и направлено с утра на вечер) Кр имеет максимум ГА изменяется при постоянной величине Е=10 мВ/м в диапазоне Kp=6±2 из-за изменения только угла Uem, достигая Кр=6 при северном направлении ММП Bz из-за роста компоненты Em Uemv Kp 0_2 2_4 4_6 6_8 8_10 10_12 0_12

ОЦЕНКА ВЛИЯНИЯ КЛЮЧЕВЫХ УГЛОВ ОСНОВНЫХ МЕХАНИЗМОВ НА Kp cos 0 cosβ Kp(cos ) Kp(cosβ) cos = (V, M)/ V M = Mx/M, where |M |=1, - угол между Ме и V; cosβ (угол V и Mе, квазивязкое взаимодействие и равноденственная модель) отвечает за ~15% изменений Kp, cosα (угол ММП B and Mе, Dangey and RM-модель) - for 35%, Uem (E and Mе) за 50%.

ГОДОВАЯ ВАРИАЦИЯ Кр and Emv ДЛЯ РАЗЛИЧНЫХ НАПРАВЛЕНИЙ ММП В ПЛОСКОСТИ XY GSE Annual variation Emv (mV/m) for different directions of IMF at the XY plane (GSE) shown by signs Bx By Annual variation of Kp for different directions of the IMF at the XY plane (GSE) shown by Bx,By

ГОДОВАЯ ВАРИАЦИЯ Кр and Emv ДЛЯ РАЗЛИЧНЫХ НАПРАВЛЕНИЙ ММП В ПЛОСКОСТИ YZ GSE Annual variation of Kp for different directions of the IMF at the ZY plane (GSE) shown by Bz, By

НОВОЕ ОБЪЯСНЕНИЕ ПОЛУ-ГОДОВОЙ ВАРИАЦИИ ГЕОМАГНИТНОЙ АКТИВНОСТИ Emv = V(MyBy + BzMz) Bz~0 при данных за г (симметрия относительно знака Bz) Emv = VMyBy Расчеты показывают: My 0 во 2-ю, значит By 0 – во 2-ю (Emv должно быть >0) Фаза определяется годовой вариацией к-ты Му магнитного момента в плоскости, л инии Солнце-Земля: экстремумы в равноденствия, 0- в солнцестоя- ния. Амплитуда модулируется компонентой Ez cолнечного ветра, плоскости эклиптики

РАСПРЕДЕЛЕНИЕ БОЛЬШИХ ВОЗМУЩЕНИЙ В ГОДУ Annual distribution of counts with values aa>100 nT N, m is number of 14-day interval for year, Annual distribution of counts with values Kp=8- 9 N, m is number of 14-day interval for year, Peaks: 4 April±7 days (m=7), 28 September (m=19), 24 July (m=14-15), (m=23-24)-18 November. Дополнительные пики в июле и ноябре, увеличенный пик в марте- апреле.

ГОДОВОЕ И UT РАСПРЕДЕЛЕНИЕ ВОЗМУЩЕНИЙ С ИНДЕКСОМ АА>200 нT Annual distribution of quantities N with аа>200 nT for the period , m is number of 14-days interval in a year. Spring peak in March (m=7). Fall peak in September (m=19). Additional peak in July (m=14) and smoothed peak in November UT distribution of quantities with аа>200 nT N for the period , h is UT Дополнительные пики в июле и ноябре. Максимум UT распределения в 18 часов

ГОДОВАЯ ВАРИАЦИЯ ВЕЛИЧИНЫ Е-ПОЛЯ СОЛНЕЧНОГО ВЕТРА В ЧЕТНЫХ-НЕЧЕТНЫХ ЦИКЛАХ АКТИВНОСТИ Годовое распределение величины Е в году:1 (четные) - N20+22, 2 (нечетные) -N21+23, В нечетных циклах Е больше, чем в четных. В циклах 2 пики (июль, ноябрь, апрель) с Т~150 д (Rieger,1984). У нас в W T=146.2д; d/146.2=5/2 1/0.9= 1/11+ 1/1 (1/г) г=3/5г – проявление 11-л. цикла на земле Величина Е в марте и октябре не зависит от четности-нечетности солнечного цикла Наибольшие отличия величин Е в четных и нечетных циклах в ноябре, июле и апреле

Годовая вариация разности величин поля Е солнечного ветра в циклах по максимуму W Максимумы разностей в июле, мае и ноябре (Е между циклами) В июле D= 20% от среднего уровня Е= 1 mV/m (D=0.2 mV/m) В мае и ноябре D=10%. Годовая вариация разности D величины поля солнечного ветра Е мВ/м в двух группах данных D= E (1)- E(2); 1 –нечетным - четным циклом (по максимумам W), 2 –между последующим четным - нечетным май и ноябрь - пики скоростей поступления электромагнитной энергии в полярные шапки (слайд 19)

ГОДОВАЯ И UT ВАРИАЦИЯ КОМПОНЕНТЫ ВЕКТОРА ПОЙТИНГА ВДОЛЬ ГЕОМАГНИТНОГО МОМЕНТА Annual variation of Pm, component of Poyting flux (P=ExH, w/m) along moment Me, for 8 of 3-hour intervals of UT Pm>0 имеет максимум в Мае (СП); Pm

20 Концепция магнитной бури с Bz

2-я научная конференция ПОЛАР-2012, мая 2012 года, ИЗМИРАН СПАСИБО ЗА ВНИМАНИЕ ! Заинтересованным могу выслать нашу последнюю статью Advances in space Research v. 47, 978–990, 2011, doi: /j.asr