Одной из актуальных задач современной геофизики является изучение плотностного строения земной коры и подстилающей ее мантии. В этом изучении гравитационное моделирование в комплексе с сейсмической томографией представляется важнейшим инструментом определения плотностных неоднородностей. Данные неоднородности могут оказывать существенное влияние на тектонику регионов, вызывая напряжения в коре и мантии [Pogorelov, Baranov, 2010]. Это наглядно показывают активные зоны сочленения окраин материков с субдуцирующими под них океаническими плитами, образующими напряжения по всей глубине тектоносферы. Возникающие при этом процессы приводят к появлению значительных горизонтальных и вертикальных напряжений, что нередко заканчивается катастрофами глобального масштаба. Все вышеперечисленное указывает на необходимость изучения глобальных тектонических процессов в тектоносфере Земли, в котором гравитационное моделирование в комплексе с сейсморазведкой играет определяющую роль.

Данные предыдущих исследователей. Впервые распределение мантийных плотностных неоднородностей в Центральной Азии было рассчитано И.Л. Нерсесовым с соавторами [Nersesov et al., 1975]. Однако в силу недостатка данных точность полученных в то время расчетов была невелика. Вместе с тем значительные результаты были получены в изучении строения земной коры. Благодаря многочисленным проектам по ее исследованию геофизическими методами, степень изученности коры достигла такого уровня, что появились региональные модели, дающие вполне адекватную информацию о глубинном расположении ее внешних и внутренних границ, а также о распределении сейсмических скоростей и плотности. Это, в свою очередь, дает нам возможность проводить изучение мантийных плотностных неоднородностей со значительно большей точностью.

В настоящей работе представлены результаты расчетов гравитационных аномалий коры и верхней мантии Центральной и Южной Азии. Этот регион хорошо подходит для решения данной задачи, так как он охватывает области всего диапазона мощности, возраста и способов образования коры [Christensen, Mooney, 1995]. Для выполнения этой задачи использовалась оригинальная программа трехмерного гравитационного моделирования, учитывающая сферическую форму поверхности Земли, разработанная В.Н. Сеначиным [Senachin, 2015a, 2015b]. Расчет гравитационных аномалий выполнялся на основе данных цифровой модели земной коры AsCrust, разработанной А.А. Барановым [Baranov, 2010].

Изучаемый регион включает в себя Альпийско-Гималайский складчатый пояс, область тройного сочленения рифтовых зон на севере Африки и область окраинных морей Юго-Восточной Азии, обрамляющихся глубоководными желобами с сопутствующими таким зонам вулканическими поясами. Рельеф земной поверхности и морского дна включает самые высокие горы (Гималаи) и одни из самых глубоких областей на нашей планете – глубоководные желоба Индонезии. В данном регионе предполагается столкновение Индийской плиты с Азиатской, которое могло привести к надвигу и смятию края Азиатской плиты и к большому увеличению мощности континентальной коры [Oreshin et al., 2011]. Этот процесс, возможно, сопровождается деламинацией (отделением корового слоя Индийской литосферной плиты от ее мантийной «подушки», механизм которого до конца не изучен) [Jiménez-Munt et al., 2008; Krystopowicz, Currie, 2013; Ueda et al., 2012] (рис. 1).

Цифровая модель земной коры AsCrust: карта глубины до Мохо. При построении модели AsCrust было проанализировано большое количество новых данных по отраженным, преломленным и поверхностным волнам от землетрясений и взрывов. Все они были интегрированы в единую модель с разрешением 1×1°. Результаты были представлены в виде десяти цифровых карт, определяющих глубину до границы Мохо, мощность верхней, средней и нижней части консолидированной коры, а также скорость продольных волн в этих слоях [Baranov, 2010]. Для расчета гравитационных аномалий в настоящей работе использовались данные о мощности слоев коры и их плотности в каждой точке градусной сетки. Плотности в слоях были рассчитаны по скоростям продольных волн с использованием формулы из [Brocher, 2005] (рис. 2).

Алгоритм расчета гравитационных аномалий. В гравитационном моделировании при изучении крупных региональных объектов возникает необходимость учитывать кривизну земной поверхности. Для этой цели разрабатываются алгоритмы вычисления гравитационного поля от тел, ограниченных сферическими поверхностями [Kosygin et al., 1996; Starostenko et al., 1986; Strakhov et al., 1989; Jones et al., 2010; Li et al., 2011; Schmidt et al., 2007; и др.]. В настоящей работе для этой цели был использован алгоритм, основанный на применении прямых формул расчета гравитационного эффекта, которые могут быть получены для точки, находящейся на полюсе сферы. Использование прямых формул вычисления во многом упрощает алгоритм, но требует постоянного перерасчета системы координат для каждой точки расчета, что несколько усложняет задачу (рис. 3).

Исходные данные и методика расчета гравитационных аномалий. Количество слоев в нашей расчетной модели равно семи: это водный слой, три слоя осадков (глубины границ и плотности осадочных слоев) из модели [Laske, Masters, 1997] и три слоя коры (глубины границ и плотности в коре, полученные из Vp скоростей) из модели AsCrust [Baranov, 2010] для покрываемой ею территории. Для окружающих регионов данные о строении и свойствах коры брались из модели CRUST 2.0 [Bassin et al., 2000] и были проинтерполированы на сетку 1´1°. Таким образом, для осадков и коры были использованы данные с разрешением 1´1°, тогда как для водного слоя (батиметрия) разрешение составило 0.1´0.1°.

В качестве гравитационного поля Земли была взята модель GGM01, полученная на основе спутниковых наблюдений в ходе выполнения проекта «GRACE» (http://www.csr.utexas.edu). На основе этой модели были рассчитаны аномалии в «свободном» воздухе на всей поверхности нашей модели, которые учитывают поправку за высоту объекта наблюдений: гравитационное поле изменяется в пределах от –250 до +260 мГал. Зона столкновения Индийской плиты с Азиатской выделяется наличием параллельных узких аномалий разного знака, достигающих 200 мГал и более: юго-западная аномальная зона, характеризующаяся отрицательными аномалиями, соответствует, по-видимому, границе сочленения двух плит, где Индийская плита пододвигается под Азиатскую согласно [He et al., 2010; Oreshin et al., 2011; и др.]. В гравитационном поле региона также видно довольно четкое разделение Тибета и примыкающей к нему с северо-востока Таримской плиты (отрицательная аномалия до –150 мГал), границы которой оконтуриваются узкими зонами положительных аномалий. Южная часть Каспийского моря также характеризуется отрицательной аномалией до –150 мГал, тогда как Тянь-Шань – узкой полосой положительных аномалий до 110 мГал. Однако на большей части исследуемого региона поле близко к нормальному, изменяясь в пределах нескольких десятков миллигал. Умеренно повышенным уровнем гравитационного поля (в пределах 40–80 мГал) выделяется остальная часть Альпийско-Гималайского складчатого пояса. В морях юго-восточной окраины Азии отмечается слабое повышение уровня поля в акваториях морей и две узкие зоны высокоамплитудных аномалий разного знака (до 200 мГал), которые порождаются изостатически нескомпенсированными системами островных дуг и желобов окраинных морей Юго-Восточной Азии (рис. 6).

Гравитационный эффект земной коры. Гравитационный эффект земной коры, рассчитанный в Азиатском регионе, показывает наличие крупных аномалий, изменяющихся в пределах 940 мГал (от –380 до +560 мГал). Максимальная положительная аномалия наблюдается в районе африканского тройного сочленения рифтовых зон, где она достигает положительного максимума в регионе +560 мГал. Положительные аномалии также наблюдаются в Таримском бассейне (до +130 мГал), Юго-Восточном Китае (до +100 мГал), в районе Иранского нагорья (до +180 мГал) и в задуговых районах субдукционных зон Индийской и Тихоокеанской плит (до +290 мГал). Крупные отрицательные аномалии соответствуют районам Каспийского и Черного морей (до –380 мГал), Гималаев (до –280 мГал) и восточной части Тибета (до –330 мГал). Восточное Средиземноморье характеризуется отрицательной аномалией до –310 мГал.

Восточная часть Аравийского полуострова и Месопотамская низменность характеризуются отрицательными аномалиями до –220 мГал. Также на карте расчетных аномалий коры выделяются подводные хребты (до +280 мГал), тянущиеся субмеридионально с юга на север, происхождение которых предположительно определяют как след «горячих точек, прожигающих проходящую над ней литосферную плиту (рис. 7).

Мантийные гравитационные аномалии. Мантийные гравитационные аномалии изучаемого региона были получены вычитанием гравитационного эффекта земной коры из наблюденного гравитационного поля. Они изменяются в пределах от –570 до +350 мГал, что примерно вдвое превышает пределы изменения этого поля. Это прямо указывает на наличие крупных плотностных неоднородностей в литосферной мантии, которые должны скомпенсировать коровые аномальные массы. Наиболее крупные положительные мантийные плотностные неоднородности в изучаемом регионе наблюдаются в районах узкой полосы Гималаев (до +330 мГал) и Восточного Тибета (до +350 мГал). В Каспийском и Черном морях они достигают +250 и +300 мГал соответственно. Восточное Средиземноморье характеризуется положительной аномалией до +280 мГал. Восточная часть Аравийского полуострова и Месопотамская низменность характеризуются положительными аномалиями до +220 мГал. Отрицательные аномалии наблюдаются в районах Таримского бассейна (до –190 мГал), над подводными хребтами в Индийском океане (до –340 мГал), в юго-восточной части Китая (до –120 мГал), в центральной части Индостана (до –80 мГал), а также в районах Гиндукуша и Каракорума (до –150 мГал). Субдукционные зоны Индийской и Тихоокеанской плит также характеризуются отрицательными аномалиями до –250 мГал. На северо-востоке Африканского континента, где имеет место тройное сочленение рифтовых зон (Красное море, Аденский залив, Африканский рифт), наблюдается область максимальных отрицательных мантийных аномалий. Здесь гравитационное поле понижается до –570 мГал (рис. 8).

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