Магнитосфера Земли не выдержит экстремальный солнечный шторм: спутниковые замеры 40 лет искажали реальную силу космической погоды

В околоземном космическом пространстве непрерывно идет взаимодействие между магнитным полем Земли и солнечным ветром — потоком заряженных частиц, который выбрасывает Солнце. Десятилетиями в геофизике считалось доказанным фактом, что реакция Земли на этот поток имеет естественный верхний предел. Ученые полагали: если солнечный шторм слаб или умерен, электрические токи в верхних слоях атмосферы растут пропорционально силе воздействия, но при экстремальных событиях этот рост останавливается, выходя на постоянный уровень. Это явление в научной литературе называют «насыщением геомагнитной активности».

Для обоснования этого предела исследователи создали десять различных физических теорий. Все они пытались объяснить, какие именно процессы в магнитосфере заставляют ее сдерживать напор солнечной плазмы.

Однако статья, опубликованная в журнале Nature, показала, что никакого физического предела не существует. Группа исследователей под руководством Нитина Сивадаса доказала: эффект насыщения представляет собой математическую погрешность наблюдений, известную как «регрессия к среднему». Магнитосфера Земли не ограничивает поглощение энергии. При сильнейших вспышках интенсивность геомагнитных возмущений растет линейно, без замедления, а реальное воздействие экстремальных бурь на энергосети и спутники в два раза превосходит величины, заложенные в современные модели безопасности.

Физика взаимодействия: от Солнца до ионосферы

Солнечный ветер представляет собой поток плазмы, состоящий преимущественно из протонов и электронов. Этот поток несет с собой собственное магнитное поле, которое называют межпланетным магнитным полем. Когда направление этого поля оказывается противоположным направлению магнитного поля Земли, происходит процесс магнитного пересоединения: силовые линии солнечного ветра соединяются с силовыми линиями магнитосферы Земли.

В результате пересоединения в околоземном пространстве возникает электрическое поле. Это поле проецируется вдоль магнитных линий вниз, в полярную ионосферу — ионизированный слой атмосферы на высотах от 80 до нескольких сотен километров.

В ионосфере возникают мощные электрические потенциалы и крупномасштабное движение плазмы. Напряжение, возникающее поперек полярной области, ученые называют кросс-полярным потенциалом. Чем выше скорость солнечного ветра и чем сильнее его магнитное поле, тем выше разность потенциалов над полюсами и тем мощнее электрические токи, текущие в ионосфере и магнитосфере.

Для оценки силы этих процессов используют два основных параметра:

1. Электрическое поле пересоединения (Em*). Это показатель интенсивности внешнего воздействия. Его вычисляют на основе измерений скорости солнечного ветра и величины межпланетного магнитного поля.
2. Индекс полярной шапки (PCI). Это показатель реакции самой Земли. Его рассчитывают по данным наземных магнитометров, установленных в Арктике и Антарктике. Индекс напрямую связан с величиной кросс-полярного потенциала и показывает, сколько энергии поступило в систему.

На протяжении многих лет анализ наблюдений приводил к парадоксальному выводу: когда поле солнечного ветра достигает значений выше 10-15 милливольт на метр, индекс полярной шапки перестает расти. График зависимости реакции Земли от силы воздействия изгибался и становился горизонтальным.

Сорок лет теоретических поисков

Начиная с 1980-х годов ученые пытались найти физические причины, из-за которых реакция ионосферы не может расти бесконечно. Появилось десять основных теоретических моделей, объяснявших это явление:

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

Каждая из этих теорий опиралась на строгие физические уравнения, но все они были созданы для объяснения одного и того же феномена: видимого изгиба на графиках наблюдений. При этом авторы теорий исходили из предпосылки, что данные спутниковых измерений точно отражают реальность.

Неопределенность дистанционных измерений

Чтобы понять причину ошибки, необходимо рассмотреть способ сбора исходных данных.

Ученые не могут измерить характеристики солнечного ветра в точке его прямого контакта с магнитосферой Земли, поскольку эта область крайне нестабильна. Вместо этого измерительные приборы размещают на автоматических космических аппаратах, работающих в первой точке Лагранжа (L1) системы Солнце — Земля. Эта точка находится на прямой линии между Солнцем и нашей планетой, на удалении около 1,5 миллиона километров от Земли.

Данные, полученные в точке L1, пересчитывают с учетом времени, которое требуется плазме, чтобы долететь до околоземного пространства. Обычно эта задержка составляет от 40 до 60 минут. Затем добавляют еще около 15-20 минут на то, чтобы возникшее возмущение спустилось по силовым линиям в полярную ионосферу.

В этой схеме измерений возникают три неизбежных источника погрешности:

1. Неточность времени движения плазмы. Поток солнечного ветра не является плоской однородной стеной. Фронт волны искривлен, а скорость движения отдельных сгустков плазмы постоянно меняется. В результате реальное время прибытия плазмы от точки L1 к Земле имеет случайный разброс в пределах десятков минут.
2. Неточность времени отклика ионосферы. Задержка между ударом плазмы о магнитосферу и формированием токов над полюсами непостоянна. Исследования показывают, что она может колебаться от нескольких минут до почти часа.
3. Физическая трансформация потока. Проходя через головную ударную волну и турбулентный слой возле магнитосферы, плазма меняет свои свойства. Направление и сила локального магнитного поля в месте пересоединения отличаются от тех, что были зарегистрированы датчиками спутника за 1,5 миллиона километров.

Прямое сравнение данных со спутников в точке L1 (Em*) с измерениями аппаратов, пролетавших непосредственно возле границы магнитосферы Земли (Em_sh), показало: относительная случайная ошибка составляет не менее 30%. При этом ошибка обладает свойством гетероскедастичности — ее абсолютный разброс увеличивается пропорционально силе самого солнечного ветра.

Механизм статистической иллюзии: регрессия к среднему

Случайная погрешность измерений в науке обычно устраняется усреднением большого массива данных: считается, что случайные отклонения в большую и меньшую сторону компенсируют друг друга. Однако это правило работает только тогда, когда измеряемые величины распределены равномерно или подчиняются обычному нормальному закону распределения.

Параметры солнечного ветра распределены иначе — они подчиняются логнормальному закону. Это означает, что слабые и умеренные значения величины поля регистрируются почти непрерывно, а экстремально высокие показатели случаются крайне редко.

Из-за этого возникает математический эффект, известный как регрессия к среднему (или регрессия к более вероятному значению):

1. Предположим, прибор на спутнике в точке L1 зафиксировал редчайшее значение электрического поля — 25 милливольт на метр.
2. С точки зрения теории вероятностей существуют два объяснения этого замера: либо мимо спутника действительно прошел уникальный сверхмощный поток плазмы, либо прошел гораздо более распространенный поток средней мощности (например, 12 милливольт на метр), но случайный шум и погрешность времени исказили сигнал в сторону резкого завышения.
3. Поскольку умеренные потоки встречаются в тысячи раз чаще экстремальных, вероятность ошибки в сторону завышения во много раз превосходит вероятность истинного экстремума. Следовательно, в подавляющем большинстве случаев реальное воздействие на Землю было существенно слабее, чем показал спутниковый датчик.
4. Магнитосфера Земли реагирует на реальное физическое поле у своей границы (Em_sh), а не на завышенную цифру в протоколе спутника. Соответственно, наземные магнитометры фиксируют умеренный отклик ионосферы.

Когда геофизики строили графики зависимости, они сопоставляли реальный умеренный отклик Земли с ошибочно завышенным показателем спутника. В результате точка на графике неизбежно смещалась вправо и вниз от теоретической прямой.

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

Математическое доказательство и выпрямление графиков

Чтобы строго доказать эту гипотезу, авторы исследования построили статистическую модель ошибок методом Монте-Карло. Они взяли за основу базовое предположение: реакция Земли всегда строго линейна по отношению к реальному магнитному полю.

В модель ввели математические законы распределения ошибок:

• Случайные задержки во времени пролета плазмы (описанные распределением Стьюдента);
• Вариации времени передачи энергии в ионосферу (описанные распределением Вейбулла);
• Случайные изменения напряженности поля в магнитослое, возрастающие с силой бури.

Результат моделирования полностью совпал с экспериментальными данными наблюдений за 25 лет (с 1995 по 2019 год). Модель с искусственно заложенной линейной физикой выдала точно такую же кривую насыщения, какую физики десятилетиями видели на реальных приборах. Искусственное насыщение возникло исключительно из-за структуры статистического шума.

Следующим шагом исследователей стала математическая калибровка реальных данных наблюдений. Авторы применили статистический метод калибровки регрессии, позволяющий вычесть расчетную погрешность из зафиксированных спутником величин входного сигнала.

Как только искажение устранили, экспериментальная кривая распрямилась. Зависимость между интенсивностью солнечного ветра и индексом полярной шапки (PCI) оказалась строго линейной во всем диапазоне надежных наблюдений вплоть до 15 милливольт на метр.

Для дополнительной проверки авторы применили тот же метод к совершенно другому, независимому показателю — индексу западного аврорального электроджета (SML), который рассчитывается на основе отдельной глобальной сети наземных станций SuperMAG. Результат повторился: график зависимости электроджета от калиброванного поля солнечного ветра также стал строго линейным.

Это доказало, что насыщение полярного потенциала не является свойством околоземного пространства. Десять физических моделей, созданных для описания этого предела, объясняли статистическую погрешность.

Практическое значение: пересмотр устойчивости инфраструктуры

Опровержение теории насыщения имеет прямое прикладное значение для технологической инфраструктуры человечества.

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

Расчеты группы Сивадаса показывают:

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

Магнитосфера не защищает Землю от избытка поступающей энергии за счет внутреннего самоограничения. Вся энергия солнечного шторма передается в верхние слои атмосферы пропорционально его мощности. Это требует пересмотра регламентов безопасности для энергетических систем и спутниковых группировок на низкой околоземной орбите.

Системная проблема анализа экстремальных явлений

Результаты работы выходят за рамки геофизики. Исследование показало фундаментальную уязвимость любого научного анализа, в котором редкие экстремальные события изучаются по неточным косвенным измерениям.

Эффект регрессии к среднему неизбежно искажает зависимости, если соблюдаются три условия:

1. Изучаемый физический параметр распределен неравномерно (экстремальные значения редки).
2. Независимая переменная измеряется с погрешностью.
3. Величина погрешности зависит от силы измеряемого сигнала или содержит временную неопределенность.

Эти условия встречаются во многих областях современной науки:

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

Работа ученых демонстрирует: случайный шум в измерениях не исчезает автоматически при увеличении объема выборки. Без строгого математического учета характеристик погрешности измерительных приборов наука рискует принимать артефакты обработки данных за фундаментальные законы природы.

Источник: Nature