Бесплатный автореферат и диссертация по наукам о земле на тему

Численное моделирование генерации акустико-гравитационных волн и ионосферных возмущений от наземных и атмосферных источников

ВАК РФ 25.00.29, Физика атмосферы и гидросферы

Автореферат диссертации по теме "Численное моделирование генерации акустико-гравитационных волн и ионосферных возмущений от наземных и атмосферных источников"

МОСКОВСКИЙ ОРДЕНА ЛЕНИНА, ОРДЕНА ОКТЯБРЬСКОЙ РЕВОЛЮЦИИ, ОРДЕНА ТРУДОВОГО КРАСНОГО ЗНАМЕНИ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ им. М.В. ЛОМОНОСОВА

ФИЗИЧЕСКИЙ ФАКУЛЬТЕТ

На правах рукописи УДК 550.388.2+551.509.33

АХМЕДОВ Раван Рамин оглы

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ГЕНЕРАЦИИ АКУСТИКО-ГРАВИТАЦИОННЫХ ВОЛН И ИОНОСФЕРНЫХ ВОЗМУЩЕНИЙ ОТ НАЗЕМНЫХ И АТМОСФЕРНЫХ ИСТОЧНИКОВ

Специальность 25.00.29 - физика атмосферы и гидросферы

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

Москва 2004 г.

Работа выполнена на кафедре физики атмосферы физического факультета Московского Государственного Университета им. М.В. Ломоносова

Научный руководитель:

доктор физико-математических наук, КУНИЦЫН Вячеслав Евгеньевич профессор

Официальные оппоненты:

доктор физико-математических наук, КУТУЗА Борис Георгиевич профессор

доктор физико-математических наук, ШУВАЛОВ Валерий Викторович

Ведущая организация: Научно-исследовательский институт физики Санкт-Петербургского Государственного Университета

Защита диссертации состоится «22» апреля 2004 г.

в_часов на заседании Диссертационного совета Д501.001.63

при МГУ им. М.В. Ломоносова по адресу: 119992 Москва, Ленинские горы, физический факультет, аудитория_

С диссертацией можно ознакомиться в библиотеке физического факультета МГУ.

Автореферат разослан «19» марта 2004 г.

Ученый секретарь

Диссертационного совета Д501.001.63 кандидат физ.-мат. наук

В.Б. СМИРНОВ

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность. В конце XIX века ученые начали обращать внимание на то, что сила тяжести и стратификация геофизических сред сильно модифицируют распространяющиеся в них звуковые волны. С тех пор начались первые теоретические исследования свойств распространения акустико-гравитационных волн (АГВ) в атмосфере. Наряду с этим, с развитием радиотехнических средств, в середине XX столетия начались наблюдения перемещающихся ионосферных возмущений (ПИВ). Впервые в 1960м году ПИВ были интерпретированы как проявления АГВ в ионосфере [1]. С тех пор с помощью различных методов ведутся интенсивные наблюдения АГВ в нижней и верхней атмосфере. В 90-х годах, с развитием дистанционных методов зондирования ионосферы с помощью спутников [2-5], начался новый «бум» в области исследования ионосферных возмущений. Радиотомографические методы позволили получить двумерные сечения ПИВ. Существующие радиотомографические системы и сеть приемников ОР8/ГЛОНАСС дают возможность исследовать данное явление непрерывно в широких временных и пространственных масштабах. В течение почти 50 лет был накоплен огромный экспериментальный материал об ионосферных возмущениях, вызванных источниками разной природы. Среди таких источников особую роль играют сильные землетрясения, ядерные испытания, промышленные взрывы, сверхзвуковые полеты ракет и т.д.

Исследование отклика ионосферы на такие антропогенные и природные источники имеет большое значение для физики атмосферы, так как это позволяет подробно изучить различные физические процессы, происходящие в атмосфере. Земная атмосфера является уникальной лабораторией для изучения многих сложных физических процессов, таких как, генерация и распространение крупномасштабных волн, ударных волн, возникновение неустойчивых образований в ионосферной плазме и т.д. Изучение данной проблемы актуально в связи с необходимостью обоснования надежных сигнальных признаков техногенных воздействий (запуски ракет, несанкционированные взрывы и подземные ядерные испытания) [5,11]. Кроме того, результаты этих работ могут быть использованы для идентификации откликов от цунами и предвестников землетрясений в ионосфере [6,7].

Несмотря на огромное количество экспериментальных данных, лишь в относительно небольшом количестве работ были сделаны попытки с теоретических позиций объяснить наблюдаемые возмущения в атмосфере и ионосфере от наземных и атмосферных источников [6-12]. Следует отметить, что работы в данном направлении сыграли существенную роль в понимании механизма связей в системе литосфера - атмосфера - ионосфера [1]. Эти исследования внесли большой вклад в развитие теории распространения АГВ в атмосфере и т.д. Несмотря на это, до настоящего времени отсутствует полная интерпретация данных наблюдений [4,6]. Как правило, в предшествующих теоретических работах акустические волны (АВ) [11], внутренние гравитационные волны (ВГВ) [9], нелинейность среды, влияния вязкости, стратификации и зонального ветра на распространение волн рассматривались обособлено друг от друга. В результате этого, ни одна из существующих моделей не может объяснить весь спектр атмосферных и ионосферных возмущений, возбуждаемых источниками разной природы. Что

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

В последнее десятилетие в связи с увеличением быстродействия компьютеров и развитием вычислительной гидродинамики зародилось новое направление в физике атмосферы - исследование распространения атмосферных волн с помощью численного решения нелинейных уравнений геофизической гидродинамики [13]. Применение таких численных методов позволяет учесть совместно все факторы, влияющие на распространения АГВ с конечной амплитудой в реальной диссипативной среде.

Цели работы:

• Разработка численного метода для моделирования распространения АГВ с конечной амплитудой на большие горизонтальные и вертикальные расстояния;

• Моделирование генерации АГВ и волнообразных ионосферных возмущений от разного типа источников: наземных импульсных и длиннопериодных источников, сверхзвуковых полетов ракет;

• Сопоставление результатов моделирования с данными наблюдений.

Новизна результатов:

• Предложена система уравнений геофизической гидродинамики и граничные условия, позволяющие провести численное моделирование генерации и распространения АГВ от источников разной природы;

• Разработан численный алгоритм для решения системы уравнений вместе с соответствующими начальными и граничными условиями с помощью конечно-разностного метода;

• Промоделирована генерация волн от разного типа источников: сильных землетрясений, подземных ядерных взрывов, поверхностных волн Рэлея, длиннопериодных наземных источников, сверхзвуковых движений ракетоносителей;

• Впервые с помощью конечно-разностного метода были получены вместе две ветви спектра акустико-гравитационных волн - акустические волны и внутренние гравитационные волны;

• Впервые были построены кривые чувствительности (отклика) нейтральной атмосферы на воздействие поверхностных источников с разными периодами;

• Впервые были промоделированы волнообразные вариации электронной плотности в ионосфере с учетом геомагнитного поля на больших временных и пространственных масштабах, генерированные вышеуказанными источниками;

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

Научная и практическая ценность работы. Используя разработанную численную модель, можно исследовать свойства АГВ, генерированных источниками разных видов, такими как движение солнечного терминатора, экваториальные и полярные токовые системы, тропические циклоны и т.д.

Полученные численные результаты могут быть использованы для выделения и идентификации отклика в ионосфере от слабых длиннопериодных колебаний земной поверхности, наблюдаемых в некоторых случаях перед землетрясениями, зарегистрировать которые крайне трудно традиционными сейсмическими методами. Развитая модель позволяет исследовать нелинейный режим распространения АГВ и распространение слабых ударных волн в атмосфере. Результаты работы могут быть использованы в будущем для разработки и постановки экспериментов по наблюдению атмосферных и ионосферных возмущений, генерированных источниками разной природы. Полученные графические и мультимедийные материалы могут применяться в учебных курсах как наглядное пособие для изучения свойств АГВ в атмосфере.

Апробация работы. Результаты работы докладывались на XII Всероссийской школе-конференции по дифракции и распространению волн (Москва, 2001), XX Всероссийской научной конференции по распространению радиоволн (Нижний Новгород, 2002), на V сессии Байкальской молодежной школе по фундаментальной физике (Иркутск, 2002), Международной конференции «Потоки и структуры в жидкостях» (Санкт-Петербург, 2003), Генеральной Ассамблее Европейского Геофизического Союза (Nice, 2003), LVIII Научной Сессии НТОРЭС им. А.С. Попова (Москва, 2003), международной конференции и школе для молодых ученых «Вычислительно-информационные технологии для наук об окружающей среде» (Томск, 2003). По теме диссертации опубликовано 12 работ в отечественных и зарубежных изданиях.

Защищаемые положения.

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

2. Численный метод для решения соответствующих систем уравнений геофизической гидродинамики.

3. Результаты моделирования пространственных и временных распределений возмущений плотности и температуры нейтральной атмосферы, скорости гидродинамических частиц, плотности волновой энергии АГВ, генерированных наземными и атмосферными источниками.

4. Результаты моделирования пространственных и временных распределений возмущений электронной концентрации в ионосферной плазме, генерированных наземными и атмосферными источниками при разных моделях нейтральной атмосферы и ориентации магнитного поля.

5. Результаты моделирования вариаций полного электронного содержания в разных направлениях в возмущенной ионосфере, для разного рода источников.

6. Кривые чувствительности (отклика) нейтральной атмосферы на наземные источники с разными периодами.

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

Структура и объем работы. Диссертация состоит из введения, четырех глав,

заключения и списка цитируемой литературы. Работа содержит_страниц текста,

_рисунков.

Краткое содержание работы.

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

В первой главе представлен обзор теоретических и экспериментальных работ по исследованиям АГВ в атмосфере и ионосфере. Глава состоит из трех параграфов.

Первый параграф посвящен общим свойствам распространения АГВ в атмосфере. Приводятся основные спектральные характеристики АВ и ВГВ, захваченных АГВ, механизмы диссипации волн и т.д. Более энергетическая часть этих волн - ВГВ играет большую роль в циркуляции атмосферы, и их параметризация имеет важное значение в мезомасштабных и глобальных моделях атмосферы [14]. Кроме того, ВГВ являются одним из основных механизмов передачи энергии из нижних слоев атмосферы в термосферу.

Второй параграф посвящен механизмам генерации АГВ. Источники генерации могут быть антропогенными или природными. Характеристики АГВ также зависят от свойств источников. Такие источники, как конвективная неустойчивость, сдвиговая неустойчивость, орография, погодные фронты и т.д. постоянно генерируют ВГВ в нижней атмосфере. Такие источники, как сильные землетрясения, мощные взрывы, запуски ракет и т.д. способны генерировать широкий спектр волн - от высокочастотных инфразвуковых до длиннопериодных ВГВ, которые распространяются на тысячи километров без заметного затухания.

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

Во второй главе рассмотрено построение численной модели.

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

дР=-#- кр0+?')«]-#■ [(р0+рм-и дР

д/ дх дх

ди ди тт ди ди ди0 1 др'

— = -и--и 0--w--w-- - А1--—

д/ дх дх дх дх р0 + р дх

дw дТ'

-и

дw

-и

дх д Т'

-и

дw дw

дх

■ - w-

дх

-А

1

0

др'

- А2Хи р

1 Р0 + Р' дх

- А3 ^ ДХ

Р0 +Р

-и0

дТ' д((0 + Т')

■ - w■

д/ дх дх

р' = (р0 Т ' + р'Т0 + рТ')/

дх

-А

ди дw

+

V дх дх

+ А.

гд2Т' д2 Т л

дх2

дх7

тп

Здесь, первое уравнение - это уравнение неразрывности, второе и третье -уравнения Навье-Стокса, четвертое - уравнение сохранения энергии, и, наконец, последнее - это уравнение состояния идеального газа, р- плотность, Т -температура, р - давление, и0 - х компонента скорости зонального ветра, т0-

относительная молекулярная масса воздуха. Индексами «0» и « ' » обозначены невозмущенные и возмущенные части термодинамических параметров, соответственно. Здесь и и w горизонтальные и вертикальные компоненты скорости движения частиц атмосферы, соответственно, А1-5- постоянные. Для таких

относительно быстрых движений несущественна сила Кориолиса, которая не учитывалась. В данной модели сила вязкости вводится как Рэлеевское трение

Рл = -аи, а х = аР - коэффициент трения, который растет с высотой. При затухании энергии таких больших волновых движений основную роль играет теплопроводность, т.е. в уравнении для температуры диссипативный член равен Qd = кАТ, где к - коэффициент теплопроводности воздуха. Как видно из (1), в

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

д 2Т

решении уравнений будут генерироваться волны из-за того, что —^ Ф 0. А это с

дх

физической точки зрения недопустимо. При моделировании использовались профили фоновой плотности и температуры атмосферы из модели М818Б-90. Значения скорости зонального ветра были взяты из модели ИШМ-93. Таким образом, после задания параметров атмосферы и соответствующих начальных и граничных условий можно решать систему уравнений (1).

Второй параграф посвящен разработке численного алгоритма для решения (1). При выборе численного метода мы учитывали следующие особенности данной задачи:

• Решение должно содержать медленные и быстрые волновые движения.

• Метод должен быть устойчивым к резким градиентам плотности атмосферы.

• Обеспечение устойчивости решения усложняется присутствием нелинейных компонент в уравнениях.

• Большой объем вычислений из-за размера области интегрирования и длительности времени распространения волн.

Учитывая все вышеперечисленное, после анализа и сопоставления разных методов был выбран явный конечно-разностный метод Flux Corrected Transport [15]. В разработанном нами численном алгоритме используются основные свойства этого метода. Суть метода состоит в том, что для обеспечения устойчивости решения на первом этапе вводится искусственная диффузия, а потом с помощью антидиффузии минимизируется численная диффузия. Так как все переменные в (1) могут быть положительными и отрицательными, мы не использовали ограничения на антидиффузионные слагаемые для обеспечения положительности. Таким образом, мы освобождаемся от таких проблем, как синхронизация потоков и большие затраты на время расчета. Для обеспечения устойчивости и высокой точности мы аппроксимировали конвективные, адвекционные и нелинейные слагаемые в (1) особыми методами [15]. Для решения двумерной и трехмерной задачи использовался метод расщепления шага по времени для координатных направлений. Граничные условия состоят из двух частей: первая часть обеспечивает вход возмущения от источника в данную модель, вторая часть - уход волны из расчетной области. Начальные условия предполагались нулевыми.

В третьем параграфе приведен вывод формулы для вычисления изменения электронной концентрации в плоскопараллельной ионосфере, вследствие столкновений с нейтральными частицами во время прохождения АГВ. Эта формула получена из уравнении неразрывности для заряженных частиц [3]:

11Г , ч dNe <•/, , j2 \ 7 ЛТ с(, 2 дм , , i дм aw j2 aw 1 ANe (x, z, t) = -—±\[bxbzu + bz w)dz - Ne\\b2x — + bxbz — + — + b2z — dr (2) dz t 11 ax ^ ^

au aw

ydz dx j

2 aw

az

0

0

где Ne- концентрация электронов в невозмущенной ионосфере, bx и bz горизонтальные и вертикальные компоненты единичного вектора магнитного поля, соответственно. В конкретных расчетах использовался профиль фоновой ионосферы, состоящий из двух параболических слоев - E и F (с максимумом на 300 км). Для наблюдения за состоянием ионосферы имеет очень большое значение определение изменения полного электронного содержания (ПЭС) в разных направлениях, т.е. между приемником и искуственным спутником земли (ИСЗ) [4] (см. Рис. 1):

satellite

ATEC = \ANdr (3)

receiver

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

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

вертикальной компоненты (горизонтальные перемещения не учитываются)

условие будет условием сшивания над твердой

скорости граничное поверхностью земли:

w

81П

т

- )

Р

■ ехр

х - V

Бх

*о < * < г0 + Р

(4)

Хс

где wm - амплитуда, ^ - момент начала возмущения, Р - период, горизонтальная координата эпицентра источника, Бх - Гауссовский масштаб, характеризующий горизонтальный размер источника. Для сильного землетрясения или подземного ядерного испытания мы брали wm= 10 м/с, Р= 1 с, Бх = 20 км. При выборе значений этих параметров, мы учитывали то обстоятельство, что на практике определить точные значения параметров импульса возмущения не так просто, и они зависят от конкретного случая. Но так как нашей главной целью является изучение общих свойств возмущений от таких источников, мы принимали некоторые приближенные значения параметров источника. Что касается наземных промышленных взрывов, то их можно учитывать в нашей модели с помощью задания эквивалентного подземного источника.

Во втором параграфе представлены основные результаты, полученные при решении уравнений (1) для наземного источника с вышеперечисленными параметрами. Таким образом, система уравнений (1) вместе с начальными и граничными условиями решалась с помощью составленного нами пакета компьютерных программ. Параметры расчетной сетки были выбраны следующие: шаг по высоте 5 км, по горизонтальной координате 10 км и 0.1 сек по времени. Размеры моделируемой области составляли 4000 км и 500 км по ширине и высоте, соответственно.

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

На Рис.2 представлена зависимость горизонтальной скорости гидродинамических частиц на высоте И=100 км, от времени и расстояния. Она рассчитана на основе июльского профиля температуры модели М818Б-90, без учета ветра для вышесказанного источника. Здесь четко видны АВ и ВГВ. Как видно из графика, период ВГВ растет от 300 сек до 1200 сек линейно по мере удаления от источника в горизонтальном направлении. Эта связано тем, что волны с большими периодами движутся наклонно к вертикали. Период акустических волн около 200 сек. Точка насыщения, т.е. высота максимальной амплитуды АГВ находиться на высоте около 200 км. Вертикальная скорость, возмущения плотности, давления и температуры в волне, также показывают подобную зависимость. Для АВ

вертикальная скорость, а для ВГВ горизонтальная скорость играют доминирующую роль.

В третьем параграфе приведены результаты вычислений ионосферных возмущений, вызванных наземными импульсными источниками. На Рис. 3-4 представлены пространственные распределения вариации электронной плотности в ионосфере через 5000 с после землетрясения или взрыва. После толчка в течении 15 мин возмущение доходит до верхней атмосферы. Учет зонального ветра нарушает симметрию отклика атмосферы и ионосферы. Этот случай продемонстрирован на Рис. 3, где ось OX направлена на восток. Здесь наблюдаются среднемасштабные ПИВ с длинами волн около 100 км, распространяющиеся от источника с горизонтальной фазовой скоростью около 300 м/с.

Проведенные исследования показали сильную качественную и количественную модификацию ионосферного отклика при изменении ориентации геомагнитного поля. Как видно из Рис. 4, в случае когда магнитное поле имеет заметную вертикальную компоненту (наклонение /=45°), получается совсем другая картина. Эти результаты могут быть также применены к интерпретации атмосферных возмущений вызванных наземными взрывами.

Что касается данных, полученных при трансионосферном зондировании верхней атмосферы сигналами спутниковой радионавигационной системы GPS, то анализ полученных результатов показывает, что характеристики этих сигналов сильно зависят от месторасположения приемника, угла места спутника и т.д. Как показывают наши результаты, диапазон частот этих возмущений может быть очень широким. Если посмотреть на Рис. 5a, мы увидим, что этот график качественно хорошо описывает высокочастотные (период «200 с) возмущения вертикального

^_____о

ПЭС (т.е. угол места спутника 0=90 ), зарегистрированные авторами [4,5] во время сильных землетрясений. Этот график соответствует случаю, когда приемник находиться над эпицентром. В Рис. 5b представлена зависимость вариации вертикального ПЭС, когда приемник находиться на расстоянии 800 км от источника. Здесь видна ВГВ с периодом приблизительно 1000 с. Длиннопериодные колебания в ионосфере, предсказанные нами, наблюдались с помощью томографических, доплеровским зондированием, и др. методами во время землетрясений и взрывов [2,3,12]. Периоды волн, горизонтальные фазовые и групповые скорости распространения волн в нашей модели и экспериментальных данных почти совпадают.

В четвертом параграфе представлены результаты моделирования возбуждения атмосферы и ионосферы от длиннопериодных наземных источников. Как известно, неглубокие землетрясения генерируют длинные поверхностные волны Релея, которые распространяются на большие расстояния. В модели в качестве источника бегущей волны Релея предполагалась одиночная синусоидальная волна с периодом 50 си длиной волны 70 км, которая распространяется со скоростью 1400 м/с по поверхности земли. Для амплитуды скорости вертикальных колебаний земной поверхности принималось значение wm= 1 мм/с. Ионосферный отклик через 2000 с после начала распространения волны отточки О в направлении оси ОХ (см. Рис.1) продемонстрирован на Рис. 6. В этом

Рис. 1. Схематическая диаграмма, иллюстрирующая, генерацию АГВ наземными и атмосферными источниками и возможности спутникового радиозондирования. Показаны зональный ветер и магнитное поле.

Рис. 2. Изменение горизонтальной скорости частиц атмосферы в зависимости от времени и горизонтального расстояния.

Рис.3. Распределение вариаций электронной концентрации в ионосфере, после наземного возмущения, где учтен ветер и горизонтальное магнитное поле.

Рис.4. Распределение вариаций электронной концентрации в ионосфере, после наземного возмущения (наклонение магнитного поля 1=45 ).

случае ПИВ распространяются со сверхзвуковой фазовой скоростью. Амплитуда возмущения электронной плотности такого же порядка, как в случае импульсного источника. По периодам атмосферные волны, генерированные волнами Рэлея, принадлежат к высокочастотным АВ. Этот результат подтверждается результатами аналитических вычислений, проведенных в работе [8].

Кроме того, в природе существуют длиннопериодные источники, такие как цунами, колебания земной поверхности и т.д. [6-8]. Длиннопериодные колебания земной поверхности могут быть основным механизмом генерации ВГВ в сейсмически активных регионах. С этой целью мы брали источник в форме (4), с амплитудой 1 мм/с, но с разными периодами. Анализ отклика нейтральной атмосферы на такие источники показывает, что роль ВГВ в отклике растет с ростом периода источника. Результаты показывают, что в зависимости от расстояний в горизонтальном и вертикальном направлениях от источника, чувствительность атмосферы существенно меняется. На Рис. 7 показаны две резонансные кривые, рассчитанные на основе сравнения амплитуд скоростей над источником и на горизонтальном расстоянии 100 км от него. Очевидно, что первая кривая имеет максимум на частоте Брента-Вяйсяля. Из второй кривой видно, что максимум отклика смешается в сторону длинных периодов с удалением от источника в горизонтальном направлении. Что касается, амплитуды возмущений, то они сравнимы с откликами в случае импульсных источников, несмотря на то, что амплитуды источников отличаются в 104 раз. Таким образом, длиннопериодные наземные источники способны генерировать интенсивные ВГВ в атмосфере в больших масштабах. Если учесть, что в сейсмически активных регионах протяженность таких источников велика, тогда очевидно, что они могут сильно возбудить верхнюю атмосферу, и их можно зарегистрировать радиофизическими методами.

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

возможности решения трехмерной задачи и обработки полученных результатов. Поэтому здесь представлены результаты, полученные на более грубых разностных сетках в меньшей пространственной области. На Рис. 8 показано пространственное распределение вертикальной скорости гидродинамических частиц в горизонтальной плоскости, где ось ОХ направлена под углом 45° к востоку. В

V-/ гл с»

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

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

Первый параграф посвящен анализу свойств УВ, излучаемых сверхзвуковыми полетами ракет. Общеизвестно, что генерация АГВ происходит в основном во время горизонтального сверхзвукового полета ракеты после запуска [12]. При сверхзвуковом обтекании ракеты на больших расстояниях от нее вызываемые ею возмущения слабы, и поэтому их можно рассматривать как цилиндрическую звуковую волну, расходящуюся от оси, проходящей через ракету параллельно направлению обтекания. Мы ввели источник возмущения в расчетную модель следующим образом: пусть вертикальная координатная плоскость - Х02 (см. Рис. 1.) перпендикулярна к траектории горизонтального полета ракеты. Предположим, что с левой границы (участок СБ на Рис. 1.) в расчетную область входит одиночный сильно нелинейный акустический импульс, который далее превращается в пакет АГВ. Цилиндрический звуковой импульс в первом приближении можно аппроксимировать с помощью функции Гаусса. Таким образом, акустический импульс вводиться в модель в такой форме, что выражение для горизонтальной компоненты колебательной скорости частиц воздуха на оси 02 определяется формулой:

и = и 81И

т

2п( - X,)

Р

■ ехр

' ^ г ^2Л

т

У

, V Ог у , V 4 у у

Х() < X < ¿0 + Р (5)

где ит - амплитуда, Х0 - момент прихода импульса, Р - период, гт - высота оси звукового цилиндра от земной поверхности, Ог - Гауссовский масштаб,

о о гл

характеризующий поперечный размер цилиндра. Значения плотности и температуры в волне можно вычислить с помощью формул, описывающих простые нелинейные акустические волны [12]. Эти граничные условия обеспечивают передачу возмущения от среды вокруг траектории ракеты в рассматриваемую нами область атмосферы. Что касается моделирования второй области атмосферы, находящейся слева от траектории ракеты, то в этом случае такие же граничные условия будут на правой границе расчетной области. Граничные условия на других стенках области будут аналогичными, как в случае для наземных источников.

Второй параграф посвящен результатам моделирования генерации АГВ от звукового цилиндра вокруг сверхзвуковой ракеты. Здесь мы принимали

приближенные значения для параметров акустического импульса [12]: um= 260 м/с, P= 26 c, Dz = 10 км. Результаты моделирования показывают, что в этом случае генерируется широкий спектр волн - АВ, ВГВ и захваченные АГВ. Основной отклик - это длиннопериодные волны. Примечательным фактом является то, что скорость распространения возмущений близка к скорости звука в термосфере, около 800 м/с. Как в предыдущих случаях горизонтальный ветер слабо влияет на характер отклика.

В третьем параграфе приведены результаты вычислений ионосферных волнообразных возмущений, генерированных во время полета ракет. На Рис. 9 изображен график изменения электронной концентрации в зависимости от координат в момент /=5000 с, после начала возмущения, где учитывалась реальная стратификация атмосферы. Угол наклонения магнитного поля 75°. Из графика очевидно, что в ионосфере генерируются крупномасштабные ПИВ с масштабами сотни километров. В области мезосферно-термосферного волновода видны захваченные волны, которые не появляются в решении, если принимать атмосферу изотермической.

Сравнение результатов с данными радиотомографических реконструкций, полученных во время запуска ракеты с космодрома Плесецк в 1991 г. [2], показывает,что данная модель качественно описывает пространственные свойства ионосферных возмущений от такого рода источников. Кроме того, низкочастотные ВГВ, предсказанные данной моделью, наблюдались также с помощью методов некогерентного рассеяния, доплеровского зондирования и т.д. во время запусков ракет с космодромов Байконур и Kennedy Space Center [12]. Что касается данных по вариации ПЭС, здесь модель также предсказывает появление N-образного возмущения по времени [4].

В Заключении сформулированы основные результаты, полученные в ходе выполнения данной работы. Была разработана численная модель распространения АГВ в двумерной и трехмерной сжимаемой атмосфере, учитывающая стратификацию атмосферы, зональный ветер, диссипативные эффекты и нелинейность. На основе двумерной модели были промоделированы генерация АГВ от разного вида источников. Были вычислены временные и пространственные распределения возмущений электронной плотности в ионосфере, генерированные во время распространения АГВ. Рассчитаны интегральные вариации электронной концентрации в разных направлениях для сопоставления с экспериментальными данными. Результаты показывают, что разработанная численная модель достаточно хорошо предсказывает наблюдаемые возмущения в ионосфере от рассматриваемых источников. Анализ полученных результатов демонстрирует, что отклики от наземных и атмосферных источников сильно отличаются друг от друга, так как вторые возбуждают верхний волновод атмосферы и находятся прямо в ионосфере. Сравнение численных и экспериментальных результатов показывает, что для углубленного изучения данной проблемы, следует привлекать, комплекс различных методов наблюдения. Среди них самым информативным и подходящим является

Рис. 6. Распределение вариации электронной концентрации в ионосфере, генерированная при распространении волны Рэлея. (наклонение магнитного поля 1=0°)

a)

б)

1

ф 0.8

-О

л

Q.

Е 0.6

га

(D

> 0.4

■HJ га

ш

L_

0.2

300 600

period, s

300 . . 600 period, s

Рис. 7. Кривые чувствительности атмосферы к воздействию источников с разными периодами: а) на высоте 50 км от источника; б) на той же высоте, но на горизонтальном расстоянии 100 км от источника

Ь= 140 кт; 1= 1800 э

0.04

0.03

« 0.02 Е

> 0.01

-0.01 1000

у, кт

1000

500

500

0 О

х, кт

Рис. 8. Распределение вертикальной скорости частиц атмосферы в горизонтальной плоскости над импульсным источником.

Рис. 9. Распределение возмущения электронной концентрации в ионосфере, во время сверхзвукового полета ракеты. Траектория ракеты перпендикулярна к координатной плоскости и находиться слева от нее.

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

Из полученных результатов следует, что для детальной интерпретации экспериментальных результатов необходимо учитывать следующие факторы:

• геометрию и временные характеристики источника;

• состояния атмосферы и ионосферы над местом наблюдения;

• ориентацию геомагнитного поля над местом наблюдения;

• особенности данного экспериментального метода.

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

Результаты диссертации опубликованы в следующих работах:

1. Ахмедов Р.Р., Куницын В.Е. Численный метод решения задачи распространения акустико-гравитационных волн в атмосфере до ионосферных высот// Вестник Московского Университета. Серия 3. Физика. Астрономия. 2003. № 3. С. 38- 42.

2. Ахмедов Р.Р., Куницын В.Е. Моделирование ионосферных возмущений, вызванных землетрясениями и взрывами// Геомагнетизм и аэрономия. 2004. Т. 44. № 1. С. 1 - 8.

3. Ахмедов Р.Р., Кадиров Ф.А., Куницын В.Е. Моделирование атмосферных возмущений, вызванных землетрясениями// Изв. НАН Азерб. Сер. Наук о Земли. 2004. № 1. с. 59 - 68.

4. Ahmadov R.R., Kunitsyn V.E. Numerical simulation of atmospheric propagation of acoustic-gravity waves caused by impulsive sources on the earth// Geophysical Research Abstracts. 2003. V. 5. 12451.

5. Ahmadov R.R., Kunitsyn V.E. The atmospheric disturbances generated by the impulsive sources on the earth// International conference. Fluxes and Structures in Fluids. Sanct Petersburg. P. 6. 2003. 23-26 June.

6. Ахмедов Р. Р. Распространение внутренних гравитационных волн в стратифицированной атмосфере// Труды XII Всероссийской школы-конференции по дифракции и распространению волн. Москва. 19-23 декабря. 2001 г.

7. Ахмедов Р.Р. Численное моделирование ионосферных возмущений, вызванных землетрясениями и взрывами// ВНКСФ-9. Сборник тезисов. Т. 2. 871-873.

8. Ахмедов Р.Р. Численное моделирование распространения акустико-гравитационных волн в атмосфере// Труды V сессии молодых ученых «Гелио- и геофизические исследования». Иркутск. 2002 г. с. 71 - 73.

9. Ахмедов Р.Р. Численное моделирование распространения акустико-гравитационных волн в атмосфере, вызванных импульсными источниками на Земле// Международная конференция и школа молодых ученых. Томск. 1-11 сентября. 2003 г.

10. Ахмедов Р.Р., Куницын В.Е. Моделирование ионосферных возмущений, генерируемых наземными импульсными источниками// Труды XX Всероссийской конференции по распространению радиоволн. Нижний Новгород. 2-4 июля. 2002 г.

11. Ахмедов Р.Р., Куницын В.Е. Моделирование ионосферных возмущений, вызванных землетрясениями и взрывами// Труды LVIII научной сессии, посвященная дню радио. Москва. 2003.

12. Ahmadov R.R., Kunitsyn V.E. Modeling of acoustic-gravity waves generation and propagation in the atmosphere during rocket launchings// International Geomagnetism and Aeronomy (in press).

ЦИТИРУЕМАЯ ЛИТЕРАТУРА:

1. Григорьев Г.И. Акустико-гравитационные волны в атмосфере Земли (обзор)// Изв. ВУЗов Радиофизика. 1999. Т. XLII. № 1. с. 3 - 25.

2. Kunitsyn V., Tereshchenko E. Ionospheric Tomography. Springer-Verlag. 2003.272p.

3. Андреева Е.С., Гохберг М.Б., Куницын В.Е., Терещенко Е.Д., Худукон Б.З., Шалимов С. Л. Радиотомографическая регистрация возмущений ионосферы от наземных взрывов// Космич. исслед. 2001. т. 39. № 1. с. 13 - 17.

4. Calais E., Minster J.B. GPS, earthquakes, the ionosphere, and the Space Shuttle// Phys. of the Earth and Planet. Inter. 1998. 105. 167 - 181.

5. Afraimovich E.L., Perevalova N.P., Plotnikov A.V., Uralov A.M. The shock-acoustic waves generated by earthquakes// Ann. Geophysicae. 2001. 19. 395 - 409.

6. Artru J., Lognonne P., Blanc E.. Normal modes modeling of post-seismic ionospheric oscillations// Geophys. Res. Let. V. 28. No. 4. P. 697 - 700. 2001.

7. Перцев Н.Н., Шалимов С.Л. Генерация атмосферных гравитационных волн в сейсмически активном регионе и их влияние на ионосферу// Геомагнетизм и аэрономия. 1996. 36. 111 - 118.

8. Голицын Г.С., Кляцкий В.И. Колебания в атмосфере, вызываемые движениями земной поверхности// Изв. АН СССР. ФАО. 1967. Т. 3. № 10. 1045 - 1052.

9. Гаврилов Н.М. Внутренние гравитационные волны и их воздействие на среднюю атмосферу и ионосферу. Дис... докт. ф.-м. наук. 04.00.22. ЛГУ. 1988.

10. Шувалов В.В. Динамические процессы в атмосфере вызванные сильными импульсными возмущениями. Дис. докт. ф.-м. наук. 04.00.23. Москва. 1999.

11. Drobzheva Ya.V., Krasnov V.M. The acoustic field in the atmosphere and ionosphere caused by a point explosion on the ground// J. Atmos. Sol.-Terr. Phys. 2003. V. 65. No. 3. 369 - 377.

12. Нагорский П.М. Модификация F-области ионосферы мощными импульсными источниками волн в нейтральном газе. Дис. докт. ф.-м. наук. 04.00.23. ТГУ. 1998. 365 с.

13. Zhang S.D., Yi F. A numerical study of propagation characteristics of gravity wave packets propagating in a dissipative atmosphere// J. Geophys. Res. 2002. V. 107. D14. 1 - 9.

14. Fritts D.C., Alexander M.J. Gravity wave dynamics and effects in the middle atmosphere// Reviews of geophysics. 2003. V. 41. N. 1. 1 - 64 .

15. Durran D. Numerical methods for wave equations in Geophysical Fluid Dynamics. Springer - Verlag New York, Inc. 1999. 465 p.

Содержание диссертации, кандидата физико-математических наук, Ахмедов Раван Рамин оглы

СПИСОК СОКРАЩЕНИЙ.

ВВЕДЕНИЕ.

ГЛАВА I. ВОЛНООБРАЗНЫЕ ДВИЖЕНИЯ В АТМОСФЕРЕ И ИОНОСФЕРЕ ЗЕМЛИ

1.1. Акустико-гравитационные волны в земной атмосфере.

1.2. Механизмы генерации акустико-гравитационных волн.

1.3. Ионосферные волновые возмущения, генерированные мощными источниками.

ГЛАВА II. ПОСТРОЕНИЕ МОДЕЛИ И РАЗРАБОТКА МЕТОДА РЕШЕНИЯ УРАВНЕНИЙ

2.1. Уравнения геофизической гидродинамики, описывающие 46 распространения волн в атмосфере.

2.2. Разработка численного алгоритма для решения уравнений.

2.3. Моделирование ионосферных возмущений, генерированных 62 прохождением атмосферных волн.

2.4. Краткое содержание главы. Основные выводы.

ГЛАВА III. МОДЕЛИРОВАНИЕ АТМОСФЕРНЫХ И ИОНОСФЕРНЫХ ВОЗМУЩЕНИЙ ОТ НАЗЕМНЫХ ИСТОЧНИКОВ

3.1. Включение наземных импульсных источников в модель.

3.2. Генерация и распространение акустико-гравитационных волн, 69 вызванных наземными импульсными источниками.

3.3. Моделирование возмущений ионосферы наземными импульсными источниками, и сравнение результатов с данными наблюдений.

3.4. Отклик атмосферы и ионосферы на длиннопериодные наземные 86 источники.

3.5. Трехмерное моделирование.

3.6. Краткое содержание главы. Основные выводы.

ГЛАВА IV. МОДЕЛИРОВАНИЕ АТМОСФЕРНЫХ И ИОНОСФЕРНЫХ ВОЗМУЩЕНИЙ ПРИ ЗАПУСКАХ РАКЕТ

4.1. Ударные волны, генерируемые при запусках ракет.

4.2. Генерация и распространение акустико-гравитационных волн 101 при сверхзвуковых полетах ракет.

4.3. Моделирование возмущений ионосферы, вызванных запусками 108 ракет и сравнение результатов с данными наблюдений.

4.4. Краткое содержание главы. Основные выводы.

Введение Диссертация по наукам о земле, на тему "Численное моделирование генерации акустико-гравитационных волн и ионосферных возмущений от наземных и атмосферных источников"

Актуальность. В конце XIX века ученые начали обращать внимание на то, что сила тяжести и стратификация геофизических сред сильно модифицируют распространяющиеся в них звуковые волны. С тех пор начались первые теоретические исследования свойств распространения АГВ в атмосфере. Наряду с этим, с развитием радиотехнических средств, в середине XX столетия начались наблюдения ПИВ. Впервые в 1960м году ПИВ были интерпретированы как проявления АГВ в ионосфере. С тех пор с помощью различных методов ведутся интенсивные наблюдения АГВ в нижней и верхней атмосфере. В 90-х годах, с развитием дистанционных методов зондирования ионосферы с помощью спутников, начался новый «бум» в области исследования ионосферных возмущений. Радиотомографические методы позволили получить двумерные сечения ПИВ. Существующие радиотомографические системы и сеть приемников GPS/ГЛОНАСС дают возможность исследовать данное явление непрерывно в широких временных и пространственных масштабах. В течение почти 50 лет был накоплен огромный экспериментальный материал об ионосферных возмущениях, вызванных источниками разной природы. Среди таких источников особую роль играют сильные землетрясения, ядерные испытания, промышленные взрывы, сверхзвуковые полеты ракет и т.д.

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

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

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

Цели работы:

• Разработка численного метода для моделирования распространения АГВ с конечной амплитудой на большие горизонтальные и вертикальные расстояния;

• Моделирование генерации АГВ и волнообразных ионосферных возмущений от разного типа источников: наземных импульсных и длиннопериодных источников, сверхзвуковых полетов ракет;

• Сопоставление результатов моделирования с данными наблюдений.

Новизна результатов:

• Предложена система уравнений геофизической гидродинамики и граничные условия, позволяющие провести численное моделирование генерации и распространения АГВ от источников разной природы;

• Разработан численный алгоритм для решения системы уравнений вместе с соответствующими начальными и граничными условиями с помощью конечно-разностного метода;

• Промоделирована генерация волн от разного типа источников: сильных землетрясений, подземных ядерных взрывов, поверхностных волн Рэлея, длиннопериодных наземных источников, сверхзвуковых движений ракетоносителей;

• Впервые с помощью конечно-разностного метода были получены вместе две ветви спектра акустико-гравитационных волн - акустические волны и внутренние гравитационные волны;

• Впервые были построены кривые чувствительности (отклика) нейтральной атмосферы на воздействие поверхностных источников с разными периодами;

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

Научная и практическая ценность работы. Используя разработанную численную модель, можно исследовать свойства АГВ, генерированных источниками разных видов, такими как движение солнечного терминатора, экваториальные и полярные токовые системы, тропические циклоны и т.д. Полученные численные результаты могут быть использованы для выделения и идентификации отклика в ионосфере от слабых длиннопериодных колебаний земной поверхности, наблюдаемых в некоторых случаях перед землетрясениями, зарегистрировать которые крайне трудно традиционными сейсмическими методами. Развитая модель позволяет исследовать нелинейный режим распространения АГВ и распространение слабых ударных волн в атмосфере. Результаты работы могут быть использованы в будущем для разработки и постановки экспериментов по наблюдению атмосферных и ионосферных возмущений, генерированных источниками разной природы. Полученные графические и мультимедийные материалы могут применяться в учебных курсах как наглядное пособие для изучения свойств АГВ в атмосфере.

Апробация работы и публикации. Результаты работы докладывались на XII Всероссийской школе-конференции по дифракции и распространению волн (Москва, 2001), XX Всероссийской научной конференции по распространению радиоволн (Нижний Новгород, 2002), на V сессии Байкальской молодежной школе по фундаментальной физике (Иркутск, 2002), Международной конференции «Потоки и структуры в жидкостях»

Санкт-Петербург, 2003), Генеральной Ассамблее Европейского Геофизического Союза (Nice, 2003), LVIII Научной Сессии НТОРЭС им. А.С. Попова (Москва, 2003), международной конференции и школе для молодых ученых «Вычислительно-информационные технологии для наук об окружающей среде» (Томск, 2003). По теме диссертации опубликовано 12 работ в отечественных и зарубежных изданиях.

Защищаемые положения.

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

2. Численный метод для решения соответствующей системы уравнений геофизической гидродинамики.

3.Результаты моделирования пространственных и временных распределений возмущений плотности и температуры нейтральной атмосферы, скорости гидродинамических частиц, плотности волновой энергии АГВ, генерированных наземными и атмосферными источниками.

4.Результаты моделирования пространственных и временных распределений возмущений электронной концентрации в ионосферной плазме, генерированных наземными и атмосферными источниками при разных моделях нейтральной атмосферы и ориентации магнитного поля.

5.Результаты моделирования вариаций полного электронного содержания в разных направлениях в возмущенной ионосфере, для разного рода источников.

6. Кривые чувствительности (отклика) нейтральной атмосферы на наземные источники с разными периодами.

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

Структура и объем работы. Диссертация состоит из введения, четырех

Заключение Диссертация по теме "Физика атмосферы и гидросферы", Ахмедов Раван Рамин оглы

ЗАКЛЮЧЕНИЕ

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

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

Условия в нейтральной атмосфере и ионосфере также влияют на свойства отклика от наземных источников. В зависимости от времени года и направления зонального ветра меняются характеристики генерированных АГВ. Но сильнее всего на формирование картины ионосфеных возмущений влияет магнитное поле Земли. В зависимости от направления силовых линий геомагнитного поля, пространственное распределение ионосферных неоднородностей может усиливаться или подавляться. Если учесть, что ионное торможение, которое сильно влияет на движения нейтралов в ВА, зависит от ориентации магнитного поля, тогда становится ясным, что параметры АГВ на больших высотах также могут меняться в зависимости от направления магнитного поля.

Моделирование генерации атмосферных волн от сверхзвукового полета ракет показывает, что в этом случае генерируется АВ, ВГВ и захваченные АГВ в верхних волноводах атмосферы. В отличие от наземных источников, атмосферные источники, типа звукового цилиндра, образованного сверхзвуковым полетом ракеты, возбуждают более крупномасштабные ПИВ в ионосферной плазме. Главное отличие атмосферных источников от наземных состоит в том, что первые очень сильно возбуждают верхние волноводы атмосферы. В этом случае могут распространяться интенсивные ВГВ с большими скоростями, чего нельзя сказать о наземных источниках.

Из полученных результатов следует, что для детальной интерпретации экспериментальных данных необходимо учитывать следующие факторы: • геометрию и временные характеристики источника;

• состояние атмосферы и ионосферы над местом наблюдения;

• ориентацию геомагнитного поля над местом наблюдения;

• особенности данного экспериментального метода.

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

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

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

Также выражаю благодарность за огромную помощь в оформлении диссертации к.ф.-м.н. Бербеневой Н.А., Кожарину М.А., и моей супруге Ахмедовой И.Ф.

Библиография Диссертация по наукам о земле, кандидата физико-математических наук, Ахмедов Раван Рамин оглы, Москва

1. Адушкин В.В., Козлов С.И., Петров А.В. Экологические проблемы и риски воздействий ракетно-космической техники на окружающую природную среду. Справочное пособие. М.: Изд-во Анкил, 2000. 640 с.

2. Аки К., Ричарде П. Количественная сейсмология. Москва «Мир». Т. 1 и 2. 1983.

3. Альперович Л.С., Вугмейстер Б.О., Гохберг М.Б. и др. Об опыте моделирования магнитосферно-ионосферных эффектов при сейсмических явлениях// Докл. АН СССР. 1983. 269. № 3. 573 578.

4. Андреева Е.С., Гохберг М.Б., Куницын В.Е., Терещенко Е.Д., Худукон Б.З., Шалимов С.Л. Радиотомографическая регистрация возмущений ионосферы от наземных взрывов// Космич. исслед. 2001. т. 39. № 1. с. 13 -17.

5. Атмосфера Стандартная. Параметры. ГОСТ 4401 -81.1981.

6. Атмосфера. Справочник. Гидрометеоиздат. 1991.

7. Афраймович Э.Л., Косогоров Е.А., Плотников А.В. Ударно-акустические волны, генерируемые при запусках ракет и землетрясениях// Космические исследование. 2002. Т. 40. № 3. 261 275.

8. Афраймович Э.Л., Косогоров Е.А., Плотников А.В., Уралов A.M. Параметры ударно-акустических волн, генерируемых при землетрясениях// Физика Земли. 2001. № 6. 1 13.

9. Афраймович Э.Л., Чернухов В.В., Кирюшкин В.В. Пространственно-временные характеристики ионосферного возмущения, обусловленного ударно-акустическими волнами, генерируемыми при запусках ракет// Радиотехника и электроника. 2001. Т. 46. № 11. 1299 1307.

10. Ахмедов P.P. Распространение внутренних гравитационных волн в стратифицированной атмосфере// Труды XII Всероссийской школыконференции по дифракции и распространению волн. Москва. 19-23 декабря. 2001 г.

11. Ахмедов P.P. Численное моделирование ионосферных возмущений, вызванных землетрясениями и взрывами// ВНКСФ-9. Сборник тезисов. Т. 2. 871-873.

12. Ахмедов P.P. Численное моделирование распространения акустико-гравитационных волн в атмосфере// Труды V сессии молодых ученых «Гелио- и геофизические исследования». Иркутск. 16-21 сентября. 2002 г.

13. Ахмедов P.P. Численное моделирование распространения акустико-гравитационных волн в атмосфере, вызванных импульсными источниками на Земле// Международная конференция и школа молодых ученых. Томск. 1-11 сентября. 2003 г.

14. Ахмедов P.P., Кадиров Ф.А., Куницын В.Е. Моделирование атмосферных возмущений, вызванных землетрясениями// Изв. НАН Азерб. Сер. Наук о Земли. 2004. № 1. с. 59 68.

15. Ахмедов P.P., Куницын В.Е. Моделирование ионосферных возмущений, генерируемых наземными импульсными источниками// Труды XX Всероссийской конференции по распространению радиоволн. Нижний Новгород. 2-4 июля. 2002 г.

16. Ахмедов P.P., Куницын В.Е. Моделирование ионосферных возмущений, вызванных землетрясениями и взрывами// Труды LVIII научной сессии, посвященная дню радио. Москва. 2003.

17. Ахмедов P.P., Куницын В.Е. Моделирование ионосферных возмущений, вызванных землетрясениями и взрывами// Геомагнетизм и аэрономия. 2004. Т. 44. № 1.С. 1-8.

18. Ахмедов P.P., Куницын В.Е. Численный метод решения задачи распространения акустико-гравитационных волн в атмосфере доионосферных высот// Вестник Московского Университета. Серия 3. Физика. Астрономия. 2003. № 3. С. 38 42.

19. Белашов В.Ю. Динамика нелинейных внутренних гравитационных волн на высотах F-области ионосферы// Геомагнетизм и аэрономия. Т. 30. № 4. С. 637-641. 1990.

20. Брюнелли Б.Е., Намгаладзе А.А. Физика ионосферы. Москва. «Наука». 1988.

21. Габов С.А. Новые задачи математической теории волн. «Наука». 1998. 448 с.

22. Гаврилов Н.М. Внутренние гравитационные волны и их воздействие на среднюю атмосферу и ионосферу. Дис. докт. ф.-м. наук. 04.00.22. Ленинград. ЛГУ. 1988.

23. Гаврилов Н.М. Распространение внутренних гравитационных волн в стратифицированной атмосфере. Изв. АН СССР. ФАО. 1985. 21. 921 -927.

24. Гаврилов Н.М., Юдин В.А. Численное исследование вертикальной структуры внутренних гравитационных волн от тропосферных источников//Изв. АН СССР. ФАО. 1986. Т. 22. № 6. 563 572.

25. Голицын Г.С., Кляцкин В.И. Колебания в атмосфере, вызываемые движениями земной поверхности// Изв. АН СССР. ФАО. 1967. Т. 3. № 10. 1045- 1052.

26. Голицын Г.С., Романова Н.Н., Чунчузов Е.П. О генерации внутренних волн в атмосфере морским волнением// Изв. АН СССР. ФАО. 1976. Т. 12. №6. 669-673.

27. Госсард Э., Хук У. Волны в атмосфере. «Мир». 1975. 532 с.

28. Григорьев Г.И. Акустико-гравитационные волны в атмосфере Земли (обзор)// Изв. ВУЗов Радиофизика. 1999. Т. XLII. № 1. с. 3 25.

29. Григорьев Г.И., Савина О.Н. Об орографическом возбуждении акустико гравитационных волн// Изв. АН СССР. ФАО. 1991. Т. 27. № 5. 545 -549.

30. Данилов А.Д., Казимировский Э.С., Вергасова Г.В., Хачикян Г.Я. Метеорологические эффекты в ионосфере. «Гидрометеоиздат». Ленинград. 1987. 268.

31. Дикий JI.A. Теория колебаний земной атмосферы. 1969.

32. Дэвис К. Радиоволны в ионосфере. «Мир». 1973.

33. Еланский Н.Ф., Кожевников В.Н., Кузнецов Г.И., Волков Б.И. Влияние орографических возмущений на перераспределение озона в атмосфере при обтекании антарктического полуострова// Изв. РАН. ФАО. 2003. 39. № 1. 105- 120.

34. Казимировский Э.С., Кокоуров В.Д. Движение в ионосфере. Наука. 1979.

35. Карлов В.Д., Козлов С.И., Ткачев Г.Н. Крупномасштабные возмущения в ионосфере, возникающие при полете ракеты с работающим двигателем (обзор)// Космич. исслед. 1980. Т. 18. Вып.2. С. 266 277.

36. Кожевников В.Н. Возмущения атмосферы при обтекании гор. Москва. 1999. 160 с.

37. Куличков С.Н. Дальнее распространение инфразвуковых волн в атмосфере. Дис. докт. ф.-м. наук. 01.04.06. Москва. 1999.

38. Куницын В.Е., Терещенко Е.Д. Томография ионосферы. М.: Мир. 1991.

39. Ландау Л.Д., Лифщиц В.М. Гидродинамика. «Наука». 1988.

40. Линьков Е.М., Петрова Л.Н., Осипов К.С. Сейсмогравитационные пульсации Земли и возмущения атмосферы как возможные предвестники сильных землетрясений//Докл. АН СССР. 1990. 313. № 5. 1095 1098.

41. Липеровский В.А., Похотелов О.А., Шалимов С.Л. Ионосферные предвестники землетрясений. «Наука». 1992. 304 с.

42. Нагорский П.М. Анализ отклика КВ-радиосигнала на возмущения ионосферной плазмы, вызванные ударно-акустическими волнами// Изв. ВУЗов. Радиофизика. Т. XLII. №1. 36-44. 1999.

43. Нагорский П.М. Модификация F-области ионосферы мощными импульсными источниками волн в нейтральном газе: Дис. . докт. физ.-мат. наук. Томск: ТГУ. 1998. 325 с.

44. Нагорский П.М. О возмущениях электронной концентрации в ионосфере, вызываемых наземными взрывами// Изв. АН СССР. ФЗ. 1985. № 11.66-71.

45. Оран Э., Борис Дж. Численное моделирование реагирующих потоков: Пер. с англ. М. Мир. 1990. 660 с.

46. Перцев Н.Н., Шалимов С.Л. Генерация атмосферных гравитационных волн в сейсмически активном регионе и их влияние на ионосферу// Геомагнетизм и аэрономия. 1996. 36. 111 118.

47. Плотников А.В. Детектирование с помощью GPS решеток ударно-акустических волн, генерируемых при запусках ракет, землетрясениях, и взрывах. Дис. канд. ф.-м. наук. 01.04.03. Иркутск. 2001.

48. Погорельцев А.И., Перцев Н.Н. Влияние фонового ветра на формирование структуры акустико гравитационных волн в термосфере// Изв. АН. ФАО. 1995. № 6. 755 - 760.

49. Поттер Д. Вычислительные методы в физике. Москва. «Мир». 1975. 392 с.

50. Уизем Дж. Линейные и нелинейные волны. 1977. М. Мир.

51. Ратклифф Дж. Введение в физику ионосферы и магнитосферы. «Мир». 1975.

52. Романова Н.Н. О нелинейном распространении акустико-гравитационных волн в изотермической атмосфере. Изв. АН СССР. ФАО. 1971. Т. VII. № 12. 1251 1262.

53. Романова Н.Н., Чунчузова Е.П. О нелинейном затухании коротких внутренних гравитационных волн в верхней атмосфере// Изв. АН СССР. ФАО. 1982. 18. №2. 191-193.

54. Роуч П. Вычислительная гидродинамика. 618 с. М.: Мир. 1980.

55. Руденко О.В., Солуян С.И. Теоретические основы нелинейной акустики. 1975. Москва. «Наука».

56. Савина О.Н. Акустико-гравитационные волны в атмосфере с реалистичным распределением температуры// Геомагнетизм и аэрономия. 1996. 36. 104- 110.

57. Самарский А.А., Попов Ю.П. Разностные методы решения задач газовой динамики. Наука. 1979. 422 с.

58. Смертин В.М., Намгаладзе А.А. Исследование зависимости характеристик внутренних гравитационных волн от параметров источника// Геомагнетизм и аэрономия. 1981. Т. XXI. 302 308.

59. Таран В.И., Подьячий Ю.И., Смирнов А.Н., Герштейн Л.Я. Возмущения ионосферы после наземного взрыва по наблюдениям методом некогерентного рассеяния// Изв. АН СССР. Физика Земли. 1985. №11. 75 79.

60. Тетевин В.Б. Справочник необходимых знаний. Москва. 2000. 768 с.

61. Трухин В.И., Показеев К.В., Куницын В.Е., Шрейдер А.А. Основы экологической геофизики. 2000. 291 с.

62. Шалимов С.Л. Ионосферные неоднородности инициированные интенсивными магнитосферными токами и атмосферными волнами. Дис. докт. ф.-м. наук. 04.00.23. Москва. 1998.

63. Шувалов В.В. Динамические процессы в атмосфере вызванные сильными импульсными возмущениями. Дис. докт. ф.-м. наук. 04.00.23. Москва. 1999.

64. Яновский Б.М. Земной магнетизм. Ленинград. 1978. 592 с.

65. Afraimovich E.L., Kosogorov Е.А., Perevalova N.P., Plotnikov A.V. The use of GPS arrays in detecting shock-acoustic waves generated during rocket launchings//J. Atmos. Sol.-Terr. Phys. 2001. 1 17.

66. Afraimovich E.L., Perevalova N.P., Plotnikov A.V. Registration of ionospheric responses to shock-acoustic waves generated by carrier rocket launches// Geomagnetism and Aeronomy. 2002. V. 42. No. 6. P. 755 762.

67. Afraimovich E.L., Perevalova N.P., Plotnikov A.V., Uralov A.M. The shock-acoustic waves generated by earthquakes// Ann. Geophysicae. 2001. 19. 395 -409.

68. Ahmadov R.R., Kunitsyn V.E. Modeling of acoustic-gravity waves generation and propagation in the atmosphere during rocket launchings// International Geomagnetism and Aeronomy. (in press).

69. Ahmadov R.R., Kunitsyn V.E. Numerical simulation of atmospheric propagation of acoustic-gravity waves caused by impulsive sources on the earth// Geophysical Research Abstracts. 2003. V. 5. 12451.

70. Ahmadov R.R, Kunitsyn V.E. The atmospheric disturbances generated by the impulsive sources on the earth// International conference. Fluxes and Structures in Fluids. Sanct Petersburg. P. 6. 2003. 23-26 June.

71. Alexander M.J., Holton J.R. A model study of zonal forcing in the equatorial stratosphere by convectively induced gravity waves// American Meteorological Society. 408-419. 1997.

72. Arendt P.R. Ionospheric undulations following "Apollo 14" launching// Nature. 1971. V.231. P.438-439.

73. Artru J., Lognonne P., Blanc E. Normal modes modeling of post-seismic ionospheric oscillations// Geophys. Res. Let. V. 28. No. 4. P. 697 700. 2001.

74. Beynon W.J. Evidence of horizontal motion in region F2 ionization. 1948. Nature. 162. 887.

75. Blanc E., Rickel D. Nonlinear wave fronts and ionospheric irregularities observed by HF sounding over a powerful acoustic source// Radio. Sci. 1989. V. 24. N3. P. 279-288.

76. Boris J.P., Landsberg A.M., Oran E.S., Gardner J.H. LCPFCT A Flux -Corrected Transport Algorithm for Solving Generalized Continuity Equations. 1993.

77. Calais E., Minster J.B. GPS, earthquakes, the ionosphere, and the Space Shuttle// Physics of Earth and Planetary Interiors. 1998. 105. 167- 181.

78. Davies J., Archambeau C. Modeling of atmospheric and ionospheric disturbances from shallow seismic sources// Phys. of the Earth and Planet. Inter. 1998. 105. 183 199.

79. Didebulidze G.G., Pataraya A.D. Ionosphere F2-region under the influence of the evolutional atmospheric gravity waves in horizontal shear flow// J. Atmos. Sol.-Terr. Phys. 1999. V. 61. 479 -489.

80. Drobyazko I.N., Gavrilov N.M. Wave energy and momentum fluxes coming to the middle and upper atmosphere from Tropospheric mesoscale turbulence// Phys. Chem. Earth. 2001. V. 26. No. 6. 449 452.

81. Drobzheva Ya.V., Krasnov V.M. The acoustic field in the atmosphere and ionosphere caused by a point explosion on the ground// J. Atmos. Sol.-Terr. Phys. 2003. V. 65. No. 3. 369 377.

82. Ducic V., Artru J., Lognonne P. Ionospheric remote sensing of the Denali Earthquake Rayleigh surface waves// Geophys. Res. Lett. 2003. V. 30. No. 18.

83. Durran D. Numerical methods for wave equations in Geophysical Fluid Dynamics. Springer Verlag New York, Inc. 1999. 465 p.

84. Fitzgerald T.J. Observations of total electron content perturbations on GPS signals caused by a ground level explosion// J. Atmos. Sol.-Terr. Phys. 1997. V. 59. № 7. 829 834.

85. Francis S.H. Acoustic-gravity modes and large-scale traveling ionospheric disturbances of a realistic, dissipative atmosphere// J. Geophys. Res. 1973. V. 78. № 13.2278-2301.

86. Francis S.H. Global propagation of atmospheric gravity waves: a review// J. Atmos. Terr. Phys. 1975. 37. 1011 1054.

87. Friedman J.P. Propagation of internal gravity waves in a thermally stratified atmosphere//J. Geophys. Res. 1966. V. 71. № 4. 1033 1053.

88. Fritts D.C., Alexander M.J. Gravity wave dynamics and effects in the middle atmosphere//Reviews of geophysics. 2003. V. 41. N. 1. 1-64.

89. Gardner C.S. Reply to Hines' comments on "Testing theories of atmospheric gravity wave saturation and dissipation"// J. Atmos. Sol.-Terr. Phys. 1998. V. 60. 663 665.

90. Gardner C.S. Testing theories of atmospheric gravity wave saturation and dissipation//J. Atmos. Terr. Phys. 1996. V. 58. No. 14. 1575 1589.

91. Gavrilov N.M. Parametrization of momentum and energy depositions from gravity waves generated by Tropospheric hydrodynamic sources// Ann. Geophysicae. 1997. 15. 1570- 1580.

92. Hines C.O. Comments on the paper by C.S. Gardner "Testing theories of atmospheric gravity wave saturation and dissipation"// J. Atmos. Sol.-Terr. Phys. 1998. V. 60. 655-662.

93. Hines C.O. Internal atmospheric gravity waves at ionospheric heights// Can. J. Phys. 1960. V. 38. 1441 1481.

94. Носке К., Schlegel К. A review of atmospheric gravity waves and traveling ionospheric disturbances: 1982 1995// Ann. Geophysicae. 1996. 14. 917 -940.

95. Holton J.R., Beres J.H., Zhou X. On the vertical scale of gravity waves excited by localized thermal forcing//Notes and correspondence. 2002. 20192023.

96. Holton J.R., Millendore G.L., Zhou X. Forcing of secondary waves by breaking of gravity waves in the mesosphere// J. Geophys. Res. V. 107. NO. D7, 10. 1029/2001JD001204, 2002.

97. Huang C.S., Sofko G.J. Numerical simulations of midlatitude ionospheric perturbations produced by gravity waves// J. Geophys. Res. 1998. V. 103. NO. A4. 6977 6989.

98. Imamura Т., Ogawa T. Radiative damping of gravity waves in the terrestrial planetary atmospheres// Geophys. Res. Lett. 1995. V. 22. NO. 3. 267 270.

99. Jacobson A.R., Carlos R.C. Observations of acoustic-gravity waves in the thermosphere following Space Shuttle ascents// J. Atmos. Terr. Phys. V. 56. No. 4. P. 525-528. 1994.

100. Jacobson A.R., Carlos R.C., Blanc E. Observations of ionospheric disturbances following a 5-kt chemical explosion. 1. Persistent oscillation in the lower thermosphere after shock passage// Radio Sci. 1988. V.23. № 5. P. 820-830.

101. Klostermeryer J. Numerical calculation of gravity wave propagation in a realistic thermosphere. J. Atmos. Terr. Phys. 1972. 34. 765 774.

102. Kohl H., Ruster R., Schlegel K. Modern ionospheric science. Lindau. 1996.

103. Koshevaya S.V., Perez-Enriquez R., Kotsarenko N.Y. The detection of electromagnetic processes in the ionosphere caused by seismic activity// Geophisica Internacional. 1997.4. 1 10.

104. Kunitsyn V., Tereshchenko E. Ionospheric Tomography. Springer-Verlag, 2003, 272 p.

105. Kunitsyn V.E., Tereshchenko E.D., Andreeva E.S. et al. Radiotomographic investigations of ionospheric structures at auroral and middle latitudes// Ann. Geophysicae. 1995. V. 13. No. 12. P. 1351-1359.

106. MacKinnon R.F. Effects of winds on atmospheric pressure waves produced by hydrogen bombs// Journal of the meteorological Society of Japan. 1968. Ser. II. V. 46. P. 45 59.

107. Mayr H.G., Harris F.A., Herrero F.A., Spencer N.W., Varosi F., Pesnell W.D. Thermospheric gravity waves: observations and interpretation using the transfer function model (ТЕМ)// Space Sci. Rev. 1990. 54. 297 375.

108. Noble S.T. A large-amplitude traveling ionospheric disturbance exited by the Space Shuttle during launch// J. Geophys. Res. 1990. V. 95. 19037 19044.

109. Piani C., Durran D., Alexander M.J., Holton J.R. A numerical study of three-dimensional gravity waves triggered by deep tropical convection and their role in the dynamics of the QBO// J. Atmos. Sci. 2000. V. 57. No. 22. 3689 -3701.

110. Pitteway M., Hines C. The viscous damping of atmospheric gravity waves// Can. J. Phys. 1963. V. 41. 1935 1948.

111. Pitteway M.L.V., Rickel D.G., Wright J.W., Al-Jarrah M.M. Modeling the ionospheric disturbance caused by an explosion on the ground// Ann. Geophysicae. V. 3. No. 6. P. 695 704. 1985.

112. Row R.V. Acoustic-gravity waves in the upper atmosphere due to a nuclear detonation and an Earthquake// // J. Geophys. Res. 1967. V. 72. No. 5. 1599 -1610.

113. Row R.V. Evidence of long-period acoustic-gravity waves launched into the F region by the Alaskan Earthquake of March 28, 1964// J. Geophys. Res. 1966. V. 71.No. 1.343-345.

114. Rudenko G.V., Uralov A.M. Calculation of ionospheric effects due to acoustic radiation from an underground nuclear explosion// J. Atmos. Sol.-Terr. Phys. V. 57. No. 3. P. 225 236. 1995.

115. Sauli P., Boska J. Tropospheric events and possible related gravity wave activity effects on the ionosphere// J. Atmos. Sol.-Terr. Phys. 2001. 63. 945 -950.

116. Tolstoy I., Montes H., Rao G., Willis E. Long-period sound waves in the thermosphere from Apollo launches// J. Geophys. Res. 1970. V. 75. № 28. P. 5621 -5625.

117. Whiteway J.A., Duck T.J. Evidence for critical level filtering of atmospheric gravity waves. Geophys. Res. Lett. 1996. 23. 145 148.

118. Yongqi W., Kolumbau H. Comparisons of numerical methods with respect to convectively dominated problems// Int. J. Numer. Meth. Fluids. 2001. 37. 721 -745.

119. Zaslavski Y., Parrot M., Blanc E. Analysis of TEC measurements above active seismic regions// Phys. of the Earth and Planet. Inter. 1998. 105. 219 -228.

120. Zhang S.D., Yi F. A numerical study of nonlinear propagation of a gravity-wave packet in compressible atmosphere// J. Geophys. Res. 1999. V. 104. D12. 14261 14270.

121. Zhang S.D., Yi F. A numerical study of propagation characteristics of gravity wave packets propagating in a dissipative atmosphere// J. Geophys. Res. 2002. V. 107. D14. 1-9.