Использование дистанционного зондирования Земли для прогнозирования землетрясений

УДК 528.584: 550.343.6 И.В. Степанов1

ИСПОЛЬЗОВАНИЕ ДИСТАНЦИОННОГО ЗОНДИРОВАНИЯ ЗЕМЛИ ДЛЯ ПРОГНОЗИРОВАНИЯ ЗЕМЛЕТРЯСЕНИЙ

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

Ключевые слова: землетрясение, данные дистанционного зондирования, прогнозирование землетрясений, облачные структуры, предвестники землетрясений.

Введение. Под угрозой землетрясений находятся обширные территории, часто густонаселенные. Наибольшая опасность землетрясений заключается в их неожиданности и неотвратимости. Исследования, направленные на объяснение причин землетрясений, с переменным успехом ведутся уже много лет. Однако научные достижения последних лет в геофизике, сейсмологии и других смежных науках приблизили решение задачи прогнозирования землетрясений [1, 14].

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

Методы прогнозирования землетрясений. В геофизике (и других науках) известно множество различных по генезису и особенностям проявления предвестников землетрясений. Перечислим некоторые из них [3, 8, 11, 19, 20]: форшоки, изменение наклона земной поверхности, деформации горных пород, аномальные вариации магнитного поля Земли, изменение электросопротивления горных пород, изменение уровня воды в колодцах и скважинах, изменение температурного режима поверхностных земных слоев, изменение химического состава вод и газов, необычное поведение животных, облачные структуры определенного типа и др.

Известно, что прогноз землетрясений сводится к определению основных параметров: даты места (ф, X — широта и долгота эпицентра) и интенсивности (или магнитуды М) землетрясения [14].

Для решения проблемы прогноза землетрясений в Научном центре оперативного мониторинга Земли (НЦ ОМЗ) проводились исследования, в которых разрабатывали новые подходы к прогнозированию

землетрясений. Речь идет о землетрясениях, возникающих до глубины 700 км. Ранее [3, 4] специалисты НЦ ОМЗ на основе данных наземного и космического мониторинга выявили закономерности в структуре солнечно-земных связей, которые позволили определить классы признаков подготовки землетрясений. Установлено, что наличие определенных по строению и физическим свойствам облачных структур, формирующихся в высоких слоях атмосферы, — еще один признак землетрясений.

В классе облачных признаков установлены следующие факты: на космических снимках с геостационарных спутников, а также со спутников «Aqua», «Terra», NOAA и других фиксируется образование облачных структур, не имеющих объяснения в рамках классических представлений об облачности [2, 3, 11, 20]. В НЦ ОМЗ такие облачные структуры назвали облачными сейсмотектоническими индикаторами (ОСТИ). В научной литературе нет устоявшегося термина, обозначающего облачный предвестник землетрясения. В частности, можно встретить такие названия, как тектонические облака, сейсмические облака, облачные сейсмоиндуцированные структуры, линейные облачные аномалии, в зарубежной литературе — earthquake clouds, seismic clouds [1, 2, 8, 11, 20, 21]. Они образуются в результате дегазации Земли, прежде всего при выделении водорода. Гипотезу водородной дегазации выдвинул В.И. Вернадский. Позднее применительно к тектонике плит и сейсмике ее развили В.Н. Ларин, В.Л. Сывороткин [5, 9, 10, 11, 15, 20] и др. Помимо водорода в состав продуктов дегазации входят радон, гелий, аэрозоли, заряженные частицы и пр.

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

1 Московский государственный университет имени М.В. Ломоносова, географический факультет, кафедра картографии и геоинформатики, аспирант, e-mail: johnms@mail.ru

цо. Выявлен минимальный порог длины облачного предвестника — приблизительно 400 км. Эта длина определена на основании методики расчета магниту-ды (М) возможного сейсмического события; значение минимальной магнитуды принято равным 6. Согласно шкале сейсмической интенсивности М8К-64 и ее соотношения с магнитудой, при М > 6 следует ожидать сильное разрушение зданий, что может привести к многочисленным жертвам среди населения, нанести экологический и экономический ущерб.

Рассматриваемые облачные структуры могут перемещаться и/или поворачиваться относительно какого-либо разрывного нарушения под воздействием общей циркуляции атмосферы и/или под действием силы Лоренца. Они могут выглядеть как четкие линии облаков, как «пушистые», размытые облачные структуры, причем облака могут состоять из нескольких линейных образований (3-пальцевые, квадратные, Т-образные и др.). На космоснимках отмечается также отсутствие облачности на фоне сплошного облачного покрова, т.е. своего рода структура «провала», которая может иметь различную геометрическую форму (клинья, трапециевидные образования и др.). Кроме того, обнаружено, что циклоны при перемещении над активными тектоническими зонами изменяют траекторию [2, 4].

Время жизни этих облачных структур колеблется от нескольких часов до 2—3 дней. Как правило, облака большей протяженности (а значит, с большей потенциальной магнитудой) живут около 10—15 часов, высота их образования составляет около 5—10 км. Определение высоты облачной структуры проводится с помощью космических снимков в тепловом диапазоне.

Установлена также [3] логарифмическая связь между облачной структурой максимальной протяженности и потенциальной магнитудой:

М = ln D,

где М — магнитуда, D — максимальная протяженность облачной структуры над активной сейсмотектонической областью, км.

Автором составлена база данных космических снимков с облачными предвестниками землетрясений с 2003 г., полученных со спутников «Terra», «Aqua», «Meteosat» и других. Наиболее часто проявление этих процессов отмечено по границам Тихоокеанского сейсмического пояса (Индонезия, Филиппины, Камчатка, Япония и др.), в рифтовых зонах (Африка, Атлантический хребет), в Тибетско-Гималайском регионе. Регулярный мониторинг регионов Земли, где возможны образования ОСТИ, позволяет своевременно обнаружить эти признаки и рассчитать зоны, где вероятно возникновение землетрясений.

Как сказано выше, прогноз землетрясения должен основываться на вычислении трех параметров — даты, места и магнитуды возможного землетрясения. Но прежде чем описать методику вычисления этих

Рис. 1. Процесс возникновения облачных сейсмотектонических индикаторов

параметров, которая используется в НЦ ОМЗ, необходимо вкратце упомянуть о механизме образования облачных сейсмоиндикаторов (рис. 1).

В рассматриваемой методике прогноза землетрясений учитывается несколько параметров. Во-первых, это аномалии параметров вращения Земли (ПВЗ), которые вызваны множеством факторов, прежде всего гравитационным влиянием Солнца и Луны (и других планет). Эти аномалии приводят к деформации геоида, которая в свою очередь влияет на процессы, происходящие на тектонических границах различного ранга [16, 17]. В зонах активных тектонических границ происходит выброс флюидов и сопутствующих газов, которые способствуют образованию ОСТИ. Такой подход к прогнозу предполагает отслеживать аномалии, происходящие в литосфере и атмосфере Земли, а также изменение параметров солнечной активности по данным наземно-космических наблюдений.

Определение даты прогнозируемого землетрясения основано на связи солнечной, геомагнитной и сейсмической активности. Установлено, что сейсмический процесс имеет магнитно-меридиональную направленность (так называемый эффект D-триггера), наблюдаемую в виде цепочек землетрясений, происходящих в среднем на 14-е или/и 21-е сутки (±2 суток) после геоэффективных явлений на Солнце [3, 4, 15]. Подобные явления на Солнце оцениваются по максимальным вариациям магнитного поля Земли в тектонических зонах, накопивших достаточный для сброса сейсмической энергии потенциал напряжений и деформаций (ПНД).

Критические значения ПНД формируются в результате накопления нестабильностей (аномалий) во вращении и обращении Земли. Это приводит к деформациям геоида, аномальным процессам на границе ядро—мантия, выносу флюидов и сопутствующих газов (Н2, Ял, Не) и в итоге к накоплению напряжений на границах разломов тектонических плит и блоков. Накопленный ПНД проявляется в виде аномалий геофизических полей в лито- и атмосфере. В частности, в тропосфере такие аномалии могут быть выражены в виде специфических облачных структур, в ионосфере — в виде концентрации электронов и сопутствующих электромагнитных аномалий.

Место возможного землетрясения определяется на основе анализа положения облачных сейсмотектонических индикаторов относительно каких-либо

разрывных нарушений, а также активных на данный момент магнитных меридианов. Последние определяются на основе данных астрономических наблюдений солнечной активности. Геоэффективность солнечных явлений определяется по значениям так называемого ^-индекса, который характеризует изменение геомагнитной активности, фиксируемое на конкретной обсерватории в 3-часовом интервале времени, начиная с 0 час. 00 мин. UT [5].

Магнитуда прогнозируемого землетрясения определяется по формуле (1), основанной на связи между максимальной протяженностью ОСТИ и потенциальной магнитудой возможного землетрясения. Методика выявления сейсмопризнаков и прогноза землетрясений опирается на созданную в НЦ ОМЗ автором статьи подсистему ГИС. В ней анализируется сейсмомагнитная обстановка, учитывается результат анализа облачных предвестников землетрясений (преимущественно по космическим снимкам со спутников «Meteosat», «MTSAT-1R», «Terra», «Aqua» и др.), а также положение магнитных меридианов и других признаков в сейсмически активных районах Земли.

В качестве средства реализации ГИС [8] используется ГИС-пакет ArcGIS 9.1, в котором проводится сбор данных, а также анализ и оценка сейсмической ситуации, мониторинг и прогнозирование землетрясений. В рамках решения задачи используются следующие пространственные данные (слои ГИС):

— цифровые карты: картографическая основа в проекции Меркатора, в качестве которой используется: а) карта мира, поставляемая с пакетом ArcGIS, на нее наносят все данные, применяемые для создания подсистемы ГИС как составной части методики прогнозирования землетрясений; б) карты границ литосферных плит различного ранга, а также карты разломов. Источником данных служат векторная карта-модель литосферных плит (Plate Boundary Model PB2002) П. Бирда, векторизованная схема глубинных разломов территории бывшего СССР (карта 1977 г.), векторизованные карты других зон Земли;

— регулярно пополняемая и поддерживаемая база данных (БД) землетрясений. В ней содержатся сведения о произошедших землетрясениях (источник: сайт Геологической службы США (USGS) и Европейского сейсмологического центра (EMSC));

— цифровые космические снимки с различных спутников («Terra», «Aqua», «Meteosat», «MTSAT-1R» и др.), по которым исследуются облачные структуры на предмет обнаружения признаков землетрясений методом визуального компьютерного дешифрирования;

— база данных сейсмомагнитных меридианов, которая используется для расчета потенциальных зон риска землетрясений. В НЦ ОМЗ регулярно создаются карты магнитных меридианов на конкретные даты. По этим данным выявляются предварительные зоны возможных землетрясений.

Рассмотрим подробнее содержание и применение отображаемых отдельных слоев ГИС.

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

Слой «Магнитные меридианы» отображается в виде линейных знаков. Все меридианы имеют 2 атрибута в базе данных — долготу и направление. Долгота определяется по точке пересечения с экватором. Меридианы имеют наклонение ±79° относительно плоскости экватора, обусловленное полем геомагнитного диполя, и образуют «кольца» вокруг Земли, графически отображая магнитные силовые трубки, в пределах которых возможны землетрясения.

Слои «Границы плит» и «Разломы» также содержат линейные объекты. В подсистеме ГИС имеются границы литосферных плит различного ранга и типы границ плит на Земле в целом, основанные на указанной модели П. Бирда, включающей 52 плиты и 7 типов границ плит, отдельно выделяются зоны орогенеза и субдукции [18].

Рассматривались также модели отечественных исследователей, в частности схема «Современные литосферные плиты и их кинематика», приведенная в монографиях В.Е. Хаина, где выделено 7 плит 1-го порядка (мегаплиты), 11 плит 2-го порядка (мезопли-ты) и плиты 3-го порядка (микроплиты), расположенные в крайней северо-восточной части Тихого океана, между Индо-Австралийской плитой и Тихоокеанскими плитами на юго-востоке Тихого океана; а также 3 типа границ (дивергентные—спрединговые, конвергентные—субдукционные и трансформные) [16, 17].

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

Решение задач мониторинга землетрясений на основе облачных предвестников представляется, по мнению автора, более продуктивным при использовании модели П. Бирда. В пользу этого можно привести ряд аргументов. Во-первых, за счет высокой подробности содержания и относительной сложности рисунка границ разломов даже визуальный анализ позволяет отметить приуроченность большинства ОСТИ именно к границам, отмеченным на карте П. Бирда. Это обстоятельство позволяет объяснить, в частности, возникновение некоторых ОСТИ в районе Корейского п-ова, так как именно там (по модели П. Бирда) проходит граница между Амурской и Евразийской плитами. На карте В.Е. Хаина эти плиты объединены в одну Китайскую плиту. Подобный пример можно привести для Камчатки и некоторых других регионов.

Рис. 2. Примеры космических снимков с ОСТИ перед землетрясением (03.10.2009) в районе о-ва Тайвань с М ~ 6,0: А и Б — облачные сейсмотектонические индикаторы (выделены овалом), возникшие: А — 21.09.2009, Б — 28.09.2009; 1 и 2 — землетрясения, произошедшие: 1 — 28.09.2009 с М = 5,8; 2 — 3.10.2009 с М = 6,0

Рис. 3. Карта ОСТИ перед землетрясением в районе о-ва Тайвань 03.10.2009 с М ~ 6,0: 1 — землетрясения, произошедшие: 1 — 28.09.2009 с М = 5,8; 2 — 03.10.2009 с М = 6,0; 2 — зона ОСТИ, возникшие А — 21.09.2009, Б — 28.09.2009; 3—6 — тектонические границы плит (3 — конвергентные, 4 — дивергентные, 5 — трансформации, 6 — зоны субдукции); плиты: Еи — Евразийская, Р8 — Филиппинская, АМ — Амурская, ОК — Окинавская, YA — Южно-Китайская

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

Космические снимки с геостационарных, а также с орбитальных спутников поступают из баз данных специализированных серверов в интернете для территорий с повышенным риском землетрясений. При отсутствии привязки космоснимков операция трансформирования производится в полуавтоматическом режиме в программном пакете ERDAS Imagine S.4 либо в ENVI 3.4. Затем работа со снимками продолжается в Arc GIS.

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

С помощью подсистемы ГИС успешно определяются дата, место и магнитуда возможного землетрясения. Так, магнитуда определяется по формуле (1) на основе измерения максимальной протяженности облачной структуры (D) на космических снимках. Из тех предвестников землетрясений, которые используются в НЦ OMЗ для прогноза землетрясений, наиболее точные по дате — данные о солнечной активности. Mесто возможного сейсмического события определяется на основании проекции магнитных меридианов (на определенную дату) и положения предвестников землетрясений (прежде всего облачных) на карты плит и разломов, а также в результате анализа различных классов признаков подготовки землетрясений (протонные признаки, вариации полюса

вращения Земли, слабые предшествующие землетрясения и др.).

К продуктам ГИС для решения задач мониторинга и прогнозирования землетрясений относятся серии карт-схем прогноза землетрясений на определенные даты. Эти карты-схемы, как правило, включают, помимо основной карты зон риска возможного землетрясения, врезки — различные графики (например, график траектории полюса вращения Земли), космические снимки облачных структур-признаков, карты границ плит и линии разломов; линий магнитных меридианов и др. Такие карты-схемы строят в различном масштабе, который зависит от размеров района исследования (как правило, от 1:10 000 000 для отдельных регионов Земли до 1:150 000 000 для карт мира) в проекциях иТМ, Меркатора (для карты мира) либо в прямоугольной географической системе координат.

В качестве примера рассмотрим фрагменты космических снимков с двумя сменяющимися ОСТИ (рис. 2), возникшими 21 и 28 сентября 2009 г., т.е. накануне землетрясения, и карту ОСТИ перед землетрясением в районе о-ва Тайвань (3 октября 2009 г.) с М = 6,0 (рис. 3). На рис. 2, А длина ОСТИ составила около 1100 км, что по формуле (1) дает М ~ 7,0. По рис. 2, Б длина ОСТИ равна приблизительно 590 км, а М ~ 6,3. Этот пример наглядно иллюстрирует динамику Окинавской микроплиты и последовательно возникших ОСТИ и землетрясений. Несоответствие расчетной магнитуды и произошедших землетрясений специалисты НЦ ОМЗ объясняют возникшим в этом же районе мощным тайфуном (существует гипотеза об «отбирании» сейсмической энергии активизированными вулканами и тайфунами/циклонами).

Аналогичный подход был успешно реализован в ряде прогнозов землетрясений, в том числе катастрофических в Суматрано-Яванской зоне в 2004—2008 гг. (например, землетрясение на Суматре 20 февраля 2008 г., прогноз которого был дан на сайте НЦ ОМЗ (http://www.ntsomz.ru/projects/earthquake/ doda_news_201207) и в Сейсмосовете России.

Заключение. К настоящему времени автором накоплен определенный опыт прогнозирования землетрясений с использованием облачных признаков за период с 2004 г. по настоящее время в зоне Тихоокеанского сейсмического пояса (Индонезия, Филиппины, восточная часть РФ — Камчатка, Сахалин, а также Япония и др.), в рифтовых зонах (Африка, Антлантический хребет), в Тибетско-Гималайском регионе и др. Указанное позволяет сделать следующие выводы:

— мониторинг облачных структур по космическим снимкам с целью выявления облачных сейс-мопризнаков наиболее целесообразно проводить над тектоническими границами (над границами плит разного ранга и разломами);

— наиболее эффективен мониторинг с использованием сканирующей аппаратуры, установленной

на геостационарных спутниках в видимом и ИК-диапазонах и имеющей пространственное разрешение 4—8 км. Уточнение и детализация сейсмопризнаков проводится с использованием аналогичной аппаратуры на спутниках «Terra» и «Aqua» с разрешением 1—2 км;

— оперативную обработку полученных космо-снимков и доведение их до тематического продукта целесообразно осуществлять с помощью ГИС-техно-логий и методик, описанных в статье;

— для наиболее сейсмоопасных регионов мира созданы банки данных типовых облачных сейсмотектонических индикаторов с выделением характерных для конкретного региона. Примеры типичных ОСТИ для Камчатки приведены на рис. 4 и 5;

— приведенные в статье примеры анализа признаков подготовки землетрясений показывают состоятельность применяемых методик прогнозирова-

Рис. 4. Космический снимок с характерным для Камчатки ОСТИ (выделен овалом). Облачный сейсмотектонический индикатор возник 10.11.2008 перед землетрясением в Охотском море (обозначено кружком), которое произошло недалеко от западного побережья Камчатки 24.11.2008 с М = 7,1. Протяженность ОСТИ составляла 1200 км, по ней была рассчитана магнитуда, равная 7,1

Рис. 5. Космический снимок ОСТИ (выделен овалом). Облачный сейсмотектонический индикатор возник 25.03.2009 перед землетрясением на одном из Курильских островов, произошедшим 07.04.2009 с М = 6,9. Протяженность ОСТИ составляла 1050 км, по ней была рассчитана магнитуда, равная 6,9

ния на основе использования геоинформационных технологий.

Автор считает необходимым создание постоянного космического мониторинга границ тектонических плит, прежде всего Индостанской, Индонезийской, Охотской и др. (в перспективе для всей Земли). Это связано с тем, что на границах плит происходит большинство землетрясений и расположены сейсмоо-пасные зоны с повышенной плотностью населения, промышленными объектами. В таких зонах находятся страны Южной, Юго-Восточной и Восточной Азии,

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

1. Геофизика на рубеже веков. М.: ОИФЗ РАН, 1999. 330 с.

2. Дода Л.Н., Натяганов В.Л., Тихонкова Ю.А. Облачные сейсмоиндуцированные структуры как следствие протонного электротеплового пробоя литосферы // Современные проблемы аэрогидродинамики. М.: Изд-во Моск. ун-та, 2007. С. 43-44.

3. Дода Л.Н., Новикова Н.Н., Пахомов Л.А.. и др. Прогнозирование и космический мониторинг предвестников землетрясений: проблемы, надежды, реалии // Проблемы прогнозирования чрезвычайных ситуаций: Тез. докл. V науч.-практ. конф. М.: МТП-инвест, 2006. С. 114-120.

4. Дода Л.Н., Новикова Н.Н., Пахомов Л.А., Степанов И.В. Космический мониторинг предвестников землетрясений // Наука в России. 2009. № 6. С. 37-43.

5. Жамалетдинов А.А., Митрофанов Ф.П., Токарев А.Д., Шевцов А.Н. Влияние лунно-солнечных приливных деформаций на электропроводимость и флюидный режим земной коры // Докл. РАН. 2000. Т. 371, № 2. С. 235-239.

6. Заболотная Н.А. Индексы геомагнитной активности: Справ. пособие. 2-е изд. М.: ЛКИ, 2007. 88 с.

7. Икея М. Землетрясения и животные. От народных примет к науке. М.: Научный мир, 2008. 320 с.

8. Капралов Е.Г., Кошкарев А.В., Тикунов B.C. и др. Основы геоинформатики. Кн. 1. М.: Академия, 2004. 480 с.

9. Ларин В.Н. Наша Земля (происхождение, состав, строение и развитие изначально гидридной Земли). М.: Агар, 2005. 248 с.

10. Ларин В.Н., Ларин Н.В. Водородная дегазация на Русской платформе, ее плюсы и минусы. URL: http:// hydrogen-future.com/rus/page-id-18.html. 2008 (01.02.2009).

I.V. Stepanov

Мексика, западная окраина США и др. Для решения задач оперативного мониторинга и прогнозирования землетрясений необходимы данные дистанционного зондирования Земли с периодичностью 30 мин. При этом съемка должна осуществляться в видимом и инфракрасном диапазонах с разрешением в 1—2 км. Наличие таких данных позволило бы своевременно выявить появление ОСТИ и тем самым обнаружить активные на данный момент времени зоны Земли, что позволит повысить точность прогнозирования землетрясений как в пространстве, так и во времени.

11. Летников Ф.А. Флюидные фации континентальной литосферы // Смирновский сборник — 99. М.: Изд-во Моск. ун-та, 1999. С. 63-98.

12. Морозова Л.И. Спутниковый мониторинг землетрясений. Владивосток: Дальнаука, 2005. 136 с.

13. Официальный сайт Института физики Земли им. О.Ю. Шмидта РАН. URL: http://seismos-u.ifz.ru (01.02.2009).

14. Соболев Г.А., Пономарев А.В. Физика землетрясений и предвестники. М.: Наука, 2003. 270 с.

15. Сывороткин В.Л. Глубинная дегазация Земли и глобальные катастрофы. М.: Геоинформцентр, 2002. 250 с.

16. Хаин В.Е., Тектоника континентов и океанов. М.: Научный мир, 2001. 606 с.

17. Хаин В.Е., Ломизе М.Г. Геотектоника с основами геодинамики. М.: КДУ Университет, 2005. 560 с.

18. Bird P. An updated digital model of plate boundaries // J. Geochemistry, Geophysics, Geosystems. An Electronic Journal of the Earth Society. 2003. 14 March. 52 р.

19. Gere J.M., Shah H.C. Terra Non Firma. Understanding and preraring for earthquakes. N.Y.: W.H. Freeman and Company, 1984. 203 р.

20. Liperovsky V.A., Meister C.-V, Doda L.N. et al. On the possible influence of radon and aerosol injection on the atmosphere and ionosphere Earth System Sciences // EGU. 2005. Vol. А-10347.

21. Shou Z. Earthquake clouds a reliable precursor // Science and Utopya. 1999. P. 53-57.

Поступила в редакцию 06.03.2009

APPLICATION OF REMOTE SENSING OF THE EARTH FOR EARTHQUAKE FORECASTING

The article deals with earthquake forecasting on the basis of space monitoring data. A method of earthquake forecasting has been developed which uses cloud seismic indicators easily revealed on space images as the earthquake precursors. A system of space monitoring of the earthquake zones is suggested.

Key words: earthquake, remote sensing data, earthquake forecasting, cloud structures, earthquake precursors.