В атмосфере Земли постоянно генерируются и развиваются различные крупномасштабные вихревые структуры: циклоны, тайфуны и торнадо. Внутри таких вихрей часто происходят многочисленные вспышки молний. Систематическое дистанционное зондирование подобных крупномасштабных вихревых структур доказывает положительную корреляцию между усилением электромагнитной активности (масштабом разделения зарядов, увеличением напряжённости электрических полей, количеством молниевых разрядов) и последующей интенсификацией крупномасштабного вихря [1–3]. Значит, плазменные вихри, составляющие электромагнитную основу структуры, могут оказывать заметное воздействие на структуру и генезис атмосферных вихрей в целом. К примеру, плазменные вихри взаимодействуют с атмосферными волнами Россби как вихри скорости частиц (т.е. на роторном «генетическом» уровне). Таким образом, определённая доля энергии крупномасштабных вихревых образований в земной атмосфере генерируется и подпитывается плазменными вихрями.
Заметную роль в процессе генерации и усиления плазменных вихрей в атмосфере Земли играет аэрозольная примесь. Аэрозоль, как единый объект, обладает заметно меньшим потенциалом ионизации по сравнению с атомами и молекулами, входящими в его структуру. Тем не менее такая примесь играет ключевую роль в возбуждении плазменных вихрей и последующем накоплении усиливающимися атмосферными вихрями энергии и массы в процессе конденсации атмосферной влаги. С возникающими электрическими полями в атмосферной аэрозольной плазме связано возбуждение и нарастание плазменных вихрей [4–7]. Немонотонная пространственная зависимость коэффициента диэлектрической проницаемости по отношению к стохастическим электростатическим возмущениям детерминированно приводит к усилению плазменных неоднородностей и появлению ячеистых структур в переменных электромагнитных полях, возникающих при разделении зарядов и электрических пробоях. Возникающее расслоение способствует мозаичному распределению аэрозольных подсистем. При неоднородном нагреве ячеистых структур плазмы в атмосфере возбуждаются вихри.
Влияние ионизующего потока энергичных космических частиц на земную атмосферу, погоду и климат описано в [8, 9]. Ионизация аэрозолей высыпающимися заряженными частицами провоцирует генерацию атмосферных вихрей [10–13].
Цель исследования – показать, что на процесс усиления вихревых структур в атмосфере влияют обратные связи. Электромагнитные волны, излучаемые при разрядах молний, при распространении в плазменных волноводах в геомагнитных силовых трубках стимулируют высыпание энергичных заряженных частиц из радиационных поясов Земли, в частности протонов внутреннего радиационного пояса с энергией до 100 МэВ. В каскадах высыпающихся частиц ионизуются аэрозоли. На градиентах давления плазмы, ортогональных геомагнитному полю, возбуждаются плазменные вихри. С нарастанием концентрации загрязнений, в частности аэрозолей, происходит усиление мощности атмосферных вихревых структур и, как следствие, раскачка климатических качелей. Усиливаются торнадо, засухи и наводнения. При образовании волноводных каналов в геомагнитных силовых трубках над атмосферной вихревой структурой, генерирующей разряды молний, вихревая структура влияет на ионизацию аэрозоля также и в магнитосопряжённом районе.
Положительные обратные связи и динамика вихрей
Грозовые разряды связаны с плазменными вихрями. Возбуждение электрических полей зависит от координат немонотонно даже при монотонном исходном распределении плотности аэрозольной плазмы. Относительная скорость движения плазменных неоднородностей влияет на спектры электрических полей, которые приобретают многомодовый характер. Одна из причин образования мощных вихрей в системе подвижных плазменных вихрей связана с проявлением резонансов в спектрах электрических полей взаимодействующих неоднородностей.
Грозовые разряды излучают пакеты электромагнитных волн, проникающие на ионосферные и магнитосферные высоты. В результате электростатические и электромагнитные возмущения в неустойчивой плазме вызывают стимулированное высыпание энергичных заряженных частиц из радиационных поясов Земли. В каскадах высыпающихся энергичных заряженных частиц ионизуются атмосферные аэрозоли. Так возникают положительные обратные связи, которые влияют на процесс усиления вихревых структур в атмосфере и на усиление грозовой активности.
В первом приближении распространение электромагнитных волн в плазменном волноводе можно представить по следующей модели. Для вычисления координатной зависимости амплитуды исследуемой волны, распространяющейся вдоль бесконечного цилиндрического волновода, воспользуемся, например, волновым уравнением для векторного потенциала:
(можно исходить из подобного же уравнения, записанного для векторов электрической или магнитной напряжённости E, H поля волны). При этом вдоль оси волновода компоненты электромагнитного поля будут: в частном случае ТМ-волны магнитное поле H = 0, а в частном случае ТЕ-волны электрическое поле E = 0. В нашем случае можно считать длину образующегося волновода порядка расстояния вдоль силовой трубки между двумя плотными слоями ионосферной плазмы, которые отражают эти распространяющиеся электромагнитные волны (т.е. между магнитосопряжёнными районами). Данная длина, как правило, имеет порядок десятков тысяч километров, т.е. существенно больше длины волны захваченного электромагнитного излучения.
Рассмотрим вначале ТЕ-волну, которая движется вдоль волновода круглого сечения, образованного геомагнитной силовой трубкой, и для простоты пренебрежём её затуханием:
(1)
(2)
где n – нормаль к поверхности волновода. Данная задача для осевого компонента магнитного поля волны H может быть решена методом разделения переменных, который приводит к следующему общему ответу:
, (3)
здесь Jn – функции Бесселя первого рода, коэффициенты αn,m = λn,m / α, где λn,m – m-ный положительный корень дифференциального уравнения dJn / dx = 0, x = rω / c.
Пусть теперь имеем ТМ-волну вдоль того же идеального волновода. Решая аналогичную уравнению (1) задачу для электрического компонента E осевого волнового поля, но вместо (2) – с другим граничным условием: 
аналитическое выражение для E также получается методом разделения переменных. Решение по структуре совпадает с (3), только вместо H пишем поле E, и здесь λn,m – m-ный положительный корень другого уравнения: Jn(x) = 0. Граничные условия определяют порядок величины волнового вектора (в направлении, ортогональном геомагнитному полю) захваченных волн: λ┬ ~ 1 / k┬ ~ a. Поскольку природные волноводы искривлены, то нужно учесть изменение поперечного сечения волновода S и групповой скорости волнового пакета вдоль оси волновода vg. Тогда получим для поля
.
В зависимости от вида функции распределения частиц в ионосферной плазме волны либо усиливаются в плазменном волноводе, либо затухают. При этом, даже при частичном захвате волновых пакетов в плазменный волновод, эти волны детектируются на больших расстояниях, так как мощность электромагнитного сигнала вдоль канала волновода убывает пропорционально первой, а не второй степени r, как в случае волн от точечного источника. В субавроральной зоне магнитное поле Земли способно захватывать поляризационные потоки, например, которые появляются в области с пониженной плазменной плотностью от затухающих электростатических колебаний. Заряженные частицы могут покидать волновод перпендикулярно геомагнитному полю (например, за счёт поляризационного дрейфа). Поперечная диффузия потоков заряженных частиц в волноводных каналах происходит при столкновениях с ионосферными частицами и при рассеянии на флуктуирующих электромагнитных полях. Изменение питч-угла θ рассеянной частицы имеет порядок:
,
где Ω – циклотронная частота, B – амплитуда магнитного поля волны, Δt – итоговое время рассеяния частиц в излучаемых электромагнитных полях (например, от молниевых разрядов). Все виды рассеяния частиц приводят к увеличению их попадания в конус потерь и высыпания в атмосферу (что проявляется, например, в виде полярных сияний).
При поляризационном дрейфе положительно и отрицательно заряженные потоки движутся с разной по модулю скоростью в направлении, перпендикулярном геомагнитному полю. В результате может появиться поперечное электрическое поле (из-за частичного разделения зарядов) и плазменный вихрь. Вихревая структура поддерживает далее существование волноводного канала, увеличивая устойчивость возникшей плазменной неоднородности к диффузионным процессам.
На спутниках детектируются электромагнитные волны от различных источников, в частности от грозовых разрядов. При этом электрические и магнитные компоненты поля электромагнитной волны можно записать в виде комбинаций:
Известны численные и аналитические решения для инкрементов (декрементов) роста волн различной поляризации. Так, в магнитоактивной плазме с учётом взаимодействий волна-частица для геликонов имеем
, (4)
где 
, eσ – собственные векторы следующей задачи Эрмита: 
; 
, 
– тензор проводимости.
Так как собственные векторы взаимно ортогональны [14], то для свистов вместо (4) инкременты можно записать в виде
, (5)
а матричный элемент имеет следующую операторную форму:
,
,
где f(v) – функция распределения частиц по скоростям.
Выберем для электронов следующую функцию v-распределения с потоком:
, (6)
здесь vzb – среднеквадратичная скорость потока вдоль геомагнитного поля (переменная z); средняя скорость направленного потока вдоль z будет Ub; 
– среднеквадратичная скорость в потоке для направления поперёк магнитного поля B0; δ – дельта-функция.
Для фоновой плазмы возьмём двухтемпературное распределение электронов:
, (7)
где α, β – тепловые скорости соответственно для продольного и поперечного направлений (по отношению к геомагнитному полю).
Для введённых в (6) и (7) функций распределения f(v), f(vb) с отношением плотности потока к плотности фоновой плазмы, равным p, из (5) компоненты тензора 
для свистов имеют следующий вид:
,
,
где циклотронная частота для электронов 
, 
, а функции 
– производные от плазменной дисперсионной функции 
, которые можно последовательно найти из дифференциального уравнения: 
– гипергеометрическая функция; 
, а Nb , No – обозначают соответственно плотности электронного потока и фоновой плазмы. Суммирование ведётся по натуральным числам (n > 0).
Известно, что для свистов частота лежит в пределах 
, где 
– плазменная частота. При этом резонансное условие 
отвечает движению волн против потока: kv<0, т.е. свисты эффективно взаимодействуют с набегающим потоком частиц, что приводит к усилению волн. Инкременты роста свистовых волн оказываются на несколько порядков меньше, чем инкременты электростатических колебаний. Из закона дисперсии вистлеров следует ограничение величиной π / 2 на угол между волновым вектором k и групповой скоростью. Для анизотропного распределения по скоростям волны свистового диапазона усиливаются в волноводном потоке электронов в ионосферной плазме.
Следует учитывать, что в области магнитных зеркальных точек продольная скорость захваченных геомагнитным полем частиц близка к нулю. В результате заряженные частицы успевают рассеиваться даже в более медленных полях электростатических возмущений. Высыпание заряженных частиц, захваченных геомагнитным полем, в конус потерь связано с изменением питч-угла частиц при столкновениях с частицами среды, при рассеянии в полях электромагнитных волн и электростатических колебаний в неустойчивой плазме. Мощные пакеты электромагнитных волн от грозовой активности усиливают электростатическую турбулентность в верхней ионосфере, провоцируя возбуждение параметрических взаимодействий волн и высыпание заряженных частиц. Распространение волн сквозь слой неоднородной плазмы можно описать также как туннелирование через потенциальный волновой барьер [15], включая безотражательное прохождение. В каскадах высыпающихся энергичных заряженных частиц радиационных поясов Земли ионизуется атмосферный аэрозоль, что способствует усилению плазменных вихрей и атмосферной вихревой структуры, излучающей электромагнитные волны при грозовых разрядах. Вихри возбуждаются на градиентах давления плазмы, ортогональных геомагнитному полю. Плазменные вихри взаимодействуют с высотными волнами Россби как вихри скорости частиц (на роторном, «генетическом» уровне), приводя к интенсификации атмосферных вихрей. Усиление мощности и устойчивости вихревых структур нелинейным образом связано с ростом концентрации загрязнений в атмосфере. Грозовая активность не только напрямую, но и через обратные связи влияет на усиление вихревых структур. Одна из причин образования мощных вихрей в системе подвижных плазменных вихрей связана с проявлением резонансов в спектрах электрических полей взаимодействующих вихрей. Спектры электростатических возмущений в структуре торнадо имеют многомодовый характер.
Крупномасштабная атмосферная вихревая структура генерирует разряды молний, которые излучают пакеты электромагнитных волн, способные проникать на ионосферные высоты. Если над такой структурой формируются волноводные каналы в геомагнитных силовых трубках, такое влияние усиливается. Электромагнитные моды, излучаемые разрядами молний, при попадании в волноводный канал, стимулируют высыпание энергичных заряженных частиц из радиационных поясов Земли. Ионизация аэрозолей в каскадах высыпающихся частиц приводит к увеличению активности плазменных вихрей в земной атмосфере. В процессе взаимодействия плазменных вихрей с волнами Россби происходит подпитка крупномасштабных вихревых структур атмосферы дополнительной энергией плазменных вихрей. Рост концентрации атмосферных загрязнений коррелирует с интенсификацией таких крупномасштабных вихревых структур, как циклоны (включая тропические), смерчи (торнадо) и антициклоны. Так, в 2021 г. резко увеличилось количество и мощность торнадо и циклонов. Усиливаются засухи и наводнения, раскачиваются климатические качели. Мощные вихревые структуры, генерирующие молнии, влияют на раскачку кризисных атмосферных процессов также и в магнитосопряжённых районах.