Геология и Геофизика Юга РоссииГеология и Геофизика Юга России2221-3198Федеральное государственное бюджетное учреждение науки Федеральный научный центр "Владикавказский научный центр Российской академии наук"112710.46698/VNC.2025.53.30.003Геотектоника и геодинамикаОсобенности проведения радонового мониторинга в режиме принудительной доставки почвенного воздуха к детектору в регионах с различными геодинамическими условиямиЮрковА.К.akyurkov@mail.ruЗаалишвилиВ.Б.vzaal@mail.ruКозловаИ.А.БирюлинС.В.МельковД.А.БутанаевЮ.В.jyra3@mail.ruРычковаК.М.ИмашевС.А.sanzhar.imashev@gmail.comЛазареваЕ.А.КаменевП.А.Институт геофизики УрО РАН, Россия, 620016, г. Екатеринбург, ул. Амундсена, 100Геофизический институт Владикавказского научного центра РАН, Россия, 362002, г. Владикавказ, ул. Маркова, 93аТувинский Институт комплексного освоения природных ресурсов СО РАН, Россия, 667007, г. Кызыл, ул. Интернациональная, 117аНаучная станция РАН в г. Бишкеке, Кыргызстан, Бишкек-49, 720049Институт морской геологии и геофизики ДВО РАН, Россия, 693022, г. Южно-Сахалинск, ул. Науки, 12403202513244Copyright (c) 2024 Геология и геофизика Юга России2024Геология и геофизика Юга России<p>Актуальность работы. Установление связи геодинамических процессов с данными геологогеофизического мониторинга остается актуальной задачей. К наиболее известным и хорошо исследованным относятся явления, связанные с изменением напряженно-деформированного состояния горных пород, которое влияет на состояние подземной гидросферы и атмосферы. Несмотря на большой объем радоновых исследований, до сих пор не выявлены надежные закономерности поведения объемной активности радона, связанные с процессом подготовки сейсмического события. Цель исследования. Целью работы является изучение особенностей изменения радонового поля (вариаций объемной активности почвенного радона) при подготовке тектонического землетрясения в сейсмически активных регионах с отличающейся геолого-тектонической обстановкой. Методы работы. Мониторинг радона проводился в режиме принудительной доставки почвенного радона к детектору. Интерпретация проводилась по методике, разработанной по результатам, полученным для условий зоны субдукции. Результаты работ. Выполнен анализ вариаций объемной активности радона (ОАР) при подготовке землетрясений, полученных при мониторинговых наблюдениях в регионах России отличающихся различной геодинамической активностью. Подтверждена причинно-следственная связь вариаций ОАР с процессом подготовки землетрясений. Установлено, что получаемые закономерности проявления в кривых объемной активности радона процесса подготовки землетрясений не противоречат ранее полученным. Изучены и обобщены ошибки при организации наблюдений с использованием режима принудительной доставки радона к детектору. Показано, что для получения достоверных результатов, позволяющих сравнивать получаемые закономерности поведения ОАР, необходимо учитывать геологические условия установки станций и соблюдать методику измерений в режиме принудительной доставки почвенного радона к детектору.</p>earthquakeradonvariationsmonitoringforced deliveryforecastземлетрясениерадонвариациимониторингпринудительная доставкапрогноз<p>Введение</p> <p>Практически каждый день по всему земному шару происходят землетрясения, которые остаются одними из самых опасных и до конца не изученных природных явлений. В настоящее время более 30 % территории России относится к сейсмоопасным зонам. В этих районах расположены крупные мегаполисы, объекты атомной промышленности и энергетики, поэтому уменьшение опасности для жизни людей и снижение экономических потерь от катастрофических явлений остается актуальной задачей. Очевидно, что определение точного времени, места и магнитуды события задача чрезвычайно сложная, но даже определение составляющих прогноза землетрясения может уменьшить человеческие потери и экономический ущерб от воздействия природных катастроф. Несмотря на достаточно многочисленные исследования, связанные с изучением предвестников землетрясений – прогноз остается нерешенной задачей. Тем не менее, развитие сейсмологических наблюдений, проведение лабораторных экспериментов и теоретических исследований, накопление большого статистического материала при натурных исследованиях позволяет полнее понять процесс подготовки землетрясения и на этой основе разрабатывать новые методы его прогноза.</p> <p>Предвестниками обычно выступают особенности поведения различных физических полей, не всегда одинаково реагирующих на процесс подготовки землетрясения. Поэтому для успешного прогноза необходимо использовать несколько предвестников, подтвердивших свою эффективность на большом статистическом материале [Завьялов, 2006]. Радон, как индикатор изменений напряженного состояния горных пород, приводящих к тектоническим землетрясениям, широко применяется во всем мире [Уткин, Юрков, 1998; Фирстов, Макаров, 2018; Hatuda, 1953; Kasahara, 1979; King et al., 1996; Teng, 1980]. Возникновение аномальных временных изменений концентраций радона изучалось рядом авторов, специализирующихся на измерениях Rn в почвенном газе [Kuo et al., 2010; Mogro-Campero et al., 1980; Planinic et al., 2001; Ramola et al., 1990; Reddy, Nagabhushanam, 2011], в грунтовых водах [Favara et al., 2001; Gregoric et al., 2008; Singh et al., 1999], в приземной атмосфере [Oh, Kim, 2015; Omori et al., 2007].</p> <p>Получение статистически значимого материала, в первую очередь, связано с количеством происходящих землетрясений. Причем одинаково плохо как малое их количество, так и большое. В последнем случае трудно выделить закономерности подготовки конкретного события. Уникальным объектом для исследований стали Южные Курилы, где частота происходящих землетрясений позволяла проследить подготовку отдельных событий. При поддержке сотрудников Института морской геологии и геофизики ДВО РАН удалось организовать мониторинговые наблюдения объемной активности радона на нескольких Южно-Курильских островах.</p> <p>В период с 2011 по 2023 гг. сотрудниками Института геофизики УрО РАН получены результаты непрерывного радонового мониторинга на островах Кунашир, Итуруп, Шикотан. На основании этих наблюдений была разработана методика интерпретации и установлены статистически обоснованные закономерности поведения ОАР при подготовке землетрясений применительно к зоне субдукции [Бирюлин, 2022]. Для решения вопроса о применимости полученных закономерностей для других регионов, в последние годы была расширена сеть радонового мониторинга с охватом Урала, Северного Кавказа, Сахалина, Саян, Северного Тянь-Шаня, Крыма. Понятно, что для получения конкретных количественных значений для каждого региона необходимо проведение длительного мониторинга ОАР при достаточном количестве происходящих землетрясений. Кроме этого, условием объективного сравнения получаемых закономерностей поведения ОАР при подготовке землетрясений в регионах, различающихся геодинамической активностью, является соблюдение методики проведения измерений в режиме принудительной доставки радона к детектору и оптимальный выбор геолого-тектонических условий для размещения станций наблюдения. Выбор информативных точек наблюдения очень важен, так как вариации ОАР определяются способностью среды к деформациям в точке наблюдения. В данной статье проведем анализ соответствия характеристик работающих станций радонового мониторинга этим требованиям и покажем особенности отражения подготовки землетрясений в поле радона.</p> <p>Объекты и методика исследований</p> <p>Объектом исследований являются вариации объемной активности почвенного радона, полученные в результате проведения радонового мониторинга в различных регионах с отличающейся геодинамической активностью. Методика исследований заключается в проведении анализа отражения сейсмических событий в аномалиях радона, полученных при непрерывных измерениях, выполненных в режиме принудительной доставки почвенного (породного) радона к детектору. Режим принудительной доставки радона из почвы осуществляется с помощью насоса из зонда, установленного на глубину не менее 0.7 м. Данная глубина является наименее градиентной, учитывая вертикальное распределение радона, и на ней существенно снижено влияние на вариации ОАР суточных перепадов температуры и атмосферного давления [Козлова, Юрков, 2005]. Тонкий зонд, перфорированный на конце, не искажает трещинно-поровое пространство и практически не влияет на напряженно-деформированное состояние среды. Применение способа вынужденной доставки радона к детектору существенно увеличивает эффективный объем горных пород, что приводит к «осреднению» их петрофизических характеристик. Как показали исследования по сравнению диффузионного режима измерений и режима вынужденной доставки, регистрируемые аномалии ОАР в последнем случае имеют более высокую амплитуду, небольшую продолжительность (несколько суток), что является важным при краткосрочном прогнозе землетрясений. Для мониторинговых измерений используется недорогая и серийно-выпускаемая аппаратура: регистраторы радона – SIRAD, Radex-107, Radex-107+, насос Sonic-108. </p> <p>При интерпретации полученных результатов анализировались сейсмические события, происходящие в области установки станции наблюдения. Оценка радиуса проявления радонового предвестника землетрясения основана на соотношении магнитуды к логарифму расстояния от эпицентра события до станции наблюдения (геодинамический критерий). По разработанной методике отражение в поле радона подготовки землетрясения возможно и при соотношении М/lgR≥2 [Бирюлин, 2022]. Данная методика позволяет оценить временной интервал проявления землетрясения в поле радона и отнести событие к «дальней» или «ближней» зоне. </p> <p>Результаты исследований и их обсуждение</p> <p>Станция Южно-Курильск (о. Кунашир)</p> <p>Расположена в зоне субдукции, зонд установлен в подвале сейсмостанции п. г. т. Южно-Курильск на глубине 0.7 м в коре выветривания пород андезито-базальтового состава, слагающих монолитный лавовый поток. Условия проведения измерений на этой станции полностью соответствуют режиму принудительной доставки радона к детектору. Выделяемые аномалии объемной активности радона (ОАР) контрастные, с хорошим превышением над фоном. При интерпретации результатов учитывались события с магнитудой более 4 и радиусом не более 500 км, для выполнения геодинамического критерия М/lgR≥2 [Бирюлин, 2022]. За период с 2011 по 2023 гг. было проанализировано 289 событий, удовлетворяющих геодинамическому критерию, и установлено, что 267 предшествовали аномалии ОАР (что составляет 96%). На рисунке 1 представлен пример отражения землетрясений «ближней» (меньше 130 км) зоны.</p> <p><img class="frame-18" src="Геол_журн__№1_2025-web-resources/image/139564.png" alt="139564.png" /> </p> <p>Рис. 1. Фрагмент записи ОАР на Станции Южно-Курильск /</p> <p>Fig. 1. A fragment of the recording of the VRA at the Yuzhno-Kurilsk Station</p> <p> </p> <p>Станция Коуровка (Средний Урал)</p> <p>Средний Урал не относится к сейсмоактивным районам. Редкие ощутимые землетрясения произошли в 1914 г. – Билимбаевское, с магнитудой около 5.5, в 2015 г. – Среднеуральское с магнитудой 4.7, в 2018 г. – Катав-Ивановское с магнитудой 5.6.</p> <p>Станция наблюдений «Коуровка» расположена на массиве монолитных известняков девонского возраста. Зонд установлен на глубине 0.7 м в подвальном помещении наблюдательного павильона. Наблюдения проводятся в режиме принудительной доставки радона к детектору. Так как события происходят редко, то их достаточно легко выделять и соотносить с аномалиями ОАР. Пример поведения ОАР перед землетрясениями приведен на рисунке 2. После прохождения аномалии через 5‒6 дней произошли землетрясения, которые являются «дальними» (расстояние от пункта наблюдений более 300 км), что соответствует ранее выявленным закономерностям. Особенностью проявления землетрясений в кривых ОАР стал геодинамический критерий М/lgR≥1.6, т.е. меньше 2.0. Это может свидетельствовать о том, что для условий Среднего Урала геодинамический критерий, при котором землетрясения отражаются в кривых ОАР, может быть меньше.</p> <p><img class="frame-19" src="Геол_журн__№1_2025-web-resources/image/139578.png" alt="139578.png" /> </p> <p>Рис. 2. Фрагмент записи ОАР на станции Коуровка [Юрков и др., 2023в] /</p> <p>Fig. 2. Fragment of the recording of the VRA at the Kourovka station [Yurkov et al., 2023c] </p> <p> </p> <p>Еще одна уральская станция наблюдений была установлена в г. Катав-Ивановск после землетрясения в 2018 г. Станция «Катав-Ивановск» расположена на вершине антиклинальной складки в коре выветривания известняков. Режим измерений соответствует принудительной доставке радона, но зонд расположен в подвале административного здания в центре города. Большое количество техногенных помех вносит отрицательный вклад в измерения почвенного радона. С октября 2018 г. по настоящее время в радиусе 200 км не произошло землетрясений с магнитудой больше 2, удовлетворяющих геодинамическому критерию больше 1,6. На рисунке 3 приведен пример записи ОАР, аномалии 16.02.24 г., 18.02.24 г. вызваны техническими помехами. Станция оказалась малоинформативна, в связи с отсутствием происходящих землетрясений.</p> <p><img class="frame-20" src="Геол_журн__№1_2025-web-resources/image/%d0%a0%d0%b8%d1%81.%203.%20%d1%81%d1%82%d1%80.%2037_fmt.jpeg" alt="%d0%a0%d0%b8%d1%81.%203.%20%d1%81%d1%82%d1%80.%2037.jpg" /> </p> <p>Рис. 3. Фрагмент записи ОАР на станции Катав-Ивановск /</p> <p>Fig. 3. A fragment of the recording of the VRA at the Katav-Ivanovsk station</p> <p>Станция «Кармадон» (Северный Кавказ)</p> <p>Станция расположена в районе развития вулканогенных отложений на склоне вулкана Казбек. Зонд помещен в кору выветривания в павильоне сейсмостанции на глубину 0.7 м. Условия измерений на станции соответствуют методике измерений в режиме принудительной доставки радона к детектору. Пример отражения землетрясения в поле радона приведен на рисунке 4. Аномалия объемной активности радона (ОАР) достаточно контрастная, с хорошим превышением над фоном, что позволяет легко ее выделить. Землетрясение произошло через 10 дней после появления аномалии и выхода ее на фоновый уровень (26.10.18), что соответствует проявлению «близкого» события. К сожалению, на этой станции за период наблюдений дальние события не проявились в вариациях ОАР. Возможная причина – экранирующее влияние системы хребтов Кавказского орогена и геодинамический критерий событий был меньше 2.0.</p> <p><img class="frame-21" src="Геол_журн__№1_2025-web-resources/image/139608.png" alt="139608.png" /> </p> <p>Рис. 4. Фрагмент записи ОАР на станции Кармадон, Северный Кавказ [Заалишвили и др., 2022]/</p> <p>Fig. 4. Fragment of the recording of the VRA at Karmadon station, North Caucasus [Zaalishvili et al., 2022]</p> <p> </p> <p>Станция Владикавказ (г. Владикавказ)</p> <p>Расположена в цокольном этаже Геофизического Института ВНЦ РАН. Зонд установлен в валунно-галечных отложениях. Режим измерений – принудительная доставка радона. На данной станции, кроме неудачного расположения в валунно-галечных отложениях, присутствует отрицательное воздействие техногенных помех, влияющих на измеряемую величину ОАР (рис. 5). Проявление процесса подготовки землетрясений в кривой ОАР проявляется нечетко. За период с 20.07.2021 г. по 04.09.2021 г. произошло 3 землетрясения. Исходя из геодинамического критерия, они относятся к «дальней» зоне, но при анализе соотнести их с радоновыми аномалиями сложно как по амплитуде, так и по положению на кривой ОАР.</p> <p> </p> <p> </p> <p><img class="frame-22" src="Геол_журн__№1_2025-web-resources/image/139625.png" alt="139625.png" /> </p> <p>Рис. 5. Фрагмент записи кривой ОАР на стации Владикавказ /</p> <p>Fig. 5. A fragment of the recording of the curve of the UAR at Vladikavkaz station</p> <p> </p> <p>Станция Кызыл (Республика Тыва)</p> <p>Расположена в пойменных и русловых отложениях реки Ка-Хем (Малый Енисей). Зонд установлен в грубообломочных пролювиальных отложениях, слабо передающих возникающие деформации. Измерения соответствуют режиму принудительной доставки радона. Пример записи приведен на рисунке 6. Видно, что характер вариаций ОАР затрудняет проведение их интерпретации. Величина аномалий незначительно превышает фоновые значения, что может быть объяснено удаленностью событий. Суточные вариации ОАР нечеткие, геодинамический критерий для дальних событий не достигает величины 2.0, при которой возможно их проявление в кривых ОАР. Практически все произошедшие в этот период землетрясения отмечаются как дальние. Близкое событие 15 сентября 2024 г. ориентировочно можно связать с незначительной аномалией 2‒4 сентября. Хотя по величине геодинамического критерия – 2.43, данное землетрясение должно было уверенно проявится в кривой ОАР. Основные причины нечетких аномалий – неудачное расположение радоновой станции и возможные технические помехи из-за ее размещения в здании Института.</p> <p><img class="frame-23" src="Геол_журн__№1_2025-web-resources/image/139651.png" alt="139651.png" /> </p> <p>Рис. 6. Фрагмент записи кривой ОАР на станции Кызыл /</p> <p>Fig. 6. A fragment of the recording of the curve of the VRA at the Kyzyl station</p> <p>Станция «Иссык-Ата» (Кыргызстан)</p> <p>Станция находится на Северном Тянь-Шане, в зоне складчато-разрывного сочленения Киргизского хребта и Чуйской впадины. Станция установлена в ~1 км от тальвега на правом борту р. Иссык-Ата в разнообломочных пролювиальных средне-верхнечетвертичных отложениях (pQ2-3), лежащих на гранодиоритах ордовика с ксеноблоками верхнепротерозойских сланцев. Породы палеозойского фундамента сильно трещиноваты и с признаками дробления из-за близости (~1.5 км) линии Шамсинского взброса с падением его плоскости сместителя под углом 40‒50° в сторону станции.</p> <p>В подвальном помещении павильона Научной станции РАН установлена квадратная труба сечением 20 мм на 20 мм длиной 1 м с перфорированным основанием, заглублена на 0.6 м, в нее установлен зонд с выходящей из него трубкой с выходом в операторское помещение, где установлен прокачивающий насос и детектор радона. Такая схема установки зонда не соответствует режиму принудительной доставки радона. Наличие дополнительного объема за счет внешней квадратной трубы приводит к разбавлению и усреднению величины ОАР. Связь открытого конца квадратной трубы с помещением обуславливает поступление в нее воздуха с низким содержанием радона и при отключении насоса, детектор регистрирует содержание радона в воздухе помещения. При такой схеме измерений регистрация короткопериодных высокоамплитудных аномалий затруднительна, а получаемые вариации неприменимы для качественной интерпретации.</p> <p>На рисунке 7 приведен фрагмент записи ОАР с 31.01.2023 г. по 19.02.2024 г. Произошедшие за этот промежуток времени землетрясения в ноябре 2023 г. и январе 2024 г., по геодинамическому критерию должны были отразиться в поле радона. Но процесс их подготовки не отразился в вариациях объемной активности почвенного радона, что, скорее всего, объясняется осреднением величины ОАР в дополнительном объеме – квадратной трубе. </p> <p>Резкий спад величины ОАР после землетрясения 22.01.2024 г. вероятнее всего, вызван техническими неисправностями работы насоса или разгерметизацией соединительных элементов от зонда к детектору. На это указывают последующие фоновые значения ОАР.</p> <p><img class="frame-24" src="Геол_журн__№1_2025-web-resources/image/139673.png" alt="139673.png" /> </p> <p>Рис. 7. Фрагмент записи кривой ОАР на станции «Иссык-Ата» /</p> <p>Fig. 7. A fragment of the recording of the VRA curve at the «Issyk-Ata» station</p> <p>Станция «Петропавловская» (о. Сахалин)</p> <p>Пункт наблюдения в с. Петропавловское, расположен в южной части острова, непосредственно в зоне регионального субмеридионального Центрально-Сахалинского разлома. В структурном плане, разлом представляет собой взбросо-надвиг, с падением сместителя на запад под углом 55°. На земной поверхности его трасса ярко выражена наползанием резко пересеченного гористого рельефа, сложенного холмской и невельской свитами, на западный фланг низменной равнины Сусунайской депрессии, заполненной Маруямскими отложениями. Непосредственно в месте установки шпура для отбора газовых проб, породы представлены аллювиально-пролювиальными желтыми и серыми глинами и суглинками современного четвертичного возраста с мощностью толщи до 4 м.</p> <p>Условия измерений соответствуют режиму принудительной доставки радона. За период наблюдений с 17.10.2024 г. по 29.12.2024 г. отразилось только одно землетрясение 23.10.2024 г. с геодинамическим критерием 2.2 на расстоянии 928 км. В кривой ОАР оно отметилось как дальнее событие, на спаде аномалии, что соответствует полученным ранее закономерностям (рис. 8). Аномалии ОАР на этой станции имеют небольшое превышение над фоновыми значениями ОАР, что затрудняет их выделение. Возможной причиной этого может быть неудачное размещение зонда, под которым находится подвальное помещение, за счёт этого возможно искажение деформационной картины при подготовке землетрясения. Наблюдаемая высокая влажность объясняется условиями депрессионной воронки, возникающей из-за подвального помещения.</p> <p><img class="frame-25" src="Геол_журн__№1_2025-web-resources/image/139689.png" alt="139689.png" /> </p> <p>Рис. 8. Фрагмент записи ОАР на станции «Петропавловская», о. Сахалин /</p> <p>Fig. 8. Fragment of the recording of the VRA at the Petropavlovsk station, Sakhalin Island</p> <p>Выводы</p> <p>Станции Южно-Курильск, Петропавловская, Коуровка и Кармадон расположены в коре выветривания монолитных горных пород, реагирующих изменением структуры трещинно-порового пространства при подготовке землетрясения, что отражается в кривой ОАР. Условия измерения полностью соответствуют режиму принудительной доставки почвенного радона к детектору. Отражения землетрясений в кривых объемной активности радона (ОАР), получаемых на этих станциях, укладываются в закономерности для «ближней» и «дальней» зон, полученные для зоны субдукции. </p> <p>Станции в Кызыле (Саяны) и Владикавказе расположены в грубообломочных породах, плохо передающих деформации, и даже при измерениях в режиме принудительной доставки радона к детектору, не удается получать качественные кривые ОАР, пригодные для выделения закономерностей отражения в них процесса подготовки землетрясений. </p> <p>Станция в Иссык-Ате не соответствует методике измерений в режиме принудительной конвекции из-за наличия дополнительного объема (объема наружной трубы) в системе доставки почвенного радона к детектору. Эта причина и учитывая то, что используемый насос для создания режима принудительной доставки радона, имеет незначительную производительность, не позволяет проследить в кривых ОАР процесс подготовки землетрясений. </p> <p>Таким образом, при выборе места установки станций радонового мониторинга, необходимо учитывать геолого-тектонические условия точки мониторинга [Маркин, 2023; Zaalishvili et al., 2024]. Обустройство станции и проведение измерений необходимо осуществлять строго в соответствии с методикой режима принудительной доставки радона к детектору. </p> <p>Правильно работающая сеть станций радонового мониторинга позволит получить количественные характеристики закономерностей отражения процесса подготовки землетрясения в вариациях объемной активности радона в различных геодинамических условиях.</p>1.Бирюлин С.В. Пространственно-временные закономерности проявления аномалий объёмной активности радона перед тектоническими землетрясениями в районе Южных Курил. Автореф. дисс. канд. геол.-минерал. наук: 25.00.10 / Бирюлин Сергей Викторович. – Екатеринбург: ИГФ УрО РАН, 2022. – 24 с.2.Заалишвили В.Б., Бирюлин С.В., Козлова И.А., Мельков Д.А., Юрков А.К. Закономерности проявления аномалий объёмной активности радона перед землетрясениями в различных сейсмоактивных регионах. // Геология и геофизика Юга России. – 2022. – Т. 12. No 1. – С. 35–50. DOI: 10.46698/VNC.2022.53.28.003.3.Завьялов А.Д. Среднесрочный прогноз землетрясений: основы, методика, реализация. М.: Наука, – 2006. 242 с.4.Козлова И.А., Юрков А.К. Методические вопросы измерения содержания радона-222 в почвенном воздухе при мониторинговых наблюдениях. // Уральский геофизический вестник. – 2005. – No 1. – С. 30–34.5.Маркин М.Ю. Оценка объемной активности радона тектонических разломов в грани-цах городской агломерации Усть-Каменогорск (Республика Казахстан). // Геология и геофизика Юга России. – 2023. – Vol. 13. No. 2. – pp. 29–39. DOI: 10.46698/VNC.2023.93.88.003.6.Уткин В.И., Юрков А.К. Динамика выделения радона из массива горных пород как краткосрочный предвестник землетрясения. // Доклады РАН. – 1998. – Т. 358. No 5. – С. 675–680.7.Фирстов П.П., Макаров Е.О. Динамика подпочвенного радона на Камчатке и сильные землетрясения. Петропавловск-Камчатский: КамГУ им. Витуса Беринга, 2018. – 148 с.8.Favara R., Grassa F., Inguaggiato S., Valenza M. Hydrogeochemistry and stable isotopes of thermal springs: earthquake-related chemical changes along Belice Fault (Western Sicily). // Applied Geochemistry. – 2001. – Vol. 16. No. 1. – pp. 1–17.9.Gregoric A., Zmazek B., Vaupotic J. Radon concentration in the thermal water as an indica-tor of seismic activity. // Collegium Antropologicum. – 2008. – No. 32. – pp. 95–98.10.Hatuda Z. Radon content and its change in soil air near the ground surface. // Memoirs of the College of Science, University of Kyoto. – 1953. – Ser. B. Vol. 20. – pp. 285–306.11.Kasahara K. Migration of crustal deformation. // Tectonophysics. – 1979. – Vol. 52. No. 1-4. – pp. 329–341.12.King C.-Y. King B.-S., William C.E., Zhanget W. Spatial radon anomalies on active faults in California. // Applied Geochemistry. – 1996. – Vol. 11. Issue 4. – pp. 497–510.13.Kuo T., Su C., Chang Ch., Lin Ch. et al. Application of recurrent radon precursors for fore-casting large earthquakes (Mw > 6.0) near Antung, Taiwan. // Radiation Measurements. – 2010. – Vol. 45. No. 9. – pp. 1049–1054.14.Mogro-Campero A., Fleischer R.L., Likes R.S. Changes in subsurface radon concentra-tion associated with earthquakes. // Journal of Geophysical Research. – 1980. – Vol. 85. No. NB6. – pp. 3053–3057.15.Oh Y., Kim G. A radon-thoron isotope pair as a reliable earthquake precursor. // Scientific Reports. – 2015. – No. 5. – Art. No. 13084. DOI: 10.1038/srep13084.16.Omori Y., Yasuoka Y., Nagahama H., Kawada Y., et al. Anomalous radon emanation linked to preseismic electromagnetic phenomena. // Natural Hazards and Earth System Sciences. – 2007. – No. 7. – pp. 629–635.17.Planinic J., Radolic V., Lazanin Z. Temporal variations of radon in soil related to earth-quakes. // Applied Radiation and Isotopes. – 2001. – Vol. 55. No. 2. – pp. 267–272.18.Ramola R.C., Singh M., Sandhu A.S., Singh Sh.V. et al. The use of radon as an earthquake precursor. // Nuclear Geophysics. – 1990. – Vol. 4. No. 2. – pp. 275–287.19.Reddy D.V., Nagabhushanam P. Groundwater electrical conductivity and soil radon gas monitoring for earthquake precursory studies in Koyna, India. // Applied Geochemistry. – 2011. – Vol. 26. No. 5. – pp. 731–737.20.Singh M., Kumar M., Jain R.K., Chatrath R.P. Radon in ground water related to seismic events. // Radiation Measurements. – 1999. – Vol. 30. No. 4. – pp. 465–469.21.Teng T. Some recent studies on groundwater radon content as an earthquake precursor. // Journal of Geophysical Research. – 1980. – Vol. 6. No. 85. – pp. 3089–3099.22.Zaalishvili V.B., Melkov D.A., Martyushev N.V. Klyuev R.V., Kukartsev V.V., Konyuk-hov V.Yu., Kononenko R.V., Gendon A.L., Oparina T.A. Radon emanation and dynamic process-es in highly dispersive media. // Geosciences (Switzerland). – 2024. – Vol. 14. No. 4. – pp. 102. DOI 10.3390/geosciences14040102.