WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Афраймович Э.Л., Астафьева Э.И., Кирюшкин В.В. Ионосферное возмущение в ближней зоне эпицентра землетрясения на о. Хоккайдо 25 сентября 2003 г.

Научная статья

 

Электронный журнал «Исследовано в России»           1955               http://zhumal.ape.relarn.ru/articles/2004/183.pdf

Ионосферное возмущение в ближней зоне эпицентра землетрясения на о. Хоккайдо 25 сентября 2003 г.

Э.Л. A(bpaiiMOBH4(afra@iszf.irk.ru), Э.И. Астафьева, В.В. Кирюшкин Институт солнечно-земной физики, СО РАН

1.   ВВЕДЕНИЕ

Исследованию ионосферного отклика на возмущения, возникающие при импульсном воздействии на земную атмосферу (землетрясения, взрывы вулканов, промышленные взрывы, запуски ракет, подземные испытания), посвящено множество работ [Row, 1967; Павлов, 1979, 1986; Li et al., 1994; Rudenko, Uralov, 1995; Calais at al., 1995, 1998, 2003; Afraimovich et al., 2001; Афраймович и др., 2001, 2002; Ахмедов, Куницын, 2003, 2004, и др]. При этом во многих случаях значительная доля энергии исходного атмосферного возмущения оказывается сосредоточенной в акустической ударной волне [Row, 1967; Павлов, 1979, 1986]. В литературе для обозначения ионосферного отклика ударной волны используют термины, отличающиеся различной физической интерпретацией, в том числе «ударно-акустическая волна (УАВ)» или более общее обозначение «ионосферное возмущение (ИВ)» [Павлов, 1979, 1986; Li et al., 1994].

В настоящее время основные параметры, характеризующие ионосферный отклик УАВ во временной области (форма, амплитуда, период), исследованы достаточно широко. Следующим шагом является детальное исследование пространственно-временных характеристик УАВ: измерение скорости и динамики распространение фронта УАВ, и на этой основе определение местоположения и время воздействия импульсного источника ИВ. Это необходимо также для селекции УАВ, возникающих при техногенных воздействиях, от ИВ естественного происхождения, таких как эффекты магнитных бурь, погодные фронты, тропосферная турбулентность и струйные течения, солнечный терминатор.

Методы анализа вариаций полного электронного содержания (ПЭС), измеренного при трансионосферном зондировании двухчастотными навигационными сигналами GPS, оказались весьма продуктивным для решения поставленной задачи.

В работах [Афраймович и др., 2001; Afraimovich et al. 2001] был разработан метод определения параметров УАВ (включая фазовую скорость, угловые характеристики волнового вектора УАВ, направление на источник и его положение) с помощью элементарной 3-компонентной GPS-решетки, расположенной в дальней зоне источника (приближение плоского фронта возмущения). В отличие от известных радиофизических методов этот метод дает оценку параметров УАВ без априорных сведений о месте и времени момента импульсного воздействия.

Однако методы [Афраймович и др., 2001; Afraimovich et al. 2001] применимы только тогда, когда фронт волны возмущения в пределах анализируемой области распределения ПЭС мало отличается от плоского, что оказывается неприемлемым для совместного анализа данных далеко разнесенных станций GPS (в том числе и в обсуждаемом эксперименте). В связи с этим, актуальна задача определения характеристик УАВ, зафиксированной в ближней зоне точечного источника, т.е. когда сферичностью фронта волны пренебречь нельзя [Афраймович и др., 2002].

Дальнейшее развитие метода GPS-детектирования УАВ основано на использовании концепции фазированной антенной решетки (ФАР). В основу положен предложенный в


Электронный журнал «Исследовано в России»           1956              http://zhurnal.ape.relarn.ru/articles/2004/183.pdf

работе [Afraimovich, 2000] принцип когерентной обработки вариаций ПЭС в ионосфере, измеренных одновременно для всей совокупности лучей на ИСЗ, на всех выбранных для анализа станциях сети GPS, составляющих неэквидистантную ФАР. Когерентная обработка полученных вариаций ПЭС использует известные алгоритмы ФАР. В зависимости от вида анализируемого возмущения выбирается вид комплексной характеристики ФАР, соответствующей пространственно-временным свойствам данного возмущения.

Целью настоящей работы является определение характеристик ионосферного отклика на главный толчок землетрясения 25 сентября 2003 года в 19:50:07 UT в 60 км юго-восточнее острова Хоккайдо с магнитудой М=8.3.


2.   ОБЩАЯ ХАРАКТЕРИСТИКА ЭКСПЕРИМЕНТА

Эпицентр землетрясения (Ф=41.8°, Л=143.85°) находился на глубине 27-33 км под водой между Курильской и Японской впадинами (рис.1). В районе эпицентра будущего

землетрясения          было

Широта (N)

зарегистрировано         10

форшоков                        с

магниту дами 3.9-5.7. Землетрясение сопровождалось цунами по всему побережью Японских островов. В последующие 5 дней было зафиксировано больше 100 афтершоков магнитудой от 3.8 до 7.4. Приведенные сведения получены по данным сайта http://www.neic.cr .usgs . gov.

32

142                   144

Долгота (Е)

Рис. 1.

136

138

140

146

150

Данные GPS с 30-

c разрешением в

стандартном формате

RINEX                       [26],

используемые в данном исследовании, представлены   на   сайте http://sopac.ucsd.

edu/cgi-bin/dbDataByPate.cgi.

Геометрия эксперимента представлена на рис.  1. Звездочкой отмечено положение

эпицентра землетрясения. Жирными точками и крупным шрифтом отмечены положение и

названия станций GPS.

3. ПЕРВИЧНАЯ ОБРАБОТКА ДАННЫХ GPS. МЕТОД ПРОСТРАНСТВЕННО-ВРЕМЕННОЙ ОБРАБОТКИ

Исходными данными для формирования пространственно-временного распределения ионосферных откликов УАВ являются временные ряды высокочастотных вариаций ПЭС dl(t) и соответствующие им ряды азимута ocs(t) и угла места 0s(t) луча зрения на ИСЗ.


Электронный журнал «Исследовано в России»                          I957                    http://zhurnal.ape.relarn.ru/articles/2004/183.pdf


Методика  обработки   данных   подробно   описана  в   работах   [Afraimovich   et   al.,   2001;

Афраймович и др., 2001, 2002]. Соответствующие временные зависимости исходных рядов

«наклонного» Ii(t) для USUD; PRN13 (слева) и MIZU; PRN24 (справа) представлены на рис.

2а и 2ж. Отфильтрованные в диапазоне 1-10 мин зависимости dli(t) для всех 10 лучей даны на

панелях б-е; ж-м рис. 2.

Ж.

Из рис. 2 видно, что на

фоне  медленных  изменений

ПЭС   уверенно   выделяются

обусловленные

распространением                 УАВ

-0.1  —

-MIZU; PRN 24

-0.2 -

0.1

MIZU; PRN 13

----- '------ 1----- >"

И.

0

-0.1 — -0.2 — -0.3

о

о

быстрые колебания N-образной формы, параметры которых лежат в следующих пределах: амплитуда А меняется от 0.1 до 0.25 TECU (TECU - Total Electron Content Units - единица ПЭС, равная 10 м"), а период возмущения AT - около 500 с.

Для              устранения

TSKB; PRN 24

0.1    —|

о  

-0.1     —

-0.2

-0.3

Г.

К.

^л^^лА/Ч

о

неопределенности

локализации ионосферного

отклика УАВ, вызванной

интегральным характером

значения ПЭС, будем

считать, что значение ПЭС

формируется         в        точке

USUD; PRN 24

0.1   —i

0 -0.1 -0.2

пересечения луча зрения на

ИСЗ с плоскостью на высоте

hd          выше          максимума

Д.

Л.

о

О

ионизации области F2

ионосферы, которая вносит

основной            вклад            в

KGNI; PRN 24

KGNI; PRN 13

-0.3 0.1

формирование ПЭС. Ниже во

всех расчетах мы принимали

hd=400          км,          наиболее

о

-0.1 — -0.2 — -0.3

М.

о

О

соответствующей ночным

условиям            эксперимента

(около            06:00            LT).

KMSV; PRN 13

---- ^Ап

KMSV-PRN 24

---- г±-

п------ <~

19.5           20

Время, UT

1----- <~

19.5           20

Время, UT

20.5

20.5

Ееографические координаты

подыоносферной точки (Ф,Л)

определяются                      как

Рис. 2.

координаты проекции на поверхность Земли точки пересечения луча зрения на

ИСЗ с плоскостью на высоте hd. Траектории подыоносферных точек приведены на рис.1

тонкими линиями.

Вычислив координаты подыоносферных точек для моментов tmin (этот момент для

MIZU; PRN13 обозначен на рис. 26), определим пространственное положение ионосферного

отклика УАВ. При этом положение каждого луча на ИСЗ определяется значениями 6s(t) и

ocs(t) для моментов времени tmin- Определив моменты времени tmhy для каждого ряда dli(t) и


Электронный журнал «Исследовано в России»           1958               http://zhumal.ape.relarn.ru/articles/2004/183.pdf

соответствующие координаты (Ф;, Л;) подыоносферных точек, получаем пространственно-временное распределение ионосферных откликов УАВ (рис.1, рис. 2).

Значения tmin;i, Ф;, Л; для каждого луча «приемник GPS-ИСЗ», по которому зафиксированы ионосферные отклики, приведены в Таблице 2. Для удобства дальнейших расчетов широта и долгота подыоносферных точек пересчитываются в горизонтальные декартовы координаты (х;, у;) топоцентрической системы координат, в которой начало совпадает с одной из подыоносферных точек - опорной точкой О (хо, уо). В этом приближении воображаемые детекторы расположены на горизонтальной плоскости; ошибка такого приближения растет при увеличении расстояния между станциями GPS.

За опорную точку принимаем положение ионосферного отклика УАВ с минимальным значением времени tm;n. Для рассматриваемого события это подыоносферная точка MIZU; PRN13 с временем tmin;o=20.058 UT. Тогда временной ряд ПЭС с удаленным трендом dIo(t) для GPS-станции MIZU; PRN13 также будем считать опорным.

Используем упрощенную модель, в которой эпицентральный излучатель УАВ при землетрясении подменяется точечным источником УАВ, расположенным на высоте hs в ионосфере. При hs=0 км он совпадает с эпицентром землетрясения. Фронт УАВ представляет собой полусферу, распространяющуюся над поверхностью Земли с постоянной скоростью. Возможные рефракционные эффекты, возникающие при распространении акустического возмущения от источника до подыоносферной точки, в первом приближении игнорируем. Кроме того, не учитываем преобразование амплитуды отклика ПЭС на прохождение УАВ, которое определяется пространственным затуханием амплитуды УАВ, ракурсными соотношениями между волновым вектором возмущения, направлением на ИСЗ и ориентацией вектора магнитного поля.

При этих условиях алгоритм пространственно-временной обработки сводится к суммированию предварительно сфазированных рядов dli(t) к опорному ряду dI0(t), в результате чего получаем полный сигнал пространственной сборки рядов ПЭС с удаленным трендом:

^(0 = ^о(0+И(^-А^)                                 (1)

г=1

где Ati=(tmax,i-tmax,о) -время сдвига i-ro ряда dl/t) относительно опорного ряда dlo(t); Р - количество суммируемых рядов.

Вследствие того, что полезные колебания ПЭС, обусловленные распространением УАВ, в каждом суммируемом ряде хорошо коррелированы между собой, а фоновые шумовые колебания некоррелированы (см. рис.2 б-е; з-м), отношение сигнал/шум в полном сигнале пространственной сборки возрастает не менее чем в л/Р раз.

Время сдвига суммируемых рядов будем определять следующим образом:

ААм ,                                                             (2)

А г.

Vk,i

где Aph,i=(ph,i-Ph,o) - разность горизонтальных удалений i-ой и опорной подыоносферных точек от источника ИВ; Vh,i=Vr/cos9 - горизонтальная составляющая скорости распространения ИВ, Vr - радиальная скорость.

Скорость Vh,i зависит от величины угла места 0 нормали к волновому фронту возмущения в iподыоносферной точке. С учетом принятой модели сферического фронта УАВ величина этого угла определяется отношением горизонтального и радиального удаления i-й подыоносферной точки от источника ИВ:


COS0 =

В свою очередь, выражение для Арьд имеет вид:


Ph,i


(3)


Электронный журнал «Исследовано в России»           1959               http://zhumal.ape.relarn.ru/articles/2004/183.pdf

Ph,i = ^(Хг-\)2+{Уг~Уи)2 '                                                        (4)

где (х;,у;) - координаты i-ой подыоносферной точки; (хи,уи) - координаты источника ИВ.

В этом случае с учетом того, что начало выбранной топоцентрической системы координат совпадает с опорной подыоносферной точкой (т.е. х0=0 и уо=0), выражение (2) принимает вид:

Ат= (V (*,- - хи У + {у, -у и)2 - У*»2 + у и У(х, -хиУ+ iy, -yJ(5)

Тогда для всей совокупности подыоносферных точек, в которых зафиксированы ионосферные отклики УАВ, имеем систему из (Р-1) уравнений:

Ат= (V fa ~хи)2+ (ух -уи)2- лк2 + у и М*1 - х« Т + (л - yJ

K^-xuf+(y-yuf+hs2

д = Wfe-i - хи У+(yF_, -yuf-^xu2+уи2 A/fe-i - хи У+(у„_. - л У        (6)

^ л/^-i - х» )2 + (у,-. - ^«)" + К2

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

С помощью формул (2-6) были вычислены значения координат Ф8 и Л8, момента времени ts и запаздывания включения источника ИВ Ъ\ относительно момента времени главного толчка to, а также радиальной скорости Vr перемещения ИВ и пространственной невязки о для различных значений высоты hs предполагаемого источника ИВ (Таблица 1).

Координаты источника, время и запаздывание включения источника ИВ, скорость перемещения ИВ для различных значений высоты hs:

Таблица   .

hs, км

5h, км

Ф,,Е°

Л,,№

to,UT

5t, с

Vr, м/с

Vv, м/с

о, км

0

400

42,6

143,8

20

594

700

673

6,8

100

300

42,8

143,8

19,989

554

700

541

6,9

200

200

43

143,6

19,95

414

650

483

6,5

300

100

43

143

19,907

259

600

386

6,2

400

0

43,2

143,6

19,842

25

550

-

6

-

-

41,8

143,85

19,835

4.   ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ

Об обнаружении методом GPS-просвечивания ионосферных возмущений, предположительно генерируемых вертикальным смещением земной поверхности непосредственно в эпицентральной зоне землетрясения, сообщается в работах [Row, 1967; Afraimovich   et   al.,   2001;   Афраймович   и   др.,   2001].   Результаты   настоящей   работы


Электронный журнал «Исследовано в России»           1960              http://zhurnal.ape.relarn.ru/articles/2004/183.pdf


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

Общепринятой считается модель ионосферного возмущения при землетрясении, в которой источником возмущения является не сам эпицентр, а область в ионосфере над эпицентром [Row, 1967; Павлов, 1979, 1986; Ахмедов, Куницын, 2003, 2004]. Согласно этой модели, на ионосферные высоты волновое акустические возмущение проникает только в узком секторе зенитных углов меньше 20°. Остальные волны захватываются в акустический волновод волны и не достигают ионосферных высот. Наш эксперимент вполне согласуется с такой моделью (см. ниже).

Геометрия эксперимента в вертикальной плоскости представлена на рис.3 а. Здесь 5h и р - высотное и радиальное удаление детектора ИВ от вторичного источника ИВ, локализованного на высоте hs (обозначено жирными точками). Предполагается, что детектор ИВ локализован в окрестности главного максимума электронной концентрации на высоте 400 км (обозначен на оси ординат закрашенным ромбом). Стрелками дано обозначение волновых векторов ИВ, распространяющегося от вторичного источника до детектора с радиальной скоростью Vr и акустического возмущения, распространяющегося от эпицентра

а. землетрясения до вторичного источника в ионосфере с вертикальной скоростью Vv.

На рис. 36 данные

Таблицы 1 приведены в виде

графика зависимости от высоты hs

времени запаздывания включения

источника          5t          (пунктир),

радиальной скорости Vr (тонкая кривая) и вертикальной скорости Vv (жирная кривая). Как видно из рис.3, наименьшее значение невязки получено для высоты hs = б-    400 км.

Полученные нами оценки скорости Vr близки к скорости звука С на высоте порядка 400 км [Li et al., 1994]. Таким образом, принятое выше значение высоты hd=400 км, на которой расположен воображаемый детектор ИВ, вполне согласуется с высотной зависимостью скорости звука C(h) [Lietal, 1994].

Рис. 3                                                              Значение                скорости

распространения V=700 м/с, определенное методом [Афраймович и др., 2002], немного меньше оценок, представленных в [Afraimovich et al., 2001; Афраймович и др., 2001]. В этих работах фазовая скорость УАВ, измеренная с помощью GPS-решеток в дальней зоне эпицентров землетрясений в Турции 17 августа и 12 ноября 1999 г., составила соответственно V=l 173 м/с и V=l 157 м/с.

Интересно сравнить результаты пространственно-временной обработки [Афраймович и др., 2002] с оценками средней скорости распространения УАВ Va, определяемой по запаздыванию At=tmin-to и известной длине пути L между эпицентром землетрясения и подыоносферной точкой. Длина пути L определяется по дуге большого круга. Соответствующие данные представлены в Таблице 2.


Электронный журнал «Исследовано в России»           1961               http://zhumal.ape.relarn.ru/articles/2004/183.pdf

Основные параметры ионосферных откликов при землетрясении на о. Хоккайдо:

Таблица 2.

Станция -PRN

ЧШП5  U А

Фт1„,°

Л •   °

AT, с

8t,4

Va, М/С

1

MIZU-13

20,067

42,654

145,847

600

0,232

428

2

USUD-13

20,1

40,498

143,528

510

0,265

352

3

KSMV-13

20,117

40,763

146,244

360

0,282

376

4

TSKB-13

20,108

40,337

145,112

570

0,273

364

5

KGNI-13

20,092

40,613

145,088

510

0,257

374

6

MIZU-24

20,158

39,723

140,804

417

0,323

398

7

USUD-24

20,3

36,759

138,349

871

0,465

483

8

KGNI-24

20,325

35,832

139,795

838

0,49

467

9

TSKB-24

20,308

36,724

139,986

838

0,473

431

10

KSMV-24

20,3

36,099

140,884

748

0,465

455

Как можно заметить, полученные таким образом значения скорости в среднем в 1.5 раза меньше, чем определенные выше значения Vr (Таблица 1). Это различие можно объяснить тем, что длина пути L определяется вдоль земной поверхности, а детектирование производится в области главного максимума электронной концентрации на высоте 400 км. Это удлиняет путь в нашем случае в 1.5 раза. Это согласуется с тем, что источником возмущения является не сам подземный толчок, а вторичный источник, расположенный на высоте hs в ионосфере, и включается он позже, чем происходит толчок.

Как можно убедиться, результаты определения координат источника ИВ (Ф8=43.0°, As=143.0°) методом пространственно-временной обработки хорошо согласуются с координатами эпицентра землетрясения (Фе=41.8°, Ле=143.85°). Точность определения положения источника ИВ этим методом составляет соответственно около 60 км. Отличие полученным значений от истинных (определенных по сейсмическим данным) можно объяснить тем, что в используемом здесь методе [Афраймович и др., 2002] использовалось приближение плоской земли, которое на расстояниях порядка 500 км дает соответствующее отклонение.

Присутствие зонального и меридионального ветров на ионосферных высотах приводит к смещению и деформации волнового фронта и, соответственно, появлению зависимости интенсивности акустической волны от направления распространения. При этом решающую роль играет градиент скорости ветра [Li et al., 1994; Ахмедов, Куницын, 2003, 2004]. Это может объяснить смещение горизонтального положения вторичного источника возмущения в ионосфере относительно эпицентра землетрясения.

5.   ЗАКЛЮЧЕНИЕ

С помощью разработанного авторами метода пространственно-временной обработки распределения вариаций полного электронного содержания, полученных по данным станций GPS, расположенных в ближней зоне землетрясения вблизи о. Хоккайдо 25 сентября 2003 г. (магнитуда 8.3 балла), определены координаты эпицентра, время главного толчка и скорость распространения ионосферного возмущения, вызванного землетрясением. Указанные параметры хорошо согласуются с истинными координатами эпицентра и временем толчка (сейсмические данные), а также с результатами, полученными ранее при анализе ионосферных возмущений во время сильных землетрясений 17.08. и 12.11 1999 г. в Турции, 4.06.2000 г. на Южной Суматре и 13.01.2001 г. вблизи побережья Центральной Америки.

Полученные в работе результаты согласуются с моделью ионосферного возмущения при землетрясении, в которой источником возмущения является не сам эпицентр, а область в ионосфере над эпицентром [Row, 1967; Павлов, 1979, 1986; Ахмедов, Куницын, 2003, 200].


Электронный журнал «Исследовано в России»           1962               http://zhumal.ape.relarn.ru/articles/2004/183.pdf

С целью более детального исследования пространственного затухания волны, дисперсионных свойств атмосферы, определения фазовой и групповой скорости распространения волнового пакета ИВ, представляет большой интерес использование предложенного метода для большего количества станций GPS [Calais et al., 2003].

Авторы благодарны СВ. Воейкову за помощь при первичной обработке данных. Эта работа выполнена при поддержке Российского фонда фундаментальных исследований (гранты 00-05-72026 и 03-05-64100), а также гранта N НШ-272.2003.5 государственной поддержки ведущих научных школ Российской Федерации. Мы благодарим сотрудников Scripps Orbit and Permanent Array Center (SOPAC) за предоставление используемых в настоящей статье первичных данных глобальной сети наземных двухчастотных приемников GPS.


Электронный журнал «Исследовано в России»           1963               http://zhurnal.ape.relarn.ru/articles/2004/183.pdf

СПИСОК ЛИТЕРАТУРЫ

1.    Афраймович Э.Л., Е.А. Косогоров, А.В. Плотников, A.M. Уралов. // Физика Земли.

2001. Т.67. С.16.

  1. Afraimovich,  E.L.,  N.P.  Perevalova,  A.V.  Plotnikov  and  A.M.  Uralov.  // Annales Geophysicae. 2001. V.19. P.395.
  2. Афраймович Э.Л., B.B. Кирюшкин, Н.П. Перевалова. // Радиотехника и электроника.

2002. Т.47. N7. С.822.

    • Afraimovich E.L.//Radio Sci. 2000. V.35(6).P.1417.
    • Ахмедов P.P.,  B.E.  Куницын.  // Вестник Московского  Университета.   Серия  3. Физика. Астрономия. 2003. N3. С.38.
    • Ахмедов P.P., B.E. Куницын. //Геомагнетизм и аэрономия. 2004. Т. 44. N1. С.1.
    • Павлов В.А. //Изв. ВУЗов. Радиофизика. 1979. Т.22. С.18.
    • Павлов В.А. //Геомагнетизм и аэрономия. 1986. Т.26. N.5. С.807.
    • Calais, Е., and J.B. Minster. // Geophys. Res. Lett. 1995. V.22. P. 1045.
    • Calais, E., J.B. Minster, M.A. Hofton, and M.A.H. Hedlin. // Geophys. J. Int. 1998. V.132. P.191.
    • Calais, E., J.S. Haase, J.B. Minster. // Geophys. Res. Letters. 2003. V.30. N.12. 1628, 10.1029/2003GL017708.
    • Li, Y.Q., A.R. Jacobson, R.C.  Carlos, R.S. Massey, Y.N.  Taranenko, and G. Wu. Geophys. Res. Lett. 1994. V.21. P.2737.
    • Row, R.V. //J. Geophys. Res. 1967. V.72. P.1599.
    • Rudenko, G.V., A.M. Uralov. // J. Atmos. Terr. Phys. 1995. V.57. P.225.
     



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.