Космические исследования, 2019, T. 57, № 6, стр. 451-460

Мелкомасштабные флуктуации плазмы солнечного ветра в быстрых и медленных потоках

М. О. Рязанцева 1*, Л. С. Рахманова 1, Г. Н. Застенкер 1, Ю. И. Ермолаев 1, И. Г. Лодкина 1, Л. С. Чесалин 1

1 Институт космических исследований РАН
Москва, Россия

* E-mail: orearm@gmail.com

Поступила в редакцию 18.04.2019
После доработки 18.04.2019
Принята к публикации 25.05.2019

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

Аннотация

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

ВВЕДЕНИЕ

Солнечный ветер (СВ) не однороден по своей структуре, в нем наблюдаются вариации различных пространственных и временных масштабов, которые в совокупности формируют каскад турбулентных флуктуаций [см. обзор 1]. При этом вариации, ассоциированные с крупномасштабными структурами в короне солнца, являются источниками энергии в каскаде (так называемые интегральные масштабы). Далее, в результате взаимодействия между структурами различных масштабов осуществляется энергетический обмен в системе на так называемых инерционных шкалах каскада (используется также термин магнитогидродинамические (МГД) масштабы). И затем на ионных и субионных шкалах начинают преобладать процессы диссипации энергии, приводящие к нагреву и ускорению частиц в плазме (так называемые диссипативные или кинетические масштабы). При этом согласно теоретическим подходам для развитой турбулентности на интегральных масштабах предсказывается спектр флуктуаций с наклоном –1 [2], на инерционных масштабах – спектр с наклоном –5/3 (так называемая колмогоровская турбулентность) [3], а на диссипативных масштабах – спектр с наклоном –7/3 (согласно теории кинетической альфвеновской турбулентности) [4]. Однако экспериментальные наблюдения показывают, что течение в СВ часто носит не стационарный характер и наблюдающиеся спектры флуктуаций параметров не всегда могут быть описаны в рамках общепринятых теоретических подходов. В частности, как правило, на диссипативных масштабах в эксперименте наблюдаются более крутые, чем это предсказывается в теориях, наклоны спектров флуктуаций ~–2.8 [см. обзор 5]. Близкие показатели наклона спектра –8/3 предсказываются в подходах, учитывающих вклад нелинейного затухания Ландау [6] или влияние перемежаемости структур в плазме [7]. Тем не менее, все наблюдающееся разнообразие спектров флуктуаций объяснить в рамках этих моделей, пока не удается. Кроме того, как правило, в упомянутых выше моделях рассматриваются только альфвеновские флуктуации и предполагается, что флуктуации параметров плазмы являются пассивно связанными с флуктуациями магнитного поля [4, 7], тогда как в эксперименте эта связь не прослеживается однозначно. Показатели наклонов спектров флуктуаций параметров СВ, наблюдающихся в экспериментах, могут сильно варьировать, и соответствуют модельным представлениям лишь в среднем [например 8, 9, 10]. На спектрах плазменных флуктуаций могут наблюдаться различные особенности в окрестности перехода между МГД и кинетическими масштабами [11, 12]. В работе [12] было показано, что спектры флуктуаций потока ионов, могут быть аппроксимированы в соответствии с традиционными теоретическими представлениями линейными зависимостями на МГД и кинетических масштабах с четким изломом между ними не более чем в 50% случаев. При этом часто могут наблюдаться спектры с уплощением в окрестности излома (~20–30%) или с нелинейным укручением на кинетических масштабах (~5–10%). Отличия спектров флуктуаций от вида предсказанного классическими теориями могут появляться вследствие высокого уровня перемежаемости и сложной пространственной конфигурации структур СВ [7, 13, 14]. Свойства флуктуаций потока на малых масштабах зависят от крупномасштабных характеристик потока СВ, в котором они наблюдаются, при этом одним из важнейших факторов считается скорость потока. Так, в работе [15] показано, что в быстрых потоках наблюдается более высокий уровень флуктуаций модуля магнитного поля, при этом спектры флуктуаций спадают быстрее, чем в медленных потоках СВ. В работе [16] исследовались изменения спектров флуктуаций потока ионов в высокоскоростных течениях за фронтами межпланетных ударных волн. Авторами было показано, что за фронтом ударной волны уровень флуктуаций заметно возрастает, а наклон спектра резко укручается. Как правило, при подобном анализе производится селекция на медленные и быстрые потоки СВ по величине направленной скорости, но не уделяется внимание их происхождению. Однако ранее в работе [17] было показано, что характеристики спектров флуктуаций могут в большей степени зависеть от крупномасштабного типа течения СВ чем от локального значения параметров. В данной работе исследуются характеристики спектров флуктуаций СВ в быстрых квазистационарных потоках, и в быстрых возмущенных потоках, ассоциированных с транзиентными явлениями солнечной короны, а также проводится сравнение с аналогичными свойствами в медленных квазистационарных потоках СВ. Анализируются мощность и наклоны спектров флуктуаций потока ионов для широкого диапазона масштабов, а также их зависимость от скорости СВ. Следует отметить, что ранее большинство исследований спектров турбулентных флуктуаций проводились лишь по измерениям межпланетного магнитного поля ввиду ограниченности временного разрешения данных по плазме для большинства космических экспериментов [1, 5]. Отдельные прямые измерения параметров плазмы с высоким временным разрешением были доступны на нескольких КА [1822], но носили не систематический характер. В связи с этим в данной работе особое внимание уделяется ионно-кинетическим масштабам, систематическое исследование которых по измерениям плазмы стало доступно по данным с КА СПЕКТР-Р.

1. ДАННЫЕ И МЕТОДИКА АНАЛИЗА

Анализ, приведенный в работе, основан на плазменных измерениях прибора БМСВ [23, 24] в рамках эксперимента ПЛАЗМА-Ф на борту КА СПЕКТР-Р [25]. Прибор позволяет проводить измерения параметров плазмы в СВ и магнитослое с временным разрешением 31 мс. Это дает возможность систематически проанализировать плазменные флуктуации в широком диапазоне масштабов в том числе на ионно-кинетических масштабах, что ранее не было доступно по данным других КА (см. Введение). В данной работе для исследования мелкомасштабных флуктуаций используются измерения потока ионов СВ доступные с максимальным временным разрешением полное время работы прибора. Полный набор плазменных параметров Np, Vp, Tp, Na/Np доступен с временным разрешением 3 с и определялся из полных энергетических распределений, измеренных прибором БМСВ. В связи с отсутствием на борту КА измерений по магнитному полю используются соответствующие измерения прибора MFI с КА WIND [26] из базы данных CDAWeb (https://cdaweb.sci.gsfc.nasa.gov), сдвинутые на время распространения плазмы между аппаратами. Для анализа были выбраны интервалы в СВ с длительностью не менее 3 ч, и с ориентацией прибора максимально ориентированной на Солнце. Отобранные интервалы были разделены на подинтервалы длительностью ≈17 мин для последующего анализа, причем каждый последующий интервал по времени перекрывал предыдущий на половину его длительности. На каждом интервале вычислялись фурье-спектры флуктуаций потока ионов СВ и определялись их основные характеристики. В работах [9, 27] было показано, что спектры флуктуаций потока ионов, измеренных прибором БМСВ, имеет смысл рассматривать на частотах не более ~10 Гц, так как на больших частотах приборный шум может вносить заметный вклад. Однако, в случаях, где наблюдается низкий уровень потока ионов, соотношение сигнал шум может снижаться, и предельные частоты будут смещаться в область более низких значений. В связи с этим, для подобных случаев в данной работе предельная частота спектра определялась автоматически по характерному уплощению спектра на высоких частотах связанному с повышением уровня шума. Случаи, где предельная частота снижается до значений 4 Гц и менее, исключались из рассмотрения.

Для вычисления наклонов спектров флуктуаций на МГД и кинетических масштабах проводилась линейная аппроксимация спектров в соответствующих диапазонах частот. Как было сказано выше во введении, спектры флуктуаций параметров плазмы в СВ не всегда могут быть аппроксимированы двумя линейными функциями на МГД и кинетических масштабах соответственно, в частности достаточно часто могут наблюдаться спектры с уплощением в окрестности излома между масштабами. В связи с этим, изначально каждый спектр аппроксимировался тремя линейными функциями в предположении наличия уплощения между МГД и кинетическим интервалами. В случае если второй (промежуточный) наклон оказывался достаточно близким к первому (разница показателей наклонов составляла не более 10%) или область уплощения не была ярко выражена (имела ширину менее 0.1 Гц по частотам) то проводилась повторная аппроксимация двумя линейными функциями. При этом область возможной линейной аппроксимации на частотном диапазоне определялась автоматически по минимальному отклонению экспериментальных точек от линейной функции, а условие минимального отклонения от линейных функций должно было соблюдаться сразу для всех трех исследуемых диапазонов частот. Примеры подобной аппроксимации можно видеть ниже на рис. 1, подробно обсуждаемым в разделе 2. По полученным линейным аппроксимациям определялись наклоны спектров на соответствующих масштабах, а также частоты пересечений соответствующих линейных функций. Сравнение мощности спектра флуктуаций PSD (Power Spectral Density) проводилось для удобства в фиксированных диапазонах частот: 0.02–0.08 и 3–4 Гц заведомо лежащих внутри МГД и кинетического масштабов cоответственно.

Рис. 1.

Сравнение спектров флуктуаций потока ионов для интервалов (снизу вверх) в медленном СВ 23.V.2014 02:35–02:52, аппроксимация сплошной толстой черной линией (1), в быстром СВ 2.VIII.2012 18:50–19:07, аппроксимация толстой пунктирной черной линией (2), и в магнитном облаке MC 17.VI.2012 02:16–02:33, аппроксимация тонкой черной линией (3). Величины наклонов спектров на МГД масштабе P1, на кинетическом масштабе P2, и в области уплощения между этими масштабами P0 приведены в табл. 1. Тонкой пунктирной линией изображен характерный наклон т.н. колмогоровского спектра.

Каждому рассмотренному интервалу ставились в соответствие значения средних параметров плазмы и модуля магнитного поля на интервале, а также крупномасштабный тип СВ, в котором он наблюдался. Селекция интервалов по типу течения проводилась с помощью каталога крупномасштабных явлений СВ (ftp://ftp.iki.rssi.ru/pub/omni/) описанного в работе [28]. В соответствие с приведенным каталогом отобранные интервалы были разделены на следующие типы: медленный СВ (квазистационарные медленные потоки из корональных стримеров), быстрый СВ (квазистационарные быстрые потоки из корональных дыр), межпланетные проявления корональных выбросов массы EJECTA и MC (магнитные облака), области компрессии SHEATH перед EJECTA и MC и области компрессии плазмы между высокоскоростными и низкоскоростными потоками CIR. В задачи данной работы входило сравнение характеристик турбулентности в наиболее быстрых потоках СВ, поэтому ниже мы будем рассматривать 3 типа обладающих максимальными средними значениями скорости в отобранных нами интервалах: быстрый СВ, MC и SHEATH перед MC (во всех трех наблюдалась скорость в среднем более 450 км/с). Для сравнения со спокойными условиями также были рассмотрены интервалы в медленном спокойном СВ. В табл. 1 представлено распределение рассмотренных в работе интервалов по крупномасштабным типам СВ, перечисленным выше, а также средние значения параметров СВ и модуля межпланетного магнитного поля в каждом типе. Значения параметров, представленные в таблице, в среднем соответствуют значениям в каталоге крупномасштабных потоков [28], однако наблюдаются и локальные отличия. Эти отличия связаны с тем, что в текущей работе средние параметры определяются на сравнительно коротких интервалах, тогда как в каталоге крупномасштабных явлений усреднения проводятся по длительным многочасовым интервалам за многие годы наблюдений. Так в работе [28] показано, что в течениях СВ ассоциированных с MC наблюдаются скорости в среднем ниже, чем в потоках быстрого СВ из корональных дыр, однако по выборке представленной в данной работе (см. табл. 1) максимальная средняя скорость наблюдается внутри MC. В итоге в рассмотренных выше типах течения СВ в сумме было выделено 1867 интервалов (количество интервалов в каждом типе СВ представлено в табл. 1).

Таблица 1.  

Средние значения параметров плазмы и магнитного поля для рассмотренных интервалов в различных крупномасштабных типах течения СВ

  N, шт. Vp❭, км/с Fp❭, 108 см–2 с–1 Np❭, см–3 Tp❭, эВ ❬|B|❭, нТл βp
Медленный СВ 1092 369 ± 1.5 4.5 ± 0.1 11.4 ± 0.2 6.4 ± 0.2 7.6 ± 0.2 0.80 ± 0.03
Быстрый СВ 295 471 ± 1.5 4.0 ± 0.1 8.1 ± 0.3 12.1 ± 0.6 7.7 ± 0.2 0.95 ± 0.08
MC 289 483 ± 2 4.7 ± 0.2 8.0 ± 0.3 4.9 ± 0.2 8.5 ± 0.3 0.57 ± 0.05
SHEATH MC 191 451 ± 5 7.7 ± 0.2 16.7 ± 0.6 11.8 ± 0.6 9.4 ± 0.2 0.98 ± 0.06

2. СРАВНЕНИЕ ЧАСТОТНО-СПЕКТРАЛЬНЫХ ХАРАКТЕРИСТИК ФЛУКТУАЦИЙ ПОТОКА ИОНОВ НА ПРИМЕРАХ ВЫСОКОСКОРОСТНЫХ ПОТОКОВ РАЗЛИЧНОЙ ПРИРОДЫ

Во введении отмечалось, что характеристики турбулентного каскада могут существенно изменяться в высокоскоростных потоках СВ. Однако ниже будет показано, что эти изменения в большей степени характерны для быстрых потоков, ассоциированных с транзиентными явлениями. На рис. 1 представлены примеры спектров флуктуаций потока ионов для интервалов в медленном спокойном СВ, быстром квазистационарном СВ (ветер из корональных дыр), и быстром потоке внутри магнитного облака MC (см. обозначения в подписи к рис. 1). Для наглядности сравнения были отобраны спектры с уплощением в окрестности излома, достаточно часто наблюдающиеся в СВ [11, 12], и, согласно работе [29], являющиеся результатом комбинации альфвеновской и кинетической альфвеновской турбулентности. В табл. 2 представлены средние значения параметров СВ для рассмотренных примеров на рис. 1, значение RSD относительного стандартного отклонения величины потока ионов на интервалах, а также значения характеристик спектров флуктуаций на соответствующих масштабах: наклоны P1 – на МГД масштабе, P2 – на кинетическом масштабе, P0 – в области уплощения спектра между МГД и кинетическим масштабами; и значения частот излома спектров F1 – частота перехода между линейной частью спектра с наклоном P1 и линейной частью спектра с наклоном P0, F2 – частота перехода между линейной частью спектра с наклоном P0 и линейной частью спектра с наклоном P2.

Таблица 2.  

Значения параметров СВ и спектральных индексов для примеров на рис. 1

  Vp❭, км/с Fp❭, 108
см–2 с–1
Np❭,
см-3
Tp❭, эВ ❬|B|❭,
нТл
❬βp RSD P1 P2 P0 F1❭, Гц F2❭, Гц
23.V.2014 372.7 ± 0.3 2.0 ± 0.05 5.2 ± ± 0.1 5.1 ± ± 0.1 10.7 ± ± 0.1 0.1 ± ± 0.01 3% –1.6 ± 0.3 –2.6 ± 0.1 –1.0 ± 0.2 0.1 ± 0.05 1.5 ± 0.2
2.VIII.2012 440.1 ± 0.1 3.9 ± 0.05 8.4 ± ± 0.1 9.5 ± ± 0.1 10.3 ± ± 0.1 0.31 ± ± 0.02 7% –1.7 ± 0.1 –3.1 ± 0.1 –1.3 ± 0.2 0.4 ± 0.2 1.9 ± 0.3
17.VI.2012 478.9 ± 0.2 13.2 ± 0.1 24.0 ± ± 0.2 6.5 ± ± 0.1 40.2 ± ± 0.2 0.04 ± ± 0.002 12% –1.6 ± 0.2 –3.4 ± 0.2 –0.9 ± 0.1 0.26 ± 0.1 3.3 ± 0.4

На рис. 1 видно, что практически на всем диапазоне масштабов мощность спектра флуктуаций PSD в быстрых потоках СВ в несколько раз превышает мощность спектра в медленном СВ, что соответствует наблюдениям в работах предшественников [15]. Однако, следует отметить, что максимальная спектральная мощность наблюдается в интервале внутри MC, при этом для кинетического масштаба мощность спектра флуктуаций внутри MC более чем на порядок превышает мощность спектра в высокоскоростном потоке из корональной дыры. При этом, как видно из табл. 2, средняя скорость потока в рассматриваемом интервале в MC (3-я строка) не значительно превышает скорость потока в быстром потоке из корональной дыры, также показанном в табл. 2 (2‑я строка). Изменение уровня флуктуаций в целом на рассмотренных интервалах также можно наглядно показать по величине параметра RSD (см. табл. 2). Уровень RSD на интервале внутри MC почти вдвое превышает уровень RSD в быстром СВ, и в 4 раза превышает уровень RSD в медленном СВ. Ниже на статистике в качестве критерия уровня флуктуации мы будем использовать величину спектральной мощности PSD1 на МГД масштабе и PSD2 на кинетическом масштабе (способ вычисления описан в разделе 1) чтобы наглядно сравнить уровень флуктуаций в различных диапазонах частот.

Для всех трех спектров флуктуаций на МГД масштабе показатели наклонов близки и в среднем соответствуют –5/3, характерных для колмогоровского подхода [3]. Далее после первого излома спектра на частоте F1 для всех приведенных примеров наблюдается уплощение спектров (показатели наклонов P0 составляют ~–1), за которым следует второй излом на частоте F2 и снова укручение спектров на кинетических масштабах. При этом наклоны спектров на кинетических масштабах P2 для высокоскоростного потока из корональной дыры немного выше по модулю, чем для медленного спокойного СВ, однако максимальные по модулю наклоны наблюдаются для потока в MC. Причем, можно отметить, что для последнего примера внутри MC наблюдается наиболее длинная область уплощения на частотах 0.2–3 Гц, а затем, на более высоких частотах, происходит резкое укручение спектра (спектр спадает быстрее, чем на использованной линейной аппроксимации). Ранее в работах [17] было показано, что подобное нелинейное укручение спектра на кинетических масштабах сравнительно часто наблюдается внутри MC и может быть связано с увеличением темпов диссипации в высокоскоростных потоках [16]. Следует отметить, что максимальная скорость потока наблюдалась для 3-его интервала в области MC (см. табл. 2). Ниже на статистике будет показано, что вне зависимости от абсолютного значения переносной скорости, в областях, ассоциированных с MC и SHEATH MC, часто наблюдается более высокий уровень флуктуаций с более крутыми наклонами спектров на кинетических масштабах.

3. СТАТИСТИЧЕСКИЕ РАСПРЕДЕЛЕНИЯ СПЕКТРАЛЬНЫХ ХАРАКТЕРИСТИК ФЛУКТУАЦИЙ ПОТОКА ИОНОВ В СОЛНЕЧНОМ ВЕТРЕ

Выше на отдельных примерах была показана разница, наблюдаемая между спектрами флуктуаций потока ионов СВ в высокоскоростном квазистационарном потоке из корональной дыры и высокоскоростном потоке внутри MC в сравнении со спектром в медленном спокойном СВ. В данном разделе представлены статистические распределения описанных выше характеристик спектров флуктуаций для отобранных интервалов наблюдаемых в медленном СВ, быстром СВ (ветер из корональных дыр), и в областях MC и SHEATH MC (см. раздел 1). На рис. 2 показаны распределения для следующих спектральных характеристик: наклоны спектров P1 и P2 на МГД и кинетических масштабах соответственно, спектральные мощности PSD1 и PSD2 на тех же масштабах, наклон P0 для области уплощения спектров между вышеуказанными масштабами (в случаях, когда оно наблюдалось) и частота F2 перехода на кинетические масштабы. Методика определения вышеуказанных параметров описана в разделах 1 и 2. Для спектров с уплощением (например, на рис. 1) частота F2 соответствует излому спектра между областью выполаживания и областью укручения на кинетических интервалах. Для спектров без уплощения, которые могут быть аппроксимированы двумя линейными функциями на МГД и кинетическом масштабе соответственно (например, рис. 1 из работы [12]) частота F2 соответствует точке излома между этими линейными областями. В табл. 3 приведены средние значения описанных выше распределений. Из приведенных распределений и таблицы видно, что не смотря на широкий разброс параметров, ситуация представленная на примерах, описанных в разделе 2 является в целом типичной для всего статистического материала в целом.

Рис. 2.

Распределения характеристик спектров флуктуаций потока ионов для интервалов в медленном СВ (черные бины), в быстром СВ (светло-серые бины), в MC (темно-серые бины), в SHEATH MC (белые бины): наклон спектра на МГД масштабах P1 (а), наклон спектра на кинетических масштабах P2 (б), наклон спектра на переходных масштабах для спектров с уплощением P0 (в), частота перехода на кинетические масштабы F2 (г), десятичный логарифм спектральной мощности на МГД масштабах PSD1 (д), десятичный логарифм спектральной мощности на кинетических масштабах PSD2 (е).

Таблица 3.  

Средние значения спектральных индексов для отобранных интервалов

  P1 P2 P0 F2, Гц PSD1 PSD2
Медленный СВ –1.57 ± 0.01 –2.83 ± 0.01 –1.12 ± 0.02 1.22 ± 0.02 2.8 ± 0.4 (7.1 ± 0.5) · 10–4
Быстрый СВ –1.62 ± 0.02 –2.81 ± 0.02 –1.17 ± 0.03 1.33 ± 0.04 2.7 ± 0.3 (7.4 ± 0.8) · 10–4
MC –1.57 ± 0.02 –3.04 ± 0.02 –1.09 ± 0.02 1.76 ± 0.04 4.5 ± 0.6 (17.4 ± 1.9) · 10–4
SHEATH MC –1.68 ± 0.02 –3.14 ± 0.02 –1.22 ± 0.02 1.89 ± 0.05 6.1 ± 0.7 (20.6 ± 1.5) · 10–4

На МГД масштабах разница между средними наклонами P1 для различных типов СВ не значительная и в среднем они соответствуют –5/3, характерных для Колмогоровского подхода (рис. 2а). Однако разброс значений в соответствующих распределениях достаточно велик. В работе [30] отмечалось, что спектры с большими по модулю показателями наклонов на МГД масштабах (P1 ~ –2) могут наблюдаться в результате анизотропности турбулентного каскада, в случае, когда направление локального магнитного поля совпадает с направлением потока. Более плоские спектры флуктуаций на МГД масштабах в сравнении с модельными спектрами могут наблюдаться в результате локального дисбаланса энергии на границах структур, например подобные уплощения наблюдаются за фронтом околоземной ударной волны [31]. Отметим, что примерно в 10% случаев наблюдаются наклоны с показателями по модулю <1.2. Возможно, это связано с тем, что область уплощения между МГД и кинетическими масштабами может занимать достаточно широкий диапазон частот (см. рис. 1), и, в ряде случаев, частота перехода к МГД масштабам F1 не может быть идентифицирована при рассмотрении сравнительно коротких интервалов. При этом, спектр автоматически будет аппроксимирован как спектр с двумя наклонами с аномально низким значением |P1|.

Распределение наклона P2 на кинетических масштабах (рис. 2б) показывает, что в быстром СВ спектры в целом более крутые чем в медленном СВ, а в MC и SHEATH MC более крутые чем в быстром СВ из корональных дыр. При этом средние значения показателей наклонов P2 (см. табл. 3) для медленного и быстрого СВ близки и соответствуют наклонам спектров флуктуаций магнитного поля, многократно описанным в литературе [см. обзор 5]. Для областей MC и SHEATH MC наблюдаются большие по модулю показатели наклонов спектров, что свидетельствует об ускорении диссипации энергии в этих случаях. Следует отметить, что средние значения |P2| для медленного СВ в данной работе выше чем аналогичные значения в работе [17]. Это связано с тем, что в текущей работе предельная частота, до которой проводилась линейная аппроксимация спектра, определялась автоматически с учетом шумового уплощения спектра на высоких частотах (см. раздел 1). В предыдущих работах предельная частота была фиксирована значением в 10 Гц, что в ряде случаев, с низкими значениями потока ионов (преимущественно в медленном СВ) приводило к занижению показателей наклонов спектра на высоких частотах. Кроме того в данной работе использовался новый алгоритм с изначальной аппроксимацией тремя линейными функциями, что позволило аккуратно определить наклоны спектров и в случаях с уплощением между масштабами, тогда как в предыдущей работе спектры флуктуаций рассматривались в приближении описания двумя линейными функциями, а спектры со значительным уплощением исключались из рассмотрения вручную. Тем не менее, это уточнение не изменяет основных выводов работы [17]. В работе было проведено сравнение наклонов спектров на различных диапазонах высоких частот, и показано что в областях MC и SHEATH MC наблюдается аномальное укручение спектров на высоких частотах, что подтверждается также и результатами, представленными в данной работе.

В распределениях для P0 – наклонов спектров в области уплощения (рис. 2в), не наблюдается статистически значимых отличий для разных крупномасштабных типов потока. Таким образом одинаковая степень уплощения может наблюдаться в любых типах СВ и вне зависимости от крутизны наклонов на МГД и кинетических интервалах спектров. Следует отметить, что спектры удалось автоматически аппроксимировать в предположении трех линейных областей (как в рис. 1) в ~30%, что совпадает с оценками, полученными на меньшей статистике в [12]).

Четкая разница между потоками различной природы видна и в распределении частоты F2 (рис. 2г). Для интервалов наблюдающихся внутри MC и SHEATH MC частота F2 значительно превышает аналогичную частоту для медленных потоков и быстрых потоков из корональных дыр, что, по всей видимости, связано с нелинейным укручением спектра на кинетических масштабах, часто наблюдающимся в этих областях [17]. Как было показано в [11], частота F2 хорошо коррелирует с гироструктурной частотой Fs = V/(2πRT), где V – направленная скорость потока СВ, а RT – тепловой гирорадиус протона. Соответственно вполне закономерно, что F2 больше для высокоскоростных потоков с сильным магнитным полем в областях MC и SHEATH MC (см. табл. 1). Мы не приводим распределения для F1, так как эта частота не может быть определена с достаточной точностью по рассматриваемым относительно коротким интервалам (незначительные отличия в определении наклона на МГД масштабе могут приводить к сильным отличиям в значениях частоты F1).

Наиболее ярко выраженные отличия между спектрами флуктуаций наблюдаемыми в потоках различной природы наблюдаются в распределениях спектральной мощности, отражающей уровень флуктуаций в тех или других диапазонах частот (рис. 2д, 2е). Спектральная мощность как на МГД так и на кинетических масштабах в среднем выше в быстрых потоках, чем в медленных. Однако, также как и в случае с наклонами на кинетических масштабах, максимальная спектральная мощность наблюдается в интервалах, ассоциированных с MC и SHEATH MC.

На рис. 3 представлены зависимости мощности (рис. 3а) и наклона (рис. 3б) спектра на кинетическом масштабе от локальных значений скорости в рассмотренных интервалах (показано точками), а также средние значения этих параметров для фиксированных диапазонов величины направленной скорости в разных крупномасштабных типах течений (линиями). Отметим, что локальные значения величины скорости на рассмотренных 17-минутных интервалах могут отличаться от средних значений внутри длинного интервала соответствующего крупномасштабному течению. Наблюдается широкий разброс значений параметров PSD2 и P2, который не позволяет выявить четкую зависимость этих параметров от величины направленной скорости. Тем не менее, при усреднении для разных крупномасштабных типов СВ в фиксированных диапазонах скоростей (см.обозначения на рис. 3) можно видеть, что внутри каждого диапазона скорости максимальная мощность спектра наблюдается для областей SHEATH MC (рис. 3а). Так, для значений скоростей 300–400 км/с мощность спектра для интервалов внутри SHEATH MC примерно в 3–4 раза больше мощности спектра для интервалов в медленном спокойном СВ, а для значений скоростей >400 км/с в 4–8 раз больше мощности спектра для быстрого СВ из корональных дыр. Мощность спектров для MC как правило также в несколько раз превышает мощность спектров для быстрого СВ, однако она меньше чем для SHEATH MC в тех же диапазонах скоростей. Аналогичная ситуация наблюдается и для наклонов спектров на кинетическом масштабе (рис. 3б). Для SHEATH MC наклон спектра по модулю больше наклонов для медленных и быстрых потоков в соответствующих диапазонах скоростей. Наклон спектра в MC больше чем в быстром СВ, но меньше чем в SHEATH MC. Следует отметить, что разделение на отдельные группы по величине скорости, приводит к уменьшению статистики для каждой группы и к возрастанию стандартных ошибок среднего значения. При этом ошибки средних приведенных в ступенчатых диаграммах не превышают 2–3% для показателей наклона спектра, но могут достигать 10–30% для значений спектральных мощностей. Тем не менее обсуждаемая выше разница между средними значениями обеих величин в различных типах СВ значительно превышает выше указанные ошибки средних значений, а приведенные результаты явно показывают, что определяющую роль играет источник происхождения рассматриваемого потока, а не непосредственно величина скорости в потоке.

Рис. 3.

Зависимость спектральной мощности PSD2 (а) и показателя наклона спектра P2 (б) на кинетическом масштабе от величины направленной, скорости на интервале. Усреднения для различных типов СВ в фиксированных диапазонах скоростей показаны ступенчатой диаграммой: черная тонкая сплошная линия – медленный СВ, серая пунктирная – быстрый СВ, черная пунктирная – MC, серая толстая сплошная – SHEATH MC.

ОБСУЖДЕНИЕ И ВЫВОДЫ

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

В работе проведено сравнение спектральных характеристик флуктуаций потока ионов по данным прибора БМСВ в высокоскоростных потоках различного происхождения. Выявлено, что средняя мощность спектров флуктуаций для MC и SHEATH MC в 2–3 раза выше чем в медленном и быстром квазистационарном потоке СВ на всем диапазоне частот, тогда как мощность спектров для медленного и быстрого СВ практически совпадает. Наклоны спектров по модулю также больше для MC и SHEATH MС на кинетических масштабах, при этом на МГД масштабах значимых отличий не обнаружено. Максимальная мощность и наиболее крутые наклоны спектров на кинетических масштабах наблюдаются в области компрессии плазмы SHEATH MC, как для медленных, так и для быстрых потоков. Вероятно, существует определенная связь между возрастанием спектральной мощности и укручением спектра на высоких частотах, как это было показано в работах предшественников по флуктуациям магнитного поля [8, 15]. Таким образом, сдвиг скорости на границах различных нестационарных структур СВ может оказывать более важную роль в формировании турбулентных характеристик потока, чем сама по себе величина скорости. Можно заключить, что, в областях MC и SHEATH MC возрастает уровень плазменных флуктуаций, что способствует возникновению нестационарности в турбулентном каскаде и ускорению диссипации энергии, тогда как быстрый квазистационарный СВ, ассоциированный с течением из корональных дыр, не меняет значимым образом свойств каскада.

Авторы благодарны своим коллегам в ИКИ РАН, в НПО им. С.А. Лавочкина и в Карловом Университете (Прага, Чехия) за помощь в разработке, отладке, калибровке, полетном контроле, а также за сбор, передачу и первичную обработку научной информации. Работа авторов поддержана грантом РНФ 16-12-10062.

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

  1. Bruno R., Carbone V. The solar wind as a turbulence laboratory // Living Reviews in Solar Physics. 2013. V. 10. № 2. https://doi.org/10.12942/lrsp-2013-2

  2. Matthaeus W.H., Goldstein M.L. Low-frequency 1/f noise in the interplanetary magnetic field // Physical Review Letters. 1986. V. 57. P. 495–498.

  3. Frisch U. Turbulence. Cambridge University Press, 1995. P. 312.

  4. Schekochihin A.A., Cowley S.C., Dorland W. et al. Astrophysical gyrokinetics: kinetic and fluid turbulent cascades in magnetized weakly collisional plasmas // Astrophysical J. Supplement Series. 2009. V. 182. P. 310–377.

  5. Alexandrova O., Chen C.H.K., Sorisso-Valvo L. et al. Solar wind turbulence and the role of ion instabilities // Space Science Reviews. 2013. V. 178. P. 101–139. https://doi.org/10.1007/s11214-013-0004-8

  6. Howes G.G., Cowley S.C., Dorland W. et al. A model of turbulence in magnetized plasmas: implications for the dissipation range in the solar wind // J. Geophysical Research. 2008. V. 113. № A12. P. 5103. https://doi.org/10.1029/2007JA012665

  7. Boldyrev S., Perez J.C. Spectrum of kinetic-Alfvén turbulence // Astrophysical J. 2012. V. 758. P. 44. https://doi.org/10.1088/2041-8205/758/2/L44

  8. Smith C.W., Hamilton K., Vasquez B.J., Leamon R.J. Dependence of the Dissipation Range Spectrum of Interplanetary Magnetic Fluctuationson the Rate of Energy Cascade // Astrophysical J. 2006. V. 645. № 1. L85–L88. https://doi.org/10.1086/506151

  9. Safrankova J., Nemecek Z., Prech L., Zastenker G. Ion kinetic scale in the solar wind observed // Physical Review Letters. 2013. V. 110. 025004. https://doi.org/10.1103/PhysRevLett

  10. Riazantseva M.O., Budaev V.P., Zelenyi L.M. et al. Dynamic properties of small scale solar wind plasma fluctuations // Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences. 2015. A. V. 373. 20140146. https://doi.org/10.1098/rsta.2014.0146

  11. Safrankova J., Nemecek Z., Nemec F. et al. Solar wind density spectra around the ion spectral break // Astrophysical J. 2015. V. 803. P. 107. https://doi.org/10.1088/0004-637X/803/2/107

  12. Рязанцева М.О., Рахманова Л.С., Застенкер Г.Н., Ермолаев Ю.И. Типы спектров флуктуаций потока ионов в солнечном ветре и магнитослое на стыке инерционного и диссипативного интервалов // Геомагнетизм и аэрономия. 2017. Т. 57. № 1. С. 3–9. https://doi.org/10.1134/S001679321701011X

  13. Зеленый Л.М., Милованов А.В. Фрактальная топология и странная кинетика: от теории перколяции к проблемам космической электродинамики // УФН. 2004. Т. 174. № 8. С. 809–852. https://doi.org/10.1070/PU2004v047n08ABEH001705

  14. Budaev V.P., Zelenyi L.M., Savin S.P. Generalized self-similarity of intermittent plasma turbulence in space and laboratory plasmas // J. Plasma Physics. 2015. V. 81. 395810602. https://doi.org/10.1017/S0022377815001099

  15. Bruno R., Trenchi L., Telloni D. Spectral Slope Variation at Proton Scales from Fast to Slow Solar Wind // Astrophysical J. Letters. 2014. V. 793. № 1. L15. P. 5. https://doi.org/10.1088/2041-8205/793/1/L15

  16. Pitna A., Safrankova J., Nemcek Z. et al. Density fluctuations upstream and downstream of interplanetary shocks // Astrophysical J. 2016. V. 819. P. 41. https://doi.org/10.3847/0004-637X/819/1/41

  17. Рязанцева М.О., Рахманова Л.С., Застенкер Г.Н. и др. Особенности спектральных характеристик плазменных флуктуаций в различных крупномасштабных потоках солнечного ветра // Геомагнетизм и аэрономия. 2019. Т. 59. № 2. С. 139–147.

  18. Unti T.W.J., Neugebauer M., Goldstein B.E. Direct measurements of solar-wind fluctuations between 0.0048 and 13.3 Hz // Astrophysical J. 1973. V. 180. P. 591–598. https://doi.org/10.1086/151987

  19. Celnikier L.M., Muschietti L., Goldman M.V. Aspects of interplanetary plasma turbulence // Astronomy and Astrophysics. 1987. V. 181. P. 138–154.

  20. Vaisberg O.L., Smirnov V.N. Further analysis of non-linear density fluctuations in the foot of quasi-perpendicular shock // Advances in Space Research. 1995. V. 15. № 8/9. P. 297.

  21. Kellogg P.J., Horbury T.S. Rapid density fluctuations in the solar wind // Annales Geophysicae. 2005. V. 23. № 12. P. 3765–3773. https://doi.org/10.5194/angeo-23-3765-2005

  22. Shevyrev N., Zastenker G. Some features of plasma flow in the magnetosheath behind the quasi-parallel and quasi-perpendicular bow shocks // Planetary and Space Science. 2005. V. 53. P. 95–102. https://doi.org/10.1016/j.pss.2004.09.033

  23. Застенкер Г.Н., Шафранкова Я., Немечек З. и др. Быстрые измерения параметров солнечного ветра с помощью прибора БМСВ // Космич. исслед. 2013. Т. 51. № 2. С. 88–99. (Cosmic Research. P. 78.) https://doi.org/10.1134/S0010952515010098

  24. Safrankova J., Nemecek Z., Prech L. et al. Fast solar wind monitor (BMSW): description and first results // Space Science Reviews. 2013. V. 175. № 1–4. P. 165–182. https://doi.org/10.1007/s11214-013-9979-4

  25. Зеленый Л.М., Застенкер Г.Н., Петрукович А.А. и др. Эксперимент ПЛАЗМА-Ф на спутнике “Спектр-Р” // Космич. исслед. 2013. Т. 51. № 2. С. 83–87. (Cosmic Researc. P. 73.)

  26. Lepping R.P., Acuna M.H., Burlaga L.F. et al. The WIND magnetic field investigation // Space Science Reviews. 1995. V. 71. P. 207. https://doi.org/10.1007/BF00751330

  27. Chen C.H.K., Sorriso-Valvo L., Safrankova J., Nemecek Z. Intermittency of solar wind density fluctuations from ion to electron scales // Astrophysical J. Letters. 2014. V. 789 (L8). P. 5. https://doi.org/10.1088/2041-8205/789/1/L8

  28. Ермолаев Ю.И., Николаева Н.С., Лодкина И.Г., Ермолаев М.Ю. Каталог крупномасштабных явлений солнечного ветра для периода 1976–2000 // Космич. исслед. 2009. Т. 47 № 2. С. 99–113. (Cosmic Research. P. 81.) https://doi.org/10.1134/S0010952509020014

  29. Chandran B.D.G., Quataert E., Howes G.G. et al. Constraining low-frequency alfvénic turbulence in the solar wind using density-fluctuation measurements // Astrophysical J. 2009. V. 707. № 2. P. 1668–1675. https://doi.org/10.1088/0004-637X/707/2/1668

  30. Horbury T.S., Forman M., Oughton S. Anisotropic scaling of magnetohydrodynamic turbulence // Physical Review Letters. 2008. V. 101. 175005. https://doi.org/10.1103/PhysRevLett.101.175005

  31. Rakhmanova L., Riazantseva M., Zastenker G., Verigin M. Kinetic-scale ion flux fluctuations behind the quasi-parallel and quasi-perpendicular bow shock // J. Geophysical Research: Space Physics. 2018. V. 123. P. 5300–5314. https://doi.org/10.1029/2018JA025179

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