CC BY
19
GEODYNAMICS & TECTONOPHYSICS
PUBLISHED BY THE INSTITUTE OF THE EARTH'S CRUST SIBERIAN BRANCH OF RUSSIAN ACADEMY OF SCIENCES
ISSN 2078-502X
2017 VOLUME 8 ISSUE 4 PAGES 1065-1082
https://doi.org/10.5800/GT-2017-8-4-0333
Specific velocity structure of the upper mantle in the Transbaikalia segment of the Mongolia-Okhotsk orogenic belt
V. M. Soloviev1, V. V. Chechelnitsky2, A. S. Salnikov3, V. S. Seleznev4, A. V. Liseikin4, N. A. Galyova1
1 Altai-Sayan Branch of Federal Research Center 'Geophysical Survey of RAS', Novosibirsk, Russia
2 Baikal Branch of Federal Research Center 'Geophysical Survey of RAS', Irkutsk, Russia
3 Siberian Research Institute of Geology, Geophysics and Mineral Resources, Novosibirsk, Russia
4 Seismological Branch of Federal Research Center 'Geophysical Survey of RAS', Novosibirsk, Russia
Abstract: The paper presents the results of deep seismic studies on Geophysical Reference Profile 1-SB (Srednear-gunsk - Ust-Karenga - Taksimo - Vitim) in East Transbaikalia, Russia. The 1200 km long profile crosses the major tectonic structures of the Central Asian fold belt: the Argun median massif, the Selenga-Stanovoy and Transbaikalia folded regions, and the Baikal rift zone. Its northwestern fragment extends into the Angara-Lena monocline of the Siberian platform. The southeastern (Transbaikalia) and northwestern (Baikal-Patom) fragments of the profile are based on the spot (differential) seismic sounding technique using explosions and 40-tonne vibrators. The southeastern (Transbaikalia) fragment shows small crustal thickness values (~40 km), an almost horizontal position of the Moho, and high velocities of longitudinal waves (~8.4 km/sec) beneath the Moho. The analysis of parallelism graphs and the dynamic expression of the wave refracted from the Moho suggests a less than 5-10 km thick layer of high velocities and low gradients below Moho. The database on the territory of Transbaikalia includes ~200 wave arrival times from large earthquakes, which were refracted at the Moho at distances of ~200-1400 km. As part of the tomographic interpretation, using additional DSS data on the Moho, the territory of Transbaikalia has been mapped to show the patterns of the threshold velocity values at the Moho. The seismic data was used to contour an area with high velocity values in the mantle in the central part of the Mongolia-Okhotsk orogenic belt and the neighboring fold structures of Transbaikalia. According to the analysis of the seismic and geologic data on the study area, the mantle layer with high velocity values in the Mongolian-Okhotsk orogenic belt may be represented by the eclogitic rock plates.
Key words: deep seismic sounding (DSS); velocity of longitudinal waves; hodographs from explosions and earthquakes; Moho; Mongolia-Okhotsk orogenic belt
RESEARCH ARTICLE Received: August 17, 2017
Revised: October 27, 2017
Recommended by E.V. Sklyarov Accepted: November 13, 2017
For citation: Soloviev V.M., Chechelnitsky V.V., Salnikov A.S., Seleznev V.S., Liseikin A.V., Galyova N.A., 2017. Specific velocity structure of the upper mantle in the Transbaikalia segment of the Mongolia-Okhotsk orogenic belt. Geodynamics & Tectonophysics 8 (4), 1065-1082. doi:10.5800/GT-2017-8-4-0333.
Для цитирования: Соловьев В.М., Чечельницкий В.В., Сальников А.С., Селезнев В.С., Лисейкин А.В., Галёва Н.А. Особенности скоростного строения верхней мантии Забайкалья на участке Монголо-Охотского орогенного пояса // Геодинамика и тектонофизика. 2017. Т. 8. № 4. С. 1065-1082. doi:10.5800/GT-2017-8-4-0333.
Особенности скоростного строения верхней мантии Забайкалья на участке Монголо-Охотского орогенного пояса
В. М. Соловьев1, В. В. Чечельницкий2, А. С. Сальников3, В. С. Селезнев4, А. В. Лисейкин4, Н. А. Галёва1
1 Алтае-Саянский филиал Федерального исследовательского центра «Единая геофизическая служба РАН», Новосибирск, Россия
2 Байкальский филиал Федерального исследовательского центра «Единая геофизическая служба РАН», Иркутск, Россия
3 Сибирский научно-исследовательский институт геологии, геофизики и минерального сырья, Новосибирск, Россия
4 Сейсмологический филиал Федерального исследовательского центра «Единая геофизическая служба РАН», Новосибирск, Россия
Аннотация: Представлены материалы глубинных сейсмических исследований в Восточном Забайкалье на опорном геофизическом профиле 1-СБ (п. Среднеаргунск - п. Усть-Каренга - г. Таксимо - п. Витим) длиной свыше 1200 км. Профиль пересекает крупные тектонические структуры Центрально-Азиатского складчатого пояса: Аргунский срединный массив, Селенгино-Становую и Забайкальскую складчатые области, Байкальскую рифтовую зону и в северо-западной части выходит в Ангаро-Ленскую моноклизу Сибирской платформы. На юго-восточном (Забайкальском) и северо-западном (Байкало-Патомском) фрагментах профиля по методике точечных дифференциальных сейсмических зондирований с использованием взрывов и мощных 40-тонных вибраторов отработаны детальные системы наблюдений глубинного сейсмического зондирования. В юго-восточной части профиля (Забайкальский фрагмент) выявлена небольшая мощность земной коры (~40 км) при практически горизонтальном залегании границы Мохоровичича и высокая скорость продольных волн (~8.4 км/с) под ней. Анализ графиков параллельности и динамической выразительности преломленной волны от границы Мохоровичича позволил предположить наличие высокоскоростного слабоградиентного слоя мощностью не менее 5-10 км ниже границы Мохоровичича. По территории Забайкалья были собраны данные о ~200 временах пробега волн от крупных землетрясений, преломленных на границе Мохо-ровичича на удалениях ~200-1400 км. В рамках томографической интерпретации, с использованием дополнительных сведений о границе Мохоровичича по данным ГСЗ, построены площадные карты распределений граничной скорости по границе Мохоровичича на территории Забайкалья. Привлечение данных сейсмологии позволило оконтурить участок площади с высокими значениями скорости в подкоровой мантии центральной части Монголо-Охотского орогенного пояса и сопредельных складчатых структур Забайкалья. На основе анализа полученных сейсмических данных и геологических сведений по региону исследований высказано предположение, что в районе Монголо-Охотского орогенного пояса высокоскоростной слой мантии может быть представлен пластинами эклогитоподобных пород.
Ключевые слова: профиль ГСЗ; скорость продольных волн; годографы от взрывов и землетрясений; граница Мохоровичича; Монголо-Охотский орогенный пояс
1. Введение
Территория Забайкалья хорошо изучена геофизическими методами, в особенности методом глу-
бинного сейсмического зондирования (ГСЗ) [Ригу-1981, 1993; МгзИеп'кт еЬ а!., 1999; Seleznev еЬ а!., 2003; Suvorov еЬ а!., 2002]. Главным образом это от-
носится к западной и северной части Забайкалья в пределах Байкальской рифтовой зоны (БРЗ). По результатам работ ГСЗ 70-80-х годов прошлого столетия были составлены детальные карты глубины залегания подошвы земной коры (поверхно-
сти Мохоровичича) и скорости продольных волн непосредственно под ней [Puzyrev, 1981; М1зНеп'кт еЬ а!., 1999; Suvorov еЬ а!., 2002]. Согласно [Puzyrev, 1981; МгзИеп'кт еЬ а!., 1999; Suvorov еЬ а!., 2002], толщина коры в среднем составляет 40 км (от 35 до 45 км), а скорость под корой изменяется от пониженных значений 7.7-7.8 км/с в пределах БРЗ до нормальных значений 8.1-8.2 км/с за ее границей. Практически неизученной до настоящего времени оставалась юго-восточная часть Забайкалья. В 2014-2015 гг. Алтае-Саянский филиал ФИЦ ЕГС РАН совместно с ООО «Геотех» и АО «СНИИГГиМС»
по заданию Федерального агентства по недропользованию (Роснедра) на опорном геофизическом профиле 1-СБ (Восточный) по маршруту п. Средне-аргунск - п. Усть-Каренга - г. Таксимо - п. Витим, протяженностью свыше 1200 км, проведены глубинные сейсмические исследования методом ГСЗ с использованием взрывов в озерах и 40-тонных вибраторов в качестве источников возбуждения. В статье представлены материалы сейсмических исследований границы Мохоровичича на юго-восточном (Забайкальском) фрагменте опорного профиля 1-СБ (восточный участок), пересекающем крупные тектонические структуры Центрально-Азиатского складчатого пояса: Аргунский срединный массив, Селенгино-Становую и Забайкальскую складчатые области (рис. 1, а). Представлены также кинематические годографы по Байкало-Патомскому фрагменту профиля 1-СБ, пересекающему в северной части Байкальскую рифтовую зону. Результаты профильных исследований методом ГСЗ в Прибайкалье и Забайкалье дополнены данными о распределении скорости непосредственно под границей Мохоровичича (граничной скорости). Данные получены при площадной томографической интерпретации времен пробега преломленной на этой границе волны от землетрясений (рис. 1, б).
2. ИЗУЧЕНИЕ ГРАНИЦЫ МОХОРОВИЧИЧА НА ПРОФИЛЕ 1-СБ (Забайкальский фрагмент)
Работы ГСЗ на 640-километровом профиле выполнялись по методике точечных дифференциальных сейсмических зондирований [Иейупзку, Kosmin-Бкауа, 1978]. Была реализована достаточно плотная система наблюдений (рис. 1), с расстояниями в среднем 15-30 км между источниками возбуждения (взрывы в озерах массой заряда 3-6 т, мощные вибраторы) и в среднем 4-5 км между регистрирующими станциями (РОСА, Байкал) с группами вертикальных приборов СВ-5 и SerselSG-5 и 20 км - с трехкомпонентными приборами GS-20DX.
Для определения граничной скорости под границей Мохоровичича (М) использовались годографы преломленных продольных волн на удалениях ~200-350 км; с отдельных годографов (ПВ 10, 12, 21, 22) были зарегистрированы преломленные волны от границы М на удалениях до 400-500 км (рис. 2, а, б). Анализ протяженных годографов преломленных волн показывает, что резких изломов годографов (к более высоким значениям скорости) не происходит даже на удалениях свыше 450 км.
На основании графиков параллельности для ПВ 2, ПВ 4 в прямом и ПВ 22, ПВ 18 во встречном направлении (рис. 2, в) установлено, что в пределах протяженной совместной области прослеживания
(~110-170 км) годографы не сходятся. При этом отмечается разброс значений разности редуцированных времен в пределах 0.10-0.15 с, что практически соответствует погрешности определения времен вступлений на удалениях 200-400 км. Отсутствие схождения годографов на таких больших расстояниях свидетельствует о наличии высокоскоростного слабоградиентного слоя, расположенного ниже границы М. Мощность такого слоя должна быть не меньше 5-10 км, что следует из динамической выразительности преломленной волны (^пр) на этой границе.
По годографам Рп^-волн, достроенным с использованием начальных точек отраженных волн от границы М (рис. 2, а, б), были определены графики значений кажущейся скорости в прямом и встречном направлении, осредненные значения кажущейся скорости и граничная скорость (с учетом сейсмического сноса) (рис. 3) [Seleznev еЬ а!., 1987; Soloviev еЬ а!., 2014].
По данным продольных отраженных волн от поверхности Мохоровичича (РоТ,-волн) на профиле 1-СБ (восточный участок, забайкальский этап) определялась средняя (эффективная) скорость продольных волн для всей толщи земной коры (Уэф). С учетом высокой плотности годографов РоТ?р-волн (рис. 4, а) определение Vэф производилось по центрам зондирований (рис. 4, б) (согласно технике вычислений, описанной в работе [Soloviev еЬ а!., 2014]) и по сводному годографу в квадратичных координатах (рис. 4, в). Осредненные значения эффективной скорости продольных волн в земной коре оказались достаточно близки по разным определениям и составили 6.4±0.1 км/с. Наибольший разброс данных при определении по центрам зондирований получен на северо-западном участке профиля, что связано с невысокой плотностью наблюдений.
Глубины до поверхности Мохоровичича на исследуемом профиле определялись с использованием данных отраженных и преломленных волн от границы М в рамках простых соотношений: для отраженных волн Н=^((У2эф^2-12)/4), где Н, Vэф, ^ 1 - соответственно глубина до границы М, эффективная скорость продольных волн в земной коре, время пробега отраженных волн от границы М и удаление; для преломленных волн Н= =(t-1/Vг)•(Vэф/2cos(i)), где Н, Vэф, t, 1, Vг - соответственно глубина до границы М, эффективная скорость продольных волн в земной коре, время пробега преломленных волн от границы М, удаление и
граничная скорость по границе М, I = arcsin (-¡э^)
[Гейу^ку, Kosminskaya, 1978]. При расчетах использовались представленные выше осредненные зна-
Y, км 800
600
400
200
б о
-200
-400
-600
-800
-1000 -800 -600 -400 -200 0 200 400 600 800 1000 1200 1400 1600 X, км
Рис. 1. Район исследований границы Мохоровичича в Забайкалье. а - схема тектонического районирования согласно [РвЬтоу, 2012] в полосе опорного профиля 1-СБ (восточный участок). Цифрами на карте представлены структурные элементы Байкальской рифтовой системы: 37 - Чарско-Чаруодинская рифтогенная впадина, 38 - Муйская рифтогенная впадина, 39 - Верхнеангарская рифтогенная впадина, 40 - Ципинская рифтогенная впадина, 41 - Байкальский рифт, 42 - Баргузинская рифтогенная впадина, 43 - Амалатское лавовое поле, 44 - Киренская рифтогенная впадина, 334 - Яросимо-Олгуйдахская гряда, 340 - Илимская впадина, 341 - Илгинская впадина; б - схема площадных сейсмологических наблюдений.
Fig. 1. The Moho study area in Transbaikalia, Russia. a - scheme of tectonic zoning (after [Petrov, 2012] in the band of Reference Profile 1-SB (Eastern section). The structural elements of the Baikal rift system are numbered: 37 - Chara-Charuodin rift valley, 38 - Muya rift valley, 39 - Upper Angara rift valley, 40 - Tsipa rift valley, 41 - Baikal rift, 42 - Barguzin rift valley, 43 - Amalat lava field, 44 - Kirensk rift valley, 334 - Yarosimo-Olguydakh ridge, 340 - Ilim depression, 341 - Ilga depression; 6 - scheme of the area covered by seismological observations.
T-L/8, С 9.0
8.5
5.5 5.0 4.5 4.0 L
100 200 300 400 500 600 700 X, км
T-L/8, с 9.0
8.5
4 00 100 200 300 400 500 X, км
220 240 260 280 300 320 340 X, км
0.2
в 0
-0.2
At, с
At, с 1.0
0.9
0.8
0.7
0 20 40 60 80 100 120 140 160 180 X, км
Рис. 2. Интерпретация преломленных волн от поверхности Мохоровичича на профиле 1-СБ (восточный участок, забайкальский этап). а, б - редуцированные годографы преломленных волн от границы Мохоровичича с начальными точками по данным отраженных волн в прямом и встречном направлении, соответственно; в - графики параллельности для ПВ 2, ПВ 4 (верхний график) в прямом и ПВ 22, ПВ 18 (нижний график) во встречном направлении (по оси Х - километраж по профилю ГСЗ, по оси Y - различие редуцированных времен).
Fig. 2. Interpretation of the refracted waves from the Moho in Profile 1-SB (Eastern section, Transbaikalia stage). а, б - reduced hodographs of refracted waves from the Moho with initial points according to the reflected waves in the forward and opposite directions, respectively; в - parallelism graphs for ПВ 2, ПВ 4 (upper graph) in the forward direction, and ПВ 22, ПВ 18 (lower graph) in the opposite direction (axis X - kilometers along the DSS profile; axis Y - reduced time difference).
V, km/C
\
\ V
д
/ 1
/ • <-—-Vrp
0 100 200 300 400 500 600 700 X, км
Рис. 3. Определение граничной скорости преломленных волн под границей Мохоровичича на профиле 1-СБ (восточный участок, забайкальский этап): графики кажущейся скорости в прямом и встречном направлении и граничной скорости под границей М (выделена черным цветом), соответственно.
Fig. 3. Determination of the threshold velocity of refraction waves beneath the Moho in Profile 1-SB (Eastern section, Transbaikalia stage): curves of the apparent velocity in the forward and opposite directions and the threshold velocity beneath the Moho (highlighted in black), respectively.
чения скоростей продольных волн в земной коре (дополнительно сглаженные окном в 50 км). Граничная скорость усреднялась на участке пробега вдоль границы Мохоровичича по данным, представленным на рис. 3.
Результаты изучения границы Мохоровичича на профиле 1-СБ (восточный участок, забайкальский этап) представлены на глубинном сейсмическом разрезе на рис. 5. Одной из особенностей глубинного сейсмического разреза на Забайкальском участке профиля 1-СБ является практически плоское горизонтальное залегание границы Мохоровичича (на глубинах ~40 км) с высокими значениями граничной скорости 8.3-8.5 км/с. Эффективная скорость распространения сейсмических волн в земной коре составляет 6.4±0.1 км/с. Наиболее высокие значения граничной скорости Р-волн (до 8.55 км/с) непосредственно под корой отмечаются на северо-западном участке Газимурского блока (Х=190-240 км на рис. 5); несколько пониженные значения граничной скорости Р-волн вдоль границы М (8.25-8.30 км/с) отмечены в северной части Витимо-Урюмского блока на северо-западе профиля (Х=465-505 км). Полученные на Забайкальском участке профиля ГСЗ высокие значения граничной скорости ранее при глубинных сейсмических исследованиях не отмечались [Puzyrev, 1981; Mishen'-Мп et al., 1999; Seleznev et al., 2003; Suvorov et al., 2002], что вызывает повышенный интерес геоло-
гов и геофизиков к этому участку, диктуя необходимость дальнейшего его изучения и, в первую очередь, площадного оконтуривания высокоскоростного блока подкоровой мантии.
3. Площадные сейсмологические исследования по
ИЗУЧЕНИЮ СКОРОСТИ продольных волн по границе М в Забайкалье
С целью локализации высокоскоростного блока под границей Мохоровичича в Забайкалье проводилось площадное изучение граничной скорости с использованием данных землетрясений. В обработке было использовано около 200 времен пробега преломленной на границе М продольной волны на удалениях от 200 до 1470 км из базы сейсмологических данных Байкальского филиала ФИЦ ЕГС РАН за период с 2002 по 2014 г. [Masalsky et al., 2007; International Seismological Centre, 2017]. Эти данные отличаются высокой точностью, поскольку были получены на цифровой регистрирующей аппаратуре от крупных землетрясений с энергетическим классом К от 12.0 до 15.9 (М~4.5-6.5); для них были определены с высокой точностью эпицентры и с меньшей точностью (в среднем не грубее чем ±5 км) - глубины землетрясений. Лучевая схема, иллюстрирующая площадное освещение изучаемого участка данными от землетрясений, представ-
1
6+
2 3
+ +
4
+
5
+
6
+
7 8
+ +
9
+
10 11 12
+ + +
13
+
14 15 16 17 18 19 20
+ + + + + + +
21
+
22
+
23 24
+ +
25
+
26
+
100
200
300
400
500
600 X, км
V, км/с 6.50
6.45
6.4ok
6.35
6.30
6.25
• •
• в
• •
• * •
• •
• •
100
200
300
400
500
600 X, км
в
T2-L2/6.42, с2 200
150
100
# • • 6.4 км/с
•---fS ■*/.. • • * •
10000 20000 30000 40000 50000 60000
L2, км2
Рис. 4. Интерпретация отраженных волн от границы Мохоровичича на профиле 1-СБ (восточный участок, забайкальский этап), а - редуцированные годографы отраженных волн от границы М; б, в - определение эффективной скорости продольных волн в земной коре (УЭф) по центрам зондирований (окно 50 км, шаг 25 км) и по сводному редуцированному годографу в квадратичных координатах, соответственно.
Fig. 4. Interpretation of the reflected waves from the Moho in Profile 1-SB (Eastern section, Transbaikalia stage), a - reduced hodographs of the reflected waves from the Moho; 6, e - determination of the effective velocity of the longitudinal waves in the crust (V3Ij,) at the seismic sounding centers (box: 50 km; spacing: 25 km) and the composite reduced hodograph in quadratic coordinates, respectively.
о
Восточно-Забайкальская складчатая область
Монголо-Охотский пояс
Н, км Заурулюнгуйский О
Газимурский
Пришил кинская ,1ЙзонаЖ
Букачача-
БорщовочныйзонаЖирекенский Сырыгичинский Юмурченский
Саяно-Байкальская складчатая область
Нерча-
Витимо- Нюкжа-Урюмский Юмурченский
-40
по данным отраженных волн по данным преломленных волн
V, км/с
■ 8.40
= 7.80
= 7.40
i 7.15
7.00
6.85
— 6.70
= 6.55
= 6.40
— 6.10
5.80
= 5.50
— 5.20
4.90
0 100 200 300
| Рис. 5. Глубинный сейсмический разрез по профилю 1-СБ (Забайкальский участок). | Fig. 5. Deep seismic section in Profile 1-SB (Transbaikalia).
400
500
600 Xm„ км
Y, км 800
600
400
200
0
-200 -400 -600 -800
-1000 -800 -600 -400 -200 0 200 400 600 800 1000 1200 1400 1600 X, км
Рис. 6. Лучевая схема сейсмологических наблюдений в Прибайкалье и Забайкалье.
Fig. 6. Radial diagram of seismic observations in Pribaikalie and Transbaikalia.
• города профиль 1 профиль 2 :| землетрясение станция
лена на рис. 6. Для расчетов томографических моделей географические координаты землетрясений и регистрирующих станций были пересчитаны в плоскую прямоугольную систему координат с ор-тографической (азимутальной) проекцией с началом координат в точке 110°Е и 54°^ Видно, что наиболее полно данные о временах пробега волн от границы М представлены для западного и северного участков Забайкалья. Из анализа материалов ГСЗ прошлых лет [Puzyrev, 1981; Mishen'kin et а1., 1999; Seleznev et а1., 2003; Suvorov et а1., 2002] и представленных выше данных по Забайкальскому участку профиля 1-СБ следует, что средние глубины до границы М составляют 40 км, а средние скорости в земной коре (V) - 6.4 км/с. С использованием этих данных и сведений о глубинах гипоцентров землетрясений по каждому отрезку (эпицентр землетрясения - сейсмостанция) в рамках простых соотношений из работы [Soloviev et а1., 2016] определялись значения граничной скорости (УТ) по границе М:
^ = ~~(2Н-Й)С05(1)' , (1)
V
где I = arcsin I, ^ Н и h - эпицентральное расстояние, время пробега от гипоцентра до сейсмо-
станции, глубина до границы М и глубина очага землетрясения, соответственно. Покрывающая среда при этом предполагается однородной, а граница - плоской и горизонтальной.
Эти же данные (средние скорости в коре, средняя глубина до границы М, данные времен от землетрясений, сведения о глубинах и координатах землетрясений) позволяют ввести поправку за глубину очага во времена вступлений продольных волн от границы М и построить сводный годограф от землетрясений для территории Прибайкалья и Забайкалья. Разность времен одной и той же головной волны от поверхностного и заглубленного источников от к-того события до приемных ^станций связана с глубиной гипоцентра (^ следующим соотношением:
^ = 7^(0, где £ = агат (2)
Ниже приведены фрагменты сводных годографов продольной преломленной волны от границы Мохоровичича, полученные как от данных регистрации землетрясений на площадной сейсмологической сети, так и от данных с Забайкальского и Байкало-Патомского участков профиля 1-СБ от взрывов и мощных вибраторов (рис. 7). Из представленного рисунка видно, что достаточно близ-
T-L/8, с 10
9
8 •
7
•
6
5
4
3 200
T-L/8, с
ф
• • • % • #
...... • • • • • * ч V» у П ' '
• « V » • • • • • • • • • • • •
• • • . • • •
•
250
300
350
400
450
500
550
600 L, км
. Ф Z. А. «/
« / «* * ♦ ♦ ♦ /н * ♦ ♦ ♦ \ ** ♦ ♦ х ♦
♦ ♦♦ ♦ ♦ ♦ ** » ♦ ♦
600 L, км
T-L/8, с
в
• t*
• •
30 25 0 ЗС )0 3£ >0 4С о 4е >0 5( Ю 5£
550 L, км
Рис. 7. Редуцированные годографы продольных преломленных волн от границы Мохоровичича от землетрясений, химических взрывов и мощных вибраторов. а - данные от землетрясений, полученные на площадной сейсмологической сети в Забайкалье (см. рис. 6); б, в - данные от взрывов и вибраторов, полученные на Байкало-Патомском и Забайкальском фрагментах профиля 1-СБ, соответственно.
Fig. 7. Reduced hodographs of the longitudinal refracted waves from the Moho from earthquakes, chemical explosions and strong vibrators. а - earthquake data obtained by the areal seismic network in Transbaikalia (see Fig. 6); б, в - data from explosions and vibrators in the Baikal-Patom and Transbaikalia fragments of Profile 1-SB, respectively.
кие годографы получены от землетрясений и по данным ГСЗ на Байкало-Патомском фрагменте профиля 1-СБ (рис. 7, а, б). Осредненные значения кажущейся скорости по ним также близки и составляют 8.0±0.2 км/с. Большой разброс значений времен на этих годографах (в особенности на годографе от землетрясений) обусловлен разнородно-
стью данных, полученных как с большой территории, так и с участков с различным строением зоны перехода кора-мантия, и, соответственно, скоростным строением (нормальная и аномальная мантия) (по [Puzyrev, 1981; Mishen'kin et al., 1999; Seleznev et al., 2003; Suvorov et al., 2002]). Значительно отличается от этих данных сводный годограф продольной
V, км/с 8.8 8.6 , 8.4 8.2 8.0 7.8 7.6 7.4
»
• * i j • j • • •
»¿¿J . * • •
wr« > V •
• \ * • •
•
0 200 400 600 800 1000 1200 1400 1600 L, км
I Рис. 8. Значения кажущейся скорости Р^-волны для разных эпицентральных расстояний в Прибайкалье и Забайкалье.
Fig. 8. The apparent velocity values of P^-wave for different epicentral distances in Pribaikalie and Transbaikalia.
преломленной волны от границы М, полученный на Забайкальском участке профиля 1-СБ (рис. 7, в). По сравнению с годографами на рис. 7, а, б, здесь, при общем понижении примерно на 1 секунду, существенно меньше разброс значений редуцированных времен, а осредненные значения скорости по годографу повышены и составляют 8.4±0.1 км/с. Эти данные еще раз указывают на обособленность по скоростному строению мантии центрального участка Забайкальского фрагмента профиля 1-СБ от прилегающих территорий.
С целью его пространственной локализации по сейсмологическим данным проведена оценка значений кажущейся скорости Рп^-волны для различных интервалов регистрации и выполнена площадная томографическая интерпретация годографов от землетрясений.
Из представленного на рис. 8 графика всей совокупности распределения кажущейся скорости (рассматриваемой, в первом приближении, как распределение граничной скорости Рп^-волны, вычисленной по формуле (1) для нарастающих интервалов регистрации), видно, что до удалений в 700 км отмечается значительный разброс данных. В интервале удалений 200-400 км значения скорости находятся в пределах 7.6-8.3 км/с, а на удалениях до 400-600 км они несколько повышены и составляют 7.8-8.4 км/с. На больших интервалах удалений (от 800 до 1200-1470 км) разброс существенно уменьшен и (за исключением пары точек) не превышает 0.1 км/с; значение скорости при этом составляет ~8.05 км/с. Представленный график указывает на два обстоятельства: существенное влияние аномальной низкоскоростной мантии на удалениях до 600 км и наличие одного или нескольких непротяженных высокоскоростных блоков по границе М со скоростью 8.2-8.6 км/с. Наличие же протяженного
высокоскоростного блока по границе М в Прибайкалье и Забайкалье существенно бы подняло значение скорости на удалениях свыше 800 км. Была выполнена оценка влияния ошибок в определении глубин землетрясений на значения кажущейся скорости Рп^-волны для разных удалений. Для единичных определений ошибка в определении времен в ±0.5 секунды (за счет средних ошибок в глубине гипоцентра ~ в ±5км) приведет к ошибкам в значениях скорости ~ в ±0.15 км/с на малых базах (200-400 км) и к менее значительным ошибкам ~ в ±0.05 км/с на больших удалениях - 800-1200 км. При осреднении данных ошибки в определении скорости Рп^-волны за счет ошибок в глубинах гипоцентров будут существенно (в Тп, где п - количество определений) меньше, в особенности для больших удалений.
Полученные значения граничной скорости на отрезках лучей в пределах площади на рис. 6 были использованы для дальнейшей томографической интерпретации. Первоначально в исходные времена и расстояния вводились поправки за сейсмический снос, чтобы полученные таким образом времена относились к участку пробега по границе М. Площадь была разбита на равные квадраты со стороной 200 км. Далее была составлена система уравнений:
^ = (3)
где £( - времена пробега волны вдоль /-го луча по границе М; а^ - длина сегмента /-го луча в у-м блоке разбиения поверхности М; = ^ - длина /-го луча по границе М; Vj - неизвестная граничная скорость на поверхности мантии в у-м блоке разбиения поверхности. После ввода замены
-200
-400
1200 X, км
Рис. 9. Схема распределения граничной скорости по поверхности Мохоровичича, полученная при томографической инверсии системы уравнений (5) методом QR-разложения.
Fig. 9. The pattern of the threshold velocity on the Moho, according to the tomographic inversion of the system of equations (5) by the QR decomposition.
Xj = —
J Vj
система (3) приобретает вид:
ti = l
(4)
(5)
Полученная система (5) является разреженной (большое число нулевых значений) и переопределенной (уравнений намного больше, чем неизвестных). Существуют разные способы решения подобных систем. В представленной работе использовано два способа: ^-разложение [ВакИуа!оу, 1988] и метод регуляризации [Ыо!е% 1987].
Метод ^-разложения используется для приближения методом наименьших квадратов. Риском
решения системы этим методом является возможное большое число обусловленности. На рис. 9 приведено решение данным методом для всех лучей.
Вторым способом решения является метод регуляризации. Данный метод - это способ предотвратить влияние экспериментальных ошибок в правом столбце на решение. Вместо (5) решается задача:
(6)
где йЬ - разница экспериментального времени и времени, рассчитанного по референтной модели среды. Значения скоростей для референтной модели основаны на схемах распределения граничной
Y, KM 800
-400
-800
-600
-400
-200
200
400
600
800
1000
1200 X, км
Рис. 10. Схема распределения граничной скорости по поверхности Мохоровичича, полученная при томографической инверсии системы уравнений (6) методом регуляризации (с использованием референтной модели).
Fig. 10. The pattern of the threshold velocity on the Moho, according to the tomographic inversion of the system of equations (6) by the regularization method (using the reference model).
скорости на поверхности мантии, полученной методом ГСЗ в 70-80-х годах прошлого столетия для территории Байкальской рифтовой зоны, юга Сибирской платформы и разновозрастных складчатых областей Прибайкалья и Забайкалья [Puzyrev, 1981, 1993; Mishen'kin et а1, 1999; Seleznev et а1., 2003; Suvorov et а1., 2002], и данных профиля 1-СБ 20142015 гг. ^о^^ et а1., 2016; Soloviev et а1., 2017].
Из (6) находится решение системы
х=[АТА + Л21У1АТ(И, (7)
где х - отличие медленности (величина, обратная скорости), искомой от референтной модели. Искомая скорость находится следующим способом:
V; «УР-х^0)2. (8)
Полученное распределение граничной скорости методом регуляризации представлено на рис. 10.
Анализ результатов расчетов граничной скорости разными способами показывает, что метод @Р-разложения дает более контрастную карту распределения граничной скорости. Различие минимальных и максимальных значений скорости на исследуемой площади составляет более 1.25 км/с, что в полтора раза выше, чем на карте граничной скорости, полученной методом регуляризации. При сравнении данных, полученных разными методами, более высокие значения граничной скорости Р-волн (до 0.33 км/с) получены методом ^Р-разло-
жения в центральной части площади (в пределах Забайкальского фрагмента профиля 1-СБ); более низкие значения скорости получены методом ^-разложения по сравнению с данными, полученными методом регуляризации (от 0.29 до 0.63 км/с) в северо-восточной части площади. Примерно в половине точек площади получены достаточно близкие значения скорости разными способами (не более ±0.1 км/с).
К чисто модельным погрешностям построенных карт распределений граничной скорости следует отнести частичное выпадение непротяженных участков площади с пониженной скоростью на границе Мохоровичича (участков аномальной мантии со скоростью 7.7-7.9 км/с и мощностью слоев в 10-15 км). По данным [Ршугеу, 1981, 1993; МзНеп'-кш еЬ а1., 1999; Быуогоу еЬ а1., 2002], в годографах сейсмических волн в Забайкалье продольная волна от слоя пониженной скорости («аномальной мантии») и небольшой мощности прослеживается до удалений 280-300 км. Далее она сменяется волной от границы М1 с нормальной скоростью в 8.0-8.2 км/с [Ршугеу, 1981; МИзНеп'кт еЬ а1., 1999]. Фактически на этих участках мы имеем дело с двухслойной моделью верхней мантии. При использовании для томографии отрезков годографов для расстояний до 1000 км и более в результатах интерпретации в большей мере будет присутствовать информация о скорости по границе нормальной мантии. В Прибайкалье (район оз. Байкал) слой аномальной мантии сильно увеличен [Ри2угеу, 1993] и длина годографов менее сказывается на результатах определения граничной скорости. В настоящей работе были выполнены расчеты с максимальными удалениями от землетрясений не более 600 км. На ряде участков БРЗ были получены более низкие значения скорости продольных волн по границе М, близкие к значениям в референтной модели, построенной по данным ГСЗ, и адекватно отражающие аномальную мантию данных участков.
4. Обсуждение результатов
Несмотря на присутствие модельных погрешностей, результаты региональных сейсмических исследований вдоль опорного профиля 1-СБ и данные определений глубин землетрясений существенно дополнили построенные ранее карты граничной скорости. В северо-западной части района исследований (в пределах Байкальской рифтовой зоны) и северо-восточной части (Забайкалье) граничные скорости имеют пониженные значения в 7.80-7.95 км/с. В юго-восточной части площади в верхней мантии выделен высокоскоростной мантийный блок с поперечными размерами около
600x600 км. В пределах Забайкальского фрагмента профиля 1-СБ это область сочленения Восточно-Забайкальской и Саяно-Байкальской складчатых областей Центрально-Азиатского складчатого пояса, разделенных Монголо-Охотским поясом (см. рис. 5). Высокоскоростной участок от профиля 1-СБ вытянут в южном и в юго-западном (за Читу) направлении. В целом же, высокоскоростная зона по границе М прослеживается через всю исследуемую площадь с юго-запада на северо-восток (см. рис. 9, 10) и совпадает с простиранием Монголо-Охотского орогенного пояса. Центральным звеном последнего является Главный Монголо-Охотский разлом [ИНапоУБку, 1996; Раг/епоу еЬ а1., 2003]. Установленные при сейсмических исследованиях высокоскоростные неоднородности в верхней части разреза (см. рис. 5) могут являться интрузиями, проникшими в верхнюю кору через указанный глубинный разлом.
В работах [Бгйепко еЬ а1., 2013; Раг/епоу еЬ а1., 2003; БНеуеНепко еЬ а1., 2011] формирование поздне-палеозойско-раннемезозойского Монголо-Охотского орогенного пояса связывается с закрытием существовавшего на его месте океанического бассейна и субдукцией океанической коры под смежные континентальные окраины. В зонах субдукции океанической плиты возможно образование глаукофа-новых сланцев и эклогитов. Последние выступают как химические эквиваленты габбро и базальтов, но превышают их по плотности на 20-25 % [БоЬ-геЬБоу еЬ а1., 1989]. Их плотность (3.3-3.7 г/см3), а значит и скорость, выше, чем у пород верхней части верхней мантии. Можно предположить, что выделенный по данным сейсмологии и ГСЗ высокоскоростной слой в верхней мантии является пластиной эклогитов (или эклогитоподобных пород) в районе Монголо-Охотского орогенного пояса.
5. ЗАКЛЮЧЕНИЕ
Территория Прибайкалья и Забайкалья в структурно-тектоническом плане является частью крупного Центрально-Азиатского орогенического пояса. Развитие современных представлений о геологическом строении и характере эволюции Центрально-Азиатского складчатого пояса требует привлечения новых данных по глубинной геодинамике, тектонике и магматизму. Чрезвычайно важны при этом сведения по глубинному строению земной коры и верхов мантии по данным геофизических, и прежде всего сейсмических, методов. В последние годы в России и за рубежом главный объем информации о строении литосферы Земли, в том числе о земной коре и верхней мантии, поступает при сейсмотомографических исследованиях.
Распространены также исследования с временными сетями широкополосных станций для получения информации о строении коры и мантии при обработке данных методом обменных волн от землетрясений (МОВЗ). При таких построениях фактически размываются опорные сейсмические границы, такие как поверхность Мохоровичича, с большими погрешностями в определении скорости упругих волн вдоль нее. В то же время в результатах ГСЗ граница М строится очень точно (±3 %) по данным отраженных и преломленных волн с распределением вдоль нее скорости упругих волн, с точностью ±0.1 км/с. На больших удалениях (400500 км и более, Западная Якутия, Алтай, Забайкалье) преломленная волна от взрывов от границы М зачастую прослеживается без значительного проницания, фактически как головная. Распределение скорости по границе М можно определять без учета рефракции. В настоящей работе сделана попытка совместной интерпретации данных ГСЗ и сейсмологии для построения площадной карты распределения скорости продольных волн по границе М территории Прибайкалья и Забайкалья. Были использованы данные о двух сотнях времен пробега преломленной волны от границы М от крупных землетрясений (энергетический класс К выше 12.5), зарегистрированные на удалениях от 200 км до 1400 км. Результаты глубинных сейсмических исследований 70-80-х гг. прошлого столетия и данные профилей ГСЗ 2014-2015 гг. по границе М при этом выступали как опорные. Хотя за счет неравномерного распределения землетрясений и
сейсмостанций по площади недостаточно освещенной по границе М оказалась юго-восточная часть Забайкалья, все же была получена новая площадная информация о распределении скоростей под границей М на значительной территории Прибайкалья и Забайкалья. Привлечение данных сейсмологии позволило оконтурить участок площади с высокими значениями скорости (до 8.4-8.5 км/с) в подкоровой мантии центральной части Монголо-Охотского орогенного пояса и сопредельных складчатых структур Забайкалья. В северо-западной части площади (в пределах Байкальской рифтовой зоны) и северо-восточной части площади исследований установлены пониженные до 7.80-7.95 км/с граничные скорости продольных волн.
Дальнейшее развитие площадных сейсмических исследований территории Прибайкалья и Забайкалья связывается авторами с увеличением плотности сейсмического освещения площади за счет привлечения данных от менее сильных землетрясений, с использованием данных поперечных волн.
6. Благодарности
Авторы выражают признательность коллегам из Байкальского филиала ФИЦ ЕГС РАН, участвующим в регистрации и создании базы сейсмологических данных по Прибайкалью и Забайкалью, а также коллегам из АСФ и СЕФ ФИЦ ЕГС РАН и АО «СНИИГГиМС» за ценные замечания при обсуждении статьи.
7. Литература / References
Bakhvalov N.S., 1988. Numerical Methods. Nauka, Moscow, 314 p. (in Russian) [Бахвалов Н.С. Численные методы. М.: Наука, 1988. 314 с.].
Didenko A.N., Efimov A.S., Nelyubov P.A., Sal'nikov A.S., Starosel'tsev V.S., Shevchenko B.F., Goroshko M.V., Gur'yanov V.A., Zamozhnyaya N.G., 2013. Structure and evolution of the Earth's crust in the region of junction of the Central Asian fold belt and the Siberian platform: Skovorodino-Tommot profile. Russian Geology and Geophysics 54 (10), 12361249. https://doi.Org/10.1016/j.rgg.2013.09.008.
DobretsovN.L., SobolevN. V., Shatsky V.S., AschepkovI.V., BakirovA.B., GabovN.F., GilbertA.E., Dobretsova L.V., Kozyreva N. V., Kolcheva K., Sklyarov E.V, 1989. Eclogites and Glaucophane Schists in Fold Areas. Siberian Branch, Nauka, Novosibirsk, 236 p. (in Russian) [Добрецов Н.Л., Соболев Н.В., Шацкий В.С., Ащепков И.В., Бакиров А.Б., Габов Н.Ф., Гильберт А.Э., Добрецова Л.В., Козырева Н.В., Колчева К., Скляров Е.В. Эклогиты и глаукофановые сланцы в складчатых областях. Новосибирск: Наука. Сибирское отделение, 1989. 236 с.].
Fedynsky V.V., Kosminskaya I.P. (Eds.), 1978. Technique and Equipment for Regional Seismic Studies in Remote Areas and Their Application in Siberia. Siberian Branch, Nauka, Novosibirsk, 206 p. (in Russian) [Методика и аппаратура для региональных сейсмических исследований в труднодоступной местности и их применение в Сибири / Ред. В.В. Федынский, И.П. Косминская. Новосибирск: Наука. Сибирское отделение, 1978. 206 с.].
International Seismological Centre, 2017. On-line Bulletin. Thatcham, United Kingdom. Available from: http://www. isc.ac.uk (last accessed December 6, 2017).
Masalsky O.K., Chechelnitsky V.V., Gileva N.A., 2007. The state-of-the-art of seismic observations in Pribaikalie. In: Modern methods of seismological data processing and interpreting. Proceedings of the Second International Seismological School. GS RAS, Obninsk, p. 140-144 (in Russian) [Масальский О.К., Чечельницкий В.В., Гилева Н.А. Современное состояние сейсмических наблюдений в Прибайкалье // Современные методы обработки и интерпретации сейсмологических данных: Материалы Второй международной сейсмологической школы. Обнинск: ГС РАН, 2007. С. 140-144].
Milanovsky E.E, 1996. Geology of Russia and the Neighboring Countries (Northern Eurasia). Moscow State University, Moscow, 448 p. (in Russian) [Милановский Е.Е. Геология России и ближнего зарубежья (Северной Евразии). М.: МГУ, 1996. 448 с.].
Mishen'kin B.P., Mishen'kina Z.R., Petrik G.V., Shelud'ko I.F., Mandel'baum M.M., Seleznev V.S., Solov'ev V.M., 1999. Deep seismic sounding of the Earth's crust and upper mantle in the Baikal rift zone. Izvestiya, Physics of the Solid Earth 35 (7-8), 594-611.
Nolet G. (Ed.), 1987. Seismic Tomography with Applications in Global Seismology and Exploration Geophysics. Kluwer, Dordrecht, 386 p. [Русский перевод: Сейсмическая томография. С приложениями в глобальной сейсмологии и разведочной геофизике / Ред. Г. Нолет. М.: Мир, 1990. 416 с.].
Parfenov L.M., Berzin N.A., Khanchuk A.I., Badarch G., Belichenko V.G., Bulgatov A.N., Dril S.I., Kirillova G.L., Kuzmin M.I., Nokleberg W.J., ProkopyevA.V., Timofeev V.F., Tomurtogoo O., Yang H., 2003. A model for the formation of orogenic belts in Central and Northeast Asia. Tikhookeanskaya Geologiya (Russian Journal of Pacific Geology) 22 (6), 7-41 (in Russian) [Парфенов Л.М., Берзин Н.А., Ханчук А.И., Бадарч Г., Беличенко В.Г., Булгатов А.Н., Дриль С.И., Кириллова Г.Л., Кузьмин М.И., Ноклеберг У., Прокопьев А.В., Тимофеев В.Ф., Томуртогоо О., Янь Х. Модель формирования орогенных поясов Центральной и Северо-Восточной Азии // Тихоокеанская геология. 2003. Т. 22. № 6. С. 7-41].
Petrov O.V. (Ed.), 2012. Geological Map of Russia and Adjoining Water Areas. Cartographic Factory of VSEGEI, Saint Petersburg (in Russian) [Геологическая карта России и прилегающих акваторий. 1:2500000 / Ред. О.В. Петров. СПб.: Картографическая фабрика ВСЕГЕИ, 2012].
Puzyrev N.N. (Ed.), 1981. Interior of Baikal from Seismic Data. Nauka, Novosibirsk, 173 p. (in Russian) [Недра Байкала (по сейсмическим данным) / Ред. Н.Н. Пузырев. Новосибирск: Наука, 1981. 105 с.].
Puzyrev N.N. (Ed.), 1993. Detailed Seismic Studies of the Lithosphere Using P and S waves. Nauka, Novosibirsk, 199 p. (in Russian) [Детальные сейсмические исследования литосферы на P и S-волнах / Ред. Н.Н. Пузырев. Новосибирск: Наука, 1993. 199 с.].
Seleznev V.S., Soloviev V.M., Suvorov V.D., Kreinin A.B., Uarov V.F., 1987. Using the transverse waves in deep seismic sounding in West Yakutia. Geologiya i Geofizika (Russian Geology and Geophysics) (1), 109-117 (in Russian) [Селезнев В.С., Соловьев В.М., Суворов В.Д., Крейнин А.Б., Уаров В.Ф. Использование поперечных волн при глубинном сейсмическом зондировании в Западной Якутии // Геология и геофизика. 1987. № 1. С. 109-117].
Seleznev V.S., Tibo G., Perchuk E., Emanov A.F., Suvorov V.D., Soloviov V.M., Tat'kov G.I., Liseikin A.V., Annenkov V.V., Mishen'kina Z.R., 2003. Using the new deep seismic technology to study the deep structure of the southwestern flank of the Baikal rift zone. In: Problems of seismology in the third millennium. Proceedings of the International Geophysical Conference (15-19 September 2003, Novosibirsk). Publishing House of SB RAS, Novosibirsk, p. 324-329 (in Russian) [Селезнев В.С., Тибо Г., ПерчукЭ., Еманов А.Ф., Суворов В.Д., Соловьев В.М., Татьков Г.И., Лисейкин А.В., Анненков В.В., Мишенькина З.Р. Использование новой технологии глубинных сейсмических исследований при изучении глубинного строения юго-западного фланга Байкальской рифтовой зоны // Проблемы сейсмологии III тысячелетия: Материалы международной геофизической конференции (15-19 сентября 2003 г., г. Новосибирск). Новосибирск: Изд-во СО РАН, 2003. С. 324-329].
Shevchenko B.F., Goroshko M.V., Didenko A.N., Gur'yanov V.A., Starosel'tsev V.S., Sal'nikov A.S., 2011. The junction of the eastern Central Asian fold belt and the Siberian platform: deep structure and Mesozoic tectonics and geodynamics. Russian Geology and Geophysics 52 (12), 1672-1679. https://doi.org/10.1016/j.rgg.2011.11.008.
Solov'ev V.M., Seleznev V.S., Sal'nikov A.S., Shibaev S.V., Timofeev V.Y., Liseikin A.V., Shenmaier A.E., 2016. Deep seismic structure of the boundary zone between the Eurasian and Okhotsk plates in Eastern Russia (along the 3DV deep seismic sounding profile). Russian Geology and Geophysics 57 (11), 1613-1625. https://doi.org/10.1016/j.rgg. 2016.10.007.
Soloviev V.M., Salnikov A.S., Seleznev V.S., Elagin S.A., Romanenko I.E., Galyova N.A., 2017. Deep seismic studies of the Baikal-Patom fragment of the Eastern section of Reference Profile 1-SB. Interexpo Geo-Siberia 2 (4), 106-112 (in Russian) [Соловьев В.М., Сальников А.С., Селезнев В.С., Елагин С.А., Романенко И.Е., Галёва Н.А. Глубинные сейсмические исследования на Байкало-Патомском фрагменте восточного участка опорного профиля 1-СБ // Интерэкспо Гео-Сибирь. 2017. Т. 2. № 4. С. 106-112].
Soloviev V.M., Seleznev V.S., SalnikovA.S., Kashubina T.V., Shenmaier A.E., 2014. Effective seismic models of media used in DSS operations on base lines in Eastern Russia. Geology and Mineral Resources of Siberia (4), 74-86 (in Russian) [Соловьев В.М., Селезнев В.С., Сальников А.С., Кашубина Т.В., Шенмайер А.Е. Использование эффективных сейсмических моделей сред при работах ГСЗ на опорных профилях в восточной части России // Геология и минерально-сырьевые ресурсы Сибири. 2014. № 4. С. 74-86].
Soloviev V.M., Seleznev V.S., Salnikov A.S., Liseikin A.V., Romanenko I.E., Elagin S.A., Shenmaier A.E., 2016. Features of seismic structures in the Transbaikalia segment of the Central Asian fold belt in Geophysical Reference Profile 1-SB. Interexpo Geo-Siberia 2 (2), 234-238 (in Russian) [Соловьев В.М., Селезнев В.С., Сальников А.С., Лисейкин А.В., Романенко И.Е., Елагин С.А., Шенмайер А.Е. Особенности сейсмического строения структур забайкальской части Центрально-Азиатского складчатого пояса в створе опорного геофизического профиля 1-СБ // Интерэкспо Гео-Сибирь. 2016. Т. 2. № 2. С. 234-238].
Suvorov V.D., Mishenkina Z.M., Petrick G.V., Sheludko I.F., Seleznev V.S., Solovyov V.M., 2002. Structure of the crust in the Baikal rift zone and adjacent areas from deep seismic sounding data. Tectonophysics 351 (1-2), 61-74. https:// doi.org/10.1016/S0040-1951(02)00125-7.
Соловьев Виктор Михайлович, канд. геол.-мин. наук, зам. директора по научной работе Алтае-Саянский филиал Федерального исследовательского центра «Единая геофизическая служба РАН» 630090, Новосибирск, просп. Академика Коптюга, 3, Россия e-mail: solov@gs.nsc.ru
Soloviev, Victor M., Candidate of Geology and Mineralogy, Deputy Director for Science Altai-Sayan Branch of Federal Research Center 'Geophysical Survey of RAS' 3 Academician Koptyug ave., Novosibirsk 630090, Russia e-mail: solov@gs.nsc.ru
Чечельницкий Владимир Васильевич, канд. геол.-мин. наук, зам. директора по научной работе Байкальский филиал Федерального исследовательского центра «Единая геофизическая служба РАН» 664033, Иркутск, ул. Лермонтова, 128, Россия Тел.: (3952)427490; e-mail: chechel@crust.irk.ru
Chechelnitsky, Vladimir V., Candidate of Geology and Mineralogy, Deputy Director for Science Baikal Branch of Federal Research Center 'Geophysical Survey of RAS' 128 Lermontov street, Irkutsk 664033, Russia Tel.: (3952)427490; e-mail: chechel@crust.irk.ru
Сальников Александр Сергеевич, докт. геол.-мин. наук
Сибирский научно-исследовательский институт геологии, геофизики и минерального сырья 630091, Новосибирск, Красный проспект, 67, Россия e-mail: seispv@sniiggims.ru
Salnikov, Aleksander S., Doctor of Geology and Mineralogy Siberian Research Institute of Geology, Geophysics and Mineral Resources 67 Krasnyy ave., Novosibirsk 630091, Russia e-mail: seispv@sniiggims.ru
Селезнев Виктор Сергеевич, докт. геол.-мин. наук, директор
Сейсмологический филиал Федерального исследовательского центра «Единая геофизическая служба РАН» 630090, Новосибирск, просп. Академика Коптюга, 3, Россия e-mail: svs0428@mail.ru
Seleznev, Victor S., Doctor of Geology and Mineralogy, Director Seismological Branch of Federal Research Center 'Geophysical Survey of RAS' 3 Academician Koptyug ave., Novosibirsk 630090, Russia e-mail: svs0428@mail.ru
Лисейкин Алексей Владимирович, канд. геол.-мин. наук
Сейсмологический филиал Федерального исследовательского центра «Единая геофизическая служба РАН» 630090, Новосибирск, просп. Академика Коптюга, 3, Россия e-mail: lexik1979@mail.ru
Liseikin, Aleksei V., Candidate of Geology and Mineralogy Seismological Branch of Federal Research Center 'Geophysical Survey of RAS' 3 Academician Koptyug ave., Novosibirsk 630090, Russia e-mail: lexik1979@mail.ru
Галёва Наталья Александровна, геофизик
Алтае-Саянский филиал Федерального исследовательского центра «Единая геофизическая служба РАН» 630090, Новосибирск, просп. Академика Коптюга, 3, Россия e-mail: tatapelya@gmail.com
Galyova, Natalia A., geophysicist
Altai-Sayan Branch of Federal Research Center 'Geophysical Survey of RAS' 3 Academician Koptyug ave., Novosibirsk 630090, Russia e-mail: tatapelya@gmail.com
