сторов вместо развития местного бизнеса.
Таким образом, стратегия устойчивого развития является инструментом для укрепления локальной территориальной идентичности, при этом для жителей побережья озера Тахо - это возврат к утраченным
традициям, для жителей Байкальского региона - их актуализация, а для жителей Тофаларии - часть повседневности.
Статья поступила14.03.2014 г.
Библиографический список
1. Батоцыренов Э., Куклина В.В. «Байкальская гавань» -особая экономическая зона туристско-рекреационного типа и ее восприятие местными жителями // Известия РГО. 2010 Т. 142. Вып. 1. С. 75-81.
2. Город после комбината: социально-экономические стратегии жителей города Байкальска / И. Корюхина, Т. Тимофеева, Т. Гребенщикова, И. Абдулова, В. Куклина, М. Рожан-ский; под ред. М. Рожанского. Иркутск: Изд-во Центра независимых социальных исследований и образования, 2012. 144 с.
3. Каганский В.Л. Россия. Провинция. Ландшафт // Отечественные записки. 2006. № 6. С. 244-257.
4. Куклина В.В. Социальное пространство постсоветского моногорода (на примере Байкальска и Усть-Илимска) // Оценка современных факторов развития городов и урбани-зационных изменений в Сибири: кол. монография; отв. ред. Л.М. Корытный, Н.В. Воробьёв. Новосибирск: Академическое изд-во «Гео», 2011. С. 128-139.
5. Мантатов В.В., Потапов Л.В. Стратегия устойчивого развития Республики Бурятия // Информационный портал «Бай-кал-1_аке» [Электронный ресурс]. 11131-: ||Нр:/Л№мм.Ьа1ка1-сетег.ги/Ьоок8/е1етет.р||р?Ю=1703 (04 янв. 2013).
6. Программа ООН по населённым пунктам (ООН-Хабитат),
URL:
2012 [Электронный ресурс].
http://www.unhabitat.ru/ru/index (30 янв. 2013).
7. Родоман Б. Экологическая специализация России // Интелрос (Интеллектуальная Россия) [Электронный ресурс]. URL: http://www.intelros.ru/subject/figures/boris-rodoman/12628-ekologicheskaya-specializaciya-rossii.html (14 фев. 2014).
8. Antipov A.N., Kravchenko V.V., Semenov Yu.M. and others. Landscape Planning: Tools and Experience in Implementation. Bonn - Irkutsk: V.B. Sochava Institute of Geography SB RAS Publishers, 2006. 149 p.
9. Medvedev D.A. Statement of H.E. Mr. Dmitry Medvedev, PM of Russia, at Rio+20. October 2, 2012 // Интелрос (Интеллектуальная Россия) [Электронный ресурс]. URL: http://sustainabledevelopment.un.org/index.php?page=view&nr= 3136&type=255&menu=35 (04 янв. 2013).
10. O'Riordan T.R., Voisey H. The Transition to Sustainability: The Politics of Agenda 21 in Europe. London: Earthscan, 1998.
11. Sustainable Development for the Next Twenty Years United States Views on Rio+20. November 1, 2011 // Sustainable development: knowledge platform [Электронный ресурс]. URL: http://www.state.gov/e/oes/sus/releases/176863.htm (07 марта 2014).
УДК 550.343
ОПЫТ ПРОГНОЗА ЗЕМЛЕТРЯСЕНИЙ В СЕЙСМИЧЕСКОЙ ЗОНЕ РИФТА ШАНЬСИ (КИТАЙ) И ЕГО ОПЕРЯЮЩИХ РАЗЛОМАХ НА ОСНОВЕ ВОЛНОВОЙ КОНЦЕПЦИИ С ИСПОЛЬЗОВАНИЕМ ОГРАНИЧЕННОЙ БАЗЫ СЕЙСМИЧЕСКИХ СОБЫТИЙ
© М.Г. Мельников1
Институт земной коры СО РАН,
664033, Россия, г. Иркутск, ул. Лермонтова, 128.
На основе ограниченной базы сейсмических событий и карты активной тектоники Китая выполнен прогноз сейсмических событий для рифтовой системы Шаньси. В работе использован метод и концепция деформационных волн С.И. Шермана. Построены графики распространения деформационных волн вдоль оси сейсмической зоны рифта, а также его оперяющих активных разломов. Прогноз сделан в двух масштабах соответственно. Также построена карта прогноза с доверительными границами областей и доверительными интервалами по времени для прогнозируемых событий. Ил. 7. Табл. 3. Библиогр. 25 назв.
Ключевые слова: деформационные волны; сейсмическая зона; активный разлом; землетрясение; рифт.
EARTHQUAKE PREDICTION EXPERIENCE IN SHANXI RIFT SEISMIC ZONE (CHINA) AND ITS FEATHERING FAULTS BASED ON WAVE CONCEPTION USING LIMITED DATABASE OF SEISMIC EVENTS M.G. Melnikov
Institute of the Earth's Crust SB RAS, 128 Lermontov St., Irkutsk, 664033, Russia.
Based on the limited database of seismic events and the map of active tectonics in China the forecast of seismic events for Shanxi rift system is made. S.I. Sherman's method and conception of deformational waves are used in this work. The plots of deformational wave propagation along the axis of the rift seismic zone and its feathering active faults are constructed. The forecast is made in two scales respectively. The map of the forecast is created with the confidence contours of areas and confidence time intervals for predicted events. 7 figures. 3 tables. 25 sources.
1Мельников Михаил Геннадьевич, аспирант, тел.: 89501101058, e-mail: [email protected] Melnikov Mikhail, Postgraduate, tel.: 89501101058, e-mail: [email protected]
Key words: deformational waves; seismic zone; active fault; earthquake; rift.
Введение
Прогноз сейсмических событий сегодня является одной из самых актуальных задач сейсмологии. Сейсмическое событие - явление сложное, и его возникновение (место, время и энергия) зависит от множества факторов. Известны различные методы прогноза землетрясений [14, 22]. Чаще всего, методы и методики связаны с авторскими концепциями. Например, концепция сейсмического процесса С.И. Шермана, основанная на понятии о деформационных волнах [6], сейсмической зоне [10, 23], а также разработанная методика и специальный программный пакет [15, 16]. Природа деформационных волн еще до конца не изучена, эта теория редко применяется на практике, хотя существуют работы, объясняющие их возможное происхождение [18-21]. Теория деформационных волн занимает сегодня важное место среди теоретически и социально значимых исследований, ей посвящено большое количество работ российских и зарубежных авторов [4; 5; 13], но этот вопрос до сих пор остается малоизученным.
На основе работ С.И. Шермана по сейсмическому процессу Байкальского региона, автор данной статьи сделал прогноз сейсмических событий для рифтовой
системы Шаньси по аналогии с выполненными ранее работами для Байкала [6; 8; 9; 16; 23].
Объект исследования
Объектом настоящего исследования была выбрана рифтовая система Шаньси (Китайская Народная Республика). Она расположена вдоль осевой зоны куполообразного поднятия нагорья Шаньси, протяженность которого составляет 1200 км (рис. 1). На севере и на юге она ограничена горами Иньшань-Яньшань и Цинлинь. Западным плечом рифта является поднятие Лулиангшань, которое лежит на востоке от Лессового плато и Ордоса и простирается на северо-восток. Восточное плечо рифта - поднятие Тай-ханшань, лежащее к западу от Северокитайской равнины [12]. Впадина, которую образует рифт, представляет собой долину Фенвей протяженностью 600 км и шириной от 30 до 90 км. Как и в Байкальской рифтовой системе, в рифте Шаньси присутствует правая сдвиговая составляющая [17]. Рельеф долины Фенвей является выражением ее разломно-блоковой структуры. В ее основе лежат многочисленные разломы, которые почти все имеют углы падения 53-85° [11].
Рис. 1. Карта активной тектоники Китая (фрагмент) [1]. Основные условные обозначения: 1 - коренные выходы пород; 2 - события М > 8; 3 - события 7 > М > 7,9; 4 - события 7 > М > 7.9; 5 - сейсмоактивные разломы;
6 - пассивные разломы
Рифт Шаньси является активной внутриплитной тектонической зоной. Имеющиеся данные о сейсмичности показывают, что в Ордосском массиве и поднятии Тайханшань (на западе и востоке рифта) крупных землетрясений не было. В противоположность этому между двумя стабильными блоками рифт Шаньси является интенсивной сейсмической зоной, где с 231 г. н. э. произошло 32 события с магнитудой М > 6 [3].
По многим своим параметрам рифтовая система Шаньси схожа с Байкальской, что неоднократно обсуждалось многими исследователями [11; 17]. В частности, подобие их структурного положения состоит в
том, что своей западной частью они примыкают к жестким массивам с утолщенной корой, а восточной -к подвижным областям с более тонкой корой, и образуют Б-образные на плане формы, созданные сочетанием разломов и примыкающими к ним кайнозойскими впадинами [9]. У обеих систем схожи параметры сейсмичности, региональное поле напряжений, а также ряд других геологических и геофизических параметров (рис. 2). Поэтому изучение параметров и процессов в рифте Шаньси представляет большой интерес для выявления общих тенденций и закономерностей для рифтовых систем в целом.
Рис. 2. Векторы главных осей напряжений в рифтовых системах: а - Шаньси [17]: 1 - векторы главных осей напряжений; 2 - тип сдвиговых смещений по сегментам рифта; 3 - активные разломы; б - Байкальской [7; 26]: 1 - векторы главных осей напряжений; 2 - тип сдвиговых смещений по сегментам рифта; 3 - оси сегментов
рифтовой системы
Фактический материал и методика его обработки
Деформационные волны возбуждают и контролируют короткопериодную активизацию разрывов в реальном времени и пространственно-временную локализацию очагов в деструктивных зонах литосферы. Активизация происходит при прохождении фронта волны через разломы сейсмической зоны, который провоцирует подвижки в некоторых из них. В зависимости от текущего состояния конкретного разлома, а также энергетических параметров волны в нем происходит (или не происходит) сейсмическое событие. По методике, описанной в работах [8; 23], где подобные исследования проводились для Байкальской рифто-вой системы и других систем Центральной Азии, и с помощью программного пакета Digital Faults [15], в котором разработаны приемы для выявления направлений возбуждения очагов землетрясений в областях динамического влияния разломов и оценки параметров деформационных волн [6; 8; 16], была изучена специфика пространственно-временных локализаций очагов землетрясений. В качестве фактического материала была использована карта активной тектоники Китая [1], с которой были взяты активные разломы, и каталог сейсмических событий (бюллетень Международного Сейсмологического Центра (ISC) [2]) (рис. 3).
Разломы различных рангов контролируют сейсмичность на разных уровнях [6; 23], в связи с этим исследования проводились двумя способами: по сегментам сейсмической зоны рифта и, в более крупном масштабе, по его оперяющим активным разломам. В каждом случае понятие области активного динамического влияния разлома использовалось в соответствии с масштабом исследования [24]. В случае с зо-
ной она составляла 70 км в обе стороны от оси, что соответствует протяженности и ширине сейсмической зоны рифта. Для оперяющих разломов - 20 км в обе стороны от разлома, что также соответствует масштабу и ранее проведенным исследованиям [25].
В рифтотвой системе Шаньси была выделена ось сейсмической зоны, которая состоит из трех сегментов: I - центральный, II - юго-западный и III - северовосточный (см. рис. 3). Такое же разделение характерно и для Байкальской рифтовой системы. Сегменты образуют собой структуру S-образной формы, фланги которой имеют субширотное простирание, а центральная часть - северо-восточное.
Для оси сейсмической зоны рифта Шаньси был построен общий график в координатах «время-место», на котором все три сегмента представлены как единая структура (рис. 4). В данном случае был использован диапазон магнитуд событий M > 4,8.
На графике прослеживаются четыре активизации за период с 1955 по 2010 годы, направленные от северо-восточного фланга рифта к юго-западному. На северо-восточном конце пунктиром показана возможная новая активизация, для которой пока нет фактических данных. Тонкие пунктирные линии на графиках показывают доверительные интервалы для времени возникновения событий с вероятностью 90%, которые рассчитывались по формуле
А = ±/ ,
п
где I - критерий Стьюдента для вероятности 0,9 и количества элементов в выборке п; а - среднеквадрати-ческое отклонение по времени от линии регрессии.
Рис. 3. Сегменты и активные разломы рифта Шаньси: 1 - сегменты сейсмической зоны (I - центральный, II - юго-западный, III - северо-восточный); 2 - активные разломы [1]
2020
2010
2000
1990
о. 1980 СО
1970
1960
у = -0,1121 х + 2680,9 *» ч ■.у = -0,1565х +2325,5
^ 1^ = 0,5414 —\ • —ч ч Ч®-^ \ '■•. Р2 = 0,4589 4 . \ ■•■ • ч % 1
V X. о ^ о.
О \ О , » О-.
V \ • П X
V и \ \ СХ .
\ X®
О
у = -0,1639х + 2009,^>>«ч у = -0,1137х + 2б3"|34/ ''I
К* = 0,9433 —Р2 = 0,759
, д , , о —Д-,-X-А_,
100
200
300
400
500
600
700
800 900
1000
Расстояние, км
Рис. 4. График в координатах «время-место» для оси сейсмической зоны рифта Шаньси как единой структуры
Параметры выявленных волн, уравнения регрессии с доверительными интервалами и коэффициенты детерминации графиков представлены в табл. 1 и 2. Полученный график позволяет рассчитать скорость и направление деформационных волн для сейсмической зоны рифта. Полученное значение скорости - 7,2 км/год (тангенс угла наклона линии регрессии), направление распространения - от северо-восточного фланга к юго-западному. Это согласуется с полученными ранее результатами для Байкальской рифтовой системы [8; 23].
Вторым способом для прогноза землетрясений был рассмотрен каждый из оперяющих разломов рифта. Для анализа было выбрано девять из них с достаточными длинами и количеством данных о сейсмичности для обеспечения статистики (рис. 5). В данном случае были также построены графики в координатах «время-место». Диапазон магнитуд событий, участвовавших в данной части анализа, М > 3,1. Направления движения волн в разломах согласуются с полученными данными для оси сейсмической зоны (с северо-востока на юго-запад).
Результаты и обсуждение
В континентальной литосфере происходят волновые процессы на трех глубинных уровнях, охватывающих всю литосферу, ее верхнюю хрупкую часть и собственно земную кору [6]. По результатам проделанной работы можно сделать заключение, что волны, полученные двумя разными способами, имеют разные параметры и относятся к категории распространяю-
щихся в верхней хрупкой части литосферы, то есть волн второго и третьего уровня (об этом свидетельствуют глубины гипоцентров) [6]. Волновые процессы третьего глубинного уровня генерируются множеством причин, в том числе наведенной сейсмичностью. Они генерируют слабые землетрясения, появление которых сложно предсказать из-за их большого количества. Для второго способа (для оперяющих разломов) в данном исследовании были использованы события диапазона магнитуд не с минимального значения, которое в используемом каталоге порядка 1, а со значения 3,1. Это можно объяснить крайней затруднительностью или невозможностью интерпретации данных по землетрясениям низкой магнитуды (ниже полученного эмпирическим путем значения 3,1), которые составляют большую часть каталога.
На основе полученных данных о волновом процессе для каждой активизации был рассчитан прогноз событий на ближайшие годы. Продлив линии регрессии графиков до 2020 года, можно отложить на них предполагаемые будущие события с интервалом, равным среднему периоду между землетрясениями (см. табл. 2). Если учесть доверительные интервалы, то на оси зоны можно выделить опасные участки с наиболее вероятным возникновением землетрясений (на рис. 6,а они выделены черными линиями). Для графиков оперяющих разрывов рифта были построены такие же графики, их примеры приведены на рис. 6,6.
Таблица 1
Параметры деформационных волн для оси сейсмической зоны рифта Шаньси
Длина оси, км Прости -рание оси, град. Максимальная глубина гипоцентров, км Количество событий, использованных в анализе; диапазон магнитуд (классов) Длина волны, км Скорость волны, км/год и направление вектора скорости Средний период волны, лет
980 35 28 34; М=4,8+6,1 (К=13+15) 230,5 7,2 СВ-ЮЗ 32
Таблица 2
Коэффициенты детерминации и уравнения регрессии графика для оси сейсмической зоны _рифта Шаньси_
Длина оси, км Количество активизаций Коэффициент детерминации Уравнение регрессии Средний период между событиями в каждой активизации, лет
980 4 R12 = 0,94 t = (-0,1639-1 + 2009,3) ± 3,3 7,1
R22 = 0,76 t = (-0,1137-1 + 2031,3) ± 5,5 5,9
Rз2 = 0,54 t = (-0,1321-1 + 2080,9) ± 7,6 5,5
R42 = 0,46 t = (-0,1565-1 + 2125,5) ± 5,5 4,6
Рис. 5. Активные разломы рифта Шаньси [1]: 1 - разломы, участвующие в анализе сейсмического процесса
и их номера; 2 - другие разломы
Рис. 6. Графики в координатах «время-место» с прогнозами для сейсмической зоны рифта (а)
и для оперяющих разрывов (б)
Каждый участок, выделенный на разломе, имеет 90%-ю вероятность возникновения на нем события с магнитудой выше значения, соответствующего масштабу исследования (М > 4,8 и М > 3,1 для оси сейсмической зоны и для оперяющих разрывов соответственно) в определенный промежуток времени (соответствующий графику).
По полученным графикам была построена карта прогноза для рифта Шаньси (рис. 7), на которой серые участки означают опасные области с наиболее вероятными землетрясениями, ширина которых соответствует областям динамического влияния, а длина -полученным из рис. 6 участкам наиболее вероятных событий. Светлые области - для событий с магнитудой М > 4,8, темные - М > 3,1. В табл. 3 показаны
опасные временные периоды для темно-серых областей, соответствующие номеру на карте.
Для данного прогноза использовался каталог событий, включающий землетрясения, произошедшие до конца 2011 года. В ряде случаев полученные доверительные интервалы по времени захватывают уже прошедшие годы, в том числе и последние годы в каталоге (это видно на графиках рис. 6). В данном исследовании ставилась задача апробирования методики и получения опыта прогноза, поэтому автор считает возможным оставить доверительные интервалы по времени без изменений, так как по формальным признакам они соответствуют логике и не противоречат задачам исследования.
Рис. 7. Карта прогноза сейсмических событий. 1 - разломы [1]; 2 - ось сейсмической зоны; 3 - выбранные для анализа оперяющие разломы рифта; 4 - предсказанные события (см. рис. 6б) с магнитудой М > 4,8; 5 - доверительные области событий для сейсмической зоны; 6 - доверительные области событий для оперяющих разломов; 7 - наиболее опасные участки оси сейсмической зоны и разломов; 8 - предполагаемые наиболее опасные участки оси сейсмической зоны и разломов
Таблица 3
Таблица прогноза для оперяющих разрывов рифта
Номер области Период, годы
1 2011-2020
2 2009-2020
3 2010-2020
4 2013-2020
5 2010-2020
6 2008-2020
7 2013-2020
8 2011-2020
9 2013-2016
10 2010-2020
11 2008-2020
12 2008-2020
Заключение
На основе ограниченной базы по сейсмичности выделены наиболее опасные участки рифта и временные промежутки для каждого из них. Прогноз сделан до 2020 года двумя способами: для сейсмической зоны в целом (события с магнитудой М > 4,8) и для оперяющих активных разломов (события с магнитудой М > 3,1). По опыту, полученному в ходе данной работы, можно сделать вывод, что наиболее четкая картина волнового процесса в разломах прослеживается при рассмотрении событий с магнитудой М > 3,1. Эта цифра может варьироваться в зависимости от масштаба рассматриваемой структуры. Она объясняется различием в масштабах и энергетических параметрах деформационных волн при переходе на более крупные структуры литосферы, а, следовательно, и на другие глубинные уровни волн, контролирующих сейсмический процесс.
Понимание базовой основы сейсмического процесса, построение модели сейсмической зоны открывают возможности для выявления закономерностей в сейсмичности и использовании волновой концепции в качестве основы для прогноза землетрясений. Для полноты картины и более точного прогноза и оценки параметров необходимо иметь расширенную базу данных о сейсмичности за прошедшие годы (как минимум за 50 лет) и разломной тектонике изучаемого региона. Для усовершенствования и доработки представленной методики необходимо более широкое ее применение на различных сейсмических зонах.
Автор выражает благодарность своему научному руководителю С.И. Шерману за постановку темы, предоставленные материалы и обсуждение результатов ра6оты.
Статья поступила 01.04.2014 г.
Библиографический список
1. Deng Qidong. Map of Active Tectonics in China (1:4000000). Beijing: Seismological Press, 2007.
2. International Seismologocal Centre bulletin: event catalogue search http://www.isc.ac.uk/iscbulletin/search/catalogue/
3. Jiankun He, Dongsheng Cai, Yanxing Li, Zaisheng Gong. Active extension of the Shanxi rift, north China: does it result from anticlockwise block rotations? // Terra Nova. 2004. V. 16. P. 38-42.
4. Kasahara K. Migration of Crustal Deformation // Tectonophys-ics. 1979. V. 52. P. 329-341.
5. Mogi K. Migration of Seismic Activity // Bull. Earthq. Res. Inst., Tokyo University. N. 46. 1968. P. 53-74.
6. Sherman S.I. Deformation waves as a trigger mechanism of seismic activity in seismic zones of the continental lithosphere // Geodynamics & Tectonophysics. 2013. V. 4. № 2. P. 83-117.
7. Sherman S.I. Faults of the Baikal Rift Zone // Tectonophysics. 1978. V. 45, Р. 31-39.
8. Sherman S.I., Gorbunova E.A. Variation and origin of fault activity of the Baikal rift system and adjacent territories in real time // Earth science frontiers. 2008. V. 15, № 3, P. 337-347.
9. Sherman S.I., Ma Jin, Dem'yanovich V.M., Guo Yanahuang. The Baikal and Shanxi rift systems: tectonophysical relationships of earthquake epicenter and hypocenter fields // Continental rifting, accompanying processes. Proceeding of the Second All-Russia symposium with international participation. Irkutsk. 2013. V. 2. P. 185-187.
10. Sherman S.I., Zlogodukhova O.G. Seismic Belts and Zones of the Earth: Formalization of Notions, Positions in the Lithosphere, and Structural Control // Geodynamics & Tectonophysics. 2011. V. 2. № 1. P. 1-34.
11. Wang Jing-Ming. The Fenwei rift and its recent periodic activity // Tectonophysics. 1987. V. 133. P. 257-275.
12. Xu Xiwei, Ma Xingyuan, Deng Qidong. Neotectonic activity along the Shanxi rift system, China // Tectonophysics. 1993. V. 219. P. 305-325.
13. Викулин А.В. Физика волнового сейсмического процесса. Петропавловск-Камчатский: Изд-во КГПИ, 2003. 150 с.
14. Гатинский Ю.Г., Владова Г.Л., Прохорова Т.В., Рундквист Д.В. Геодинамика Центральной Азии и прогноз катастрофических землетрясений // Пространство и время. 2011. № 3.
С. 124-134.
15. Горбунова Е.А., Шерман С.И. Геоинформационная система фиксирования деформационных волн в сейсмоактивных зонах литосферы. Свидетельство о государственной регистрации программы для ЭВМ № 2013612772 от 13 марта 2013 г.
16. Горбунова Е.А., Шерман С.И. Медленные деформационные волны в литосфере: фиксирование, параметры, геодинамический анализ (Центральная Азия) // Тихоокеанская геология. 2012. Т. 31. № 1. С. 18-25.
17. Милановский Е.Е. Основные этапы рифтогенеза на территории Китая. М.: Междуведомственный геофизический комитет при Президиуме АН СССР, 1991. 148 с.
18. Николаевский В.Н. Тектонические волны Земли // Природа. 1991. № 8. С. 17-23.
19. Николаевский В.Н., Рамазанов Т.К. Генерация и распространение волн вдоль глубинных разломов // Известия АН СССР, Физика Земли. 1986. № 10. С. 3-13.
20. Николаевский В.Н., Рамазанов Т.К. О волнах взаимодействия литосферы с астеносферой // Гидрогеодинамические предвестники землетрясений. М.: Наука, 1984. С. 120-128.
21. Николаевский В.Н., Рамазанов Т.К. Теория быстрых тектонических волн // Прикладная математика и механика. 1985. Т. 49. № 3. С. 426-469.
22. Уломов В.И. Волны сейсмогеодинамической активизации и долгосрочный прогноз землетрясений // Физика Земли. 1993. С. 43-53.
23. Шерман С.И. Тектонофизическая модель сейсмической зоны: опыт разработки на примере Байкальской рифтовой системы // Физика Земли. 2009. № 11. С. 8-21.
24. Шерман С.И., Борняков С.А., Буддо В.Ю. Области динамического влияния разломов (результаты моделирования). Новосибирск: Наука, 1983. 112 с.
25. Шерман С.И., Горбунова Е.А. Волновая природа активизации разломов Центральной Азии на базе сейсмического мониторинга // Физическая мезомеханика. 2008. Т. 11. № 1. С. 115-122.
26. Шерман С.И., Днепровский Ю.И. Поля напряжений земной коры и геолого-структурные методы их изучения. Новосибирск: Наука, 1989. 158 с.