CC BY
29Journal of Siberian Federal University. Engineering & Technologies, 2017, 10(6), 783-793
УДК 551.507; 551.524; 550.34
Atmospheric Temperature Anomalies During Seismic Activity
Leonid G. Sverdlik*a and Sanjar A. Imashevb
aKyrgyz-Russian Slavic University 44 Kievskaya Str., Bishkek, 720000, Kyrgyzstan bResearch Station of the Russian Academy of Sciences in Bishkek
Bishkek, 720049, Kyrgyzstan
Received 19.12.2016, received in revised form 16.05.2017, accepted 22.08.2017
Presented results were obtained by means of methods of temperature anomalies diagnosing and on the basis of retrospective analysis of temperature profiles transformation in border of troposphere and stratosphere over epicenters of earthquakes with energy class of K>13.0, occurred in Northern and Central Tien-Shan during 2003-2015. Abnormal temperature perturbations observed in about 70% of cases in 1-5 days before the main seismic event.
Keywords: earthquake, satellite data, temperature profiles, tropopause, anomalies, outgoing longwave radiation.
Citation: Sverdlik L.G., Imashev S.A. Atmospheric temperature anomalies during seismic activity, J. Sib. Fed. Univ. Eng.technol., 2017, 10(6), 783-793. DOI: 10.17516/1999-494X-2017-10-6-783-793.
© Siberian Federal University. All rights reserved
* Corresponding author E-mail address: l.sverdlik@mail.ru
Аномалии температуры атмосферы в периоды сейсмической активности
Л.Г. Свердлика, С.А. Имашевб
аКыргызско-Российский Славянский университет Кыргызстан, 720000, Бишкек, ул. Киевская, 44 бНаучная станция Российской академии наук в г. Бишкек
Кыргызстан, 720049, Бишкек
Представлены результаты, полученные с использованием разработанной методики диагностирования аномалий температуры и основанные на ретроспективном анализе процессов трансформации профилей температуры в зоне раздела тропосферы и стратосферы над эпицентрами землетрясений энергетического класса К>13.0, произошедших в Северном и Центральном Тянь-Шане с 2003 по 2015 г. Аномальные возмущения температуры наблюдались примерно в 70 % случаев за 1-5 дней до основного сейсмического события.
Ключевые слова: землетрясение, спутниковые данные, профили температуры, тропопауза, аномалии, уходящая длинноволновая радиация.
Введение
Взаимодействие двух крупных горных систем, Тянь-Шаня и Памира, определяют сложный характер геодинамических процессов и высокий уровень сейсмической активности в Центральной Азии, а соответственно, проведение многочисленных геофизических и геологических исследований с помощью наземных наблюдений и средств космической геодезии (GPS). В то же время диагностике связанных с землетрясениями возмущений в атмосфере уделяется недостаточное внимание, как и использованию спутникового дистанционного зондирования, к достоинствам которого можно отнести глобальность охвата, включая горные районы, и регулярность получения данных о пространственно-временной структуре различных параметров атмосферы [1]. Одним из направлений применения результатов спутниковых измерений в последние годы является установление связи между динамикой литосферных и атмосферных процессов, например, в период катастрофического землетрясения в Японии в 2011 г. [2, 3], а также сильных сейсмических событий на территории России [1]. Если обратиться к этим данным, наиболее сильная отрицательная корреляционная связь между изменениями температу ры в области тропопаузы совпадала с периодом высокой сейсмической активности; предполагается, что возможной причиной наблюдаемых вариаций температуры могли быть землетрясения. Мы представляем методику и результаты ретроспективного анализа спутниковых измерений, которые позволили выявить аномальные вариации температуры в верхней тропосфере как вероятные предвестники сильных сейсмических событий в регионе Северного и Центрального Тянь-Шаня.
Используемые данные
Всего в период с 2003 по 2015 г. на территории, ограниченной координатами 39^43 °N, 7Н-79 °E, произошло 5075 сейсмических событий, имеющих энергетический класс выше 6.0.
Преобладающий вклад вносили землетрясения с 7.0 по 9.0 классы. Исследование потенциальных возможностей спутникового мониторинга применительно к выявлению аномалий температуры атмосферы проведено на примере 15 наиболее сильных землетрясений, которые характеризовались энергией сейсмических волн (К) от 13.0 до 14.8. Приведенные в табл. сведения об этих сейсмических событиях взяты из каталога землетрясений Института сейсмологии НАН КР и Научной станции РАН в г. Бишкеке, составленного по данным сейсмологической сети KNET [4].
Расположения эпицентров землетрясений (рис. 1) большей частью были сконцентрированы в зоне пересечения активных региональных разломов.
Для анализа временных и высотных вариаций температуры атмосферы использовали данные спутниковых измерений (MERRA 3D) на 23-х изобарических уровнях (pk) от 800 до 10 гПа [5]. Вертикальные профили температуры T(p) были восстановлены над каждым участком размером ±1° от эпицентра землетрясения для временных периодов, составляющих 45-60 дней до и 30-45 дней после каждого из рассматриваемых событий. Выбранное время регистрации (00:00 и 12:00 UTC) и дискретность измерений (At=12 ч) позволили ослабить суточную составляющую и сохранить особенности временной изменчивости температуры.
Карты пространственного распределения уходящей длинноволновой радиации (OLR) в пределах области, охватывающей всю территорию Кыргызстана (35.0-45.0 °N; 65.0-85.0 °E), построены по данными NOAA/ESRL [6].
Таблица. Характеристики землетрясений энергетического класса К>13.0 (2003-2015 гг.) Table. Characteristics of earthquakes of energy class K>13.0 (2003-2015)
№ Date Time Latitude N Longitude E Dept, km К n
EQ01 22.05.2003 18:11:18 42-59.32 72-48.82 7.0 14.26 18
EQ02 EQ03 16.01.2004 02.06.2004 09:06:09 17:15:11 42-32.96 42-16.55 75-17.82 74-54.85 13.95 17.88 13.68 13.25 14 2
EQ04 EQ05 EQ06 EQ07 EQ08 EQ09 EQ10 EQ11 EQ12 EQ13 08.11.2006 25.12.2006 06.06.2007 02:21:27 20:00:58 11 0975 42-33.91 42-06.74 42 34 14 75-21.35 76-01.93 75 23 99 18.56 0.12 12 96 13.37 14.83 13 25 5 20 7 2 22 2 2 5 13 1
05.10.2008 02.03.2010 15:52:41 01 55 36 38-53.75 42 2 5 95 73-23.70 75 39 66 50.18 19 30 13.25 13 34
18.03.2011 09 04 9011 09:36:27 12 30 25 43-00.33 42 03 07 74-55.22 74 50 01 17.08 6 27 13.00 13 20
05.02.2012 23.11.2013 07:10:15 094206 41-23.77 42 2 5 62 74-45.67 75 41 22 13.42 9 92 13.21 13 46
14.11.2014 01:24:17 42-11.35 77-13.75 10.08 13.98
EQ14 17.11.2015 17:29:37 40-25.54 73-11.23 3.06 13.53 7
EQ15 07.12.2015 08:30:53 41-39.62 74-41.68 15.91 13.60 3
* К - энергетический класс; D - глубина очага землетрясения (км); п - количество сейсмических толчков в течение указанных суток.
к (N=15, values=[13.0 -> 14.83])
Рис. ¡.Расположение эпицентров землетрясений энергетичеекогокласса от13.0до 14.8. Обозначения сейсмическихсобытий соответствуюттаблице
Fig. 1. The loca6iaeof ^З^е epipenBersoK theoartBqcahcHoftKe exprHOclasHercmlO.Oto 14.0.Denote aeismic events correspondto the table
Методика анализа временных рядов температуры атмосферы
На начальном этапе были сформированы непрерывные временные ряды, содержащие значения температурына всехвысотныхуровнях атмосферь1ооторыезагемпривеДены к безраз-мернот^виду вттеабетг^ а исключеоиеха неливгйносви высотаого расл^деоеоия ыемпарта3рте. Для этого традационнаш образом рассчзтеевалааы ефемечнты далимики аныхелий темперыту-ры (И7Я,какотакынонле ееиащезо мыаченид темчг оаглды ыи средоемьсявного дфовня, рованное несксднекиадкатичиое дтмллненип, аименно:
глтг . з теxi,УJ, «Л Рк)^(Т<ех1,У],Рк)) АТфХ > У] > « Рк)
ш ех, у,р)=--—=--,
стт( в„уЧ аД х.
где п - день измерения; t - времяизмерения; х - широта, у - долгота; i,j - точки сеткикоорди-нат; Т(х, у, t, рк) - текущее значение температуры на каждом уровне (рк); (Тфз, уе, рв)) - среднее овичечие аемпературы]сТ ох^С-краине кеичратичаесоеоииаонение.
ивалиаеразл ичные варта ноаолредст аввеимвер емесвыа рядояаномали йтем-перату^г^Г сак енакиаа. ]е(^1^с^1^^аавквнам:виае;И1^1ч^гса1^ьеоа^пг^с^г^;«одь]^^1е^ислением средних значений температурных аномрлий в алое 600-300 гПа, который примерно соответствует свободной тропосфере (@7ро) и в области, расиороженной выше тропопаузы (&ТТРр: Ы ри гаом ьыЧор выс о тнвьо л^^^г^^г^с^^а над тропепрузой, ^1ео тот ом лроисводится усреднение значе ни й, я вля лся ситуативным и в каждом случае строился на визу альном и корреляционном сравнении поведения аномалий температуры на уровнях в интервале от 200 до 70 гПа [7], исходя из максимального подобия временной изменчивости ©Т. Как результат таких вычислений, на рис. 2 приведены 14-дневные фрагменты графиков, которые демонстрируют в том или ином виде наличие перед сильными землетрясениями аномалий температуры, превышающих 2с и характеризующихся противофазными из- 786 -
600-350 h Pa 200-150 hPa EQ (Кз13.0) EQ СКэЮ.О)
600-300 h Pa 1S0-L00 h Pa
EQ(K=>IS.05
EQ (KslO.OJ
b __
600-300 h PS 200-150 h PS EQ 0O13.O) EQ (KslO.O]
7 9 11 13 15 17 january 2004
A
24 26 2a 30 2 4 6 september j october 200s
11 13 15 17 march 2011
Рис. 2.Вариации аномалийтемпературы НТРТи ©ТТР паред зрмлатртсьниями К>13.О: 16 января 2004 г. (а); ОИ октября 2008 г. (Мы); )8тарта 2011 г. (с).Треугальяымимо]жерами показсны1моменры срльяык зпышетряменсй
Fig. 2.M^amt^c^nsy)^i^r^]^^ratua^ai^c^r^i^lpsC5TFTanMa3MTpbefyreeyt^(^i^^kyiMM03.P: JanpacilP, 2004 (ар 05 October2008(b);March18, 2011 c).Triangularmarkersshowsthemomentsofstrongearthquakes
менениямивследсовие различияповеденияградиента томпор/туры в слоях/р^иМО/А/о-иО;
уТР=АТ/Ар>0).
Что бы объективно характеризовать наблюдаемые особенности изменч ивости рядов температурнлк аномилсйы 15мополнкловь иреобразованис авеменныл радо а («Ноги ВСЯвР о использоеанеем мооода ытользищово «окна», оотвебш основыш ни предстаятснии колыдиня члена ряръ я нЕясаас дисперсим нткаворого числа (но) пб/лыеяидах знаиений. Размер «опны», как правило, прини ма лся равным т=6, что позволило анализировать присутствующие в спектре колебаний ка роткопериодные вариации температуры [3]. Обусловленные литос-фррвыма е1))(Тцеело1С1в енноамоолвне быстрые —3 днеЯ, ятомааыник )АТ(/-, м м, и р] > > аТ(х ,, уО)ъ ыбнъаии»асычно1е (Б®'1 > 1 и О^ > 1) изменения температуры определялись как произведение значений дисперсии в двух выбранных высотных диапазонах атмосферы:
Б = Б(©ТРТ) • Б(©Ттр) = БРрТ • Б™ .
Паламетр Б зн[ачителы^о преиз1зШ1ает единицу при больших значенигх Баеш нт алзщяя-ственно меньше единицы в отсутствие флуктуаций температуры, когда естественная изменчивость не выходит за пределы среднемесячных значений: АТ(х, у, п, t, р¡) < оТ(х, у), или при воздействии случайоых факооров, когда один из сомножителей мензше адил[ицы ^ 1 или < 1). Поведение этого параметра (Б) определяет положение на временной шкале аномальных вариаций температуры в обоих слоях. Графики, иллюстрирующие все этапы применения методики к временным рядам температуры атмосферы в марте-апреле 2011 г., представлены на рис. 3. Показаны ряды средних значений &ТРТ и &ТТР в слоях 600-300 и 200-150 гПа (рис. 3а), скользящие дисперсии (рис. 3Ь) и их произведение (Б) (рис. 3с), которое коррелирует с сейсмическими данными (рис. 3ф и свидетельствует о резких изменениях температуры накануне двух сильных сейсмических событий. Наряду
Phc. /oicnepciiH(b)
Bcnoflx 600-300 h 200-150 rlla, napaMeTpa ,D(c),Baihacjichhoroc hcimoht30BarnieM c^/jheep KBa/ipara otkjtopphph (Me)hod 0M,pa3mapa (MeKhodl) h kPsM-meeota b apnai^i^n (Metho d °K ce.oKhheckaiiaKthbho c o BMappeoaopejie 2ene.aa0
Fig. 3. Temporal dynamics of anomalies ofaeeragetemperature(a)andm oving variance (b) in the layers 600-300 and 200-150 hPa, the parameter D (c) calculated using average square deviation (Method 1), range (Method 2) and the coefficient of variation(Method3); seismicactivityin March-April2011(d)
с дисперсией (Method 1) рассматривались другие статистические показатели: размах ReT (Method 2) и коэффициент вариации УвТ (Method 3), но сопоставление данных показало преимущество применения дисперсии (D&T) с более надежным выделением аномалий температуры.
Результаты и обсуждение
Диагностирование аномалий температуры атмосферы
Несмотря на специфические особенности динамики развития температурных аномалий накануне различных сейсмических событий, предлагаемая методика позволила определить основные закономерности изменения температуры, характерные для большей части произошедших сильных землетрясений. Результаты диагностирования аномалий температуры атмосферы (параметр D) перед сильными сейсмическими событиями (^>13.0) в разные годы, сведенные вместе и приведенные к одному моменту землетрясения, показаны на рис. 4. Важно отметить, что во всех случаях области высоких значений параметра D либо совпадали с сейсмическими ударами (как в 2013 и 2015 гг.), либо предваряли их (в остальных эпизодах). И не было ни одного случая, когда они проявлялись после землетрясений. Поскольку все точки графиков параметра D фактически прогностические, они, вероятно, характеризуют фазу подготовки землетрясения как период быстрых изменений температуры. Такая синхронизация по времени температурных аномалий и сильных сейсмических событий может означать некоторую общность протекающих при этом физических процессов.
3.5 -3.0 -
2.5 -2.0 -
Q
1.5 -1.0 -
0.5 -0.0 -
Рис. 4. Результаты диагностирования аномалий температуры атмосферы (параметр D) перед сейсмическим! с о бытиятю К>13.0 вразные годы,приведенные к одномумоменту землетрясения
Fig. 4. Results of diagnosis of temperature anomalies of atmosphere (D parameter) before seismic events K>13.0 for various yeaisH shifteCOor thespee enent time
Амнмамии внОнастытиопоодмев1еад эпицентраноеюйеапсзтын продсто-
ящего сильного землетрясения показали большое сходство в плане синхронности возникновения во времени с типичными, в пределах 1-5 дней, интервалами до основного сейс-милеедормтолчкан нрнртедалмеь оЮ случаяхиз14 пдосналенеровемных землевровдееш, что сн^'^^^.^яет е^рммерно 7Л Ы. Крыме para, а^инянв томпзненурь^о наблюаешисе дажер M^J^iBM^^MOCTy^^aepôoj^eehHicoгoднпoсeтhмoскoгooмaоoaеумлeере^е^сомК fK>e0.H0ThfOIм оМнтзом^номтаьсое г^оевмеммет^о^мпературы - достаточно устойчиво проявляющая тенденция.
Вариации температуры и давления атмосферы на уровне тропопаузы
В результате анализа многолетних спутниковых данных было обнаружено, что наиболее сильным землетрясениям (&>10.0) часто предшествовало резкое изменение температуры на уровне тропопаузы, фактически означающее перепады атмосферного давления. Совместный анализ временных рядов аномалий температуры на изобарическом уровне 150 гПа, давления на уровне тропопаузы (pTP) и сейсмической активности в период с 28 ноября 2011 г. по 25 февраля 2012 г. указывает на согласованность изменений процессов в земной коре и атмосфере (рис. 5а, b). Как видно, пять сейсмических событий выше 10.0 энергетического класса (включая K>13.0) происходили после быстрого (в течение 1-3 сут) перепада давления, составляющего Ap=30-120 гПа. Однако при этом не просматривалась связь между величиной изменения давления Ар и энергетическим классом землетрясения.
Такие резкие колебания температуры, а соответственно, и давления являются следствием развития циклонических процессов и распространения атмосферных фронтов, перемещение которых может вызывать многократные квазипериодические деформации земной поверхности, способные спровоцировать землетрясения, назревающие в земной коре независимо от атмосферных процессов [8]. На рис. 5с показана карта пространственного распределения перепада давления на уровне тропопаузы между двумя датами: 03 февраля 2012 г. и днем землетрясе-
- 789 -
P2(2012-02-05)-Pl(2012-02-03)
2 6 10 14 13 22 26 30 3 7 li l5 19 23 27 31 4 3 12 l6 20 24
б lo l4 la 22 26 зо 3 7 ll l5 l9 23 27 3l 4 8 l2 16 20 24 DECEC2ER 2ÍlllJoPJUrR/-)EtRUC[ic 201)
Ire
Рис. 5. Ioариуцои шюмалий темпаратурог уО Л50 hPaí), давленин ica урвете тртшпарзы frropoPress) (a) a ceftcMiieeeKe) пкгги^ностп (i)) в периом с ТВ ноябре У001 т. по МО Невравр 2ti П a. рРРл сунниту да). Пространственноерагпроделениеразностн дтеурнру ртумавну тrюпoцayзы(Мlo^)мPreed) между 03 февуaar п даом земестрясениоОа (аврала Ctlf п (с)
Fig. y.—ariotiansпftemgaracоrа оноп^£^^1ют (©т /15Г |iPn/), tire pressure 1смм1 oatice tacpopoase cBïrBoPress) (a)oTpcain^io acSiaity (T/in Ohe period fromO/Novembcr 20 11im/iicOFebruHry TOll рМ - торай^Те), The
tbe 1е\^Ё^^Н(Ь|с1гарурнр^1г(Т1^с^)^мт^г^^8) Po^wei^r^ 03 РеЬглауа aral Clie day of t/o sarthquako, Fübruaoa CO, ОП12 (cl
z
AT 'T250hPa"T70hPa'
10 -^ 86 -4 -2 -O
44 °N
d 42°N
5 o
40°N
1 4 7 10 13 16 19 22 25 2в 31 december 2006
72°E 74 °E 76°E 78°E
Рис. 6. Вариации разности температуры между уровнями 250 и 70 гПа и сейсмичность в декабре 2006 г. (а)и пространственноераспределение AT 25декабря2006г. (b)
Fig. 6. Variation of the temperature difference between levels 250 and 70 hPa and seismicity in December 2006 (a) and spatial distribution of ДТ on 25 December, 2006 (b)
ния 05 февраля 2012 г. Явно выраженные зоны положительных (Ар=150 гПа) и отрицательных (Ар= -120 гПа) барических градиентов на широтах 50 °N и 30 соответственно, располагались на смежных блоках земной коры симметрично и равноудаленно относительно активных разломов и эпицентра землетрясения.
Диагностирование аномалий пространственного распределения
температуры атмосферы и уходящей длинноволновой радиации
Для локализации области подготовки землетрясений использовались следующие подходы. Во-первых, по профилям температуры Т(р) были идентифицированы случаи наличия двойной тропопаузы [4] и определена разность температур (АТ) между изобарическими поверхностями, соответствующими уровням полярной и тропической тропопаузы. На рис. 6 приведены результаты для самого сильного за последние 24 года сейсмического события класса А=14.8 (магнитудой М=6.5), произошедшего 25 декабря 2006 г. Явно выраженный минимум на временном ряду разности температур был отмечен за двое суток до землетрясения (рис. 6а), что являлось следствием увеличения вертикальных градиентов над и под полярной тропопаузой или положительного тренда контрастности тропопаузы. В то же время поведение температуры на уровне 70 гПа Сшизкосовпадалосвариациями содержания озона вэтойобласти высотсза-паздыванием на ~1 сут.
Изменения теплового поля иакартсх пространственносораспределения темпсратуиыосо-ответствующихэтимуяовням (ЭР^ряЯТр^а^предепялись, как
Ал = \ЛтолкРа (Xк У^к*к П)-лбй№а (X,., у^,г, п)\, где п - дни; t - время измерения; х - широта, у - долгота; I, ] - точки сетки координат исследуемой области с разрешением 1.25° х 1.25°.
В двумерном представлении пространственное расположение минимума разности температур (АТ), соответствующего области максимального перекрытия двух тропопауз, показало изменение положения в виде дуги вокруг эпицентра землетрясения (рис. 6Ь). Эта картина наблюдалась примерно за сутки до землетрясения.
Далее, с учетом большого количества фактического материала, подтверждающего высокую эффективность космического мониторинга теплового излучения Земли и атмосферы при исследовании землетрясений [9], для каждого пиксела (2.5°х2.5°) определялась зональная аномалия уходящей длинноволновой радиации как отклонение текущего значения от среднемесячного уровня за 6-летний период (индекс ОЬН). Вариации среднесуточных значений ОЬЯ ^Е над эпицентральной областью в декабре 2006 г. (рис. 7а) показали колебания с возрастающей амплитудой, которые начались примерно за неделю, а локальный (абсолютный) минимум совпадал с днем землетрясения. Для выделенных на графике овалом дней максимума построены карты пространственного распределения уходящей длинноволновой радиации (рис. 7Ь). Аномалии ОЬЯ, превышающие уровень 2ст, наблюдались вблизи эпицентров за 2 дня до землетрясений. Несмотря на сложный рельеф в регионе, тепловые аномалии были зарегистрированы на сравнительно небольшой площади.
Заключение
Анализ спутниковых данных показал, что возникновению большинства наиболее сильных сейсмических событий (К>13.0), произошедших в Северном и Центральном Тянь-Шане в 2003-2015 гг., предшествовало резкое изменение температуры в зоне раздела тропосферы и стратосферы. Разработанная методика и средства программного обеспечения позволили исследовать особенности трансформации вертикальной термической структуры атмосферы, вы- 791 -
OLR ANOMALY 23 Dec 2006
december 2006 72 °E 76 °E 80 °E
Рис.7.Аномалииуходящейдлинноволновойрадиациивдекабре 2006 г.(а) и распределение аномалии OLR 23 иодабря2006 г. еОлизд эпицентров зодлдтресений, произошедших ОД декаДдд ДО. 9 °N и
74.1ОЕ)и25 декабря(006 г. (К=°4.8;42Д )N е76.0°й)(Ь(
Fig. 7. Anomalies of outgoing long-wave radiationinDecember 2006 (a) and distribution of OLR anomaly on 23 December 2006, close to the epicenters of earthquakes that occurred on 24 December (K=9.3; 40.9 °N and 74.1 °E) and 25 MOR°m3er, 2006(K=14.8; PRO °N andER.0°E) (b)
яв етьоиом абиитбмпе рсту редл у ходящойнлдннололнивоИрадиации,кьтд4ыо, ввроятнл,могут служитьхлдыкдерньшпхизыакомподготовкисейсмичесдлел ообытид ииовсматриелтлоякак краткосрочный предвестник возможного сильного землетрясения. В то же время очевидна необходимость реализации комплексного подхода с привлечением дополнительной информации, других методов и расширением представительности данных.
Авторы выражают благодарность руководителям и сотрудникам научных проектов MERRA-bDM АП^М^лдванхусилия едефлбодке нездльтатол измepвнийнbДeлиорeниe св о-бодного доступа к данным, использующимся в настоящем исследовании.
Отдельная благодарность заведующему лабораторией комплексных исследований НС РАН в г. Бишкек В.Д. Брагину за предоставление цифровых записей сейсмических событий.
Исследование частично осуществлено в рамках выполнения государственного задания НС РАН, тема № 7.
Список литературы
[1] Кашкин В.Б., Романов А.А., Григорьев А.С., Баскова А.А. Тропосферные эффекты землетрясений в Туве, наблюдаемые с искусственных спутников Земли. Журнал СФУ. Техника и технологии, 2012, 5(2), 220-228 [Kashkin V.B., Romanov A.A., Grigoriev A.S., Baskova A.A. Troposphere Effects of Tuva Earthquakes Detected with Space Technology. J. Sib. Fed. Univ. Eng. Technol, 2012, 5(2), 220-228 (in Russian)]
[2] Кашкин В.Б. Внутренние гравитационные волны в тропосфере. Оптика атмосферы и океана, 2013, 26(10), 908-916 [Kashkin V.B. Inner gravity waves in the troposphere. Atmospheric and Oceanic Optics, 2013, 26(10), 908-916 (in Russian)]
[3] Имашев С.А., Свердлик Л.Г. Вариации температуры атмосферы в период высокой сейсмической активности в Японии в 2011 г. Наука, новые технологии и инновации, 2015, 1, 15-19
[Imashev S.A., Sverdlik L.G. Temperature variations of atmosphere during high seismic activity in Japan in 2011. Science, new technologies and innovations, 2015, 1, 15-19 (in Russian)]
[4] Свердлик Л.Г., Имашев С.А. Изменение температурного режима атмосферы под влиянием сейсмической активности. Наука, новые технологии и инновации, 2015, 4, 9-14 [Sverdlik L.G., Imashev S.A. Atmospheric temperature regime change influenced by seismic activity. Science, new technologies and innovations, 2015, 4, 9-14 (in Russian)]
[5] URL: http://giovanni.gsfc.nasa.gov/
[6] URL: http://www.esrl.noaa.gov/psd/cgi-bin/
[7] Свердлик Л.Г., Имашев С.А. Динамика температуры тропопаузы в период сейсмической активности в Центральном Тянь-Шане. Наука, новые технологии и инновации, 2015, 1, 23-27 [Sverdlik L.G., Imashev S.A. Tropopause temperature dynamics during a seismic activity in Central Tien Shan. Science, new technologies and innovations, 2015, 1, 23-27 (in Russian)]
[8] Адушкин В.В., Локтев Д.Н., Спивак А.А. Влияние барических возмущений атмосферы на микросейсмические процессы в земной коре. Физика земли, 2008, 6, 77-85. [Adushkin V.V., Loktev D.N., Spivak A.A. Influence of barometric atmospheric perturbations on microseismic processes in the Earth's crust. Izvestiya, Physics of the Solid Earth, 2008, 6, 77-85 (in Russian)]
[9] Venkatanathan N., Natyaganov V. Outgoing longwave radiations as pre-earthquake signals: preliminary results of 24 September 2013 earthquake (M 7.7). Current Science, 2014, 106(9), 12911297.
