RU2812095C1 - Способ краткосрочного прогнозирования землетрясений - Google Patents
Способ краткосрочного прогнозирования землетрясений Download PDFInfo
- Publication number
- RU2812095C1 RU2812095C1 RU2023119854A RU2023119854A RU2812095C1 RU 2812095 C1 RU2812095 C1 RU 2812095C1 RU 2023119854 A RU2023119854 A RU 2023119854A RU 2023119854 A RU2023119854 A RU 2023119854A RU 2812095 C1 RU2812095 C1 RU 2812095C1
- Authority
- RU
- Russia
- Prior art keywords
- seismic
- vector
- rocking
- length
- expected
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 239000013598 vector Substances 0.000 claims abstract description 29
- 230000001133 acceleration Effects 0.000 claims abstract description 14
- 230000035939 shock Effects 0.000 claims abstract description 13
- 238000012544 monitoring process Methods 0.000 claims abstract description 4
- 238000013213 extrapolation Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 abstract 1
- 239000002243 precursor Substances 0.000 description 11
- 230000010355 oscillation Effects 0.000 description 6
- 230000005686 electrostatic field Effects 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000007872 degassing Methods 0.000 description 3
- 239000001257 hydrogen Substances 0.000 description 3
- 229910052739 hydrogen Inorganic materials 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 150000002431 hydrogen Chemical class 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000003534 oscillatory effect Effects 0.000 description 2
- 229910052704 radon Inorganic materials 0.000 description 2
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000005483 Hooke's law Effects 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 239000001307 helium Substances 0.000 description 1
- 229910052734 helium Inorganic materials 0.000 description 1
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000011089 mechanical engineering Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Images
Abstract
Изобретение относится к области сейсмологии и может быть использовано для мониторинга природных сред. Предложен способ краткосрочного прогнозирования землетрясений. Согласно заявленному способу размещают в сейсмоопасном регионе две пары приемных станций глобальной навигационной системы позиционирования во взаимно ортогональных плоскостях (x, y), разнесенных между собой на измерительных базах до 50 км с высокой точностью центрирования. Отслеживают момент изменения длин баз и отождествляют его с началом сейсмического процесса. Рассчитывают длину вектора раскачки очага через его проекции на осях координат
и
определяют плоскость раскачки очага через косинус направляющих
Численным методом рассчитывают относительное ускорение длины вектора раскачки. Прогнозируют момент удара по достижению предельных напряжений земной коры ≤10-2. Фиксируют координаты и время точки достижения предельных напряжений. Прогнозируют магнитуду (М) ожидаемого сейсмического удара из соотношения
оповещают службы МЧС о прогнозируемом сейсмическом ударе. Технический результат: возможность упреждающего оповещения, достоверность, документальность, оперативность. 6 ил.
Description
Изобретение относится к области сейсмологии и может найти применение в национальных системах наблюдения и обработки данных геофизических измерений для прогнозирования землетрясений.
Существует множество статических признаков подготавливаемого землетрясения, регистрируемых системами наземных наблюдений. Систематизированный перечень известных признаков-предвестников [см., например, «Данные о предвестниках», в книге Т. Рикитаке «Предсказание землетрясений», перевод с англ., МИР, М., 1979 г., табл. 13.15, стр. 314-333]. Перечисленные признаки-предвестники имеют долговременный интервал существования и диапазон изменения, но не позволяют точно предсказать момент сейсмического удара.
Другой класс – динамические признаки-предвестники. Они появляются за несколько часов (суток) до удара, но в силу своей скрытости, инерционности во времени (период изменения 10-3…10-4 сек)не могут быть достоверно зарегистрированы известными техническими средствами.
Достоверно установлено [см., например, «Краткосрочный прогноз катастрофических землетрясений с помощью радиофизических наземно-космических методов», Доклады конференции, РАН, ОИФЗ им. О.Ю. Шмида, М., 1998 г., стр. 30], что за несколько суток до сильного землетрясения, в области его подготовки, в атмосфере над поверхностью очага, наблюдается сильное вертикальное электростатическое поле в виде купола, величиной до нескольких кв/м. Размер области на Земле, над которой это поле существует составляет ~ 120…140 км в диаметре.
Одной из причин возникновения электростатического поля считается активная эманация газов: водорода, гелия, радона и других из зоны сжатия земной коры. Дегазация радона (период полураспады 3,8 суток) сопровождаются γ-излучением, а дегазация водорода приводит к падению атмосферного давления над зоной подготавливаемого землетрясения.
Известен «»Измеритель предвестников землетрясений» Патент Ru №2.310894, 2007 г., G.01V,9/00 – аналог, в котором для повышения достоверности прогноза используют одновременное измерение и концентрация водорода (снижение атмосферного давления Δp на 20-40 мм р.т.) и увеличение фонового уровня радиации на один-два порядка (фоновый уровень ~20 .
Устройство-аналога содержит: измерительный канал из генератора оптического излучения, чувствительного элемента и фотоприемника, подключенных последовательно, а также порогового устройства, связанного с устройством обработки результирующего сигнала и синхронизации работы элементов измерителя, отличающийся тем, что чувствительный элемент выполнен в виде электрооптического модулятора на отрезке дырчатого волокна, при это длительность фотоимпульсов электрооптического модулятора должна быть пропорциональна концентрации легких газов в атмосфере, а генератор оптического излучения подключен к делителю из параллельно включенных датчиков радиоактивного фона и переменного сопротивления, при этом между фотоприемником и пороговым устройством включен пиковый детектор.
Недостатками аналога являются:
- не определяются координаты сейсмического удара;
- поскольку наземный измеритель под воздействием сейсмических волн синфазно смещается вместе с куполом электростатического поля над очагом, то значения измеряемых параметров фактически остаются «const», то есть чувствительность измерителя невысокая.
Ближайшим аналогом к заявленному техническому решению является «Способ краткосрочного предсказания землетрясений», Патент Ru №2.227.311, 2004, G.01V,9/00. Способ ближайшего аналога включает преобразование измеряемой величины в электрический сигнал, взятие дискретных отсчетов амплитуды сигнала в разнесенных точках пространства, формирование регистрограмм измерений зависимости дискретных отсчетов амплитуд от пространственных координат, обработку регистрограмм, отличающийся тем, что образуют из датчиков-измерителей, разнесенных с постоянным шагом по координатам х, у, линейные крестообразные группы, расставляют их в контролируемом регионе ориентированно по сторонам света на расстоянии, соизмеримом с размерами очага, регистрируют датчиками-измерителями амплитуду вертикального электростатического поля у поверхности земли Е, осуществляя их последовательный циклический опрос, формируют регистрограммы каждой из групп раздельно по координатам E (x, t) и E (y, t)рассчитывают гипоцентр очага как точку пересечения векторов, направляющие косинусов которых вычисляют через производные от регистрограмм крестообразной группы:
определяют период сигнала T0сигнала регистрограмм, прогнозируют магнитуду и время удара по регрессионным зависимостям М ≈ 110/T0 2 (ч); ty≈ 4,7T0.
К недостаткам ближайшего аналога следует отнести:
- невысокая чувствительность наземных датчиков измерителей вертикального электростатического поля, поскольку под действием сейсмических волн раскачки очага они синфазно смещаются вместе с куполом возникающего электростатического поля, то есть градиент изменения измеряемого сигнала отсутствует;
- определяемый фазовый центр сейсмических волн, как правило, не совпадает с местом разрыва земной коры, то есть неточность прогноза координат сейсмического удара может составлять десятки километров.
Задача, решаемая заявленным способом состоит в обнаружении и измерении сейсмических волн раскачки очага землетрясения средствами глобальной системы позиционирования, определение плоскости колебаний вектора раскачки очага, численный программный расчет длины вектора и координат превышения его ускорения предельных значений деформации .
Поставленная задача решается тем, что способ краткосрочного прогнозирования землетрясений включает размещение в сейсмоопасном регионе двух пар приемных станций космической навигационной системы позиционирования во взаимно ортогональных плоскостях (x, y), разнесенных между собой внутри плоскостей на измерительных базах до 50 км с центрированием длин баз с высокой точностью, непрерывное отслеживание изменения длин баз и отождествление изменения с началом сейсмического процесса (t нач ), определение плоскости вектора раскачки очага расчетом косинус направляющих cosα, cosβчерез его проекции (П px , П py )на осях координат соответственно:
; длина вектора ;
численный программный расчет относительного ускорения длины вектора (длR) и сличение его с предельным напряжением земной коры:
определение предельной длины вектора раскачки (длR max ) и соответствующей ей координат точки в плоскости раскачки как места ожидаемого сейсмического удара, экстраполяция времени ожидаемого удара t y / при maxдлR, расчет времени упреждающего оповещения как разницы , прогнозирование магнитуды (М) ожидаемого сейсмического удара из соотношения:
Изобретение подтверждается чертежами, где:
фиг.1 – зарегистрированные средствами GPS сейсмические волны раскачки очага состоявшегося землетрясения;
фиг.2 – определение плоскости вектора раскачки очага землетрясения через его проекции на осях координат;
фиг.3 – динамика изменения длины вектора раскачки очага;
фиг.4 – динамика изменения относительного ускорения вектора раскачки очага;
фиг.5 –гистограмма изменения вероятности сейсмического удара в функции изменения относительного ускорения вектора раскачки;
фиг.6 – функциональная схема устройства, реализующая способ.
Техническая сущность изобретения состоит в следующем. Известно «Явление раскачки очага землетрясения перед сейсмическим ударом» [см., например, Научное открытие №365, 2018 г., Бюллетень, Научные открытия, 2009 г., РАЕН, М-Санкт-Петербург, стр. 66-68]. Частота механических колебаний(ω)определяется соотношением [см., например, Вибрации в технике, Справочник под редакцией К.В. Фролова, Т.6, М., Машиностроение, 1981 г.]:
где, c – жесткость виброизолятора;
f – частота, Гц;
m – масса, кг.
Поскольку жесткость земной коры ограничена, а колебательная масса очага землетрясения велика (млн. тонн), то частота колебаний при раскачке составляет 10-3…10-4[Гц]. Такие сверхнизкие частоты невозможно измерить, к тому же если измеритель находиться внутри инерциальной системы [см., например, Геофизические методы мониторинга природных сред, Научный сборник, под редакцией В.Н. Сорокина, АН СССР, Институт общей физики, 1991 г.].
Для измерения сверхнизких колебаний измеритель необходимо вынести за пределы инерциальной системы земной коры. Таким измерителем, вынесенным за пределы земной коры является Глобальная Позиционирующая система «Глонасс».
В процессе раскачки очага землетрясения изменяется консолидированная масса земной коры вовлекаемой в колебание. Скорость изменения консолидированной массы зависит от структуры земной коры, скорости дегазации литосферы, характеристик разломной зоны, в общем является величиной не измеряемой.
Ввиду изменения массы и неопределенности начальных условий, колебательный процесс аналитически не может быть представлен, общим решением дифференциального уравнения второго порядка [см., например, Пискунов Н.С. «Дифференциальное и интегральное исчисление для ВТУЗов», учебник, пятое издание, Наука М., 1964 г., стр. 524-536].
Поэтому для выявления скрытых закономерностей использованы численные методы программной обработки регистрируемых сигналов с последующей экстраполяцией выявленных закономерностей.
На рисунке фиг. 1 представлен графиксейсмического процесса состоявшегося 21.12.1995 г. землетрясения, зарегистрированного на полигоне МЧС Кавказские Минеральные воды [см., например, Давыдов В.Ф. «Землетрясения. Телеметрия предвестников», Монография, Издательство МГУЛ, Мытищи, 2001 г., стр. 19-20]. Графики фиг. 1 фактически являются решением дифференциального уравнения второй степени при начальных условиях земной коры упомянутого полигона при изменяемой массе и периоде колебаний. В зарегистрированных литосферных волнах раскачки очага землетрясений содержится вся информация о сейсмическом процессе, времени существования признака-предвестника, координатах сейсмического удара и его магнитуде.
Задача состоит в том, чтобы по начальным параметрам регистрируемого процесса экстраполировать скрытые закономерности для получения интервала времени на упреждающее оповещение населения.
Известно, что движущийся объект сохраняет плоскость своего вектора. Вектор движения задается своими проекциями на осях координат. Фазовый центр литосферных волн совпадает с гипоцентром очага и лежит в плоскости вектора движения. Плоскость вектора раскачки очага определяется косинус направляющими cosα, cosβ через его проекции (П px , П py ) на осях координат:
; длина вектора .
На рисунке фиг. 2 иллюстрируется определение фазового центра (гипоцентра очага) через его проекции на осях координат. При раскачке изменяется амплитуда колебаний, то есть увеличивается длина вектора (длR). Одна из реализаций изменения длины и периода колебаний иллюстрируется графиком фиг. 3.
Земная кора терпит разрыв при относительных деформациях порядка 10-3[см., например, Т. Рикитаке «Предсказание землетрясений», перевод с англ., М., МИР, 1979 г.].
Из Закона Гука «в пределах упругости величина деформации пропорциональна деформирующей силе» []и Закона Ньютона «сила равна массе умноженной на ускорение» [] следует []. При потере упругости происходит разрыв земной коры, граничные условия . Изменение относительного ускорения при изменении длины вектора [(длR), фиг. 3] определялась численным методом. Предварительно, посредством сканера высокого разрешения функция фиг. 1 преобразовывалась в цифровую форму. На рисунке фиг. 4 представлен график изменения относительного ускорения длины вектора раскачки. Из рисунка фиг. 4 следует, что величина изменяется по знаку и величине. Знак (±) содержит информацию в какой полуплоскости относительно гипоцентра ожидается разрыв земной коры, то есть координаты точки сейсмического удара, а величина содержит информацию о времени существования признака-предвестника, то есть об интервале времени упреждающего оповещения. Программа численного расчета относительного ускорения представлена ниже, в примере реализации.
Пример реализации способа
Заявленный способ может быть реализован по схеме фиг. 6. Функциональная схема устройства фиг. 6 содержит приемные станции (1, 2, 3, 4)размещенные в сейсмоопасном регионе (5), осуществляющие прием эфемеридной информации от группировки космических аппаратов (6), образующих попарно (1, 2), (3, 4) ортогональную систему координат (x, y)разнесенных на измерительных базах (7, 8) порядка 50 километров. При возникновении сейсмического процесса изменяется расстояние (длина баз) между станциями, которое фиксируется станциями. Результаты измерений через порты станций (9, 10, 11, 12) по каналам связи (13) передаются в аналитический центр (14) в составе устройства вывода (15), комплекса средств обработки: центрального процессора (16), оперативного запоминающего устройства (17), винчестера (18), клавиатуры (19) и средств отображения дисплея (20), принтера (21).
В аналитическом центре (14) ведется база данных о всех зарегистрированных землетрясениях и расчет параметров отлеживаемого процесса. Обеспечивается возможность считывания файлов регистрируемых сигналов через сервер (22) в Интернет (23) и систему оповещения (24).
Функционирование комплекса измерительных средств состоит в следующем. В аналитическом центре (14) организуют круглосуточное дежурство. После обнаружения сейсмического процесса (изменение расстояния между станциями) осуществляют непрерывную регистрацию сигналов и по приведенным выше соотношениям определяют гипоцентр очага и плоскость вектора раскачки. Установленные параметры очерчивают на карте местности. Регистрируемый сигнал (фиг. 4) характеризуется изменением амплитуды и изменением периода колебаний.
Известно, что сама функция и скорость ее изменения связаны дифференциальным уравнением первого порядка, общим решением которого является экспонента. Для прогнозирования землетрясения необходимо непрерывно измерять вторую производную, то есть ускорение вектора раскачки. Поскольку аналитическое выражение второй производной получить не представляется возможным, использованы численные методы программной обработки. Для экстраполяции достижения предельных ускорений вектора раскачки использован программный канонический метод разложения функции в ряд Маклорена с ограничением трех членов разложения.
Программа определения ускорения вектора раскачки очага
programSignal;
// Из файла f считывается время и величина сдвига
// В файл g записывается время и величина скорости
// В файл h записывается время и величина ускорения
var
f,g,h:text;
a,b,c,x,x1,x2,x3,y,y1,y2,y3:real;
begin
assign(f,'data.txt');
reset(f);
assign(g,'rezult_v.txt');
rewrite(g);
assign(h,'rezult_a.txt');
rewrite(h);
read(f,x1);readln(f,y1);
read(f,x2);readln(f,y2);
read(f,x3);readln(f,y3);
while not eof(f) do
begin
a:=((y1-y2)*(x2-x3)-(y2-y3)*(x1-x2))/
((sqr(x1)-sqr(x2))*(x2-x3)-(sqr(x2)-sqr(x3))*(x1-x2));
b:=((y1-y2)-a*(sqr(x1)-sqr(x2)))/(x1-x2);
// c:=y1-b*x1-a*sqr(x1);
x:=x2; y:=2*a*x2+b;
writeln(g,x,' ',y);
y:=2*a;
writeln(h,x,' ',y);
x3:=x2; y3:=y2;
x2:=x1; y2:=y1;
if(not eof(f)) then read(f,x1);
if(not eof(f)) then readln(f,y1)
end;
close(f);
close(g);
close(h)
end.
Результат программного расчета иллюстрируется графиком фиг. 4. Расчет позволяет однозначно зафиксировать время удара t y , в левой полуплоскости раскачки очага (фиг. 1, нижняя часть).
Известны соотношения Гутенбенрга-Рихтера, связывающие время существования предвестника t y [сутки] с магнитудой ожидаемого сейсмического удара [см., например, Краткосрочный прогноз катастрофических землетрясений с помощью радиофизических наземно-космических методов, Доклад конференции ОИФЗ им. О.Ю. Шмида, РАН, М., 1998 г., стр. 10-13]:
Для размеров (R) зоны накапливаемой деформации:
Из графика фиг. 4 t y соответствует по времени регистрации сигнала (8 час).
Поскольку средства GPS имеют интервал нечувствительности (порядка 10 см разрешение) время существования признака-предвестника существенно больше времени регистрации сигнала (фиг. 3, 8 час).
Время существования предвестника рассчитывают из соотношения:
где Т 0 – предел, к которому стремится экспоненциальная зависимость [см., например, Патент Ru№2.426155, 2011 г.].
Экспоненциальная зависимость обладает тем свойством, что по трем ее дискретным отсчетам может быть восстановлена вся функция и определен предел Т 0 , к которому стремится экспонента [см., например, Патент Ru№2.426155, 2011 г.]:
В частности для значений T 1 , T 2 , T 3 фиг. 3 расчетное значение Т 0 составило 4 час. Оттуда:
Ожидаемая магнитуда удара:
Размер очаговой зоны подготавливаемого землетрясения R=126 км, Место удара левая полуплоскость ≈ 20 км от гипоцентра.
Способ может быть реализован на существующей технической базе. Эффективность способа характеризуется достоверностью, документальностью, регистрируемого сейсмического процесса, возможностью упреждающего оповещения о сейсмическом ударе за 2-3 часа и его координатах.
Claims (5)
- Способ краткосрочного прогнозирования землетрясений включает размещение в сейсмоопасном регионе двух пар приемных станций космической навигационной системы позиционирования во взаимно ортогональных плоскостях (x, y), разнесенных между собой внутри плоскостей на измерительных базах до 50 км с центрированием длин баз с высокой точностью, непрерывное отслеживание изменения длин баз и отождествление изменения с началом сейсмического процесса (tнач), определение плоскости вектора раскачки очага расчетом косинус направляющих cosα, cosβ через его проекции (Пpx, Пpy) на осях координат соответственно: ; ,
- определение длины вектора ,
- численный программный расчет относительного ускорения длины вектора (длR) и сличение его с предельным напряжением земной коры: ,
- определение предельной длины вектора раскачки (длR max) и соответствующей ей координат точки в плоскости раскачки как места ожидаемого сейсмического удара, экстраполяция времени ожидаемого удара ty при длR max,
- расчет времени упреждающего оповещения как разницы , прогнозирование магнитуды (М) ожидаемого сейсмического удара из соотношения: .
Publications (1)
Publication Number | Publication Date |
---|---|
RU2812095C1 true RU2812095C1 (ru) | 2024-01-22 |
Family
ID=
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2130195C1 (ru) * | 1998-04-30 | 1999-05-10 | Московский государственный университет леса | Способ предсказания землетрясений |
RU2227311C2 (ru) * | 2002-01-04 | 2004-04-20 | Никитин Альберт Николаевич | Способ краткосрочного предсказания землетрясений |
RU2438147C1 (ru) * | 2010-03-31 | 2011-12-27 | Государственное образовательное учреждение высшего профессионального образования "Московский Государственный Университет Леса" (ГОУ ВПО "МГУЛ") | Способ краткосрочного прогнозирования землетрясений |
RU2458362C1 (ru) * | 2011-02-08 | 2012-08-10 | Федеральное государственное бюджетное научное учреждение "Научно-исследовательский институт аэрокосмического мониторинга "АЭРОКОСМОС" | Способ краткосрочного прогнозирования землетрясений |
RU2497158C1 (ru) * | 2012-05-17 | 2013-10-27 | Федеральное государственное бюджетное научное учреждение "Научно-исследовательский институт аэрокосмического мониторинга "АЭРОКОСМОС" | Способ краткосрочного прогнозирования землетрясений |
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2130195C1 (ru) * | 1998-04-30 | 1999-05-10 | Московский государственный университет леса | Способ предсказания землетрясений |
RU2227311C2 (ru) * | 2002-01-04 | 2004-04-20 | Никитин Альберт Николаевич | Способ краткосрочного предсказания землетрясений |
RU2438147C1 (ru) * | 2010-03-31 | 2011-12-27 | Государственное образовательное учреждение высшего профессионального образования "Московский Государственный Университет Леса" (ГОУ ВПО "МГУЛ") | Способ краткосрочного прогнозирования землетрясений |
RU2458362C1 (ru) * | 2011-02-08 | 2012-08-10 | Федеральное государственное бюджетное научное учреждение "Научно-исследовательский институт аэрокосмического мониторинга "АЭРОКОСМОС" | Способ краткосрочного прогнозирования землетрясений |
RU2497158C1 (ru) * | 2012-05-17 | 2013-10-27 | Федеральное государственное бюджетное научное учреждение "Научно-исследовательский институт аэрокосмического мониторинга "АЭРОКОСМОС" | Способ краткосрочного прогнозирования землетрясений |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Geng et al. | A new seismogeodetic approach applied to GPS and accelerometer observations of the 2012 Brawley seismic swarm: Implications for earthquake early warning | |
Bernauer et al. | BlueSeis3A: Full characterization of a 3C broadband rotational seismometer | |
Allen et al. | Application of real‐time GPS to earthquake early warning | |
Bernauer et al. | Rotational sensors—A comparison of different sensor types | |
Lee et al. | Introduction to the special issue on rotational seismology and engineering applications | |
Cochard et al. | Rotational motions in seismology: theory, observation, simulation | |
Li et al. | Heterogeneous interseismic coupling along the Xianshuihe‐Xiaojiang fault system, eastern Tibet | |
Thompson et al. | Observations of premonitory acoustic emission and slip nucleation during a stick slip experiment in smooth faulted Westerly granite | |
Psimoulis et al. | Detection of ground motions using high-rate GPS time-series | |
D’Amico et al. | Fling effects from near‐source strong‐motion records: Insights from the 2016 Mw 6.5 Norcia, central Italy, earthquake | |
Fukuyama | Dynamic faulting on a conjugate fault system detected by near-fault tilt measurements | |
Farghal et al. | The potential of using dynamic strains in earthquake early warning applications | |
Bernard et al. | Onland and offshore extrinsic Fabry–Pérot optical seismometer at the end of a long fiber | |
Psimoulis et al. | Real‐time magnitude characterization of large earthquakes using the predominant period derived from 1 Hz GPS data | |
Fang et al. | GPS/BDS precise point positioning with B2b products for high-rate seismogeodesy: application to the 2021 M w 7.4 Maduo earthquake | |
RU2812095C1 (ru) | Способ краткосрочного прогнозирования землетрясений | |
Larocca et al. | First steps using two GPS satellites for monitoring the dynamic behavior of a small concrete highway bridge | |
RU2436134C1 (ru) | Способ оперативного исследования атмосферы, земной поверхности и океана | |
Kislov et al. | Rotational seismology: Review of achievements and outlooks | |
Xu et al. | Sequential modelling of the 2016 Central Italy earthquake cluster using multisource satellite observations and quantitative assessment of Coulomb stress change | |
Provost et al. | Automatic approach for increasing the location accuracy of slow-moving landslide endogenous seismicity: the APOLoc method | |
Li et al. | Using earthquake doublets to study inner core rotation and seismicity catalog precision | |
RU2463631C1 (ru) | Способ обнаружения очагов землетрясений сетью сейсмостанций | |
JP4818010B2 (ja) | 地震時の建造物変形量に基づく地震マグニチュードの早期予測方法及び地震マグニチュードの早期予測プログラム | |
RU2337382C1 (ru) | Способ краткосрочного предсказания землетрясений |