Институт физики Земли им. О.Ю. Шмидта РАН
Москва, Россия
Институт физики Земли им. О.Ю. Шмидта Российской академии наук
Москва, Россия
Москва, Россия
Москва, Россия
Институт физики Земли им. О.Ю. Шмидта РАН
Москва, Россия
Москва, Россия
Москва, Россия
УДК 55 Геология. Геологические и геофизические науки
Рассмотрена задача ранней диагностики геомагнитных бурь на основе применения моделей координат движений центров солнечных корональных выбросов массы (КВМ) и наблюдений их угловых положений, получаемых от измерительных систем космического мониторинга. Предложен метод ранней диагностики геомагнитных бурь, введена функция прогноза расстояния между Землей и центрами КВМ, сформирована процедура принятия решений. Разработан алгоритм диагностики, основанный на минимизации функционала обобщенной триангуляции. Реализованы пример вычисления функции прогноза расстояния и процедуры принятия диагностического решения на основе моделей координат и модельных наблюдений угловых положений центров солнечных КВМ. Определена эффективность процедуры принятия решений для алгоритма ранней диагностики геомагнитных бурь.
корональные выбросы масс, геомагнитные бури, космический мониторинг, диагностические решения, функционал триангуляции
ВВЕДЕНИЕ
Солнечные корональные выбросы массы (КВМ) оказывают существенное воздействие на магнитосферу Земли, вызывают магнитные бури и приводят к возможным большим проблемам при функционировании технических и биологических объектов. Ранняя диагностика геомагнитных бурь является актуальной научной задачей [Handbook of Cosmic …, 2015; Solar Eruptions and …, 2006]. По данной проблематике имеются две категории публикаций по моделям, связанным с диагностикой и предсказанием времени прибытия КВМ. К первой категории относятся публикации с описаниями моделей, в которых реализуются подходы, основанные на уравнениях магнитной гидродинамики плазменных образований, ко второй — публикации, в которых рассматриваются различные варианты нелинейных моделей движения КВМ.
В cтатье [Odstrcil, 2003] описывается модель ENLIL, используемая NASA для прогнозов в гелиосфере, которые осуществляются путем решения уравнений магнитной гидродинамики. Основной особенностью данной модели является отсутствие в ней механизма разогрева солнечной короны, который описывается феноменологически. Точность прогнозов прибытия КВМ к Земле по этой модели составляет около восьми часов. В статье [Jin et al., 2017] описывается модель BATSRUS, также основанная на уравнениях магнитной гидродинамики, но включающая механизм разогрева солнечной короны за счет альфвеновской турбулентности. Благодаря этому модель последовательно описывает не только гелиосферу, но и солнечную корону. Описан также инструмент EEGGL, позволяющий по оптическим и магнитным наблюдениям задавать начальные условия для КВМ, дальнейшее движение которого через корону и гелиосферу может быть рассчитано согласно уравнениям модели. На сайте [http://helio-weather.net/archive/2008/05] представлена модель WSA-ENLIL, являющаяся усовершенствованной версией [Odstrcil, 2003]. Подходы, использующие решения уравнений магнитной гидродинамики и основанные на детальном рассмотрении физики процессов в гелиосфере, потенциально могут обеспечить эффективное решение задачи ранней диагностики.
В статье [Owens, Cargill, 2004] рассмотрены три феноменологические модели распространения КВМ, позволяющие вычислять время прибытия. В первой модели считается, что выброс движется с постоянным ускорением на всем пути от Солнца до орбиты Земли. Во второй предполагается, что до некоторого расстояния, меньшего радиуса земной орбиты, выброс движется равноускоренно, а далее равномерно. В третьей модели считается, что ускорение выброса пропорционально разности скоростей выброса и окружающего солнечного ветра, коэффициент пропорциональности зависит от расстояния до Солнца по степенному закону. Первая модель содержит два параметра (начальную скорость и ускорение), вторая — три параметра (начальную скорость, ускорение и расстояние, на котором ускорение прекращается). Третья модель содержит четыре параметра (начальную скорость, степень и множитель в коэффициенте увлечения и асимптотическую скорость солнечного ветра). В [Gopalswamy et al., 2000] на основе 28 межпланетных событий КВМ разработана эмпирическая формула времени прибытия КВМ, работа которой исследована: показано, что ее предсказательные характеристики лучше для быстрых КВМ, чем для медленных. В работе [Mittal, Narain, 2015] для вычисления времени прибытия используется простейшая модельная линейная зависимость между скоростью на вылете КВМ от Солнца и временем распространения до Земли; предлагаются к рассмотрению случаи медленных (<500 км/с) и быстрых (>500 км/с) КВМ. Статья [Michalick et al., 2004] содержит описание подхода к получению скорости КВМ, основанного на решении специальной оптимизационной задачи.
И, наконец, материал [Gopalswamy, 2016] является обзором истории изучения КВМ, в основном сконцентрированным на данных наблюдений и их систематизации. Вопросы диагностики и соответствующих моделей распространения КВМ подробно не обсуждаются.
К рассматриваемым здесь измерительным системам космического мониторинга относятся спутниковая система SOHO [Brueckner et al., 1995; Delaboudiniere et al., 1995] (проект NASA и ESA (European Space Agency)); спутниковая система STEREO [Howard et al., 2008] (проект NASА) и мюонный годоскоп (МГ) «Ураган» [Barbashina et al., 2008; Yashin et al., 2015] (проект НИЯУ МИФИ (РФ)).
1. Barbashina N.S., Kokoulin P.А., Kompaniets K.G., et al. The URAGAN wide-aperture large-area muon hodoscope // Instruments and Experimental Techniques. 2008. V. 51, N 2. P. 180-186.
2. Brueckner G.E., Howard R.A., Koomen M.J., et al. The large angle spectroscopic coronagraph (LASCO) // Solar Phys. 1995. V. 162, iss. 1-2. P. 357-402.
3. Delaboudiniere J.-P., Artzner G.E., Brunaud, J., et al. EIT: Extreme-Ultraviolet Imaging Telescope for the SOHO Mission // Solar Phys. 1995. V. 162, iss. 1-2. P. 291-312.
4. Gopalswamy N. History and development of coronal mass ejections as a key player in solar terrestrial relationship // Geosci. Lett. 2016. 3:8, DOI:https://doi.org/10.1186/s40562-016-0039-2.
5. Gopalswamy N., Lara A., Lepping R.P., et al. Interplanetary acceleration of coronal mass ejection // Geophys. Res. Lett. 2000. V. 27. P. 145.
6. Handbook of Cosmic Hazards and Planetary Defense / Ed. J.N. Pelton, F. Allahdadi. Springer International Publishing, 2015. 1127 р.
7. Howard R.A., Moses J.D., Socker D.G., Dere K.P., Cook J.W. Secchi Consortium. Sun Earth connection coronal and heliospheric investigation // Space Sci. Rev. 2008. V. 136. P. 67-115.
8. Jin M., Manchester W.B., van der Holst B., et al. Data-constrained coronal mass ejections in a global magnetohydrodinamics model // Astrophys. J. 2017, V. 834. Iss. 2, article id. 173, 9 pp. DOI:https://doi.org/10.3847/1538-4357/834/2/173.
9. Michalick G., Gopalswany N., Lara A., Manoham P.K. Arrival time of halo coronal mass ejection in the vicinity of the Earth // Astron. Astrophys. 2004. V. 423. P. 729.
10. Mittal N., Narain U.D. On the arrival of halo coronal mass election in the vicinity of the Earth // J. Astron. Geophys. 2015. V. 4. P. 100-105.
11. Odstrcil D. Modeling 3-D solar wind structure // Adv. Space Res. 2003. V. 32, N 4. P. 497-506.
12. Owens M., Cargill P. Predictions of the arrival time of Coronal Mass Ejections at 1 AU: an analysis of the causes of errors // Ann. Geophys. 2004. V. 22. P. 661.
13. Singiresu S. Rao. Engineering Optimization // Theory and Practice. John Wiley & Sons, 2009. 813 p.
14. Solar Eruptions and Energetic Particles / Ed. N. Gopalswamy, R. Mewaldt, J. Torsti. Geophysical Monograph Ser. V. 165. American Geophysical Union, 2006. 385 p.
15. Xue X.H., Wang C.B., Dou X.K. An ice cream cone model for coronal mass ejections // J. Geophys. Res. 2005. V. 110. P. A08103.
16. Yashin I.I., Astapov I.I., Barbashina N.S., et al. Real-time data of muon hodoscope URAGAN // Adv. Space Res. 2015. V. 56, iss. 12. P. 2693-2705.
17. URL: http://helio-weather.net/archive/2008/05 (дата обращения 28 мая 2018).
18. URL: https://stereo.gsfc.nasa.gov/where (дата обращения 28 мая 2018).
19. URL: https://cdaw.gsfc.nasa.gov/CME (дата обращения 28 мая 2018).
20. URL: https://stereo.gsfc.nasa.gov/cgi (дата обращения 28 мая 2018).
