Астрономический вестник, 2020, T. 54, № 5, стр. 468-474

Радиоизлучение метеорных тел в диапазоне сверхнизких частот

А. Д. Филоненко *

Луганский государственный университет им. Владимира Даля
Луганск, @

* E-mail: uy5lo@mail.ru

Поступила в редакцию 03.04.2018
После доработки 28.01.2020
Принята к публикации 13.02.2020

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

Аннотация

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

Ключевые слова: метеор, болид, радиоэмиссия, спектр, сродство к электрону, VLF

ВВЕДЕНИЕ

Проникновение в земную атмосферу метеорного вещества, космической пыли, а также космических лучей (в основном ядер различных атомов) вызывает множество явлений, в том числе и радиоэмиссию в широком диапазоне частот. Например, механизм радиоизлучения заряженных частиц, движущихся с около световыми скоростями известен относительно давно. Около 60 лет назад опубликована работа (Аскарьян, 1961), в которой впервые был описан механизм излучения широкого атмосферного ливня, суть которого связана с когерентным черенковским излучением избыточных электронов. В дальнейшем эта идея была дополнена и расширена во множестве работ применительно к реальным ливням (Филоненко, 2015).

Электрические явления, вызванные быстро движущимися метеорами, далеко не всегда имеют определенное объяснение. Прямыми измерениями было установлено, что прохождение метеорных тел (МТ) через атмосферу Земли сопровождается электромагнитным излучением. Кроме известного всем излучения в оптическом диапазоне, МТ излучает и в радиоволновом диапазоне. Относительно недавно опубликована работа (Obenberger и др., 1405.6772), в которой сообщается об обнаружении собственного излучения метеоров декаметровом радиотелескопом LWA1. Причем для многих событий установлена корреляция с оптическим излучением, которое регистрировалось с помощью видеокамер сети кругового обзора небесной сферы (AllSkyFireballNetwork). Тщательный анализ радиосигналов показал, что они не могут быть отражением антропогенных радиоимпульсов от ионизованного следа метеора, или какими-либо другими видами радиопомех, а являются только результатом собственного излучения МТ (Obenberger и др., 1503.05150; 1510.04147; 1607.03002; 1609.04701). Модель собственного излучения метеоров была предложена в (Филоненко, 2018). Рассчитанные на основе ее характеристики излучения хорошо совпадают с экспериментальными данными (Obenberger и др., 1503.05150; 1510.04147; 1607.03002; 1609.04701).

Не менее интересное явление было описано в работе (Мусатенко и др., 2006). С помощью радиометра высокого разрешения в спектрах флуктуации радиошумов ионосферной плазмы на длине волны 2 метра обнаружена линия, которая проявляется во время высокоскоростных метеорных потоков. Частота линии колеблется в пределах 12–60 Гц в зависимости от активности источника ионизации ионосферы, локальных характеристик метеорного потока, времени суток и пр. Математическая модель этого эффекта предложена в работе (Kopnin, Popel, 2009), в которой предполагается, что в районе ионосферы во время метеорного потока образуется т.н. “металлический слой” толщиной несколько километров, состоящий из микрометровых частиц, образованных в процессе конденсации паров металлов, входящих в состав метеоров. Эти частицы вследствие различных причин приобретают электрический заряд. При прохождении радиоизлучения ионосферы через такой слой пылевой плазмы возникает взаимодействие, приводящее к модуляции ионосферного радиошума на характерных частотах близких к пылевым акустическим волнам.

Давно было замечено, что падение крупных метеоритов оказывает сильное воздействие на телефонные сети и даже на городские сети электрического освещения. Известны также случаи, когда падающие яркие болиды изменяли на несколько секунд естественную ориентацию магнитной стрелки компаса (Астапович, 1958). Из этих наблюдений вытекало предположение об образовании статического заряда и, как следствие, магнитного поля при движении МТ в атмосфере Земли. Классическое толкование этого явления приводит к результатам, несовместимым с реальными возможностями. Магнитную индукцию, наводимую таким зарядом в точке наблюдения, можно считать порядка поля земного магнетизма, т.е. ~10−5 Тл. Для расстояния r от наблюдателя до болида порядка нескольких десятков километров и характерной скорости v ~ 105, м/с, элементарная оценка индукции поля B ≈ (μ0/4π)qv/r2 дает величину заряда q ~ 106, Кл. Это соответствует напряженности электрического поля Е на поверхности болида E ~ 1018, В/м, что в принципе, не осуществимо и свидетельствует о более сложных процессах, чем элементарные представления.

Не менее загадочное явление – это электрофонные болиды. Суть его заключается в том, что наблюдатель практически одновременно видит свет и слышит звук, производимый летящим болидом (см., например, обзоры Vinkovic и др., 0211203; Cepleha и др., 1998). Поскольку характерное расстояние до болида составляет несколько десятков километров и более, то, очевидно, что звук можно услышать только через полторы-две минуты, когда фактически болида уже нет. Крупные метеориты – это достаточно редкое явление. Поэтому даже на сегодняшний день нет достаточно полноценной статистики синхронного измерений радио, оптического и акустического излучения и существование электрофонных болидов вызывает естественное сомнение у ряда исследователей.

Не найдено пока объяснение и процессам, приводящим к излучению в килогерцовом диапазоне частот (VLF). Этот диапазон насыщен сильными радиопомехами естественного и антропогенного происхождения и для определенных выводов недостаточно регистрации только одних радиоимпульсов. В настоящее время известен ряд работ, в которых одновременного записаны радиосигналы и видеонаблюдения.

Например, в работе (Beech, Brown и др., 1995) описаны результаты 80-часового визуального/видео наблюдения и одновременно с этим VLF излучения, проведенного во время метеорных потоков Лириды, Персеиды, Ориониды, Леониды и Геминиды. Только один болид из записанных событий, сообщают авторы, произвел ясный и убедительный VLF сигнал. Это был яркий болид, с визуальной величиной М = −10m(±1), у которого заключительная стадия полета сопровождалась оптической вспышкой, и одновременно с этим одним из авторов и двумя другими наблюдателями было отмечено “жужжание”, исходящее из аудиовыхода VLF приемника. В этой работе приводится осциллограмма радиосигнала, и частотный спектр, из которого следует, что характерная частота спектра составляет ~2 кГц. Исходя из величины ЭДС индукции Ui ≈ 1 B, наведенной в контуре с 360 витками и площадью S = 1 м2, они находят, что индукция магнитного поля равна В = = 8 × 10−6 Тл. Откуда авторы получают, что напряженность электрического поля в точке наблюдения имеет величину 2000 В/м.

В ноябре 1998 г. была проведена экспедиция хорватского физического общества в Монголию (Garaj, Vinkovi и др., 1999). Цель состояла в измерении электрофонных звуков и радиоизлучения очень низких частот от метеоров во время ожидаемого метеорного потока Леониды. Для визуальной записи метеоров были использованы ПЗС-камеры высокой чувствительности с покрытием неба на 55%. Для низкочастотных радиоприемников использовались вертикальные штыревые антенны. В целом экспериментальная установка могла записывать звуки, изображения и ELF/VLF радиосигналы синхронно. Точность синхронизации между различными каналами (аудио, видео, радио) была не менее ±0.02 с. Для уменьшения техногенных радиопомех наблюдения проводились в ненаселенной области в 20 км к юго-востоку от Улан-Батора.

В ходе эксперимента были получены записи, не вызывающие сомнения в существование временной корреляцией между видео и радио ELF/VLF сигналами (см. рис. 3 в (Garaj, Vinkovi и др., 1999)). Однако электрофонные звуки в этих событиях не были обнаружены. Характерный радиосигнал представлял собой группу коротких импульсов, хаотически распределенных в интервале нескольких десятых долей секунды, т.е. фактически такого же вида, как и в (Beech, Brown и др., 1995). Наиболее густое скопление импульсов соответствовало максимуму яркости метеора. Авторы (Garaj, Vinkovi и др., 1999) нисколько не сомневаются в том, что их источником является одно и то же явление.

Через год после этой экспедиции группа физиков из университета Тель-Авива провела исследования с целью проверки гипотезы об электрофонных звуках (Price, Blum, 2000). Измерения низкочастотных электромагнитных волн были выполнены во время метеорного потока Леониды в ночь на 18 ноября 1999 года. Эксперимент проводился на полигоне для наблюдения ELF/VLF сигналов в Исследовательском институте университета Бен-Гуриона. Петлевые антенны ELF/VLF высотой 10 метров были выполнены в виде двух ортогональных треугольников с геометрической площадью около 81 м2 у каждой петли. Одна из них установлена в направлении север–юг, а другая восток-запад. Чувствительность системы в широкополосном диапазоне (0.1–50, кГц) составляла 6, мкВ/м. Динамический диапазон канала антенна−предусилитель составлял 100 дБ.

Измерения показали, что метеорные спектры имеют минимум вблизи 5 кГц и большой максимум в диапазоне ELF (0.3–1.5, кГц), а также дополнительный более слабый максимум около 2, кГц. Из представленных данных в работе (Price, Blum, 2000) следует, что напряженность электрического поля в месте наблюдения не превышает нескольких сотен микровольт на метр. Результаты эксперимента, по мнению авторов (Price, Blum, 2000), дают убедительное доказательство того, что метеоры производят электромагнитное излучение при вхождении в атмосферу, что может объяснить звуки, слышимые во время наблюдений за большими метеорными телами.

Имея в виду единственное к тому времени теоретическое объяснение того, как производятся эти радиоволны (Keay, 1993; 1995), авторы (Price, Blum, 2000) приходят к выводу, что существующая теория не объясняет каким образом, едва видимые и мелкие метеоры производят радиосигналы. Кроме того, из работ (Keay, 1993; 1995) следует, что длительность радиоимпульсов должна составлять несколько секунд, а полученные ими экспериментальные значения не превышают 10 миллисекунд, причем весь этот сигнал представляет собой совокупность очень коротких (доли миллисекунды) импульсов.

Для описания электромагнитных явлений, сопровождающих вторжение метеороида в атмосферу Земли, необходимо чтобы физическая модель генерации VLF-излучения могла воспроизвести спектр импульса, его характер и спектральную интенсивность. С этой целью авторы (Beech, Foschini, 1999) привлекают для своей гипотезы явление ударной волны, возникающей перед телом, движущимся в атмосфере с гиперзвуковыми скоростями. Наличия большого градиента давления, температуры и других параметров в ударной волне, способствует сильной диффузии электронного газа по отношению к ионному газу и возникновению объемных зарядов. Чтобы оценить порядок величины электрического поля, генерируемого ударной волной, авторы (Beech, Foschini, 1999) ссылаются на результаты, опубликованные в (Зельдович, Райзер, 1996), однако не учитывают тот факт, что высокая напряженность электрического поля существует только в области двойного электрического слоя с толщиной порядка величины дебаевского радиуса. Поэтому в точке нахождения наблюдателя это поле будет чрезвычайно мало.

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

Хорошо известно, что некоторые метеоры после достижения максимума свечения разрывается на мелкие фрагменты с визуально наблюдаемой сильной вспышкой. Такое завершение события не является исключительным явлением и подтверждается многими наблюдениями. Это может случиться вследствие разогрева тела и его саморазрушения под действием механических сил со стороны набегающего потока. Такое явление характерно для высокоскоростных метеоров. Вспышка в терминальной точке может сопровождаться взрывоми распадом на многие мелкие фрагменты (Ceplehaetal, 1998). Например. в работах (Егорова, Лохин, 2015; Егорова, Лохин 2017) авторы показали, что разрушение метеороида на множество мелких фрагментов подобно разрушению твердого тела при взрыве и описали модель, заключительная стадия которой описывает наблюдаемый эффект “теплового взрыва” метеороида.

ЭЛЕКТРИЗАЦИЯ МЕТЕОРОВ В АТМОСФЕРЕ ЗЕМЛИ

Заряжение метеорного тела в атмосфере земли происходит в процессе соударения набегающих молекул воздуха с поверхностью МТ и последующего захвата электрона атомом атмосферного кислорода. Такой процесс происходит с выделением энергии, которую называют сродством к электрону. Отрицательный ион кислорода остается в атмосфере, а МТ заряжается положительно. Этот механизм заряжения лежит в основе явления “электризации трением” (см., например, (Лёб, 1963)), при котором тесное соприкосновение двух различных тел приводит к переходу электронов с одного тела на другое. Возникший на границе двух тел так называемый двойной электрический слой практически не создает внешнего поля, однако если тела раздвинуть, то на каждом из них останутся заряды противоположного знака. Разность потенциалов между этими телами может достигать несколько киловольт.

На больших высотах (60–70 км) длина свободного пробега молекул атмосферы близка к размерам самого фрагмента метеороида, поэтому при соударении c ним ион кислорода не остается на поверхности, а отразившись от нее, может приобрести скорость, близкую к скорости самого фрагмента.

Фактически элементарный процесс излучения при таком взаимодействии начинается с момента отрыва электрона. Однако молекула кислорода быстро останавливается при столкновениях с молекулами окружающего воздуха и практически не излучает, а продолжающееся замедленное движение осколка и его последующее заряжение будут определять низкочастотную часть спектра. Для этого нужно чтобы его размеры были такой величины, при которой время торможений было порядка миллисекунды. Тогда основная часть спектра будет лежать в килогерцовом диапазоне частот.

Наиболее важный факт в этом процессе заключается в том, что излучение фрагментов метеороида будет квазикогерентным. Фактически это означает, что если бы осколки разлетались в разные стороны, то их излучение не является строго когерентным, а приведенные далее оценки являются оценками сверху. Более того, если представить себе, что сферически симметрично взорвался покоящийся заряженный метеороид, то, как хорошо известно, излучение в этом случае отсутствовало бы. Однако скорость разлета фрагментов намного меньше скорости метеороида (${{v}_{0}}\sim 5 \times {{10}^{4}}\,\,{{\text{м}} \mathord{\left/ {\vphantom {{\text{м}} {\text{с}}}} \right. \kern-0em} {\text{с}}}$), поэтому направление движения осколков будет в основном совпадать с направлением движения метеороида и это означает (как было сказано выше), что излучение будет квазикогерентным. Разумеется, для этого вывода важен тот факт, что время разрушения метеороида при взрыве порядка миллисекунды. В конечном итоге это означает, что верхняя оценка мощности излучения будет пропорциональна квадрату суммарного заряда фрагментов.

Для метеороида с массой несколько килограммов режим свободномолекулярного обтекания переходит в режим сплошной среды на высоте около 80–70, км. На высотах 70–50, км и ниже вследствие образования ударной волны следует ожидать “теплового взрыва” (Егорова, Лохин, 2015; Егорова, Лохин, 2017), при котором метеороид разрушится на мелкие осколки. Однако и для них на этой высоте может сохранится режим свободномолекулярного обтекания, если их характерный размер близок к длине свободного пробега молекул атмосферы (λ = 10–3 м на высоте 70 км).

Рассмотрим эту часть модели более подробно. Для упрощения оценок будем считать, что оторвавшийся осколок имеет форму, близкую к сферической. Под действием силы торможения $~F = \Gamma S{{\rho }_{a}}{{v}^{2}}$ он движется замедленно. Здесь v – скорость фрагмента, Γ ~ 1 – аэродинамический коэффициент сопротивления, ρа – плотность атмосферы, S – площадь поперечного сечения, и m – его масса. В этом случае состояния осколка можно описать уравнением движения $m\frac{{dv}}{{dt}} = - \Gamma S{{\rho }_{a}}{{v}^{2}}$ и уравнением потери массы $\frac{{dm}}{{dt}} = - \Lambda \frac{{S{{\rho }_{a}}{{v}^{3}}}}{{2Q}}$, где Λ – коэффициент теплопередачи равный или меньший единице, так как энергия, идущая на абляцию фрагмента, не превышает кинетической энергии набегающего потока молекул и Q – удельная теплота плавления или испарения материала фрагмента. Отношение Λ/Г2Q обозначают обычно, как σ, характерное значение которого принимают равным $\sigma = 2 \times {{10}^{{ - 8}}}$, в единицах с22. Эти уравнения дают решения, которые можно записать в виде (см., например, Бронштэн, 1981):

1. зависимость характерного размера фрагмента от начальной v0 и конечной v скоростей

(1)
$r\left( v \right) = {{r}_{0}}{{e}^{{\frac{\sigma }{6}\left( {{{v}^{2}} - v_{0}^{2}} \right)}}},$

2. интервал времени, затраченный для этого

(2)
$t\left( v \right) = - \frac{{4{{r}_{0}}{{\rho }_{m}}}}{{3{{\rho }_{a}}\Gamma }}\int\limits_{{{v}_{0}}}^v {\frac{{{{e}^{{\frac{\sigma }{6}\left( {{{x}^{2}} - x_{0}^{2}} \right)}}}}}{{{{x}^{2}}}}} dx,$

3. пройденный путь за это время

(3)
$z\left( v \right) = - \frac{{4{{r}_{0}}{{\rho }_{m}}}}{{3{{\rho }_{a}}\Gamma }}\int\limits_{{{v}_{0}}}^v {\frac{{{{e}^{{\frac{\sigma }{6}\left( {{{x}^{2}} - x_{0}^{2}} \right)}}}}}{x}} dx,$

4. величина мгновенного ускорения, как функция скорости

(4)
$a\left( v \right) = - \frac{{3{{\rho }_{a}}}}{{4{{r}_{0}}{{\rho }_{m}}}}{{v}^{2}}{{e}^{{ - ~\frac{\sigma }{6}\left( {{{v}^{2}} - v_{0}^{2}} \right)}}}.$

На рис. 1 графически представлены все эти зависимости для выбранных значений следующих величин: ${{v}_{0}} = 5 \times {{10}^{4}}~\,\,{{\text{м}} \mathord{\left/ {\vphantom {{\text{м}} {\text{с}}}} \right. \kern-0em} {\text{с}}}$ – начальная скорость фрагмента, плотность атмосферы на высоте 70 км ${{\rho }_{a}} = 0.8 \times {{10}^{{ - 4}}}~\,\,{{{\text{кг}}} \mathord{\left/ {\vphantom {{{\text{кг}}} {{{{\text{м}}}^{3}}}}} \right. \kern-0em} {{{{\text{м}}}^{3}}}}$, ${{r}_{0}} = {{10}^{{ - 4}}}\,\,{\text{м}}$, Γ = 1 – характерная величина коэффициента сопротивления. Размер фрагмента r(v) изменяется очень быстро (1). Например, в момент достижения максимального ускорения (4) радиус фрагмента равен приблизительно $r = 6.4 \times {{10}^{{ - 8}}}{\text{\;м}}$, что соответствует скорости $v = 1.76 \times {{10}^{4}}\,\,~{{\text{м}} \mathord{\left/ {\vphantom {{\text{м}} {\text{с}}}} \right. \kern-0em} {\text{с}}}$ в момент времени (2) $t = 0.77 \times {{10}^{{ - 3}}}~\,\,{\text{с}}$.

Рис. 1.

Зависимость характерных величин фрагмента с начальным размером ${{r}_{0}} = {{10}^{{ - 4}}}\,\,~{\text{м}}$ от скорости.

Средняя длина свободного пробега молекулы на характерной высоте ~60–70 км равна приблизительно λ = 10–4 метров и сравнима с размером предполагаемого осколка метеороида. По мере увеличения заряда Q(v) этого фрагмента при соударении с атомом кислорода его электрическое поле препятствует выходу иона из зоны притяжения. Поэтому предельная величина заряда Q(v) определится условием равенства энергии ${{qQ\left( v \right)} \mathord{\left/ {\vphantom {{qQ\left( v \right)} {4\pi {{\varepsilon }_{0}}r\left( v \right)}}} \right. \kern-0em} {4\pi {{\varepsilon }_{0}}r\left( v \right)}}$ электрического взаимодействия иона с зарядом qи величины его кинетической энергии, которая в предельном случае может быть порядка ${{{{m}_{i}}{{v}^{2}}} \mathord{\left/ {\vphantom {{{{m}_{i}}{{v}^{2}}} 2}} \right. \kern-0em} 2}$, т.е., приравнивая, получим

(5)
$Q\left( v \right) = \frac{{2\pi {{m}_{i}}{{\varepsilon }_{0}}}}{q}{{v}^{2}}r\left( v \right).$

Выражение (5) фактически означает, что величина заряда фрагмента в процессе движения зависит от его скорости и размера. Рассмотрим в качестве примера движение тела с формой шара и неизменными радиус и скорость. Например, пусть скорость равна $v = 5 \times {{10}^{4}}\,\,{{\text{м}} \mathord{\left/ {\vphantom {{\text{м}} {\text{с}}}} \right. \kern-0em} {\text{с}}}$ и радиус $r = {{10}^{{ - 3}}}\,\,{\text{м}}$. Тогда выражение (5) соответствует заряду $Q \approx {{10}^{{ - 10}}}{\text{\;}}\,\,{\text{Кл}}$. Теперь оценим минимальное расстояние L, необходимое для сообщения такого заряда фрагменту из равенства $L\pi {{r}^{2}}{{N}_{a}} = {Q \mathord{\left/ {\vphantom {Q q}} \right. \kern-0em} q}$, где ${{N}_{a}} = 2 \times {{10}^{{21}}}$ – концентрация молекул атмосферы на высоте около 70 км. Это дает величину пути $L \approx {{10}^{{ - 7}}}\,\,{\text{м}}$. Если предположить, что из количества соударения равном 107 хотя бы одно из них заканчивается отрывом электрона, то даже в этом случае на пути в 1 метр будет передан максимально возможный заряд. Этот пример показывает, что максимальная величина заряда отслеживается выражением (5) практически мгновенно.

СПЕКТР ИЗЛУЧЕНИЯ

Для оценки напряженности электрического поля волны нужно учитывать тот факт, что наблюдатель находится от МТ на расстоянии нескольких десятков километров и больше, а длина электромагнитной волны λ = с/ν для частоты ν ≈1 кГц равна 300 км. Это значит, что в оценку напряженности поля с помощью формул дипольного приближения будет внесена существенная погрешность, более того, в ближней зоне вклад кулоновского поля может быть очень существенным. Это, прежде всего, означает, что измерения поля классическими антеннами, рассчитанными для дальней зоны, вряд ли целесообразно.

Далее оценим поле, вызванное одним фрагментом с зарядом Q(5) с учетом вышесказанного. Поскольку размеры фрагмента намного меньше длины излучаемых волн и расстояния до наблюдателя, то объемную плотность тока точечного заряда (т.е. фрагмента МТ) можно выразить посредством дельта-функции (см., например, Ландау, Лифшиц, 1967)

(6)
$j\left( {r{\kern 1pt} ',t} \right) = \overrightarrow {{{e}_{{z~}}}} Q\left( t \right)v\delta \left[ {z{\kern 1pt} '\,\, - z\left( t \right)} \right]\delta \left( {x{\kern 1pt} '} \right)\delta \left( {y{\kern 1pt} '} \right).$

Здесь выбрана система отсчета, в которой заряд Q движется вдоль оси Z, а наблюдатель находится на оси X. Переменные (x', y', z') соответствуют координатам произвольного элементарного объема dV ' = dx 'dy'dz' пространства. Поскольку (6) выражено явно от времени, то для нахождения спектральной плотности тока

$\begin{gathered} {\mathbf{j}}\left( {r{\kern 1pt} ',\omega } \right) = {{{\mathbf{e}}}_{z}}\int\limits_0^\tau {{{e}^{{i\omega t}}}} Q\left( t \right)v\left( t \right) \times \\ \times \,\,\delta \left[ {z{\kern 1pt} '\,\, - z\left( t \right)} \right]\delta \left( {x{\kern 1pt} '} \right)\delta \left( {y{\kern 1pt} '} \right)dt. \\ \end{gathered} $

Следует перейти к переменной v = v(t) с помощью выражений (1)–(4)

$\begin{gathered} {\mathbf{j}}\left( {r{\kern 1pt} ',\omega } \right) = {{{\mathbf{e}}}_{z}}\int\limits_0^\tau {{{e}^{{i\omega t\left( v \right)}}}Q\left( v \right)} v \times \\ \times \,\,\delta \left[ {z{\kern 1pt} '\,\, - z\left( v \right)} \right]\delta \left( {x{\kern 1pt} '} \right)\delta \left( {y{\kern 1pt} '} \right)dt\left( v \right). \\ \end{gathered} $

Тогда индукцию электромагнитного поля выразится, как ${\mathbf{B}} = {\text{rot}}{\mathbf{A}}\left( {{\omega },{\text{r}}} \right)$, причем, А – векторный потенциал электромагнитного поля, определяемый как:

${\mathbf{A}}\left( {\omega ,r} \right) = \frac{{{{\mu }_{0}}{{e}^{{ikr}}}}}{{4\pi r}}\int {j\left( {r{\kern 1pt} ',\omega } \right)} {{e}^{{ - ikr{\kern 1pt} '}}}dV{\kern 1pt} ',$
где ${\mathbf{k}} \cdot {\mathbf{r}}{\kern 1pt} '$ = k(x'sinθ + z 'cosθ) скалярное произведение волнового вектора $\left| {\mathbf{k}} \right| = {\omega \mathord{\left/ {\vphantom {\omega c}} \right. \kern-0em} c}$ на радиус-вектор r' произвольного элемента пространства dV ', r − расстояние до наблюдателя ${{{\mu }}_{0}} \approx 4{\pi }{{10}^{{ - 7}}}{\text{\;}}{{{\text{Гн}}} \mathord{\left/ {\vphantom {{{\text{Гн}}} {\text{м}}}} \right. \kern-0em} {\text{м}}}$ − магнитная постоянная, θ − угол наблюдения.

Таким образом, имеем индукцию электромагнитного поля для количества осколков равное M/m, где М – масса метеороида.

(7)
$\begin{gathered} {\mathbf{B}} = \frac{{{{\mu }_{0}}}}{{4\pi }}{\text{rot}}\left( {{{{\mathbf{e}}}_{z}}\frac{{{{e}^{{ikr}}}}}{r}} \right)\frac{M}{m}\int {\int\limits_0^\tau {{{e}^{{i\omega t\left( v \right)}}}} } Q\left( v \right) \times \\ \times \,\,\left. {v\delta \left[ {z{\kern 1pt} '\,\, - z\left( v \right)} \right]\delta \left( {x{\kern 1pt} '} \right)\delta \left( {y{\kern 1pt} '} \right){{e}^{{ - ikr{\kern 1pt} '}}}\frac{{dv}}{{a\left( v \right)}}dV{\kern 1pt} '} \right]. \\ \end{gathered} $

Сначала следует проинтегрировать по пространству V':

${\mathbf{B}} = \frac{{{{\mu }_{0}}}}{{4\pi }}{\text{rot}}\frac{{{{e}^{{ikr}}}}}{r}\frac{M}{m}\left[ {\int\limits_{{{v}_{0}}}^{{{v}_{1}}} {{{e}^{{i\omega t\left( v \right)}}}} Q\left( v \right)v{{e}^{{ - ikZ\left( v \right)}}}\frac{{dv}}{{a\left( v \right)}}} \right],$
где ${\text{rot}}\frac{{{{e}^{{ikr}}}}}{r} = {{e}_{y}}\left( {ik - \frac{1}{r}} \right)\frac{2}{r}{{e}^{{ikr}}}\cos \theta $. Энергия, излучаемая за секунду на площадь в 1 м2 в единичном интервале частот равна $\Pi \left( {\omega ,r} \right) = \frac{c}{{{{\mu }_{0}}}}{{\left| {{\mathbf{B}}\left( {\omega ,r} \right)} \right|}^{2}}$. На рис. 2 графически показана зависимость удельной интенсивности (модуль вектора Умова–Пойнтинга) от частоты и размера фрагмента.

Рис. 2.

Спектральная интенсивность излучения, сопровождающего разрыв метеороида на мелкие фрагменты. Для размера фрагмента r0 = 10−4 м шкала интенсивности расположена в правой части рисунка.

Максимальное значение интенсивности растет с уменьшением размера фрагмента. Одновременно с этим расширяется спектральный диапазон. Даже для размера фрагмента $r = {{10}^{{ - 4}}}\,\,{\text{м}}$ значений спектральной интенсивности ${\text{П}}\sim 5 \times {{10}^{7}}{\text{\;Ян}}$, соответствующее спектральной напряженности $E\left( \omega \right) \approx 1.5 \times {{10}^{{ - 9}}}{{~В} \mathord{\left/ {\vphantom {{~В} {\left( {{\text{м}}\,{\text{Гц}}} \right)}}} \right. \kern-0em} {\left( {{\text{м}}\,{\text{Гц}}} \right)}}$, что, вполне достаточно для уверенной регистрации радиоимпульса. Например, для простейшей антенны с действующей высотой h ~ 10 метров ЭДС на клеммах составит ${{U}_{i}} = E\left( \omega \right)\omega h = 15~\mu V$. Этого вполне доступно даже для простейших радиоприемников.

Характерной особенностью этого спектра является максимум при частоте ω = 0. Это связано с тем, что в выражении (7) для ротора экспоненты вмножителе (ik – 1/r) оставлено слагаемое 1/r, так как фактически часто наблюдения проводятся в ближней зоне, где кулоновское поле дает значительный вклад. В экспериментальных наблюдениях в диапазоне VLF эта составляющая спектра, чаще всего, обрезается преселектором радиоприемника. Поэтому эта часть спектра не регистрируется.

ЗАКЛЮЧЕНИЕ

Генерация электромагнитных полей в диапазоне ELF/VLF, вызванная конечной стадией движения метеорных тел, являются экспериментально доказанным фактом. Однако в отношении механизма этого явления в настоящее время существует несколько гипотез, не подкрепленных численными расчетами, которые не противоречили бы экспериментальным фактам. В предлагаемой модели без введения экзотических гипотез на основе хорошо известных явленийпоказано, каким образом могут реализоваться электромагнитные поля на заключительном этапе движения метеороида. Весьма вероятно, что часто наблюдаемое “жужжание”, исходящее из радиоприемника VLF (Beech, Brown, 1995), связано с разлетом мелких фрагментов, которые излучают в килогерцовом диапазоне. Численно полученные оценки не претендуют на высокую точность вследствие отсутствия надежных экспериментальных данных. Тем не менее, вызванные этим фактом ошибки вряд ли превышают один порядок. Слабым местом в этой модели является допущение о возможном дроблении метеороида на мелкие осколки, однако теоретические работы (Егорова, Лохин, 2015; 2017) подкрепляют такую возможность.

На высоте, где происходит разрушение в виде “теплового взрыва”, возникает ударная волна и скорость молекул кислорода достигших метеороида не столь велика, чтобы вызвать существенную электризацию. Однако после взрыва размеры большинства разлетающихся частиц сравнимы с длиной свободного пробега молекул (λ = 3 × 10–4–10−3 м на высоте 60–70 км (см. Атмосфера стандартная, 2004)) и обтекание будет иметь характер свободного молекулярного потока. Поэтому использование выше записанных формул не противоречит принятой модели электризации.

Из содержания модели VLF-излучения следует, что его спектральный состав должен зависит от размера фрагментов и плотности метеороида, а это, в свою очередь, зависит от происхождения метеорного потока. Другими словами, характер взрывной стадия метеора и распределение осколков по размерам зависят от агрегатного состояния вещества метеора и его химического состава. Поэтому следует ожидать, что характерные особенности спектров известных метеорных потоков будут различаться. Этот вывод подтверждается в работе (Guha, Kumar, 2009).

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

  1. Аскарьян Г.А. Избыточный отрицательный заряд электрон-фотонного ливня и когерентное излучение от него // ЖЭТФ. 1961. Т. 41. № 8. С. 616–618.

  2. Астапович И.С. Метеорные явления в атмосфере Земли. М.: Физ-мат лит. 1958. 640 с.

  3. Атмосфера стандартная. Параметры. Москва, ИПК изд-во стандартов, 2004, с. 180.

  4. Бронштэн В.А. Физика метеорных явлений. М.: Наука, 1981. 416 с.

  5. Егорова Л.А., Лохин В.В. О взрывоподобном разрушении метеороида // Экологический вестник научных центров Черноморского экономического сотрудничества (ЧЭС). 2017. № 4. С. 1–5.

  6. Егорова Л.А., Лохин В.В. О двухстадийном разрушении метеороида с концевой вспышкой // Вестн. Московского университета. 2015. Сер. 1: Математика. Механика, Изд-во Моск. ун-та (М.). № 4. С. 43–47.

  7. Зельдович Я.Б., Райзер Ю.П. Физика ударных волн и высокотемпературных гидродинамических явлений. М.: Наука, 1996. 688 с.

  8. Ландау Л.Д., Лифшиц Е.М. Теория поля. М.: Наука, 1967. 460 с.

  9. Лёб Л. Статическая электризация. Госэнергоиздат. 1963. 408 с.

  10. Мусатенко C.И., Мусатенко Ю.С., Курочка Е.В., Ласточкин А.В., Чолий B.Я., Максименко О.И., Слипченко А.С. Пылевая плазма в среднеширотной ионосфере в периоды метеорных потоков // Геомагнетизм и аэрономия. 2006. Т. 46. № 2. С. 182–192.

  11. Филоненко А.Д. Радиоизлучение метеорных тел в диапазоне декаметровых волн // Геомагнетизм и аэрономия. 2018. Т. 58. № 5. С. 720–727.

  12. Филоненко А.Д. Радиоизлучение широких атмосферных ливней // УФН. 2015. Т. 185. № 7. С. 673−716.

  13. Beech M., Brown P., Jones J. VLF Detection of fireballs // Earth, Moon and Planets. 1995. V. 68. P. 181–188.

  14. Beech M., Foschini L. A space charge model for electrophonic bursters // Astron. and Astrophys. 1999. V. 345. P. L27–L31.

  15. Cepleha Z., Borovicka J., Eleford W. et al. Meteor phenomena and bodies // Space Sci. Rev. 1998. V. 84. P. 327–471.

  16. Garaj S., Vinkovi D., Zgrablic G. et al. Observational Detection of meteor-produced VLF electromagnetic radiation // Fizika A (Zagreb). 1999. V. 8. № 3. P. 91–98.

  17. Guha A., Kumar B., Rakesh R. Possible Detection of GEMINID 2007 Meteor Shower During Day-Time from VLF Radiation Spectra // Earth, Moon, and Planets. 2009. V. 105. № 1. P. 31–40.

  18. Keay C.S.L. Progress in Explaining the Mysterious Sounds Produced by Very Large Meteor Fireballs // Scientific Exploration. 1993. V. 7. № 4. P. 337–354.

  19. Keay C.S.L. Continued Progress in Electrophonic Fireball Investigations // Earth, Moon and Planets. 1995. V. 68. P. 361–368.

  20. Kopnin S.I., Popel S.I., M. Y. Yu. Phenomena associated with complex (dusty) plasmas in the ionosphere during high-speed meteor showers // Physics of Plasmas. 2009. V. 16. P. 063705.

  21. Obenberger K.S., Taylor G.B., Hartman J.M. et al. Detection of Radio Emission from Fireballs // arXiv: 1405.6772

  22. Obenberger K.S., Taylor G.B., Hartman J.M. et al. Monitoring the Sky with the Prototype All-Sky on the LWA1 // arXiv: 1503.05150

  23. Obenberger K.S., Taylor G.B., Lin C.S. et al. Dynamic Radio Spectra from two Fireballs // arXiv: 1510.04147

  24. Obenberger K.S., Dowell J.D., Hancock P.J. et al. Rates, Flux Densities, and Spectral Indices of Meteor Radio Afterglows // arXiv: 1607.03002

  25. Obenberger K.S., Holmes J.M., Dowell J.D. et al. Altitudinal dependence of meteor radio afterglows measured via optical counterparts // arXiv: 1609.04701

  26. Price C., Blum M. ELF/VLF radiation produced by the 1999 Leonid meteors // Earth, Moon and Planets. 2000. V. 82. P. 545–554.

  27. Vinkovic I.D., Garaj S., Lim P.L. et al. Global Electrophonic Fireball Survey: a review of witness reports // arXiv:astro-ph/0211203

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