Москва, Россия
Международная справочная модель ионосферы IRI использует в качестве показателя уровня солнечной активности эффективный ионосферный индекс IG12, выведенный по измерениям критической частоты foF2 сетью ионозондов. В модели ионосферы и плазмосферы IRI-Plas в качестве эффективного показателя уровня солнечной активности используется глобальное электронное содержание (GEC, global electron content), сглаженное скользящим 12-месячным фильтром, GEC12. Расчеты GEC основаны на глобальных картах полного электронного содержания (ТЕС, total electron content) с 1998 г., а реконструкция его за предыдущие годы и прогноз на будущее выполнены по эмпирической модели линейной зависимости GEC от солнечной активности. В настоящее время возникла необходимость переоценки солнечных и ионосферных индексов в ионосферных моделях в связи с переходом с 1 июля 2015 г. на обновленный расчет ряда солнечных пятен — SSN2 [Clette et al., 2014]. Использование SSN2 вместо прежнего ряда SSN1 приводит к увеличению погрешности прогнозов по ионосферной модели. Предлагается формула для преобразования сглаженного 12-месячным скользящим фильтром ряда SSN212 в эквивалент ранее принятого базового индекса SSN112=R12, используемого для долгосрочного прогнозирования в моделях IRI и IRI-Plas. Выведены регрессионные соотношения между GEC12, индексом солнечных пятен R12 и индексом потока солнечного радиоизлучения на волне 10.7 см F10.712. Сравнение расчетов по моделям IRI-Plas и IRI с данными наблюдений и прогноза по Москве в 23-м и 24-м циклах солнечной активности показало преимущество использования индекса GEC12 с моделью IRI-Plas.
Глобальное электронное содержание, число солнечных пятен, поток солнечного радиоизлучения, ионосфера, плазмосфера, модель, IRI, IRI-Plas
ВВЕДЕНИЕ
Трансионосферное распространение радиоволн обусловливает возможность передачи сигналов навигационных спутников при наземно-космической и межспутниковой связи. Для обеспечения надежности передачи и приема сигнала при долговременном планировании космических экспериментов используется прогноз состояния ионосферы с помощью известных моделей, таких как эмпирическая Международная справочная модель ионосферы IRI (International Reference Ionosphere) и ее расширение на высоты плазмосферы IRI-Plas [Bilitza et al., 2011; Gulyaeva, Bilitza, 2012], Российская стандартная модель ионосферы SMI [Chasovitin et al., 1987], эмпирическая модель для расчетов трансионосферного распространения радиоволн NeQuick [Nava et al., 2008]. Результаты расчетов зависят от задания солнечных и ионосферных управляющих параметров моделей. В настоящее время возникла необходимость переоценки солнечных и ионосферных управляющих параметров ионосферных моделей в связи с недавним обновлением многолетнего ряда чисел солнечных пятен за период с 1818 г. по настоящее время [Clette et al., 2014]. Второе поколение ряда чисел солнечных пятен - SSN2 значительно отличается от исходного многолетнего ряда SSN1. Значения SSN2 вблизи максимума солнечной активности зачастую выше значений индекса потока солнечного радиоизлучения на волне 10.7 см F10.7, которые, в свою очередь, в среднем на 60 единиц превышают значения SSN1 [Ahluwalia, 2016]. Прямое использование SSN2 вместо SSN1 в прикладных задачах, в том числе в прогнозах по ионосферным моделям, может привести к заметным ошибкам в результатах модельных расчетов.
Отдельные блоки ионосферных моделей IRI, IRI-Plas, SMI, NeQuick основаны на различных индексах солнечной активности или их ионосферных эквивалентах. В этих моделях трехмерное представление электронной плотности базируется на привязке вертикального распределения плотности электронов к максимальной электронной плотности и высоте слоя F2 ионосферы по планетарным картам, разработанным международным радиосоюзом ITU-R [CCIR Atlas, 1983]. Планетарные карты критической частоты foF2 и фактора распространения радиоволн M3000F2, связанного с высотой максимума ионизации hmF2, построены на основе коэффициентов разложения в ряд с помощью сферического гармонического анализа указанных параметров по измерениям ионозондов. Коэффициенты ряда получены в зависимости от широты, долготы, мирового времени (UT) от 0 до 24 ч с шагом в 1 ч, сезонных изменений для двенадцати месяцев года и четырех уровней солнечной активности R12=0, 50, 100, 150, с соответствующей интерполяцией для промежуточных значений этих параметров [Jones, Gallett, 1965].
Блок расчета критической частоты слоя F2 foF2 (или связанной с ней максимальной плотности электронов NmF2) в модели IRI при вводе карт ITU-R использует в качестве показателя уровня солнечной активности эффективный ионосферный индекс IG12, выведенный по измерениям критической частоты foF2 сетью ионозондов [Liu et al., 1983]. При расчете высоты максимума hmF2 в моделях IRI и IRI-Plas на основе карт M3000F2 (ITU-R) используются критические частоты foF2 и foE и 12-месячное сглаженное значение R12 на основе первичного ряда чисел солнечных пятен SSN1 [Bilitza et al., 1979]. Ионосферный индекс солнечной активности IG [Liu et al., 1983] и его сглаженный 12-месячным скользящим фильтром ряд IG12 основаны на полуденных значениях критической частоты foF2 по измерениям 13 ионозондов, приведенных к масштабу чисел солнечных пятен R12. Сравнение индексов IG12 и R12 за последние шесть солнечных циклов [Bilitza et al., 2012] показывает их различия, наиболее заметные в максимуме и минимуме солнечной активности.
В литературе приводятся различные комбинации индексов солнечной активности для использования в ионосферных моделях. В работе [Деминов, 2016] показано, что при долгосрочном прогнозировании критической частоты foF2 индекс потока солнечного радиоизлучения F12 более эффективен, чем R12. Согласно работе [Ратовский и др., 2015], локальная модель параметров максимума слоя F2 линейно зависит от индекса F10.7, нормированного к значению F10.7=100 е.и. (единица индекса 1 е.и. = 10-22 Вт·м-2·Гц-1). Для модели карт полного электронного содержания (ТЕС, total electron content) наилучшую точность показало применение комбинированного индекса потока солнечного радиоизлучения F10.7, а также числа солнечных пятен SSN1, осредненных по сумме трех компонентов: скользящего среднего за 3 дня, среднего за 27 предыдущих дней и среднего за 81 предыдущий день [Maruyama, 2010]. Регрессионная зависимость между рядом чисел солнечных пятен SSN1 и индексом F10.7, используемая в моделях IRI и IRI-Plas, стабильно сохранялась до 2000 г., но после 2001 г. эта связь изменилась [Lukianova, Mursula, 2011]. В настоящий момент все прежние оценки эффективности числа солнечных пятен в метрике SSN1 нуждаются в пересмотре, поскольку с 1 июля 2015 г. этот индекс не измеряется, а вместо него проводится мониторинг солнечной активности по индексу SSN2 [Clette et al., 2014].
В ионосферно-плазмосферной модели IRI-Plas для расчета критической частоты foF2 и высоты максимума hmF2 по картам ITU-R в качестве индикатора солнечной активности используется глобальное электронное содержание в ионосфере и плазмосфере Земли (GEC, global electron content), определяемое по картам полного электронного содержания [Афраймович, Перевалова, 2006; Afraimovich et al., 2008; Gulyaeva, Veselovsky, 2014]. Преимущество модельной системы IRI-Plas по сравнению с IRI заключается в наличии модели плазмосферы, дополняющей расчет высотных профилей электронной плотности и температуры до высоты 20200 км над Землей (орбита навигационных спутников системы GPS), тогда как IRI позволяет делать расчеты только до высоты 2000 км в ионосфере.
Исследование параметра GEC в качестве индикатора солнечной активности и сравнение его с дру-гими солнечными и ионосферными индексами проводится в данной работе с учетом модернизации ряда чисел солнечных пятен. Целью исследования является обновление набора управляющих параметров модели IRI-Plas, в том числе при расчете по исходным картам максимума электронной плотности и высоты слоя F2, и предостережение пользователей моделей от возможных ошибок, связанных с применением обновленного ряда чисел солнечных пятен для управления ионосферными моделями.
1. Афраймович Э.Л., Перевалова Н.П. GPS-мониторинг верхней атмосферы Земли. Иркутск: ИСЗФ СО РАН, 2006. 480 с.
2. Деминов М.Г. Индекс солнечной активности для долгосрочного прогноза ионосферы // Космические исследования. 2016. Т. 54, № 1. С. 3-9. DOI:https://doi.org/10.7868/S0023420616010064.
3. Ратовский К.Г., Ойнац А.В., Медведев А.В. Сходства и различия регулярных вариаций параметров F2-слоя полярной и среднеширотной ионосферы в восточно-сибирском секторе // Солнечно-земная физика. 2015. Т. 1, № 2. С. 70-79. DOI:https://doi.org/10.12737/7832.
4. Afraimovich E.L., Astafyeva E.I., Oinats A.V., et al. Global electron content: a new conception to track solar activity // Ann. Geophys. 2008. V. 26, N 2. P. 335-344.
5. Ahluwalia H.S. The descent of the solar cycle 24 and future space weather // Adv. Space Res. 2016. V. 57, iss. 2. P. 710-714, DOI:https://doi.org/10.1016/j.asr.2015.11.015.
6. Bilitza D., Sheikh N.M., Eyfrig R. A global model for the height of the F2 peak using M3000 values from the CCIR // Telecomm. J. 1979. V. 46. P. 549-553.
7. Bilitza D., McKinnell L.A., Reinisch B., Fuller-Rowell T. The international reference ionosphere today and in the future // J. Geodesy. 2011. V. 85. P. 909-920. DOI:https://doi.org/10.1007/s00190-010-0427-x.
8. Bilitza D., Brown S.A., Wang M.Y., et al. Measurements and IRI model predictions during the recent solar minimum // J. Atmos. Solar-Terr. Phys. 2012. V. 86. P. 99-106. DOI:https://doi.org/10.1016/j.jastp.2012.06.010.
9. Bilitza D. The International Reference Ionosphere - Status 2013 // Adv. Space Res. 2015. V. 55, iss. 8. P. 1914-1927. DOI: 10.1016.j.asr.2014.07.032.
10. CCIR Atlas of Ionospheric Characteristics // Comite Consultatif International des Radio Communications Rep. 340. Geneve: International Telecommunication Union, 1983.
11. Chasovitin Yu.K., Shirochkov A.V., Besprozvannaya A.S., et al. An empirical model for the global distribution of density, temperature and effective collision frequency of electrons in the ionosphere // Adv. Space Res. 1987. V. 7, iss. 6. P. 49-52.
12. Clette F., Svalgaard L., Vaquero J.M., Cliver E.W. Revisiting the sunspot number: a 400-year perspective on the solar cycle // Space Sci. Rev. 2014. V. 186, iss. 1. P. 35-103.
13. Gulyaeva T.L., Bilitza D. Towards ISO Standard Earth Ionosphere and Plasmasphere Model // New Developments in the Standard Model. Nova Science Publishers Inc., 2012. P. 1-39. Available at https://www.novapublishers.com/catalog/pro-duct_info.php?products_id=35812. (дата обращения 1 сентября 2016 г.).
14. Gulyaeva T.L., Veselovsky I.S. Imaging Global Electron Content backwards in time more than 160 years ago // Adv. Space Res. 2014. V. 53, iss. 3. P. 403-411. DOI:https://doi.org/10.1016/j.asr. 2013.11.036.
15. Hao Y.Q., Shi H., Xiao Z., Zhang D.H. Weak ionization of the global ionosphere in solar cycle 24 // Ann. Geophys. 2014. V. 32. P. 809-816. DOI:https://doi.org/10.5194/angeo-32-809-2014.
16. Jones W.B., Gallet R.M. Representation of diurnal and geographical variations of ionospheric data by numerical method // Telecomm. J. 1962. V. 29. P. 129; 1965. V. 32. P. 18.
17. Liu R., Smith P., King J. A new solar index which leads to improved foF2 prediction using the CCIR atlas // Telecomm. J. 1983. V. 50. P. 408-414.
18. Lukianova R., Mursula K. Changed relation between sunspot numbers, solar UV/EUV radiation and TSI during the declining phase of solar cycle 23 // J. Atmos. Solar-Terr. Phys. 2011. V. 73, iss. 2-3, P. 235-240. DOI:https://doi.org/10.1016/j.jastp.2010.04.002.
19. Maruyama T. Solar proxies pertaining to empirical ionospheric total electron content models // J. Geophys. Res. 2010. V. 115. A04306. DOI:https://doi.org/10.1029/2009JA014890.
20. Nava B., Coisson P., Radicella S.M. A new version of the NeQuick ionosphere electron density model // J. Atmos. Solar-Terr. Phys. 2008. V. 70, iss. 15. P. 1856-1862. DOI:https://doi.org/10.1016/j. jastp.2008.01.015.
21. URL: http://sidc.oma.be/silso/ (дата обращения 1 сентября 2016 г.).
22. URL: ftp://ftp.geolab.nrcan.gc.ca/data/solar_flux/daily_flux_ values/ (дата обращения 1 сентября 2016 г.).
23. URL: http://www.izmiran.ru/ services/iweather/ (дата обращения 1 сентября 2016 г.).
24. URL: ftp://sideshow.jpl.nasa.gov/pub/iono_daily/ (дата обращения 1 сентября 2016 г.).
25. URL: http://irimodel.org/indices/ (дата обращения 1 сентября 2016 г.).
26. URL: http://ftp.izmiran.ru/pub/ izmiran/SPIM/ (дата обращения 1 сентября 2016 г.).
27. URL: http://www.ionolab. org/ (дата обращения 1 сентября 2016 г.).
