Вариации силы тяжести и современная геодинамика юго-западной части Байкальского региона

GEODYNAMICS & TECTONOPHYSICS

PUBLISHED BY THE INSTITUTE OF THE EARTH’S CRUST SIBERIAN BRANCH OF RUSSIAN ACADEMY OF SCIENCES

2013 VOLUME 4 ISSUE 2 PAGES 119-134

ISSN 2078-502X

http://dx.doi.org/10.5800/GT-2013-4-2-0094

Gravity variations and recent geodynamics of the

SOUTH-WESTERN PART OF THE BAIKAL REGION

V. Yu. Timofeev1, Ye. N. Kalish2, Yu. F. Stus'2, D. G. Ardyukov1,

G. P. Arnautov2, M. G. Smirnov2, A. V. Timofeev1, D. A. Nosov2,

I. S. Sizikov2, E. V. Boiko1, E. I. Gribanova3

1 A.A. Trofimuk Institute of Petroleum Geology and Geophysics SB RAS, Novosibirsk, Russia

2 Institute of Automation and Electrometry SB RAS, Novosibirsk, Russia

3 Siberian Division of Geophysic Survey SB RAS, Novosibirsk, Russia

Abstract: Modern methods for determination of gravity values make it possible to obtain measurements with the accuracy up to 10-9 from g0 of the normal value (up to 1 microgal = 10 m/sec2). While all the systematic and periodic effects are excluded, a question is raised about stability of the gravity field of the Earth over time. Changes of the altitude (the Earth's radius) with time can be estimated with an accuracy of 0.1 mm by modern space geodetic techniques, such as VLBI method. Our experiments for evaluation of stability of the gravity values over the past decades are based on the data obtained by Russian and foreign observatories using absolute ballistic laser gravimeters. The results put a limit of 10-10 per year to changes of the Earth’s radius. These estimations can be useful for testing hypotheses in tectonics.

Measurements of non-tidal variations of gravity (Ag), which were obtained from 1992 to 2012 at the Talaya seismic station (located in the south-western part of the Baikal region), are interpreted together with GPS observation data. At the Talaya seismic station, the linear component of gravity variations corresponds to changes in the elevation of this site. The correlation coefficient is close to the normal value of the vertical gradient of gravity. At this site, coseismic gravity variations at the time of the Kultuk earthquake (27 August 2008, Mw=6.3) were caused by a combined effect of the change of the site’s elevation and deformation of the crust. Our estimations of the coseismic effects are consistent with results obtained by modeling based on the available seismic data.

Key words: monitoring of gravity variations, changes in the Earth's radius, rift zone, GPS-monitoring, earthquake.

Citation: Timofeev V.Yu., Kalish Ye.N., Stus' Yu.F., Ardyukov D.G., Arnautov G.P., Smirnov M.G., Timofeev A.V., Nosov D.A., Sizikov I.S., Boiko E.V., Gribanova E.I. 2013. Gravity variations and recent geodynamics of the south-western part of Baikal region. Geodynamics & Tectonophysics 4 (2), 119-134. doi:10.5800/GT-2013-4-2-0094.

Вариации силы тяжести и современная геодинамика

ЮГО-ЗАПАДНОЙ ЧАСТИ БАЙКАЛЬСКОГО РЕГИОНА

В. Ю. Тимофеев1, Е. Н. Калиш2, Ю. Ф. Стусь2, Д. Г. Ардюков1,

Г. П. Арнаутов2, М. Г. Смирнов2, А. В. Тимофеев1, Д. А. Носов2,

И. С. Сизиков2, Е. В. Бойко1, Е. И. Грибанова

1 Институт нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН, Новосибирск, Россия

2 Институт автоматики и электрометрии СО РАН, Новосибирск, Россия

3 Сибирский филиал Геофизической службы СО РАН, Новосибирск, Россия

119

Recent Geodynamics

V.Yu. Timofeev et al.: Gravity variations and recent geodynamics...

Аннотация: Современные методы определения значения силы тяжести позволяют проводить измерения с точностью до 10-9 от g0 нормального значения (до 1 микрогала = 10 нм/с2). При этом исключаются все систематические и периодические эффекты и возникает вопрос о стабильности поля силы тяжести Земли во времени. Оценить изменения высоты (радиуса Земли) во времени с точностью до 0.1 мм позволяют современные методы космической геодезии (VLBI метод). Экспериментальные оценки стабильности значения силы тяжести за последние десятилетия сделаны по материалам отечественных и зарубежных обсерваторий, использующих абсолютные лазерные баллистические гравиметры. Полученные результаты ограничивают изменение радиуса Земли значением 10-10 в год. Эти оценки можно использовать для тестирования тектонических гипотез.

Результаты измерений неприливных вариаций ускорения силы тяжести Ag, проведенных в 1992-2012 гг. на сейсмостанции «Талая» (юго-западная часть Байкальского региона), интерпретируются совместно с данными GPS-наблюдений. Линейная составляющая вариации силы тяжести на станции Талая соответствует изменениям высоты пункта. Коэффициент корреляции близок к нормальному значению вертикального градиента силы тяжести. Косейс-мические вариации силы тяжести на этом пункте в эпоху Култукского землетрясения (27.08.2008 г., Mw=6.3) вызваны комплексным эффектом изменения высоты пункта и деформации земной коры. Оценки косейсмических эффектов соответствуют результатам моделирования на основе сейсмологических данных.

Ключевые слова: мониторинг вариаций силы тяжести, изменение радиуса Земли, рифтовая зона, GPS-наблюдения, землетрясения.

1. Введение

Гравиметрические наблюдения являются важной частью комплексных геофизических измерений. Высокоточные измерения позволяют оценить стабильность значения силы тяжести и значения высоты пункта во времени - рассмотреть возможность изменения радиуса Земли, т.е. тестировать тектонические гипотезы о пульсирующей или расширяющейся Земле [Elhiny et al., 1978; Kukol, 1987]. Определение скорости постледниковых поднятий методами геодезии и гравиметрии позволяет дать оценки вязким параметрам геологической среды. Вместе с тем физически очевидная связь вертикальных смещений и изменений силы тяжести до последнего времени не находила ясного отражения в экспериментальных результатах, полученных в зонах высокой сейсмичности. В последние десятилетия с повышением точности гравиметрических измерений и появлением методов космической геодезии появилась реальная возможность прояснить этот вопрос. Важным дополнением этих работ являются измерения локальной деформации поверхности. Наиболее интересные результаты могут быть получены при проведении многолетних наблюдений, включая периоды сильных землетрясений со значительными косейсмическими изменениями. При этом, помимо данных о долговременных региональных процессах, они дают дополнительную информацию о положении и глубине разрыва при землетрясении, о мощности земной коры, охваченной деформациями. Исследования проводились в Байкальской рифтовой зоне (БРЗ), где наблюдается активный сейсмический процесс и регистрируется более 2000 землетрясений в год, при этом в центре БРЗ в основном регистрируются землетрясения со значительной вертикальной компонентой смещений. Измерения проводились отечественным аб-

солютным гравиметром ГАБЛ [Arnautov et al., 1994], геодезическими двухчастотными приемниками TRIMBLE-4700, лазерными и штанговыми деформографами и приливными наклономерами.

2. Оценка стабильности силы тяжести Земли

Для описания измерений силы тяжести g во времени можно воспользоваться простой моделью:

g

где r - вектор положения, t - время и Z - погрешность измерений.

Как известно, сезонная компонента вносит значимый вклад в изменения вертикальных смещений и силы тяжести на пункте. Изучение с помощью GPS-приемников, SG гравиметра и абсолютного гравиметра FG5 сезонной составляющей и тренда на итальянской обсерватории (Medicina, 1992-2002) показало, что отношение вертикальных смещений и g дает значение коэффициента -3 мм/мкгал, т.е. соответствует нормальному вертикальному градиенту у силы тяжести (рис. 1).

Сезонные вариации, связанные с изменением уровня грунтовых вод по данным французской обсерватории (Strasbourg, гравиметр FG5, 1997-2004) достигают 5 микрогал, а в целом по сети европейских станций со стационарными криогенными гравиметрами амплитуда сезонных вариаций составляет 3-8 микрогал [Hin-derer, Crossley, 2006]. Влияние гидрологических факторов сложно рассчитать, но их можно минимизировать, при размещении постамента для измерений на коренных породах [Tsuboi, 1982].

Следующие зависимые от времени факторы - это

120

Geodynamics & Tectonophysics 2013 Volume 4 Issue 2 Pages 119-134

a

b

Рис. 1. Результаты измерений вертикальных смещений (в метрах) с помощью GPS-приемников (а) и силы тяжести (измерения SG гравиметром на эффекте сверхпроводимости - сплошная линия и абсолютным гравиметром FG5) за 6 лет на итальянской обсерватории Medicina (05.1996-01.2002) в микрогалах (b).

Система СИ: ускорение силы тяжести, нормальное значение g0~9.8 м/с , 1-10 м/с =1 нм/с . Внесистемные единицы: 1 см/с =1 гал; 10_3 гал = 1 миллигал (mGal); 10-6 гал = 1 микрогал (pGal); 10_9 гал = 1 наногал; 1 микрогал = 10 нм/с2 ~ 10_9 g0, т.е. наноуровень эффектов.

Fig. 1. Vertical displacement measurements (in meters) by GPS receivers (a) and gravity measurements by an SG gravity meter based on the effect of superconductivity (shown as the solid line) and an FG5 absolute gravimeter (in microgal) for a period of six years, from May 1996 to January 2002; the surveys were conducted at the Medicina observatory, Italy.

SI system: gravity acceleration, the normal value g0 ~ 9.8 mps2, 1-10 9 mps2 ligal (mGal); 10-6 gal = 1 microgal (pGal); 10_9 gal = 1 nanogal. 1 microgal

: 1 nmps2. Units out of the system: 1 cmps2 = 1 gal, 10 3 gal = 1 mil-10 nmps2 ~ 10_9 g0, i.e. nanolevel of effects.

121

V.Yu. Timofeev et al.: Gravity variations and recent geodynamics...

| Рис. 2. Положение комплексных геофизических обсерваторий. I Fig- 2. Locations of integrated geophysical observatories.

эффекты земного прилива и полюса. Они рассчитываются достаточно точно. Так, для определения приливных параметров и выбора модели были проведены измерения с гравиметрами ЛаКоста-Ромберга по широтному профилю, простирающемуся от севера Франции до острова Сахалин [Ducarme et al., 2008]. Поправки за движения полюса Земли вводятся с использованием данных международного центра по вращению Земли, полученных экспериментально.

Изучение вариаций g во времени требует оценки стабильности значения силы тяжести Земли и ее высот. С учетом наличия скорости при линейной связи получаем следующее соотношение:

g

где g’ - скорость изменений и t0 - начальное время.

Связь вариаций g и высоты Z без учета деформации среды можно описать как:

Ag

где у - вертикальный градиент g и AZ(t) - изменение высоты пункта.

Для высот сложно говорить о глобальных вариациях, так как мировая система измерений находится в развитии, практически каждые пять лет с появлением новых станций происходит небольшое изменение глобальной референц-системы [Altamimi et al., 2007] и

уточнение геофизических моделей, используемых при построении и развитии глобальной отсчетной основы [Kaftan, Tsyba, 2009]. По оценкам отдельных авторов, вариации незначительно превышают ошибку; например, из обобщения по сети VLBI получена оценка 0.11±0.05 мм/год [Kolomiets et al., 2010].

Какова величина скорости g’ в соотношении (2)? Рассмотрим экспериментальные результаты, полученные на обсерваториях комплексных геофизических измерений - методами GPS-геодезии и абсолютной гравиметрии. Мониторинговые измерения с гравиметром ГАБЛ начаты на обсерватории Ключи (Новосибирск) в 1975-1977 гг. (рис. 2). Абсолютный лазерный баллистический гравиметр ГАБЛ неоднократно участвовал в международных сверках приборов этого класса [Robertson et al., 2001].

GPS-наблюдения на обсерватории Ключи начаты в 2000 г. Полученные за период 1977-2009 гг. на обсерватории Ключи значения силы тяжести приведены на рис. 3. Вариации на периоде наблюдений не превысили ошибку повторных измерений (0.5-2.0 мкгал). Стабильность значения силы тяжести за период 19772009 гг. составила 0.06 мкгал/ год.

Измерения с гравиметром ГАБЛ на астроплощадке в г. Иркутске проводились с 1995 г. с момента организации GPS наблюдений на пункте (рис. 3). Результаты измерений (микрогальная часть значения) с 1995 г. по 2006 г. приведены на рис. 4. Отклонение от среднего значения не превышает 3 мкгал. Стабильность

Geodynamics & Tectonophysics 2013 Volume 4 Issue 2 Pages 119-134

Рис. 3. Астрономическая обсерватория Иркутск, данные GPS-измерений в двух точках за период 1995-2010 гг. и 20022010 гг.

Геоцентрическое решение, скорости смещений С-Ю, В-З и вертикальные в сантиметрах. Сезонная компонента на вертикальной компоненте связана с климатическими вариациями, скорости +0.370±0.011 и -0.357±0.021 мм/год.

Fig. 3. The Irkutsk Astronomical Observatory; GPS measurements at two points for periods from 1995 to 2010 and from 2002 to 2010.

The geocentric solution; S-N, E-W and vertical displacements are given in centimetres. The seasonal component to the vertical component is associated with climatic variations; velocities of +0.370±0.011 and -0.357±0.021 mm per year.

значения за период 1995-2006 гг. составила 0.27 мкгал/год.

В Европе по данным германской обсерватории (Bad Homburg, 1994-2004 гг. использована группа абсолютных гравиметров FG5) стабильность значения g лежит в пределах ошибки измерений 2 мкгал [Hinderer, Crossley, 2006], т.е. уровень вариаций, или стабильность, составляет 0.2 мкгал/год (рис. 5).

В итоге приходим к выводу, что по экспериментальным данным стабильность значения силы тяжести Земли следует оценить как 0.2 мкгал. На современном этапе возможный эффект глобального изменения g и высоты можно ограничить 10-10 в год от нормального значения. Даже для десятилетней серии измерений эффект от возможных глобальных вариаций не превысит погрешности измерений.

3. Мониторинговые измерения в обсерватории Талая

Первые измерения с гравиметром ГАБЛ в БРЗ проведены осенью 1992 г. на сейсмостанции Талая (TLY, 51.68°N, 103.64°Е), расположенной в 7 км от югозападного берега озера Байкал и в 3 км южнее Главного Саянского разлома, являющегося границей Сибирской платформы и БРЗ (рис. 6). Здание станции расположено в узкой горной долине, дно ее заполнено моренными отложениями мощностью от 0 до 3 м, лежащими на коренных породах. Для проведения гравиметрических и сейсмологических измерений в 1987 г. в здании построен подвал с двумя постаментами, основания постаментов практически достигают коренных пород долины. В 30 м от здания расположена 90-

123

V.Yu. Timofeev et al.: Gravity variations and recent geodynamics...

504.5 504

503.5 503

502.5 502

501.5 501

500.5 500

< 1

1 I

1 1 .

1| еднее значение 502.50 мкгал [ ™ 1

1 1

4 *

Время, Годы

814

812

810

t# г»4 гФ ,-$ rS?

800

Среднеезначеннее11.Э4мкгап

Efeewi, Годэ!

2004 2006

b

a

Рис. 4. Результаты гравиметрических измерений на сибирских станциях (последние три цифры - абсолютное значение в микрогалах).

Время - годы, значения в микрогалах. а - обсерватория Ключи, стабильность значения силы тяжести за период 1977-2009 гг., 1.8 мкгал/32 года = 0.06 мкгал/год; b - астрономическая площадка (г. Иркутск), стабильность значения за период 1995-2006 гг. - 3 мкгал/11 лет = 0.27 мкгал/ год.

Fig. 4. Results of gravimetric measurements at stations in Siberia, Russia (the last three digits of the absolute value are given in microgals).

Time is given in years; values are given in microgals. a - the Klyuchi observatory, the stability of gravity values for the period from 1977 to 2009 - 1.8 microgals/32 years = 0.06 microgals per year; b - the astronomical site in Irkutsk, the stability of values for the period from 1995 to 2006 -3 microgals/11 years = 0.27 microgals per year.

метровая штольня, оборудованная для проведения сейсмологических и геофизических измерений. На другой стороне долины в 100 м от здания в коренных породах пробурена 100-метровая скважина, на которой с 1993 г. проводятся замеры уровня воды. По условиям видимости спутников для установки GPS-пункта выбрана открытая площадка в 300 м выше по долине, постамент со штангой установлен в 1999 г., а измерения проводятся с весны 2000 г. Гравиметрические измерения проводятся ежегодно, а в отдельные периоды -несколько раз в год.

На основании результатов гравиметрических наблюдений в Новосибирской области (ст. Ключи) и в г. Иркутске можно заключить, что наблюдаемые на пункте «Талая» вариации Ag являются региональными и связаны с сейсмическими процессами в БРЗ. При этом основными причинами изменений Ag могут быть две: изменения плотности земной коры в районе наблюдений и вариации высоты пункта наблюдений.

Результаты гравиметрического мониторинга на сейсмостанции «Талая» и в г. Иркутске показаны на рис. 7. В эти результаты внесены поправки (согласно рекомендациям Международной гравиметрической комиссии), учитывающие приливные воздействия Лу-

ны и Солнца, а также влияние смещения полюса Земли и притяжения изменяющихся атмосферных масс [Arnautov, 2005]. Вертикальными штриховыми линиями на рис. 7 обозначены моменты близлежащих к пункту наблюдений землетрясений (29.06.1995 г., М=5.6, L=67 км от эпицентра до пункта «Талая»; 25.02.1999 г., М=5.8, L=86 км; 17.09.2003 г., М=5.3, L=155 км; 23.02.2005 г., М=5.3, L=160 км; 27.08.2008 г., М=6.3, L=25 км).

За двадцатилетний (1992-2012 гг.) период наблюдений на пункте «Талая» зарегистрированы обратимые вариации силы тяжести. Амплитуда отклонений от начального значения достигала 25 мкгал при ошибке измерений от 0.5 до 2 мкгал. Дополнительные осложнения при проведении измерений вызваны периодической модернизацией приборного комплекса и изменением эффективной высоты гравиметра. Линейная составляющая вариаций на периодах 1992-2008 гг. и 2000-2008 гг. равна 1 мкгал/год. Зарегистрировано значимое косейсмическое уменьшение значения силы тяжести при Култукском землетрясении 27.08. 2008 г.

Приливный анализ вариаций водного уровня в скважине, пробуренной в скальных породах, показал присутствие зоны линейной трещиноватости по доли-

124

Geodynamics & Tectonophysics 2013 Volume 4 Issue 2 Pages 119-134

Absolute Gravity Measurements at Bad Homburg AA@125cm (4/93 to 1/04) FG5-101, FG5-301, FG5-202, FG5-220, FG5-215 and FG5-206 Epochs of Nordic Absolute Gravity Campaigns

Рис. 5. Геофизическая обсерватория (Bad Homburg), Германия.

Вариации значения силы тяжести в микрогалах за период 1993-2004 гг. (результаты измерений абсолютными гравиметрами типа FG-Фарелла).

Fig. 5. The Geophysical Observatory in Bad Homburg, Germany.

Variations in gravity values for the period from 1993 to 2004 (microgals per year) (measurements by FG-Farell absolute gravimeters).

51,583 I Скважина, водн.уровень, H = 552 м |—^

Гравит.пункт, Н = 553 м

GPS пункт. Штольня, Н = 563 м 'х Н  578 м

40SN

103,635 103,636 103,637 103,630 103,633 103,64 103,641 103,642 103,643 103,644 103,645

| Рис. 6. Положение обсерватории Талая и размещение пунктов измерений на обсерватории. I Fig. 6. Locations of the Talaya Observatory and measurement sites in the observatory’s area.

125

V.Yu. Timofeev et al.: Gravity variations and recent geodynamics...

Рис. 7. Вариации силы тяжести в микрогалах, время - годы.

Вверху - изменения на пунктах «Талая» (1992-2012 гг.) и «Иркутск» (1995-2006 гг.), желтая линия - линейный тренд 1992-2007 гг., зеленая линия - линейный тренд 1999-2007 гг. Внизу - вариации на пункте «Талая» - линейный тренд на периодах 1992-2005 и 19922009, изменения значения силы тяжести при снятом линейном тренде (1 мкгал/год), стрелками показаны моменты землетрясений.

Fig. 7. Variations of gravity (microgal) with time (year).

Top - changes at the Talya (1992-2012) and Irkutsk sites (1995-2006); the yellow line shows the linear trend for the period from 1992 to 2007; the green line shows the linear trend for the period from 1999 to 2007. Bottom - variations at the Talaya site - the linear trend for periods from 1992 to 2005 and from 1992 to 2009; changes of gravity values for the linear trend (1 microgal per year); arrows show time of earthquakes.

126

Geodynamics & Tectonophysics 2013 Volume 4 Issue 2 Pages 119-134

0

TALY plate North velocity, 2004-2010, 5752952.554 m

E n

~ 0 -10 -20

-M-I____

5 i

~i_______

0,12 0,08 * 0,04 0,00 -0,04

0

2002 2004

IRKT-TLY mean, 2000-2008

-u-

-*T-

T-T-T-CNCNOOOTtTt^lOlO<D<D<DNNCOCOCO

OOOOOOOOOOOOOOOOOOOOO

T-<OT-COOOT-IOOCMh.CMTtO>CM<OT-Tl-OOT-IOO

в

2,5400 2,5300 2,5200 2,5100 E 2,5000 2,4900 2,4800 2,4700 2,4600

5 5

8 8 8 8 1 1 8 8 8 8

в N $ N й ?! в я 8 я

5 § 5 § 5 5 5 5 5 5

8,9000

8,8900

8,8800

8,8700

8,8600

8,8500

8,8400

Date

TALY plate East velocity, 2004-2010, 7153497.145 m

TALY plate height velocity, 2004-2010, 578.929 m

5555555555555

rth

4023

8 8962

8.8907

8,8717

H

8.8610

Date

Рис. 8. Скорости смещения пункта «Талая» (на каждом графике сверху вниз: С-Ю, В-З, вертикальная; + на север, на восток и вверх).

(а) и (б) относительно пункта «Иркутск» (рис. 3). а - результаты обработки данных по программе GAMIT/GLOBK эпоха 2000-2006 гг.: Vn = -0.32±0.16 мм, Ve = +0.82±0.14 мм и Vh = -2.16±0.65 мм/год. Шкала в мм. б - ± результаты обработки данных по программе GPSur-vey, эпоха весна 2000 г. - весна 2008 г.: Vn = -0.1 мм в год, Ve = +1.5 мм в год и Vh = -2.8 мм в год. Шкала в метрах. в - косейсмический скачок смещений, период осень 2004 г. - весна 2010 г.: Vn = -2±2 мм (на юг), Ve = 10±2 мм на восток и опускание 20±5 мм. Решение геоцентрическое - относительно центра масс Земли. Шкала в метрах.

Fig. 8. Velocities of displacement at the Talaya site (in each curve, from top to bottom: N-S, E-W, vertical, + to the north, to the east and upward).

(a) and (б) with reference to Irkutsk (Fig. 3). a - results of data processing by GAMIT/GLOBK software for the period from 2000 to 2006: Vn = -0.32±0.16 mm per year, Ve = +0.82±0.14 mm per year, and Vh = -2.16±0.65 mm per year. The scale is given in millimetres. б - results of data processing by GPSurvey software for the period from the spring of 2000 to the spring of 2008: Vn = -0.1 mm per year, Ve = +1.5 mm per year, and Vh = -2.8 mm per year. The scale is given in metres. в - Coseismic displacement discontinuity for the period from the autumn of 2004 to the spring of 2010: Vn = -2±2 mm (to the south), Ve = 10±2 mm in the east and the lower 20±5 mm. The geocentric solution is given with reference to the Earth’s mass centre; the scale is given in metres.

не. Ежедневный опрос значения водного уровня в скважине № 1608 на станции «Талая» проводился с 1993 г. по 2003 г., цифровая запись уровня начата в 2005 г. [Timofeev et al., 2012а].

Отметим, что уровень находится около поверхности (высота патрубка 0.8 м), а его вариации связаны с летними осадками и сезонными эффектами. В то же время здесь отсутствует резкое весеннее изменение уровня, возможно, вследствие большой глубины про-

мерзания грунтов. Из анализа многолетних изменений можно уверенно говорить лишь о сезонных вариациях, вызванных обильными летними дождями (июль, август) и о периодических 10-12-летних вариациях, связанных с климатическими циклами [Timofeev et al., 2012а].

Анализ показал отсутствие корреляции полученных изменений силы тяжести с вариациями уровня воды в скважине, а в результате специальных исследований

127

V.Yu. Timofeev et al.: Gravity variations and recent geodynamics...

I Рис. 9. Результаты моделирования косейсмических смещений при Култукском землетрясении 2008 г. (М=6.3) [Toda et al., 2009].

I Fig. 9. Coseismic displacement modeling results for the Kultuk earthquake (2008, M=6.3) [Toda et al., 2009].

(2001-2003 гг.) определена амплитуда сезонных изменений g, она достигает 3 мкгал. Снятие линейной составляющей с графика силы тяжести проявило эффекты в эпохи близких землетрясений (М=5.6, 29.06.1995 г., L=67 км; М=5.8, 25.02.1999 г., L=86 км; М=6.3, 27.08.2008 г., L=25 км), при этом наиболее сильный эффект зарегистрирован при Култукском землетрясении 2008 г. Значимые изменения g зарегистрированы в периоды землетрясений с магнитудой больше 5.5 на расстояниях менее 90 км от станции (рис. 7).

Вариации 3D смещений пункта за период 20002010 гг., полученные с использованием программ GPSurvey и GAMIT/GLOBK [King, Bock, 2006; Herring et al., 2006], показаны на рис. 8.

Горизонтальные скорости согласуются с результатами, полученными по данным Байкальской сети за период 1994-2007 гг. [San’kov et al., 2009; Lukhnev et al., 2010]. Для периода 2000-2008 гг. перед Култук-ским землетрясением приведены решения для скоростей относительно пункта «Иркутск», расположенного на Сибирской платформе. Значения по компонентам составили: Vn = - 0.1 ... -0.3 мм/год, Ve = +0.8 ... +1.5

мм/год, Vv = -2.1 ... -2.8 мм/год. Составляющая вертикальной скорости может быть связана с рифтовыми процессами региона. Для периода 2004-2010 гг. приведено геоцентрическое решение. Здесь хорошо проявились косейсмические смещения для всех трех компонент. На периоде весна 2008 г. - весна 2010 г. изменения составили соответственно: N = -2±2 мм, E = = +10+2 мм, V = -20±5 мм. Моделирование скачка смещений проводилось с использованием сейсмологических данных, полученных Байкальским филиалом Геофизической службы СО РАН. На рис. 9 приведены результаты моделирования. В целом модели, построенные по данным сейсмологии и космической геодезии, соответствуют друг другу.

Рассматривая данные измерений в период 20092012 гг., следует отметить, что вертикальные смещения продолжали изменяться (рис. 10). Вариации продолжались и в значениях силы тяжести (см. рис. 7). Обратимся к комплексной интерпретации данных гравиметрии и космической геодезии.

Рассмотрим соотношения, связывающие изменения силы тяжести, вариаций высоты и плотности горных

128

II

Geodynamics & Tectonophysics 2013 Volume 4 Issue 2 Pages 119-134

TALY Up Offset 578.929 m

rate(mm/yr)= -4.21 + 0.35 nrms= 3.88 wrms= 25.5 mm # 25

Рис. 10. Вертикальное смещение пункта «Талая», мм, в период 2000-2012 гг.

Результаты обработки данных по программе GAMIT/GLOBK, геоцентрическое решение, шкала в мм. Скорость смещения до Култукско-го землетрясения -2.8 мм/год. Далее косейсмическое смещение в 20 мм, в конце 2010 г. зарегистрирован постсейсмический скачок до 40 мм, и в конце 2012 г. начинается обратное движение.

Fig. 10. Vertical displacement (in millimetres) of the Talaya site in the period from 2000 to 2012.

Results of data processing by GAMIT/GLOBK software; the geocentric solution; the scale is given in millimetres. Prior to the Kultuk earthquake, the displacement velocity amounted to 2.8 mm per year. Afterwards the coseismic displacement of 20 mm occurred; in late 2010, the post-seismic discontinuity up to 40 mm occurred, and at the end of 2012, the reverse movement commenced.

пород. Для воздействий, охватывающих часть земной коры с линейными размерами в 100-200 км, идеализируя, их можно свести к трем положениям.

1. При плотности p=const.

Связь вариаций g и высоты Z без учета деформации среды можно описать как

Ag

где у - вертикальный градиент g и AZ(t) - изменение высоты пункта.

2. При высоте Z=const.

Связь вариаций g и изменений плотности р, вызванных деформацией среды. Можно использовать формулу для притяжения бесконечного плоского слоя [Tsuboi, 1982]:

Ag = 2п к •Ар• d =2п к •Ае •р• d,

где k=6.67-10-8 см3/гс2 - гравитационная постоянная, d - толщина деформированного слоя земной коры, р=2.67 г/см3 - средняя плотность земной коры, Ар -изменение р из-за объемной деформации Ае земной коры.

3. Связь Ag и изменений смещений и деформации записывается как [Sasagawa, Zumberge, 1991]:

Ag = [у + 2nkp

где k - гравитационная постоянная, р - плотность, v -коэффициент Пуассона, Az - изменение высоты, у -вертикальный градиент.

Нормальное значение вертикального градиента составляет 308 мкгал на метр. Экспериментальное значение, полученное по данным измерений гравиметром Лакоста-Ромберга (LCR G 402) на пункте «Талая» 180 мкгал на метр. Такое значение обусловлено узкой горной долиной и расположением постамента в подвале на 2.5 м ниже земной поверхности. При нормальном градиенте получаем эффект +3.0 мкгал при опускании поверхности на 1 см (+1 мкгал ~ -3 мм), используя экспериментальное определение, оцениваем эффект в +1.8 мкгал при опускании поверхности на 1 см.

Рассматривая соотношения вертикального опускания пункта «Талая» в период 2000-2008 гг. и значение силы тяжести, получаем качественное соответствие. Скорость опускания до -2.8 мм/год и увеличение значения силы тяжести +1 мкгал/год ближе к нормальному значению вертикального градиента.

Нормальные значения параметров составляют для средней плотности земной коры р = 2.67 г/см3, для коэффициента Пуассона v = 0.25. Горные породы - архейские мраморы, которыми сложена долина р. Талая,

129

V.Yu. Timofeev et al.: Gravity variations and recent geodynamics...

a

b

Изменение объёмной деформации (х106), период 1989-2012 гг.

Фазовое запаздывание в градусах, период 1995-2012 гг.

Рис. 11. Результаты обработки данных деформографических измерений.

а - вариации фазового приливного запаздывания в градусах, период наблюдений 1995-2012 гг.; b - изменение значения объемной деформации за период 1989-2012 гг. Стрелками показаны периоды землетрясений [Timofeev et al., 2012b].

Fig. 11. Results of deformographic measurement data processing.

a - variations of phased tidal lag in degrees for the observation period from 1995 to 2012; b - changes of the volume deformation for the period from 1989 to 2012. Arrows show time of earthquakes [Timofeev et al., 2012b].

были исследованы лабораторно. На образце керна, извлеченном из 100-метровой скважины, получены значения коэффициента Пуассона v=0.17 и плотности р=2.87 г-см-3. По соотношению (5) получаем сходные оценки. Так, при изменении деформации Де=1-10~6 и мощности слоя в 1 км эффект составит 0.12 мкгал, при мощности 3 км - 0.36 мкгал, при мощности 10 км - 1.2 мкгал, при мощности 30 км - 3.6 мкгал. Эффект положителен при сжатии, при растяжении отрицателен.

Эти соотношения (5) соответствуют периодам локальных землетрясений (см. рис. 7). Путем исключения линейной составляющей получены аномальные значения до и после землетрясений, случившихся в окрестностях станции (до 90 км), они достигают 5-10 мкгал. Согласно зарегистрированным деформографами данным (рис. 11), на поверхности деформации достигают Де=2-10-6-3-10-6 [Timofeev et al., 2012b]. При мощности слоя земной коры 25 км эффект, согласно соотношению (5), составит 4.8-7.2 мкгал.

Для комплексного эффекта смещений и деформаций, используя соотношение (6), получаем при нормальных значениях вертикального градиента оценку +0.27 мкгал на +1 см вертикального смещения (опускание поверхности «-»). По данным, полученным экспериментально, оценка эффекта составит +4.07 мкгал при смещении +1 см. Последние оценки с учетом ошибок измерений соответствуют косейсмическому эффекту Култукского землетрясения (рис. 12) и результатам моделирования по данным сейсмологии (см.

рис. 9). Процесс опускания земной поверхности при Култукском землетрясении продолжался с 2008 г. по 2012 г., проходил в несколько этапов и зарегистрирован методами космической геодезии и абсолютной гравиметрии. Особенности процесса и его продолжительность, возможно, свидетельствуют о включении флюидов в постсейсмический процесс. Этот вопрос требует дополнительных многолетних наблюдений.

4. Заключение

Результаты работы можно свести к нескольким выводам. Экспериментальные оценки стабильности значения силы тяжести и высоты по данным трех обсерваторий ограничивают изменение радиуса Земли значением 10-10 в год от нормального значения. Эти оценки можно использовать для тестирования тектонических гипотез.

Линейная составляющая вариации силы тяжести на станции Талая соответствует изменениям высоты пункта. Коэффициент корреляции близок к нормальному значению вертикального градиента силы тяжести.

Косейсмические вариации силы тяжести на этом пункте в эпоху Култукского землетрясения 2008 г. вызваны комплексным эффектом изменения высоты пункта и деформации земной коры.

Процесс опускания земной поверхности при Кул-

130

Geodynamics & Tectonophysics 2013 Volume 4 Issue 2 Pages 119-134

I Рис. 12. Вариации силы тяжести, зарегистрированные абсолютным гравиметром и пересчитанные из вертикальных смещений.

| Fig. 12. Gravity variations recorded by an absolute gravity meter and those calculated from vertical displacements.

тукском землетрясении продолжался три года и проходил в несколько этапов. Природу этого процесса, влияние флюидов в земной коре [Fialko, 2004] и асейсмического скольжения можно рассматривать через несколько лет дополнительных наблюдений.

Вариации силы тяжести в периоды региональных землетрясений с магнитудой больше 5.5, произошедших на расстояниях менее 90 км от станции, обуслов-

лены изменением деформируемого состояния земной коры. Для объяснения гравитационного эффекта изменения деформации должны составить величину Де=240-6-3Л0-6 и охватывают слой мощностью 20 км.

Работа была выполнена в рамках проектов Президиума РАН № 4.1, МИП СО РАН № 76, ОНЗ РАН № 6-2.

5. Литература

Altamimi Z., Collilieux X., Legrand J., Garayt B., Boucher C., 2007. ITRF2005: A new release of the International Terrestrial Reference Frame based on time series of station positions and Earth Orientation Parameters. Journal of Geophysical Research 112 (B9), B09401. http://dx.doi.org/10.1029/2007JB004949.

Arnautov G.P., 2005. Results of international metrological comparison of absolute laser ballistic gravimeters. Avtometria 41 (1), 126-136 (in Russian) [Арнаутов Г.П. Результаты международных метрологических сравнений абсолютных лазерных баллистических гравиметров // Автометрия. 2005. Т. 41. № 1. С. 126-136].

Arnautov G.P., Kalish E.N., Smirnov M.G., Stus’ Yu.F., Tarasyuk V.G., 1994. Laser ballistic gravimeter GABL-M and gravity observation results. Avtometria 3, 3-11 (in Russian) [Арнаутов Г.П., Калиш Е.Н., Смирнов М.Г., Стусь Ю.Ф., Тарасюк В.Г. Лазерный баллистический гравиметр ГАБЛ-М и результаты наблюдений силы тяжести // Автометрия. 1994. № 3. С. 3-11].

Ducarme B., Timofeev V.Yu., Everaerts M., Gornov P.Y., Parovishnii V.A., van Ruymbeke M., 2008. A Trans-Siberian tidal gravity profile (TSP) for the validation of the ocean tides loading corrections. Journal of Geodynamics 45 (2-3), 73-82. http://dx.doi.org/10.1016/j.jog.2007.07.001.

Elhiny M.W., Taylor S.R., Stevenson D.J., 1978. Limits to the expansion of Earth, Moon, Mars and Mercury and to change in the gravitational constant. Nature 271 (5643), 316-321. http://dx.doi.org/10.1038/271316a0.

Fialko Y., 2004. Evidence of fluid-filled upper crust from observations of postseismic deformation due to the 1992 Mw7.3 Landers earthquake. Journal of Geophysical Research 109 (B8), B08401. http://dx.doi.org/10.1029/2004JB002985.

Herring T.A., King R.W., McClusky S.C., 2006. Globan Kalman Filter VLBI and GPS Analysis Program. Release 10.3. Department of Earth, Atmospheric, and Planetary Sciences, Massachussetts Institute of Technology, 87 p.

Hinderer J., Crossley D., 2006. NEWSLETTER 17, 5 June 2006, 2-12.

Kaftan V.I., Tsyba E.N., 2009. Assessment of changes in semi-axes of the Earth's geometric ellipsoid from results of satellite

131

V.Yu. Timofeev et al.: Gravity variations and recent geodynamics...

observations of the global geodetic networks. Izvestia VUSov, Geodeziya i aerofotos’emka 1, 33-40 (in Russian) [Кафтан В.И., Цыба Е.Н. Оценка изменений полуосей земного геометрического эллипсоида по результатам спутниковых наблюдений в глобальной геодезической сети // Известия ВУЗов, Геодезия и аэрофотосъемка. 2009. № 1. С. 33-40].

King R.W., Bock Y., 2006. Documentation for the GAMIT GPS. Analysis software. Release 10.3. Cambridge, Massachussetts Institute of Technology, 182 p.

Kolomiets A.G., Gerasimenko M.D., Il’nitskaya A.V., 2010. Problems of seismicity and modern geodynamics of the Far East and East Siberia. In: Reports of the Scientific Symposium, 1-4 June 2010, Khabarovsk. ITIG, Khabarovsk, p. 87-88 (in Russian) [Коломиец А.Г., Герасименко М.Д., Ильницкая А.В. Проблемы сейсмичности и современной геодинамики Дальнего Востока и Восточной Сибири // Доклады научного симпозиума 1-4 июня 2010 г., Хабаровск. Хабаровск: ИТИГ, С. 87-88.

Kukol Z., 1987. Velocities of Geological Processes. Mir, Moscow, 246 p. (in Russian) [Кукол З. Скорость геологических процессов. М.: Мир, 1987. 246 с.].

Lukhnev A.V., San'kov V.A., Miroshnichenko A.I., Ashurkov S.V., Calais E., 2010. GPS rotation and strain rates in the Baikal-Mongolia region. Russian Geology and Geophysics 51 (7), 1006-1017. http://dx.doi.org/10.1016Zj.rgg.2010.06.006.

Robertson L., Francis O., van Dam T.M., Faller J., Ruess D., Delinte J.-M., Vitushkin L., Liard J., Gagnon C., Guang G.Y., Da Lun H., Yuan F.Y., Yi X.J., Jeffries G., Hopevell H., Edge R., Robinson I., Kibbe B., Makinen J., Hinderer J., Amalvict V., Luck B., Wilmes H., Rehren F., Schmidt K., Schnull M., Cerutti G., Germmak A., Zabek Z., Pachuta A., Arnautov G., Kalish E., Stus Yu., Stizza D., Friderich J., Chartier J.-M., Marson I., 2001. Results from the fifth international comparison of absolute gravimeters, ICAG'97. Metrologia 38 (1), 71-78. http://dx.doi.org/10.1088/0026-1394/38/1Z6.

San’kov V.A., Lukhnev A.V., Miroshnichenko A.I., Ashurkov S.V., Byzov L.M., Dembelov M.G., Calais E., Deverchere J., 2009. Extension in the Baikal rift: Present-Day kinematics of passive rifting. Doklady Earth Sciences 425 (2), 205-209. http://dx.doi.org/10.1134/S1028334X09020056.

Sasagawa S., Zumberge M.A., 1991. Absolute gravity measurements in California, 1984-1989. Journal of Geophysical Research 96 (B2), 2501-2513. http://dx.doi.org/10.1029/90JB02283.

Timofeev V.Yu., Ardyukov D.G., Boiko E.V., Gribanova E.I., Semibalamut V.M., Timofeev A.V., Yaroshevich A.V., 2012b. Strain and displacement rates during a large earthquake in the South Baikal region. Russian Geology and Geophysics 53 (8), 798-816. http://dx.doi.org/10.1016/j.rgg.2012.06.007.

Timofeev V.Yu., Ardyukov D.G., Boiko E.V., Timofeev A.V., Yaroshevich A.V., Gornov P.Yu., 2012a. Estimation of coseismic deformation, poroelasticity, and fracturing of rocks from the data on water level in a borehole. Izvestiya, Physics of the Solid Earth 48 (7-8), 640-652. http://dx.doi.org/10.1134/S1069351312060067.

Toda S., Ross S., Sevilgen V., Lin J., 2009. Coulomb 3. Graphic-rich deformation & stress-change software for earthquake, Tectonic, and Volcano Research and Teaching - User Guide. Available from: http://pubs.usgs.gov/of/2011/1060 (last accessed 30.04.2013).

Tsuboi T., 1982. The Earth’s Gravity Field. Mir, Moscow, 286 p. (in Russian) [Цубои Т. Гравитационное поле Земли. М.: Мир, 1982. 286 c.].

Тимофеев Владимир Юрьевич, докт. физ.-мат. наук, заведующий лабораторией Институт нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН 630090, Новосибирск, пр. Академика Коптюга, 3, Россия Тел.: (383)3356442, И e-mail: timofeevvy@ipgg.nsc.ru

Timofeev, Vladimir Yu., Doctor of Physics and Mathematics, Head of Laboratory A.A. Trofimuk Institute of Petroleum Geology and Geophysics, Siberian Branch of RAS 3 Akademika Koptyuga Prosp., Novosibirsk 630090, Russia Tel.: (383)3356442; И e-mail: timofeevvy@ipgg.nsc.ru

Калиш Евгений Николаевич, канд. физ.-мат. наук, с.н.с.

Институт автоматики и электрометрии СО РАН 630090, Новосибирск, пр. Академика Коптюга, 1, Россия Тел.: (383)3332596, e-mail: kalish@iae.nsk.su

Kalish, Evgeniy N., Candidate of Physics and Mathematics, Senior Researcher Institute of Automation and Electrometry, Siberian Branch of RAS 1 Akademika Koptyuga Prosp., Novosibirsk 630090, Russia Tel.: (383)3332596; e-mail: kalish@iae.nsk.su

132

Geodynamics & Tectonophysics 2013 Volume 4 Issue 2 Pages 119-134

Стусь Юрий Фёдорович, канд. техн. наук, с.н.с.

Институт автоматики и электрометрии СО РАН 630090, Новосибирск, пр. Академика Коптюга, 1, Россия Тел.: (383)3332596, e-mail: stus@iae.nsk.su

Stus’, Yuri F., Candidate of Technical Sciences, Senior Researcher Institute of Automation and Electrometry, Siberian Branch of RAS 1 Akademika Koptyuga Prosp., Novosibirsk 630090, Russia Tel.: (383)3332596; e-mail: stus@iae.nsk.su

Ардюков Дмитрий Геннадьевич, с.н.с.

Институт нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН 630090, Новосибирск, пр. Академика Коптюга, 3, Россия Тел.: (383)3356442, e-mail: ardyukovdg@ipgg.nsc.ru

Ardyukov, Dmitriy G., Senior Researcher

A.A. Trofimuk Institute of Petroleum Geology and Geophysics, Siberian Branch of RAS 3 Akademika Koptyuga Prosp., Novosibirsk 630090, Russia Tel.: (383)3356442; e-mail: ardyukovdg@ipgg.nsc.ru

Арнаутов Геннадий Петрович, канд. техн. наук, в.н.с.

Институт автоматики и электрометрии СО РАН 630090, Новосибирск, пр. Академика Коптюга, 1, Россия Тел.: (383)3307931; e-mail: arnautov@iae.nsk.su

Arnautov, Gennadiy P., Candidate of Technical Sciences, Lead Researcher Institute of Automation and Electrometry, Siberian Branch of RAS 1 Akademika Koptyuga Prosp., Novosibirsk 630090, Russia Tel.: (383)3307931; e-mail: arnautov@iae.nsk.su

Смирнов Михаил Георгиевич, ведущий инженер-конструктор Институт автоматики и электрометрии СО РАН 630090, Новосибирск, пр. Академика Коптюга, 1, Россия Тел.: (383)3332596; e-mail: smirnov.m.g@iae.nsk.su

Smirnov Mikhail G., Lead Design Engineer

Institute of Automation and Electrometry, Siberian Branch of RAS 1 Akademika Koptyuga Prosp., Novosibirsk 630090, Russia Tel.: (383)3332596; e-mail: smirnov.m.g@iae.nsk.su

Тимофеев Антон Владимирович, инженер

Институт нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН 630090, Новосибирск, пр. Академика Коптюга, 3, Россия Тел.: (383)3356442; e-mail: timofeevav@ipgg.nsc.ru

Timofeev, Anton V., Engineer

A.A. Trofimuk Institute of Petroleum Geology and Geophysics, Siberian Branch of RAS 3 Akademika Koptyuga Prosp., Novosibirsk 630090, Russia Tel.: (383)3356442; e-mail: timofeevav@ipgg.nsc.ru

Носов Дмитрий Алексеевич, инженер-программист Институт автоматики и электрометрии СО РАН 630090, Новосибирск, пр. Академика Коптюга, 1, Россия Тел.: (383)3332596, e-mail: danossov@ngs.ru

Nosov, Dmitriy A., Software Engineer

Institute of Automation and Electrometry, Siberian Branch of RAS 1 Akademika Koptyuga Prosp., Novosibirsk 630090, Russia Tel.: (383)3332596; e-mail: danossov@ngs.ru

Сизиков Игорь Сергеевич, инженер-конструктор Институт автоматики и электрометрии СО РАН 630090, Новосибирск, пр. Академика Коптюга, 1, Россия Тел.: (383)3332596, e-mail: sizikov.i.s@gmail.com

Sizikov, Igor S., Design Engineer

Institute of Automation and Electrometry, Siberian Branch of RAS 1 Akademika Koptyuga Prosp., Novosibirsk 630090, Russia Tel.: (383)3332596; e-mail: sizikov.i.s@gmail.com

133

V.Yu. Timofeev et al.: Gravity variations and recent geodynamics...

Бойко Елена Валерьевна, н.с.

Институт нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН 630090, Новосибирск, пр. Академика Коптюга, 3, Россия Тел.: (383)3356442, e-mail: boykoev@ipgg.nsc.ru

Boiko, Elena V., Researcher

A.A. Trofimuk Institute of Petroleum Geology and Geophysics, Siberian Branch of RAS 3 Akademika Koptyuga Prosp., Novosibirsk 630090, Russia Tel.: (383)3356442; e-mail: boykoev@ipgg.nsc.ru

Грибанова Екатерина Ивановна, н.с.

Сибирский филиал Геофизической службы СО РАН 630090, Новосибирск, пр. Академика Лаврентьева, 13-3

Gribanova, Ekaterina I., Researcher

Siberian Division of Geophysic Survey, Siberian Branch of RAS 13 Akademika Lavrentieva Prosp., Novosibirsk 630090, Russia

134