CC BY
156
Сейсмология
Вестник ДВО РАН. 2013. № 3
УДК 550.3
Б.В. ЛЕВИН, Е.В. САСОРОВА, А.В. ДОМАНСКИЙ
Свойства «критических широт», вариации вращения и сейсмичность Земли
Проведен краткий анализ причин возникновения «критических широт» на медленно вращающейся планете, получена фундаментальная связь величины сжатия планеты с ее угловой скоростью вращения и даны основные понятия о слабых пульсациях сфероида. Показана бимодальность распределения сейсмичности Земли, которая указывает на несомненную связь сейсмической активности планеты с эффектом «критических широт».
Ключевые слова: неравномерность вращения Земли, бимодальное широтное распределение землетрясений, вариации угловой скорости вращения, критические широты, кинетическая энергия вращения.
Features of "critical latitudes", variations of rotation and seismicity of the Earth. B.W. LEVIN, E.V. SASOROVA, A.V. DOMANSKIY (Institute of Marine Geology and Geophysics, FEB RAS, Yuzhno-Sakhalinsk).
Short analysis of "critical latitudes" occurrence on slowly rotating planet was carried out. The basic connection of planet ellipticity with its angular rate of rotation was received, and the main concepts about spheroid weak pulsations were given. Bimodal distribution of the Earth seismicity which indicates clear connection of the planet seismic activity with the effect of "critical latitudes" was shown.
Key words: irregularity in rate of rotation of the Earth, bimodal latitudinal distribution of earthquakes, variations of angular speed of rotation, critical latitudes, kinetic energy of rotation.
В 2012 г. исполнилось 100 лет со дня опубликования в Париже трактата Алекса Веронне, в котором излагается точная теория вращения Земли [15]. Блестящие классические работы Маклорена, Клеро, Ньютона, Гюйгенса, посвященные выявлению основных характеристик вращающихся небесных тел, были проанализированы Пуанкаре [11], который тогда не мог отметить особенностей «критических широт». Эти широты, которые получили название критических в начале ХХ в., впервые отмечены французским математиком А. Веронне. Он же впервые привел вычисленное значение этих широт, равное ± 35°15'52", и указал на существование в этих стационарных зонах перехода от сжатия к растяжению пород Земли.
Обсуждению вопросов планетарной геологии и особенностям «критических широт» была посвящена переписка между В.И. Вернадским и Б.Л. Личковым в середине ХХ в. [4]. Влияние средних широт на геолого-геофизические процессы в твердой Земле отмечалось многими геологами и геофизиками и выделено в монографии В.Е. Хаина и Н.В. Коронов-ского [7]. Анализируя проблемы геофизики, О.Ю. Шмидт отметил, что вопрос, когда и почему стали происходить землетрясения, в конечном счете упирается в решение вопроса о происхождении Земли [8].
В настоящей работе приведен краткий анализ особенностей «критических широт» во вращающемся небесном теле, представлена связь величины сжатия тела с угловой
*ЛЕВИН Борис Вульфович - член-корреспондент РАН, директор (Институт морской геологии и геофизики ДВО РАН, Южно-Сахалинск), САСОРОВА Елена Васильевна - доктор физико-математических наук, главный научный сотрудник (Институт океанологии им. П.П. Ширшова РАН, Москва), ДОМАНСКИЙ Андрей Владимирович - доктор физико-математических наук, главный научный сотрудник (Институт морской геологии и геофизики ДВО РАН, Южно-Сахалинск). *Е-таП: [email protected]
Статья подготовлена по материалам доклада, прочитанного Б.В. Левиным на заседании Президиума ДВО РАН 25 октября 2012 г.
скоростью вращения, сформулированы понятия о слабых пульсациях геоида и продемонстрированы бимодальные распределения основных параметров сейсмичности Земли, указывающие на закономерную связь усилений сейсмоактивности планеты с эффектом критических широт.
Фигура Земли и «критические широты»
Фигура Земли под влиянием гравитационного и центробежного потенциала сформировалась как сфероид (геоид), обладающий сжатием: е = (Я - Н) / Я.
Здесь Я, Н - экваториальный и полярный радиусы планеты. Радиус-вектор г сфероида при малых сжатиях определяется формулой общего вида
г = Я0 [1 + е (1/3 - бш»], в которой Я0 - радиус сферы, эквивалентной по объему сфероиду, ф - геоцентрическая широта [13]. Из приведенного равенства видно, что радиус-вектор г = Я0 совпадает со средним радиусом планеты и в первом приближении не зависит от ее сжатия на «критических широтах», определяемых равенством Бшф = 1/3, ф1 = ± 35°15'52".
На рис. 1 дана схема эволюции сфероида и образования «критических широт».
2
2
2
л о В
Рис. 1. Схема вращающегося небесного тела: три фрагмента для трех различных величин угловой скорости (ю). а - ю0 ~ 0; б - ю> ю0; в - ю2> ю> ю0. Штриховые прямые линии обозначают критические широты для трех фрагментов вращающегося тела. Штриховые черные круги на фрагментах б и в определяют сферу, объем которой равен объему вращающихся эллипсоидов
Анализ фигур вращения [11] показал, что для всех рассмотренных подходов к определению формы поверхности планеты: эллипсоида вращения (Маклорен), самогравитирую-щего тела - постоянного геопотенциала (Ньютон) и упруго-деформируемой сферы [3] -сжатие не зависит от инерционных характеристик планеты и существует универсальная связь между сжатием и угловой скоростью вращения планеты е = С(Я0)ю2,
где С(Я0) - некоторая постоянная, зависящая от среднего радиуса планеты.
Полная (гауссова) кривизна поверхности эллипсоида вращения при малых сжатиях в первом приближении определяется уравнением
к = К1 [1 - 8еР2(БШф)/3], где к1 = 1/Я20 - полная кривизна сферы и Р^шф) = 1,5бш2ф - 0,5 - зональная сферическая функция второго порядка. Зависимости полной кривизны сфероида от широты при различных сжатиях е даны на рис. 2. Значения полной кривизны приведены к кривизне сферы. Видно, что общие точки пересечения пучка гауссовых кривых при различных малых сжатиях сфероида находятся на «критических широтах», в первом приближении не зависят от сжатия, а полная кривизна равна кривизне эквивалентной по объему сферы.
1,01 1,005
0,995 0,99 0,985
-90 -80 -70 -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60 70 80 90
Широта
Рис. 2. Зависимость полной (гауссовой) кривизны поверхности сфероида от широты для различных величин сжатия: 1. е = 0 - (сфера), 2. е = 0,00135; 3. е = 0,00335 (Земля); 4. е = 0,00535
Энергия вращения геоида и широтная зависимость объемной плотности кинетической энергии
Для выявления особенностей широтного распределения кинетической энергии введем функции:
1(8, ф) = ^Ш2ф + х12С0Б2ф)1/2, N(8, ф) = ^Ш2ф + Х6С082ф)1/2, X = (1 - е)1/3.
Согласно работе [5] можно показать, что элементы площадей поверхности эллипсоида вращения с и сферы равного объема ай относятся как
с / с0 = х4Ь / Щ4.
Точная формула для радиус-вектора эллипсоида вращения имеет вид
г = Я0 х2 / N.
Из двух последних формул следует формула для изменения элемента объема эллипсоида вращения при изменении радиус-вектора
АУ = сД |1 - х6 Ь / Щ / 3.
Кинетическая энергия вращения эллипсоида массы М равна
Е = МЯ02т2 / 5х2, АЕ = Е - Е0 = МЯ2 т2 (1 - х2) / 5х2 = Е0 (1 - х2) / х2.
Из двух последних формул следует, что
АЕ / АУ = 3Е0 (1 - х2) ■ Щ5 / (с0 Д0 х2 | Щ5 - хЬ|),
и относительная объемная плотность кинетической энергии равна
У0 ■ АЕ / (Е0 ■ АУ) = (1 - х2) Щ5 / (х2 | Щ5 - х6Ь|), У0 = с0Я0 / 3.
Зависимость относительной объемной плотности кинетической энергии вращения эллипсоида от широты приведена на рис. 3.
Отметим, что две критические зоны, в которых фиксируются пиковые значения плотности энергии вращения геоида, отделяют полярные области эллипсоида, где ярко выражена деформация сжатия среды, от экваториальной области, в которой превалируют деформации растяжения.
60 50 40 30 20 10 о
-90 -80 -70 -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60 70 80 90
Широта
Рис. 3. Зависимость относительной объемной плотности энергии вращения эллипсоида от широты
) V
Изменение угловой скорости вращения Земли и его вклад в вариацию энергии вращения
Наблюдаемая зависимость величины угловой скорости вращения от времени (т.е. нестабильность вращения) должна оказывать влияние на общую энергию вращения тела*. При выполнении условия ю = ю(£), для вариации кинетической энергии вращения тела получаем
ёЕ = 2Е ■ ёю / ю + У ■ ю2.
Таким образом, энергия тела в результате нестабильности вращения должна возрастать в результате прибавления величины угловой скорости и величины момента инерции тела. Наиболее существенным здесь является рост относительной величины угловой скорости вращения, которая оценивается, по наблюдениям [6], как ёю / ю ~ 10-8. В результате оценка величины вариации энергии вращающегося тела составит
ёЕ = 2Е ■ ёю / ю ~ 1029 ■ 10-8 = 1021 Дж в год.
Отметим, что эта величина по порядку совпадает с величиной ежегодно выделяемой энергии землетрясений на всей Земле. Таким образом, при ускорении вращения тела появляется дополнительная энергия, которая выделяется неравномерно по широтным поясам (рис. 3).
На рис. 4 приведен график изменения скорости вращения Земли по данным 1ЕЯ8 [9], из которого следует, что вариация угловой скорости в среднем составляет величину порядка 10-8. Рис. 5 показывает результаты спектрального анализа колебаний. В спектре наряду с годовыми и полугодовыми колебаниями присутствуют месячные и полумесячные вариации, связанные с движениями системы Земля-Луна.
Рис. 4. Среднемесячные значения величины относительной угловой скорости вращения Земли (V = 43200ю / п - 1) за последние 50 лет. На вертикальной оси показаны значения V • 1010
Рис. 5. Спектр ежемесячных средних значений величины V в двойном логарифмическом масштабе. Линии со стрелками определяют представительные пики спектра. По горизонтальной оси - количество циклов в год
Фридман А.М. Частное сообщение. 2007 г.
Материалы каталогов землетрясений и широтная зависимость сейсмичности Земли
Анализ глобальных закономерностей в широтном распределении землетрясений показал, что даже в эпоху формирования сейсмологии как науки отмечалась неоднородность в распределении эпицентров землетрясений на Земле. Концепция зависимости сейсмичности от широтного расположения событий сформировалась за последние два десятилетия [1, 2, 12, 14]. Для статистической обработки событий мы использовали данные из каталога ISC UK [10] с предварительной стандартизацией магнитуд и очисткой списка событий от афтершоков. Анализу подвергался массив событий с магнитудой M > 4,0 за период с 1964 по 2011 г. (всего более 250 000 событий). Сопоставление количества землетрясений в различных широтных поясах и плотности землетрясений (количество событий, отнесенное к единице длины границы литосферных плит в каждом широтном поясе) показало, что сейсмическая активность планеты практически отсутствует на полюсах и в полярных шапках, имеет четко выраженные пики в средних широтах Северного и Южного полушарий (40-50° N и 20-35° S) и стабильный локальный минимум возле экватора (10-20° N). Аналогичную форму широтного распределения демонстрирует плотность высвобожденной энергии событий (рис. 6). Совпадение особенностей распределения сейсмичности Земли и объемной плотности энергии вращения (рис. 3) представляется очевидным. Это свидетельствует в пользу наличия некоей физической закономерности, которая связывает особенности вращения Земли с проявлением сейсмичности планеты.
4<=M<4.5 = 4.5<=M<5 —5.0<=M<5.5 х 5.5<=M<6.0 ж 6.0<=M<6.5 —6.5<=M
СЛ, (Л, СЛ, СЛ, СЛ, СЛ, (Л, (Л, (/) _
I I I I I I I г II I I I I I I I I ООООООООоо'ОООООООО
---- . .. -----u М О Tf Ю
OOOI^iOlO'irOCN
(О МО О)
OOOOOOOOoo°°°tlOOOO
3.0E+20 2.5E+20 2.0E+20 1.5E+20 1.0E+20 5.0E+19 0.0E+00
(Л(/)1Л(Л(Л(Л(Л(Л(Л22г
J о1 о1 о1 о1 о1 о1 о1 о1 о1 О1 о1
o>cor^<0Ln^ro<4;i;=c4ro
ООООООООооООО Olvlfllfl^MN^ t-Nn
о'
LOU5 Г-»00 ОТ
ОО ОО О ^LO 00
Рис. 6. Распределение относительной плотности сейсмических событий для 6 магнитудных диапазонов (а) и плотности энергии землетрясений по широтным поясам (б). Черные штриховые линии соответствуют максимальным значениям наблюдаемых величин, серые - минимальным
Выводы
Наблюдаемые вариации скорости вращения Земли должны приводить к попеременным изменениям величины сжатия упругой оболочки геоида. Рост скорости вращения планеты вызывает увеличение общего сжатия (увеличение экваториального радиуса), а уменьшение скорости - снижение сжатия (и уменьшение экваториального радиуса). Учитывая особенности «критических широт», мы должны осознавать, что сжатие планеты сопровождается сжатием пород в полярных шапках и растяжением - в экваториальной зоне. Соответственно, уменьшение общего сжатия (при замедлении вращения) вызывает уменьшение сжатия пород в полярных шапках и увеличение сжатия в области экваториальной зоны. В результате должна возникать попеременная пульсация геоида относительно стационарных плоскостей, расположенных в области «критических широт». Очевидно, пульсации геоида должны сопровождаться изменениями площади поверхности фигуры планеты и, следовательно, вызывать деформации и напряжения в породах упругой оболочки геоида. Этот процесс по мере накопления поврежденности в среде должен завершаться землетрясениями.
Проверка сформулированной гипотезы возможна путем постановки спутниковых наблюдений. Серия ежесуточных измерений величины радиуса Земли в пяти точках (на экваторе, на полюсах и на широте ± 35°) с точностью порядка 200 м позволит получить необходимые материалы наблюдений.
ЛИТЕРАТУРА
1. Левин Б.В., Сасорова Е.В. Бимодальный характер широтных распределений землетрясений в Тихоокеанском регионе как проявление глобальной сейсмичности // ДАН. 2009. Т. 424, № 4. С. 538-542.
2. Левин Б.В., Сасорова Е.В. Сейсмичность Тихоокеанского региона: выявление глобальных закономерностей. M.: Янус-К, 2012. 308 с.
3. Лейбензон Л.С. Собрание трудов. Т. 4. M.: Изд-во АН СССР, 1955. 399 с.
4. Личков Б.Л. Природные воды Земли и литосфера. M.; Л.: Изд-во АН СССР, 1960. 163 с.
5. Рашевский П.К. Курс дифференциальной геометрии. 3-е изд. M.: ГИТТЛ, 1950. 428 с.
6. Сидоренков Н.С. Физика нестабильностей вращения Земли. M.: Наука, Физматлит, 2002. 384 с.
7. Хаин В.Е., Короновский Н.В. Планета Земля от ядра до ионосферы. M.: КДУ, 2007. 244 с.
8. Шмидт О.Ю. Четыре лекции о теории происхождения Земли и планет. M.: Изд-во АН СССР, 1957. 140 с.
9. IERS Annual Report. Intern. Earth Rot. Service. Observ. de Paris. 2000. 144 p. - www.iers.org
10. International Seismological Catalogue (ISC). - http://www.isc.ac.uk
11. Poincaré H. Figures d'équilibre: d'une masse fluide. Paris: C. Naud, 1902. 210 p.
12. Riguzzi F., Panza G., Varga P., Doglioni C. Can Earth's rotation and tidal despinning drive plate tectonics? // Tectonophysics. 2010. Vol. 484. P. 60-73.
13. Stacey F.D., Davis P.M. Physics of the Earth. N.Y.: Cambridge Univ. Press, 2008. 552 p.
14. Sun W. Seismic energy distribution in latitude and a possible tidal stress // Physics of the Earth and Planetary Interiors. 1992. Vol. 71. P. 205-216.
15. Véronnet A. Rotation de l'ellipsoïde hétérogène et figure exacte de la Terre // J. mathématiques pures et appliquées. 1912. Série 6e, t. 8. P. 331-463.
