Одномерная референтная модель Земли.
(1D REM).

Одномерная скоростная модель Земли – это кусочно непрерывная функция скорости Р или S волны, зависящей от глубины от поверхности Земли или радиуса исчисляемого от центра Земли к ее поверхности. Функция состоит из непрерывных частей, которые отделяются друг от друга разрывами и переходными зонами. Набор непрерывных и разрывных участков функции определяется детальностью проводимых исследований и применяемыми технологиями. Сейсмология изучает сейсмические волны во вполне определенном диапазоне частот. Это и определяет конструкцию функции, представляющей из себя Референтную модель Земли. На модели выделяются три типа особенностей. Непрерывные области, представляющие из себя области с постоянными значениями геофизического параметра. Это внешнее и внутренней ядро Земли для скорости поперечных волн и внутренне ядро для скорости продольных волн. Относительно слабо градиентные зоны. Это внешнее ядро для скорости поперечных волн и Нижняя мантия для скоростей P и S волн. Высоко градиентная зона и для продольных, и для поперечных волн – это переходная зона в низах Верхней мантии. Верхняя мантия над переходной зоной – это относительно высоко градиентная зона. Разрывы функции – это границы, отделяющие Земную кору от Верхней мантии, переходную зону от Верхней и Нижней мантии, Нижнюю мантию от внешнего ядра и внешнее ядро от внутреннего ядра. Все выше названные особенности были обнаружены в наблюдаемых геофизических данных. Мы уже говорили о том, что уменьшение скорости приводит к разрыву годографа рефрагированной волны, который часто совпадает с годографом первых вступлений. Разрыв виден на годографе первых вступлений как участок с пониженным значением кажущейся скорости. Положительная градиентная зона в скоростном разрезе приводит к появлению петель на годографе рефрагированной. А положительные скачки скорости приводят к излому годографа первых вступлений.
Хронологически, одной из наиболее ранних современных моделей Земли является модель Джеффриса – Буллена [Jeffreys, H., and K.E. Bullen, 1940 ]. Эта модель послужила надежным фундаментом для всех самых современных моделей. В процессе накопления данных об объемных телесейсмических волнах, выявились регулярные отличия во временах вступления этих волн от теоретических времен, предсказываемых моделью Джеффриса – Буллена. В результате были уточнены годографы телесейсмических P и S волн, и, соответственно одномерная референтная модель Земли. А сами стандартные годографы Р и S волн получили название Стандартные годографы Херрина [Herrin, 1968]. Для работы с телесейсмическими объемными волнами эти годографы не утратили актуальности до сегодняшних дней. Годограф Херрина представляет из себя детальную таблицу численных значений производных годографов P и S объемных волн в диапазоне эпицентральных расстояний от 00 до 1000 , с шагом 0.50 , для очагов, расположенных на глубинах от 0 до 800 километров (0 км, 15 км, 40 км, 50 км, 75 км, 100 км, 125 км, 150 – 800 км с шагом 50 км.). Эта форма представления стандартного годографа удобна для проведения расчетов таких промежуточных параметров, как "Невязка" и "Разность", что будет показано ниже.
В 1981 году была опубликована одномерная референтная модель PREM (Preliminary Reference Earth Model) [Dziewonski, A. M., and D. L. Anderson, 1981]. Модель была разработана по заказу Международного Союза по Геодезии и Геофизике (International Union of Geodesy and Geophysics). Союз заказал сферически симметричную модель, которую можно было бы использовать в геодезических и геофизических исследованиях. Модель учитывает все достижения в измерениях и интерпретации свободных колебаний Земли на момент ее создания. Эта модель – еще одна попытка лучше понять и уточнить одномерную референтную модели Земли. Было привлечено большое количество новых данных по кинематическим и динамическим характеристикам объемных и поверхностных сейсмических волн. Для лучшего согласования данных было учтено неупругое рассеяние. Введение анизотропии скоростей (2-4%) в верхних 220 километрах Верхней мантии позволило получить согласованную с исходными данными одномерную референтную модель Земли без обязательного, как предполагалось ранее, слоя с пониженным значением скорости в верхней мантии Земли, то есть глобального волновода. Этот вывод, как нам представляется, был наиболее существенным с точки зрения понимания внутреннего устройства Земли в свете новой глобальной тектоники плит. В соответствии с базовыми гипотезами, лежащими в основе Тектоники плит, глобальные плиты плавали и перемещались на относительно размягченном субстрате, располагающемся в верхней части Верхней мантии Земли. Этот сферический слой и должен был проявить себя, как слой с пониженной скоростью. Но структурно - сейсмологические исследования Верхней мантии Земли не выявляли этого, предположительно повсеместно развитого волновода. Во многом создание модель PREM было призвано объяснить причину подобного противоречия. Модель послужила эталоном для ряда трехмерных моделей Земли. Последующие измерения собственных частот Земли дают систематические расхождения с моделью. Имеются и другие несогласия с вновь получаемыми данными. Например 220 километровая граница не нашла подтверждения при анализе отраженных волн [ Shearer P.M., 1991].
В 1991 году была представлена одномерная скоростная модель IASP91 [ Kennett, B. L. N., E. R. Engdahl E. R., 1991]. Модель разрабатывалась в течение трех лет специальной подкомиссией по землетрясениям Ассоциации по Сейсмологии и Физике Земных Недр (IASPEI). Основная цель работ – создание новых глобальных таблиц годографов сейсмических фаз, которые обновят стандартный годограф Джеффриса – Буллена (1940) и модель PREM(1981). .
Результаты еще однй попытки обобщения данных с целью построения референтной модели Земли были представлены в 1993 году [Morelli, A., Dziewonski A. M., 1993]. Модель получила индекс SP6. Следуя руководству, разработанному рабочей группой IASP91 были изучены годографы объемных волн для нескольких фаз сейсмических волн, наиболее часто приводимых в бюллетенях МСЦ и была предложена модель скоростей Р и S волн, соответствующая этим данным. Данные по примерно 16000 хорошо записанных мелкофокусных землетрясений за 24 года (1964 – 1987) были обращены в одномерную скоростную модель с учетом поправок за горизонтальную скоростную неоднородность. Используя полученную скоростную модель, были рассчитаны обобщенные годографы P, PcP-P, PKPAB, PKPBC, PKIKP и S волн для каждого диапазона эпицентральных расстояний. Ограничения скорость для верхов внешнего ядра получены по данным SKS и SKKS-SKS [Hales A. L., Roberts J. L., 1970, 1971].
Для Верхней мантии результаты модели SP6 отличаются от модели IASP91 лишь немного более повышенным значением скоростей Р и S волн между 410 и 660 километрами. В Нижней мантии было получено более низкое значение градиента скорости. Что совпадает с моделью PREM. В верхах внешнего ядра скорость Р-волн меньше, чем в PREM. Скачок скорости на границе внутреннего ядра уменьшен до 0.62 км/с. Новый радиус внутреннего ядра – 1215.00 км. Рассчитанные по модели SP6 годографы заметно отличаются от IASP91, хотя в целом менее чем на 1.0 секунду. Для SKS и SKKS различия больше.
В 1995 году была представлена скоростная модель AK135 [Kennett, B. L. N., Engdahl, E. R., Buland R., 1995]. Эта модель обеспечила значительно лучшее приближение для большого количества сейсмических фаз, чем это было для моделей SP6 и IASP91. Различия между моделью AK135 и моделями SP6 и IASP91 в целом незначительные, кроме границы Внутреннего ядра Земли. Для этой границы был уменьшен градиент скорости, что привело к достижению удовлетворительных результатов для дифференциала годографа РКР фазы. При построении таблиц AK135 использованы экспериментальные годографы от землетрясений, для которых были переопределено положение очагов с использованием таблиц IASP91. То есть, таблицы AK135 можно рассматривать как следующую итерацию (после IASP91) расчета теоретических годографов. Скоростная модель AK135, как и модель IASP91 радиально стратифицированная. Теоретические годографы Р – волн очень близки к таковым, рассчитанным для модели IASP91. Более значимы отличия для годографа S волн и в частности для фаз, проникших в ядро Земли.
В свою очередь, таблицы AK135 были использованы для систематического пересчета очагов землетрясений [ Engdahl, E.R.,van der Hilst, R., Buland, R., 1998] и последующего обновления каталога этих землетрясений.
В настоящее время ведутся работы по уточнению одномерной референтной модели Земли. Некоторые важные особенности сейсмического поля, выявленные в результате обработки сейсмограмм землетрясений и больших взрывов, в настоящее время находятся в стадии осмысления. Анализируется достоверность таких особенностей, как граница 220, граница 410, граница 520, граница 660. Где находятся границы 410 и 660, насколько они изменчивы. Является граница 520 повсеместной, глобальной. Возможно ли, чтобы природа границы 220 была обусловлена анизотропией сейсмических скоростей, а не с изменением скорости с глубиной. Имеются ли глобальные границы в Нижней мантии. Еще предстоит решить нелегкие вопросы, связанные с учетом сферичности Земли при выборе ее осредненной структуры, степень детальности модели; как учесть анизотропию скорости; совместное, согласованное использование данных по Р и S волнам, пересчет данных структурной сейсмологии в плотностную модель Земли. Насколько допустимо использовать единую одномерную референтную модель Земли при реконструкции внутренней структуры Земли под континентами и под океанами, платформами и складчатыми областями и другими крупнейшими геотектоническими образованиями. Какую модель использовать при реконструкции. Как согласовывать построения, выполненные на базе региональных одномерных референтных моделей.
Литература.
Jeffreys, H., and K.E. Bullen, 1940. Seismological Tables. London: British Association for the Advancement of Science.
Hales, A. L., Roberts, J. L., 1970. The travel times of S and SKS. Bull. Seism. Soc. Am., v. 60, p. 461
Hales A. L., Roberts J. L., 1971. The velocities in the outer core. Bulletin of the Seismological Society of America; August 1971; v. 61; no. 4; p. 1051-1059
Herrin, E., 1968. Seismological tables for P-phases, Bull. Seism. Soc. Am., 60, 461-489. Dziewonski, A. M., and D. L. Anderson, 1981. Preliminary reference Earth model, Phys. Earth Planet. Inter., 25, 297–356.
Shearer P.M., 1991. Constrains on Upper Mantle Discontinuities From Observation of Long-Period Reflected and Converted Phases. Journal of Geophysical Research, vol.96, No B11, p. 18,147 – 18,182, October 10, 1991.
Kennett, B. L. N., E. R. Engdahl E. R., 1991. Traveltimes for global earthquake location and phase identification. Geophysical Journal International 122, 429–465.
Morelli, A., Dziewonski A. M., 1993. Body-wave traveltimes and a spherically symmetric P- and S-wave velocity model. Geophysical Journal International 112, 178–184.
Kennett, B. L. N., Engdahl, E. R., Buland R., 1995. Constraints on seismic velocities in the Earth from traveltimes. Geophysical Journal International 122, 108–124.
Engdahl, E.R.,van der Hilst, R., Buland, R., 1998. Global teleseismic earthquake relocation with improved travel times and procedures for depth determination. Bull. Seismol. Soc. Am., v. 88, 722-743.