中国地震  2023, Vol. 39 Issue (2): 243-259
同位素地球化学在地震监测预测研究中的应用进展
周启超1,2, 蒋雨函2, 高小其2, 汪倩2, 彭莱1,2     
1. 中国科学院大学, 应急管理科学与工程学院, 北京 100049;
2. 应急管理部国家自然灾害防治研究院, 地壳动力学重点实验室, 北京 100085
摘要:随着新观测技术和理论的进一步发展, 同位素地球化学方法在地震监测预测研究中发挥了越来越重要的作用。通过同位素地球化学方法确定地下流体来源, 研究地下流体循环特征, 分析地震前兆异常的成因, 评估地质构造活动的程度, 开展地震预测研究。同位素示踪技术还可以结合深源流体监测和地球物理方法, 揭示地震孕育、流体与震源之间的关系。此外, 同位素地球化学还可以构建断裂带流体地球化学背景特征, 用于地震监测点映震效能的评估, 提高地震监测预测的准确性, 为地震新监测点的布设和震情跟踪提供技术支撑。通过对现今同位素在地震监测预测中所使用的方法、技术及国内外应用情况的总结分析, 力图全面认识同位素地球化学在地震监测预测应用中的现状及发展趋势。
关键词地震    地下流体    地震监测    地震预测    同位素地球化学    
Recent Development of Isotope Geochemical Application in Earthquake Monitoring and Forecasting
Zhou Qichao1,2, Jiang Yuhan2, Gao Xiaoqi2, Wang Qian2, Peng Lai1,2     
1. School of Emergency Management Science and Engineering, University of Chinese Academy of Sciences, Beijing 100049, China;
2. Key Laboratory of Crustal Dynamics, National Institute of Natural Hazards, Ministry of Emergency Management of China, Beijing 100085, China
Abstract: With the development of new observation techniques and theories, isotope geochemical method is playing an increasingly important role in earthquake monitoring and forecasting studies, in which the method can be used to determine the source of subsurface fluids, to study the characteristics of subsurface fluid circulation, to analyze the causes of earthquake precursor anomalies, to assess the geological tectonic activity, and to conduct earthquake forecasting studies. Isotope tracing technology can also be combined with deep-source fluid monitoring and geophysical surveys to reveal the relationship among earthquake genesis, fluids and earthquake sources. In addition, isotope geochemistry can construct fluid geochemical background characteristics of fracture zones, which can be used for the evaluation of seismic monitoring sites to reflect seismic effectiveness, improve the accuracy of earthquake monitoring and prediction, and provide technical support for the deployment of new earthquake monitoring sites and earthquake activity situation tracking. By summarizing and analyzing the methods, techniques and applications of isotopic analysis in seismic monitoring and forecasting at home and abroad, we try to get a relatively comprehensive understanding of the current situation and development trend of isotope geochemistry in seismic monitoring and forecasting applications.
Key words: Earthquake     Underground fluid     Seismic monitoring     Earthquake prediction     Isotope geochemistry    
0 引言

同位素在地震研究领域的应用始于20世纪50年代初。地震的孕育和发生过程必然伴随着元素和同位素在流体中的迁移、分化和演化。由于地下流体携带了深部地球化学信息,利用同位素监测地下流体对地震监测具有较好的效果(Soldati et al,2020杜建国等,2000)。自1966年邢台地震以来,我国开始对地下流体进行监测。近30年来,随着科学技术的发展,同位素地球化学方法在地震监测和预测中的应用越来越广泛(Ramola et al,2008孙小龙等,2020)。尤其是随着分析化学仪器技术的快速发展,同位素示踪技术成为地震监测的重要手段,广泛应用于解释孕震介质演化机理(Apollaro et al,2020Caracausi et al,2022)和地震监测预测理论和方法的研究(Bernard,2001杜建国等,2003)。虽然地震孕育过程中伴随着地球化学异常变化,但地下流体的产生、循环等过程较为复杂,且存在多种与地震前兆相关的机理。因此,获取可靠的地震前兆异常信息和地震监测方法并排除观测环境变化所引起的干扰相对较为困难。

同位素地球化学在地震监测中的应用能够使地下流体监测数据在地震研究中发挥更重要的作用。国内外学者开展了一系列同位素示踪技术的地震监测预测研究,应用同位素示踪技术来研究地震活动性关系(Wakita et al,1980O’Neil et al,1981)。由于地震孕育过程中应力应变作用下,不同同位素组成的地下水混合交换形成同位素分馏,导致地下流体同位素组成变化(Skelton et al,2019)。通过确定地下水补给来源和循环,可以研究地震前兆形成机理,判识地震前兆异常,为地震预测提供科学判据(Claesson et al,2004Chen et al,2014Zhang et al,2021)。同位素示踪技术还可以用来研究断层气的来源、运移和赋存(Bräuer et al,2003Lombardi et al,2010Sciarra et al,2017)。结果表明,深源流体的释放与地震活动性相关,对判识地震前兆信息、预测地震具有重要意义(王基华等,1998)。以1995年日本神户地震为例,位于震中附近的放射性同位素研究所观察到主震发生前1.5~2个月大气中的氡浓度异常增加(Yasuoka et al,1997)。此外,利用氢氧稳定同位素示踪技术对断裂带温泉循环深度进行研究,结果显示地震活动性与温泉水循环深度相关,循环深度与强震孕育可能性呈负相关关系(王博等,2022)。同时,放射性同位素示踪技术也可以用于地下水或沉积物测年,研究地下水补给、径流特征、可更新能力,或建立古地震模型研究地震断层在足够多的地震间隔中的表面破裂行为,以确定地震复发、断层破裂行为和地震危险的特征(Zhang et al,2014Streig et al,2020)。研究表明,地下水年龄与地下水映震效能有关,可判断地震前兆类型,为地震监测点的选择提供依据(Sano et al,1998Wang et al,2005)。

在综合国内外相关领域的资料和研究成果的基础上,本文概述了同位素地球化学理论和方法,并阐述了前人如何运用同位素示踪技术进行研究,包括地下流体来源和循环、利用地下水年龄和水文地球化学特征来评估地震监测点的映震效能、判断地震活动性以及识别地震流体前兆异常。此外,本文还探讨了同位素地下流体监测在地震监测和预测中的应用,并对其未来发展方向进行了展望。

1 地震研究领域中的同位素地球化学方法 1.1 氢氧稳定同位素

氢氧稳定同位素广泛用于分析地下水成因和补给来源。地下水的氢氧同位素特征因成因不同而异(车用太等,2006)。成因为大气降水渗入的地下水,其δD与δ18O值遵循全球大气降水线(GMWL),区域大气降水线(LMWL)是了解当地地下水同位素相对于当地大气降水的重要依据(Craig,1961Clark et al,2013)。受同位素分馏、高程、纬度、陆地和季节效应的影响,大气降水中的氢氧同位素比值在不同条件下存在明显差异,与大气降水线偏离,因此可以通过氢氧稳定同位素的异常变化幅度来评估水-岩反应的强度、地下水的混合作用和地下流体的补给来源等(Pang et al,2017张磊等,2016)(图 1)。

图 1 用于研究水文和水-岩相互作用过程的同位素地理指示器(据Pang等(2017))
1.2 碳、氦的稳定同位素

地下气体具有较强的映震能力,加强对其监测有助于地震监测预测的研究。其中,CO2是深源流体中分布最广、含量最丰富的化学组分之一(方震等,2012Pierotti et al,2015)。深部CO2释放至地表时,其原始碳同位素组成得以保留,因此可以利用不同成因中的δ13C值来判定其来源与释放机制(上官志冠,1995)。断裂带CO2主要来源于幔源、碳酸岩变质、有机物分解和地表生物活动,其δ13C值与地下流体来源对应见表 1(Lombardi et al,2010申宝剑等,2010)。深部CO2可来源于无机物质、生物有机物质或与空气混合,形成三端元混合模型。利用同位素质量平衡原理,可以半定量计算各端元所占比例(刘洁等,2022),即

$ f_A+f_B+f_C=100 \% $ (1)
表 1 气体来源地球化学判识指标

其中,fA为空气所占比例,fB为深源所占比例,fC为生物所占比例。

$ V_A \times f_A+V_B \times f_B+V_C \times f_C=V_{\text {sample }} $ (2)

其中,Vsample为CO2气体样品体积浓度(%),VAVBVC分别为空气、深源、生物的体积浓度(%)。

$ \left(\delta{ }^{13} C_A \times f_A \times V_A+\delta{ }^{13} C_B \times f_B \times V_B+\delta{ }^{13} C_C \times f_C \times V_C\right) / V_{\text {sample }}=\delta^{13} C_{\text {sample }} $ (3)

其中,δ13Csample(‰)为气体样品碳同位素组成,δ13CAδ13CBδ13CC分别为空气、深源、生物的碳同位素组成(‰)。

联立式(1)、(2)、(3)可得三个端元所占比例。

氦由于其稳定的化学性质,能够准确反映地壳内部物质运移。较轻的3He保留在地球内部,通过火山活动和构造活动释放出来;而U与Th的放射性衰变在地壳中产生较重的4He,分布在整个岩石圈中。因此,通常使用R/Ra(样品3He/4He与大气3He/4He之比)表示氦同位素组成。利用3He/4He在地球不同圈层中的差异性(表 1),可以判识幔源流体,研究地下流体来源等问题(Sano et al,1995Marty et al,1999Xu et al,2013Zhao et al,2021)。氦同位素组成源区百分比计算公式如下

$ \left({ }^3 \mathrm{He} /{ }^4 \mathrm{He}\right)_{\text {sample }}=\left({ }^3 \mathrm{He} /{ }^4 \mathrm{He}\right)_{\mathrm{m}} \times M+\left({ }^3 \mathrm{He} /{ }^4 \mathrm{He}\right)_{\mathrm{c}} \times C+\left({ }^3 \mathrm{He} /{ }^4 \mathrm{He}\right)_{\mathrm{a}} \times A $ (4)
$ 1 /\left({ }^4 \mathrm{He} /{ }^{20} \mathrm{Ne}\right)_{\text {sample }}=M /\left({ }^4 \mathrm{He} /{ }^{20} \mathrm{Ne}\right)_{\mathrm{m}}+C /\left({ }^4 \mathrm{He} /{ }^{20} \mathrm{Ne}\right)_{\mathrm{c}}+A /\left({ }^4 \mathrm{He} /{ }^{20} \mathrm{Ne}\right)_{\mathrm{a}} $ (5)
$ M+C+A=1 $ (6)

式中,MCA分别表示地幔、地壳、大气的贡献百分比。

取(3He/4He)a=1Ra,(4He/ 20Ne)a=0.318;(3He/4He)m=6.5Ra,(4He/ 20Ne)m=1000;(3He/4He)c=0.02Ra,(4He/ 20Ne)c=1000。联立式(4)、(5)、(6)得到

$ A=0.318 /\left({ }^4 \mathrm{He} /{ }^{20} \mathrm{Ne}\right)_{\text {sample }} $ (7)
$ C=1-0.2279 /\left({ }^4 \mathrm{He} /{ }^{20} \mathrm{Ne}\right)_{\text {sample }}-\left({ }^3 \mathrm{He} /{ }^4 \mathrm{He}\right)_{\text {sample }} / 1.0985 $ (8)
$ M=1-(A+C) $ (9)

其中,(3He/4He)sample的单位为10-5

1.3 放射性同位素

根据放射性同位素衰变原理,可以确定地下水的年龄,而地下水年龄可以用于解决地震地下流体监测预测中的一些科学问题。放射性同位素示踪法主要有3H法、3H/3He比值法、14C法、85Kr法、39Ar法和36Cl法等,其中应用较多的是3H法、3H/3He比值法和14C法。单纯使用3H测定地下水年龄很困难,因此产生3H/3He比值法,用于测定50a内补给的地下水年龄(Bayari et al,2006Liu et al,2014)。刚渗入地下的大气降水14C含量与大气保持平衡,随着时间推移,地下水中的14C含量不断衰变减少,因此可利用地下水中的14C含量换算出地下水年龄,基于放射性衰变原理测定地下水年龄的表达式为

$ t=\frac{1}{\lambda} \ln \frac{P}{P_0} $ (10)

式中,λ为同位素衰变常数;P0为补给时放射性同位素浓度;P为计算时放射性同位素浓度。

在地震监测中,222Rn是主要观测对象,因为其具有很强的扩散和迁移能力,当存在压力差、温度差和梯度差时,222Rn会从下往上运移(Sukanya et al,2022)。由于岩层岩性和构造环境的不同,氡和汞气在地下形成不同的聚集,而氡更易被用于地震监测,因为其可以溶于水和有机物,且容易被活性炭和黏土吸附,在气态形式和溶解态下均能运移(Morales-Simfors et al,2020),氡的化学稳定性使其能够准确反映地壳内部物质的运移情况,从而有助于识别地震前兆(Chaudhuri et al,2013Muto et al,2021)。但222Rn本身并不能用来预测地震,因为其在生命周期内会因扩散性输入而增加,因此有学者提出利用石灰岩洞空气中的钍(220Rn,半衰期55.6s,平均寿命80s)记录地震信号而不受环境因素干扰,在2011年日本9.0级东北大地震前观察到222Rn和220Rn活动的异常峰值被认为是地震前兆(Oh et al,2015)。

1.4 其他同位素

Sr同位素具有不易分馏的特点,因此其变化可反映物质来源及迁移过程,被广泛应用于研究地表与地下水系统的沉积物物质来源、水-岩反应、水体混合以及地震活动性等问题,主要表现为87Sr/86Sr变化(陈晨等,2020McHugh et al,2020)。另外,汞同位素示踪技术可用于研究断裂带中汞的来源、汞与深部流体活动的关系,探讨断裂带深部汞及同位素的迁移路径,以及研究断裂带介质中汞异常的来源等问题,为相关研究提供了新的方法(Sherman et al,2010Gehrke et al,2011)。

2 地震活动性判定

同位素地球化学和流体元素通常在大地震前后出现异常(Ingebritsen et al,2014Chiodini et al,2020),在此基础上进行地震流体地球化学研究,探讨地震流体地球化学特征与地震活动之间的作用关系(Skelton et al,2014)。同位素组成对地震孕育和发生的动力学响应特征,使流体地球化学预测地震成为可能(李营等,2022)。在印度西部柯依那地区的研究中,研究人员对地下水进行了水化学特征监测和地震前兆研究,发现2015年第一季度出现了水化学异常,由于地震前应力增加导致不同含水层水的混合,δ18O值随时间递减,水化学参数随时间递增,综合分析认为这是地震前兆异常,并预测该地区即将发生5级以上地震(Reddy et al,2017)。在日本地区的研究中,研究人员发现2016年熊本7.3级地震前3~4个月,当地δ18O值升高了0.38‰,震后上升至0.51‰,而δD值在期间没有明显变化,认为这是因为地震促使含水层岩石产生破裂加强了水-岩作用,进而导致水和岩石之间的氧交换,从而在δD值不变的情况下使得δ18O值上升(Sano et al,2020b)。类似的现象也出现在2016年日本西南部鸟取6.6级地震前后。对含水层岩石进行破碎实验、岩石氦脱气和δ18O值位移表明,氧同位素异常与应力变化直接相关(Onda et al,2018)。以意大利中部的苏尔莫纳和马泰塞地区为2个监测节点进行水文地球化学分析,研究发现,虽然在2017—2020年期间发生了3次地震(震级分别为MW3.9、MW4.1和MW4.4),但除MW4.1事件外,这些地震并没有出现明显的水文地球化学异常。然而,MW4.1地震发生在地壳深处(约17km)的一条大型正断层上,而且在地震前出现了微弱但可探测的硼和δ11B水文地球化学异常,这些异常与2016年中部亚平宁山脉MW6.0地震前的异常相似但较为温和。因此,研究结果表明,B和δ11B具有地震敏感性,震前水文地球化学异常主要出现在MW≥5.0的地震前,而长而深的断层等特殊条件下也可能会出现微弱的异常(Franchini et al,2021)。

地热流体中的幔源氦同位素值常常在大地震前后发生变化(Klemperer et al,2013),一些研究表明,深层地下水中氦同位素异常变化与断层体积应变变化存在定量关系,可以在无钻孔应变仪的区域监测到地震前的应变变化(Sano et al,2016)。Zhou等(2017)研究了理塘MS5.1地震前的理塘断裂带沿线温泉气体地球化学时空变化,发现温泉中3He/4He比值为0.06~0.93Ra,断裂带中含有高达11.1%的幔源氦。禾尼泉、德达泉、措拉泉和毛垭地震台的3He/4He比值和δ13CCO2值在地震前发生了明显的变化(图 2),表明上涌的地幔流体进入地壳对青藏高原东南部的走滑断层活动起着重要作用(Zhou et al,2017),显示该断层作为深部地幔流体通道。因此,监测地热流体中的氦同位素可以提供有价值的信息,进而预测大地震的发生。

图 2 Rc/Ra值(a)、δ13CCO2值(b)时间变化与地震关系图(c)(据Zhou等(2017))
3 地下流体中的同位素与构造活动性

地下流体同位素特征是分析活动断层中流体耦合机制和流体对地震影响的有效工具(Buttitta et al,2020)。水-岩反应作用会削弱岩石抗压和抗剪强度以及断裂摩擦力(Wästeby et al,2014)。在不同的构造和水文地质背景下,与地震活动相关的水化学变化在世界范围内均有报道,以意大利中部亚平宁山脉为背景,2009年拉奎拉序列期间,观测到热泉中溶解气体成分变化及断裂带附近氡异常(Bonfanti et al,2012Quattrocchi et al,2012),研究者发现放射性地壳气体Ar和He的强烈增加,表明地震区深层长期储存着高压气体(Chiodini et al,2011),超压气体储层可能引发地壳深处的地震(Chiodini et al,2020)。2016—2017年亚平宁中部地震序列中,研究者将碳同位素数据与Na/Li地温计数据结合,发现在地震期间呈现出碳同位素响应的泉水具有类似北亚平宁山脉深部钻孔的CO2分压,这表明这些泉水可能受到了与地震有关的深部CO2混合效应的影响(Barbieri et al,2020)。地下流体被认为与断层弱化和孔隙压力增大有关(Kennedy et al,1997Claesson et al,2007),且可能在已经承受严重压力的脆性地壳中触发地震(Sano et al,2020a)。在圣安德烈斯断层区,Lee等(2016)提出断层削弱与富含3He的地幔流体有关。许多研究表明地震活动性与地幔流体有关,在纽波特—英格尔伍德断层报告了深层地震和高3He/4He比率之间的关联性(Inbal et al,2016)。因此,地下流体同位素特征对于描述断裂内部流体循环和地热构造特征的作用至关重要。研究发现秦岭北缘断裂带的温泉水循环深度与构造活动性密切相关,在甘肃段与青海段,温泉水循环深度浅,水的弱化作用小,容易积累应力成为强震危险区;相反,陕西段的温泉水循环深度较大,在强弱化作用的影响下不会积累大量应力难以孕育强震(汪万红等,2008),因此,温泉水循环深度是影响地震危险性的关键因素之一。

大地热流是约束地质构造与构造活动的重要指标,能够记录深部热学信息(Tardani et al,2016)。相关研究发现智利安第斯山脉南部火山区(SVZ)的热流体的特点是有广泛的氦、氮和碳同位素值,断裂系统对氦、氮和碳同位素特征有着控制作用,同时也促进了流体的深层循环(Tardani et al,2016)。在地热梯度带区域,由于两侧地层物性的差异,容易发生高频地震(唐晗晗等,2020),而局部熔融体和上地壳底部的高温水或热液可能是导致地震的直接热力源(张汝惠等,1991)。氦同位素研究有助于识别深源流体(地幔/地壳)逸出气(Burnard et al,2013),评估由岩石圈热结构不均而产生的热应力和流体对岩石的弱化作用对强震的孕育过程影响(王云等,2018Umeda et al,2013)。大地热流值包含了上地幔和地壳热的贡献,而地下流体中的氦同位素比值来自于上地幔和地壳的放射性成因。因此,氦同位素比值与大地热流值之间存在着密切的关系(Polyak et al,1979),可表示为

$ q=18.231 \times \lg R+181.82 $ (11)

式中,R为测试样品的氦同位素比值,q代表大地热流值(mW/m2)。

土耳其安纳托利亚地区长期的地震活动源于其构造环境和地热系统。地热流体来源于大气降水,在断层系统中循环积累能量,其年龄大于60年(Baba et al,2022)。氢氧稳定同位素数据表明,地热流体具有对流上升和将上地壳热量传导到地表的特征,而氚同位素数据表明热水具有深部循环的特征(Aydin et al,2020Pasvanoǧlu,2020)。西部地区的同位素数据显示,地热流体主要是由海水的局部补给和深循环形成的。交错断层对地热流体的上升和沉积非常有利(Uzelli et al,2017)。研究发现深部岩浆室周围热区的流体和气体输送机制受到侵入岩浆的深度和上层地壳断层的渗透性的影响(Karao lu et al,2019)。西南部德尼兹利盆地的地幔挥发物可能上涌到软流圈,并沿着构造不连续面上升(Alçiçek et al,2019)。因此,分析构造环境和地热系统之间的关系,评估该地区地热流体的水文地球化学特征,对于了解该地区长期的地震活动具有重要作用。

研究断裂带活动性的重要指标之一是断裂带气体(He、CO2、CH4)浓度变化(Chiodini et al,2011)及其同位素特征(3He/4He、δ13CCO2δ13CCH4)(Kulongoski et al,2013孙小龙等,2016Sciarra et al,2017),这些特征能够刻画气体在迁移过程气体特征的改变程度(Bräuer et al,2016)。有研究表明幔源氦贡献与气体CO2浓度正相关(Cornides et al,1986),δ13C值与3He/4He比值在断裂带呈正相关关系,且相当一部分逸出的CO2来自岩石圈地幔(Palcsu et al,2014),其可以归因于地幔流体的高压对断层的削弱作用(Lee et al,2019)。圣安德列斯断层系统(SAFS)中的纽波特—英格伍德断裂带(NIFZ)深部的氦同位素比值高达5.3Ra,有强烈的地幔氦信号,研究结果表明NIFZ是唯一与深部地幔相连的断层,是该地区地幔流体运移与热量输送的主要通道(Boles et al,2015)。Sano等(2014)在研究东日本大地震后的地震活动性时发现,震后一个月断裂带附近底层海水的3He/4He高于背景值,推测是受地幔影响的流体和高氦比流体在构造活动中被释放出来,这些高压流体可能促进了地震的产生。断裂带的温泉气体地球化学特征与构造活动性有紧密的联系,其稳定同位素(δ13C和3He/4He)在地震前后会发生变化(Chen et al,2015)。刘洁等(2022)则研究了陇县—宝鸡断裂带的CO2气体浓度和碳同位素特征,结合三端元混合模型发现南段岐山—马召断裂气体渗透性与地下构造连通性强,而北段固关—虢镇断裂气体活动性弱、断层闭锁程度高,未来发生地震的危险性高。研究发现同位素值的时间和空间变化与构造活动性相关(Du et al,2006)。在对金沙江—红河断裂带的温泉气体地球化学特征研究中,该区域的H2浓度出现峰值(图 3)。由于红河断裂南段深部流体上涌对断裂的弱化作用,该区域的构造活动性与校正氦同位素值呈上升趋势,说明幔源物质上涌与断裂活动性密切相关,并且深部流体对地震孕育具有促进作用(周晓成等,2020)。

图 3 金沙江—红河断裂地震活动性与幔源氦比例的关系(据周晓成等(2020))
4 地震过程分析中的非构造影响

在地震发生过程中,地壳表层的水、气体、土壤等介质受到应力应变作用的影响,其元素含量和同位素组成会出现明显的异常变化,例如断裂带土壤气体(CO2、CH4、Rn、Hg等)、断裂带附近温泉水、井水和地下水的水化学离子(Na+、K+、SO42-、Cl-、Na+、TDS等)、温泉气体(H2、He、Rn、CO2等)的浓度以及温泉气体中稳定同位素的比值(3He/4He和δC13CO2)等。判识地震因素和环境因素引起的地下水化学异常变化、排除干扰,是地震监测和预测研究中的关键问题(张国盟等,2015)。水化学组分异常通常与地下流体运移过程中的水-岩反应有关。Sato等(2020)研究了2016年熊本地震后泉水中硝酸盐浓度变化异常现象,结果表明主震后土壤中的硝酸盐通过强降雨渗入含水层,且强震影响植物吸收的硝酸盐含量,二者结合导致泉水中硝酸盐浓度增加(Sato et al,2020)。因此,结合不同的特征可以识别异常现象,并进一步研究其与地震活动的关系。水温异常反映地下深部热流活动的状态。以2022年云南宁蒗MS5.5地震为例,该地震前迪庆流体观测井水位和水温同时出现上升异常(图 4),利用氢氧同位素特征和离子浓度特征综合分析表明此次异常反映了该区域深部构造活动(杨黎等,2022)。

图 4 2011—2022年迪庆井水位、水温观测曲线(据杨黎等(2022))

较大地震后可能会出现新异常,需要判别是震后效应还是新地震的前兆异常。震后效应是地震发生后应力场向新的平衡态过渡中出现的与已发生地震有关的现象,包括已发生地震导致的异常和余震区应力调整;新地震的前兆异常是地震发生后,原有异常出现新变化或新异常,这些异常可能与后续的强震有关。2014年的南纳帕地震引起之前干涸的泉眼重新涌现水流,研究表明泉水中的δD和δ18O值显著低于地震前水平,排除浅层水补给的可能性,研究认为这是地震引起的渗透性增强所导致的同震响应,综合分析认为泉水来自于震后释放的附近山区地下水(Wang et al,2015)。2016年熊本7.0级地震后,许多研究加深了对地震如何在区域范围内影响地下流体系统的理解,震后区域内地下水最大异常水位上升11m,研究发现震后山区含水层水被释放并与断裂带含水层的地下水混合,山脚泉水δD和δ18O值更低,意味着有来自高海拔水补给(Hosono et al,2020b);Ide等(2020)使用氟氯碳化合物(CFCs)对地下水年龄进行追踪,推断观察到的变化是由于当地山脉的地震断裂导致的渗透性增强,从而释放出附近山脉的水。有研究者结合地球化学和同位素示踪技术对2016—2017年亚平宁中部地震序列对卡皮亚诺溪流域的地下水系统造成的影响进行分析,结果表明地震序列会暂时改变正常断层的行为,使得含水层孔隙压力变化渗透性增强,地下水流显著地垂直于断层走向分流(Mammoliti et al,2022),大多数离子和δD和δ18O值没有变化,表明是含水层内变化造成对地下水系统的影响(Rosen et al,2018)。综合上述研究中所观察到的水文地质变化,发现地震期间,山地水的释放输送会影响区域的地下水系统,而同位素有可能成为识别同震异常的理想追踪器。

5 地震地下流体监测点映震效能评估

将同位素地球化学方法应用于断裂带和地震监测点附近的水文地球化学背景信息调查和地震流体异常资料完善(Rufino et al,2021),对于评估地下流体的地震前兆信息、流体在孕震过程中的作用以及井孔的映震效能等方面具有重要参考价值。地震会导致含水层破裂,深部热液流体注入,不同层间流体混合,水-岩反应加剧,造成水化学变化(Barberio et al,2017Hosono et al,2018)。因此,与大气起源的地下水相比,地壳深部和深断裂有关的原生水或变质水具有更强的映震灵敏度。对于交替缓慢、封闭性强的地下水,其映震灵敏度强于离补给区近、交替强、开放性的地下水。地下流体观测点的映震能力可能与携带深部构造活动信息的地下水年龄有关。因此,地下水年龄是评估泉点映震能力的重要指标(潘树新等,2001)。地下水的年龄越大,在地下贮留的时间就越长,不能及时地反映出地壳内部地球化学环境的变化;相反,年龄过短的地下水由于水-岩反应不充分,也不能完全反映地壳内部地球化学环境的变化,因此映震能力会降低。而年龄适中的地下水中的某些化学组分可以更完整地显示异常形态,这有利于地震前兆的监测(刘耀炜等,2009)。

地下水补给来源、循环深度、水-岩反应程度对水位、水氡等映震效能的影响是亟待解决的问题之一。研究人员分析了山西省4口井(泉)的水氡观测资料,结合氢氧同位素特征、水-岩平衡状态和地下水循环深度,通过已发生的地震对观测资料进行效能检验,发现温泉井中夏县井和奇村井的映震特征均优于临猗井(穆慧敏等,2019)。分析计算无震时期地震监测点的水质类型、补给混合特征、循环深度与来源等水文地球化学特征,可以为地震监测点流体数据分析以及异常落实与短临跟踪工作奠定基础,并为地震监测点的映震效能判定提供依据。多次6级以上中强地震的震例资料表明,地下流体水位和水化学观测资料的异常与地震活动对应较好,特别是短临异常映震效果更好(汪成民,1990)。了解观测井孔地下水补给来源对于利用观测资料识别前兆异常及准确把握映震效能十分重要,而水化学实验分析则可以为观测井水的补给来源及水循环过程提供科学参考。胡小静等(2018)通过水化学特征结合区域地下水相关性分析和井下电视探视,对云南江川渔村井地下水循环状态与补给来源进行研究,确定了其与周边大型水体星云湖等无明显水力联系(胡小静等,2018)。

在地震监测方面,断层气变化特征可为断裂带异常评价提供背景数据和相关依据(Zhou et al,2021),同位素气体地球化学方法优越性显著。高断层气含量的部位其断裂带断层活动越强,是由于在构造应力的作用下,深部气体沿着发展的活动板块和活动块体边界以及其他活动性断裂等地壳薄弱地带向上运移,导致断层气在断裂带附近富集(吴华平等,2009)。Yuce等(2017)对死海断层(DSF)土壤气体与断层分布之间的关系进行研究,发现深部地壳气体组分表现出高Rn和CO2浓度,且高浓度异常沿该区域内的南北和西北-东南趋势的断层分布,与DSF和阿米克盆地的盆地边界断层以及卡拉苏断层等相符,结果表明,氡浓度和CO2(通量和浓度)是追踪该地区深部断层活动的良好指标。孙小龙等(2016)利用土壤气浓度、水化学组分等多种手段对海原断裂带的构造活动性进行分析,研究结果显示,该区水-岩作用强,水-岩平衡状态高,循环深度大,水泉离补给区远,水循环周期长,新老水体交换慢,同时含幔源气体。综合分析表明,海原断裂带闭锁程度较弱,但与地壳深部的连通性较强,这有利于地下流体的扩散和运移。已有相关的研究表明,地震可能是在岩浆气体优先迁移的断裂带内触发(Aizawa et al,2016)。

6 总结与展望 6.1 总结

地下流体同位素地球化学方法在地震监测预测研究方面有助于深入认识深部流体的活动特征,分析地震的孕育、发生以及地震前兆的成因,从而使得地下流体观测资料在地震科学研究中发挥更加重要的作用。通过同位素地球化学方法确定地下流体来源的原理,可以研究地下流体循环特征,进而评估地质构造活动的程度。而利用同位素地球化学方法能够认识流体异常成因(Hosono et al,2020a),识别地震前兆异常,开展地震预测研究。此外,深部流体同位素示踪技术结合深源流体监测和地球物理方法,有可能揭示地震孕育、流体与震源之间的作用关系(Umeda et al,2015Luetkemeyer et al,2016Sanz et al,2020)。同位素地球化学方法不仅可以构建断裂带流体地球化学背景特征,还可以用于地震监测点映震效能的评估。因此,同位素地球化学是一种具有很高潜力的地震监测预测技术。首先,同位素技术可以在微观尺度上对地震前的地质环境进行深入研究,从而获得地震发生前的信息,提高地震预测的准确性。其次,同位素技术可以尝试对地震孕育过程进行分析,从而对地震机制、地震动力学等方面进行深入研究,提升对断层系统演化和地震周期的理解,为地震预测和防御提供科学依据(Gold et al,2015Rockwell et al,2018)。此外,同位素地球化学还可以用于研究地震震源机制,以便更好地了解地震产生的原因(Okuyama et al,2016)。

近年来,地球化学分析测试技术得到长足发展,其中稀有气体同位素示踪、非传统稳定同位素示踪、氪和碳等同位素定年技术得到不断完善。这些进展使得同位素示踪和定年技术日益成熟,为地震监测预测的研究提供了更为精确的工具和方法。然而,同位素地球化学在地震监测预测方面的应用尚存诸多困难和局限。传统的地下流体同位素监测面临着严峻的挑战。首先,同位素技术的研究尚处初步阶段,在一些复杂的地质环境中缺乏系统科学的解释。其次,同位素技术在应用中受到多方面的限制,如样品数量和质量、采集方式、分析技术等。因此,在未来的研究中,有必要进一步深入探究同位素技术,以期在实际应用中更好地发挥其作用。

6.2 展望

未来同位素地球化学在地震监测预测方面的研究仍将有很大的潜力。首先,需要加强同位素地球化学基础理论的研究,以确定更多同位素比例变化与地震活动之间的关系。其次,随着观测技术的不断提高和数据处理方法的进一步完善,地震前地球内部的化学环境变化将得到更加细致和准确的描述,从而为地震监测预测提供更多技术支撑。

此外,中国地震科学实验场项目将在重点区域广泛布设观测点,与主要断裂带密集观测相结合,并应用多种地球化学观测手段对同位素进行监测分析。这将进一步加大同位素地球化学观测与研究,为认识构造地震的成因机理和预测方法建立提供新的依据与支撑。

另一方面,随着全球地震监测网络的不断完善,同位素地球化学在地震监测预测方面的应用将得到更广泛的推广和应用。为此,《中国地球物理站网(地下流体)规划(2023—2030)》已经编制完成,该规划最显著的特色是实现规模化观测。通过扩大站网规模、增加站点和测项数量,提升地震监测能力。新规划还提出了全国地震中心站配备氢氧同位素仪等仪器设备的措施,以中心站辐射基准站与基本站,完成各区域内同位素的普查和跟踪测试分析研究。这一规划的实施将大幅度提升同位素地球化学方法在地震监测预测研究中的应用。

最后,在未来的研究中,同位素地球化学与其他地球科学技术(如传统地震地下流体学、地球物理学和地质学)的结合和应用也将提高地震监测预测的精确度。通过加大对地震前兆的研究力度,同时结合其他监测预测手段(如卫星遥感),形成一个综合的监测预测体系。同位素地球化学还将与人工智能、机器学习等技术结合,以便更快速、准确地进行地震监测预测,不断提高同位素地球化学在地震监测预测方面的应用效果。同位素地球化学在地震监测预测研究方面仍有很大的发展空间,将在未来的研究中发挥重要作用,为地震监测预测提供更多有力的技术支持。

参考文献
车用太、鱼金子, 2006, 地震地下流体学, 北京: 气象出版社.
陈晨、金章东、徐阳等, 2020, 紫坪铺水库沉积物Sr同位素及元素组成变化对2008年汶川大地震的响应, 湖泊科学, 32(6): 1882-1890.
杜建国、康春丽, 2000, 地震地下流体发展概述, 地震, 20(增刊Ⅰ): 107-114.
杜建国、刘丛强, 2003, 同位素地球化学在地震研究方面的作用, 地震, 23(2): 99-107. DOI:10.3969/j.issn.1000-3274.2003.02.012
方震、刘耀炜、杨选辉等, 2012, 地震断裂带中气体来源及运移机制研究进展, 地球物理学进展, 27(2): 483-495. DOI:10.6038/j.issn.1004-2903.2012.02.011
胡小静、付虹、李利波等, 2018, 云南江川渔村井地下水补给来源分析, 地震研究, 41(4): 544-550. DOI:10.3969/j.issn.1000-0666.2018.04.008
李营、陈志、胡乐等, 2022, 流体地球化学进展及其在地震预测研究中的应用, 科学通报, 67(13): 1404-1420.
刘洁、李营、陈志等, 2022, 陇县—宝鸡断裂带CO2气体地球化学特征及成因, 地震研究, 45(2): 217-228.
刘耀炜、任宏微、王博, 2009, 环境同位素及其示踪技术在地震预测研究中的应用前景, 地学前缘, 16(2): 369-377. DOI:10.3321/j.issn:1005-2321.2009.02.028
穆慧敏、黄春玲、常姣等, 2019, 山西省4口井(泉)水化学特征与水氡映震效能研究, 中国地震, 35(4): 654-665. DOI:10.3969/j.issn.1001-4683.2019.04.006
潘树新、高安泰, 2001, 地下水年龄与泉点映震能力的关系及其在地震预报中应用的初步探讨, 西北地震学报, 23(2): 189-193.
上官志冠, 1995, 深源二氧化碳预报地震研究, 地震地质, 17(3): 214-217.
申宝剑、秦建中、胡文瑄等, 2010, 幔源CO2演化及CO2气成藏实例分析, 中国地质, 37(2): 495-507. DOI:10.3969/j.issn.1000-3657.2010.02.024
孙小龙、刘耀炜、付虹等, 2020, 我国地震地下流体学科分析预报研究进展回顾, 地震研究, 43(2): 216-231. DOI:10.3969/j.issn.1000-0666.2020.02.003
孙小龙、王广才、邵志刚等, 2016, 海原断裂带土壤气与地下水地球化学特征研究, 地学前缘, 23(3): 140-150. DOI:10.13745/j.esf.2016.03.018
唐晗晗、郭良辉、方圆, 2020, 青藏高原东南缘热流估算及与地震活动相关性分析, 地球物理学报, 63(3): 1056-1069.
王博、周永胜、钟骏等, 2022, 滇西北构造地热特征及对地震活动的影响, 地球物理学报, 65(9): 3419-3433.
汪成民, 1990, 中国地震地下水动态观测网, 北京: 地震出版社.
王基华、林元武、高松升, 1998, 怀来断层气CO2监测及张北-尚义地震的短临预报, 地震地质, 20(2): 113-116.
汪万红、张慧、苏鹤军, 2008, 秦岭北缘断裂带温泉水循环深度与地震活动性的关系研究, 西北地震学报, 30(1): 36-41.
王云、李其林、冉华等, 2018, 青藏高原东南缘地热与地震活动: 来自氦同位素的约束, 矿物岩石地球化学通报, 37(4): 652-662.
吴华平、郭良田、常郁等, 2009, 氡断层气测量在佛山西淋岗活断层探测中的应用研究, 华南地震, 29(4): 108-113. DOI:10.3969/j.issn.1001-8662.2009.04.013
杨黎、张立、胡小静等, 2022, 2022年1月2日宁蒗MS5.5地震前迪庆流体观测井异常分析, 地震地磁观测与研究, 43(2): 180-187.
张国盟、刘耀炜、张磊等, 2015, 地下水测年方法及其在地震监测中的应用展望, 中国地震, 31(1): 1-10.
张磊、刘耀炜、任宏微等, 2016, 氢氧稳定同位素在地下水异常核实中的应用, 地震地质, 38(3): 721-731.
张汝惠、武冀新、张晚霞, 1991, 辽宁南部的大地热流及岩石圈热结构, 中国地震, 7(3): 11-24.
周晓成、王万丽、李立武等, 2020, 金沙江—红河断裂带温泉气体地球化学特征, 岩石学报, 36(7): 2197-2214.
Aizawa K, Sumino H, Uyeshima M, et al, 2016, Gas pathways and remotely triggered earthquakes beneath Mount Fuji, Japan, Geology, 44(2): 127-130. DOI:10.1130/G37313.1
Alçiçek H, Bülbül A, Yavuzer Ï, et al, 2019, Origin and evolution of the thermal waters from the Pamukkale Geothermal Field(Denizli Basin, SW Anatolia, Turkey): insights from hydrogeochemistry and geothermometry, J Volcanol Geotherm Res, 372: 48-70. DOI:10.1016/j.jvolgeores.2018.09.011
Apollaro C, Caracausi A, Paternoster M, et al, 2020, Fluid geochemistry in a low-enthalpy geothermal field along a sector of southern Apennines chain(Italy), J Geochem Exp, 219: 106618. DOI:10.1016/j.gexplo.2020.106618
Aydin H, Karakuș H, Mutlu H, 2020, Hydrogeochemistry of geothermal waters in eastern Turkey: geochemical and isotopic constraints on water-rock interaction, J Volcanol Geotherm Res, 390: 106708. DOI:10.1016/j.jvolgeores.2019.106708
Baba A, Sözbilir H, Sayık T, et al, 2022, Hydrogeology and hydrogeochemistry of the geothermal systems and its direct use application: balçova-Narlıdere geothermal system, İzmir, Turkey, Geothermics, 104: 102461. DOI:10.1016/j.geothermics.2022.102461
Barberio M D, Barbieri M, Billi A, et al, 2017, Hydrogeochemical changes before and during the 2016 Amatrice-Norcia seismic sequence(central Italy), Sci Rep, 7(1): 11735. DOI:10.1038/s41598-017-11990-8
Barbieri M, Boschetti T, Barberio M D, et al, 2020, Tracing deep fluid source contribution to groundwater in an active seismic area(central Italy): A combined geothermometric and isotopic(δ13C)perspective, J Hydrol, 582: 124495. DOI:10.1016/j.jhydrol.2019.124495
Bayari S, Ozyurt N N, Hatipoglu Z, et al. 2006. Groundwater age: A vital information in protecting the groundwater dependent ecosystem. In: NATO Advanced Research Workshop on Groundwater and Ecosystems. Dordrecht: Springer, 33~46.
Bernard P, 2001, From the search of 'precursors' to the research on 'crustal transients', Tectonophysics, 338(3~4): 225-232.
Boles J R, Garven G, Camacho H, et al, 2015, Mantle helium along the newport-inglewood fault zone, Los Angeles Basin, California: A leaking paleo-subduction zone, Geochem Geophys Geosyst, 16(7): 2364-2381. DOI:10.1002/2015GC005951
Bonfanti P, Genzano N, Heinicke J, et al, 2012, Evidence of CO2-gas emission variations in the central Apennines(Italy)during the L′Aquila seismic sequence(March-April 2009), Boll Geofis Teor Appl, 53(1): 147-168.
Bräuer K, Geissler W H, Kämpf H, et al, 2016, Helium and carbon isotope signatures of gas exhalations in the westernmost part of the Pannonian Basin(SE Austria/NE Slovenia): evidence for active lithospheric mantle degassing, Chem Geol, 422: 60-70. DOI:10.1016/j.chemgeo.2015.12.016
Bräuer K, Kämpf H, Strauch G, et al, 2003, Isotopic evidence(3He/4He, 13CCO2)of fluid-triggered intraplate seismicity, J Geophys Res Solid Earth, 108(B2): 2070.
Burnard P, Zimmermann L, Sano Y. 2013. The noble gases as geochemical tracers: history and background. In: Burnard P. The Noble Gases as Geochemical Tracers. Berlin: Springer, 1~15.
Buttitta D, Caracausi A, Chiaraluce L, et al, 2020, Continental degassing of helium in an active tectonic setting(northern Italy): the role of seismicity, Sci Rep, 10(1): 162. DOI:10.1038/s41598-019-55678-7
Caracausi A, Buttitta D, Picozzi M, et al, 2022, Earthquakes control the impulsive nature of crustal helium degassing to the atmosphere, Commun Earth Environ, 3(1): 224. DOI:10.1038/s43247-022-00549-9
Chaudhuri H, Barman C, Iyengar A N S, et al, 2013, Network of seismo-geochemical monitoring observatories for earthquake prediction research in India, Acta Geophys, 61(4): 1000-1025.
Chen Z, Du J G, Zhou X C, et al, 2014, Hydrochemistry of the hot springs in western Sichuan Province related to the Wenchuan MS8.0 earthquake, Sci World J, 2014: 901432.
Chen Z, Zhou X, Du J, et al, 2015, Hydrochemical characteristics of hot spring waters in the Kangding district related to the Lushan MS=7.0 earthquake in Sichuan, China, Nat Hazards Earth Syst Sci, 15(6): 1149-1156. DOI:10.5194/nhess-15-1149-2015
Chiodini G, Caliro S, Cardellini C, et al, 2011, Geochemical evidence for and characterization of CO2 rich gas sources in the epicentral area of the Abruzzo 2009 earthquakes, Earth Planet Sci Lett, 304(3~4): 389-398.
Chiodini G, Cardellini C, Di Luccio F, et al, 2020, Correlation between tectonic CO2 Earth degassing and seismicity is revealed by a 10-year record in the Apennines, Italy, Sci Adv, 6(35): eabc2938. DOI:10.1126/sciadv.abc2938
Claesson L, Skelton A, Graham C, et al, 2004, Hydrogeochemical changes before and after a major earthquake, Geology, 32(8): 641-644. DOI:10.1130/G20542.1
Claesson L, Skelton A, Graham C, et al, 2007, The timescale and mechanisms of fault sealing and water-rock interaction after an earthquake, Geofluids, 7(4): 427-440. DOI:10.1111/j.1468-8123.2007.00197.x
Clark I D, Fritz P, 2013, Environmental Isotopes in Hydrogeology, Boca Raton: CRC Press.
Cornides I, Takaoka N, Nagao K, et al, 1986, Contribution of mantle-derived gases to subsurface gases in a tectonically quiescent area, the Carpathian Basin, Hungary revealed by noble gas measurements, Geochem J, 20(3): 119-125. DOI:10.2343/geochemj.20.119
Craig H, 1961, Isotopic variations in meteoric waters, Science, 133(3465): 1702-1703. DOI:10.1126/science.133.3465.1702
Du J G, Cheng W Z, Zhang Y L, et al, 2006, Helium and carbon isotopic compositions of thermal springs in the earthquake zone of Sichuan, Southwestern China, J Asian Earth Sci, 26(5): 533-539. DOI:10.1016/j.jseaes.2004.11.006
Franchini S, Agostini S, Barberio M D, et al, 2021, HydroQuakes, central apennines, Italy: towards a hydrogeochemical monitoring network for seismic precursors and the hydro-seismo-sensitivity of boron, J Hydrol, 598: 125754. DOI:10.1016/j.jhydrol.2020.125754
Gehrke G E, Blum J D, Marvin-DiPasquale M, 2011, Sources of mercury to San Francisco Bay surface sediment as revealed by mercury stable isotopes, Geochim Cosmochim Acta, 75(3): 691-705. DOI:10.1016/j.gca.2010.11.012
Gold P O, Behr W M, Rood D, et al, 2015, Holocene geologic slip rate for the Banning strand of the southern San Andreas Fault, southern California, J Geophys Res Solid Earth, 120(8): 5639-5663. DOI:10.1002/2015JB012004
Hosono T, Hartmann J, Louvat P, et al, 2018, Earthquake-induced structural deformations enhance long-term solute fluxes from active volcanic systems, Sci Rep, 8(1): 14809. DOI:10.1038/s41598-018-32735-1
Hosono T, Saltalippi C, Jean J S, 2020a, Coseismic hydro-environmental changes: insights from recent earthquakes, J Hydrol, 585: 124799. DOI:10.1016/j.jhydrol.2020.124799
Hosono T, Yamada C, Manga M, et al, 2020b, Stable isotopes show that earthquakes enhance permeability and release water from mountains, Nat Commun, 11(1): 2776. DOI:10.1038/s41467-020-16604-y
Ide K, Hosono T, Kagabu M, et al, 2020, Changes of groundwater flow systems after the 2016 MW7.0 Kumamoto earthquake deduced by stable isotopic and CFC-12 compositions of natural springs, J Hydrol, 583: 124551. DOI:10.1016/j.jhydrol.2020.124551
Inbal A, Ampuero J P, Clayton R W, 2016, Localized seismic deformation in the upper mantle revealed by dense seismic arrays, Science, 354(6308): 88-92. DOI:10.1126/science.aaf1370
Ingebritsen S E, Manga M, 2014, Hydrogeochemical precursors, Nat Geosci, 7(10): 697-698. DOI:10.1038/ngeo2261
Karaoǧlu Ö, Bazargan M, Baba A, et al, 2019, Thermal fluid circulation around the Karliova triple junction: Geochemical features and volcano-tectonic implications(Eastern Turkey), Geothermics, 81: 168-184. DOI:10.1016/j.geothermics.2019.05.003
Kennedy B M, Kharaka Y K, Evans W C, et al, 1997, Mantle fluids in the San Andreas Fault system, California, Science, 278(5341): 1278-1281. DOI:10.1126/science.278.5341.1278
Klemperer S L, Kennedy B M, Sastry S R, et al, 2013, Mantle fluids in the Karakoram fault: Helium isotope evidence, Earth Planet Sci Lett, 366: 59-70. DOI:10.1016/j.epsl.2013.01.013
Kulongoski J T, Hilton D R, Barry P H, et al, 2013, Volatile fluxes through the Big Bend section of the San Andreas Fault, California: helium and carbon-dioxide systematics, Chem Geol, 339: 92-102. DOI:10.1016/j.chemgeo.2012.09.007
Lee H, Kim H, Kagoshima T, et al, 2019, Mantle degassing along strike-slip faults in the Southeastern Korean Peninsula, Sci Rep, 9(1): 15334. DOI:10.1038/s41598-019-51719-3
Lee H, Muirhead J D, Fischer T P, et al, 2016, Massive and prolonged deep carbon emissions associated with continental rifting, Nat Geosci, 9(2): 145-149. DOI:10.1038/ngeo2622
Liu J, Chen Z Y, Wei W, et al, 2014, Using chlorofluorocarbons(CFCs)and tritium(3H)to estimate groundwater age and flow velocity in Hohhot Basin, China, Hydrol Process, 28(3): 1372-1382. DOI:10.1002/hyp.9659
Lombardi S, Voltattorni N, 2010, Rn, He and CO2 soil gas geochemistry for the study of active and inactive faults, Appl Geochem, 25(8): 1206-1220. DOI:10.1016/j.apgeochem.2010.05.006
Luetkemeyer P B, Kirschner D L, Huntington K W, et al, 2016, Constraints on paleofluid sources using the clumped-isotope thermometry of carbonate veins from the SAFOD(San Andreas Fault Observatory at Depth)borehole, Tectonophysics, 690: 174-189. DOI:10.1016/j.tecto.2016.05.024
Mammoliti E, Fronzi D, Cambi C, et al, 2022, A holistic approach to study groundwater-surface water modifications induced by strong earthquakes: the case of campiano catchment(central Italy), Hydrology, 9(6): 97. DOI:10.3390/hydrology9060097
Marty B, Zimmermann L, 1999, Volatiles(He, C, N, Ar)in mid-ocean ridge basalts: Assesment of shallow-level fractionation and characterization of source composition, Geochim Cosmochim Acta, 63(21): 3619-3633. DOI:10.1016/S0016-7037(99)00169-6
Mchugh C M, Seeber L, Rasbury T, et al, 2020, Isotopic and sedimentary signature of megathrust ruptures along the Japan subduction margin, Mar Geol, 428: 106283.
Morales-Simfors N, Wyss R A, Bundschuh J, 2020, Recent progress in radon-based monitoring as seismic and volcanic precursor: a critical review, Crit Rev Environ Sci Technol, 50(10): 979-1012.
Muto J, Yasuoka Y, Miura N, et al, 2021, Preseismic atmospheric radon anomaly associated with 2018 Northern Osaka earthquake, Sci Rep, 11(1): 7451.
Oh Y H, Kim G, 2015, A radon-thoron isotope pair as a reliable earthquake precursor, Sci Rep, 5: 13084.
Okuyama Y, Funatsu T, Fujii T, et al, 2016, Mid-crustal fluid related to the Matsushiro earthquake swarm(1965-1967)in northern Central Japan: geochemical reproduction, Tectonophysics, 679: 61-72.
Onda S, Sano Y, Takahata N, et al, 2018, Groundwater oxygen isotope anomaly before the M6.6 Tottori earthquake in Southwest Japan, Sci Rep, 8(1): 4800.
O'Neil J R, King C Y, 1981, Variations in stable-isotope ratios of ground waters in seismically active regions of california, Geophys Res Lett, 8(5): 429-432.
Palcsu L, Vetö I, Futó I, et al, 2014, In-reservoir mixing of mantle-derived CO2 and metasedimentary CH4-N2 fluids-Noble gas and stable isotope study of two multistacked fields(Pannonian Basin System, W-Hungary), Mar Pet Geol, 54: 216-227.
Pang Z H, Kong Y L, Li J, et al, 2017, An isotopic geoindicator in the hydrological cycle, Procedia Earth Planet Sci, 17: 534-537.
Pasvanoǧlu S, 2020, Geochemistry and conceptual model of thermal waters from Erciș-Zilan Valley, Eastern Turkey, Geothermics, 86: 101803.
Pierotti L, Botti F, D′intinosante V, et al, 2015, Anomalous CO2 content in the Gallicano thermo-mineral spring(Serchio Valley, Italy)before the 21 June 2013, Alpi Apuane earthquake(M=5.2), Phys Chem Earth, 85~86: 131-140.
Polyak B G, Tolstikhin I I, Yakutseni V P, 1979, The isotopic composition of helium and heat-flow-geochemical and geophysical aspects of tectogenesis, Geotectonics, 13(5): 339-351.
Quattrocchi F, Pizzi A, Gori S, et al, 2012, The contribution of fluid geochemistry to define the structural patterns of the 2009 L'Aquila seismic source, Italian J Geosci, 131(3): 448-458.
Ramola R C, Prasad Y, Prasad G, et al, 2008, Soil-gas radon as seismotectonic indicator in Garhwal Himalaya, Appl Radiat Isot, 66(10): 1523-1530.
Reddy D V, Kumar D, Rao N P, 2017, Long-term hydrochemical earthquake precursor studies at the Koyna-Warna reservoir site in western India, J Geol Soc India, 90(6): 720-727.
Rockwell T K, Meltzner A J, Haaker E C, 2018, Dates of the two most recent surface ruptures on the southernmost San Andreas Fault recalculated by precise dating of lake cahuilla dry periods, Bull Seismol Soc Am, 108(5A): 2634-2649.
Rosen M R, Binda G, Archer C, et al, 2018, Mechanisms of earthquake-induced chemical and fluid transport to carbonate groundwater springs after earthquakes, Water Resour Res, 54(8): 5225-5244.
Rufino F, Cuoco E, Busico G, et al, 2021, Deep carbon degassing in the Matese massif chain(Southern Italy)inferred by geochemical and isotopic data, Environ Sci Pollut Res, 28(34): 46614-46626.
Sano Y, Hara T, Takahata N, et al. 2014. Helium anomalies suggest a fluid pathway from mantle to trench during the 2011 Tohoku-Oki earthquake, Nature Communications, 5.
Sano Y, Kagoshima T, Takahata N, et al, 2020a, Groundwater anomaly related to CCS-CO2 injection and the 2018 Hokkaido eastern Iburi earthquake in Japan, Front Earth Sci, 8: 611010.
Sano Y, Marty B, 1995, Origin of carbon in fumarolic gas from island arcs, Chem Geol, 119(1~4): 265-274.
Sano Y, Onda S, Kagoshima T, et al, 2020b, Groundwater oxygen anomaly related to the 2016 Kumamoto earthquake in Southwest Japan, Proc Jpnan Acad Ser B Phys Biol Sci, 96(7): 322-334.
Sano Y, Takahata N, Igarashi G, et al, 1998, Helium degassing related to the Kobe earthquake, Chem Geol, 150(1~2): 171-179.
Sano Y, Takahata N, Kagoshima T, et al, 2016, Groundwater helium anomaly reflects strain change during the 2016 Kumamoto earthquake in Southwest Japan, Sci Rep, 6: 37939.
Sanz E, Menéndez Pidal I, Escavy J I, et al, 2020, Hydrogeological changes along a fault zone caused by earthquakes in the Moncayo massif(Iberian Chain, Spain), Sustainability, 12(21): 9034.
Sato T, Takahashi H A, Kawabata K, et al, 2020, Changes in the nitrate concentration of spring water after the 2016 Kumamoto earthquake, J Hydrol, 580: 124310.
Sciarra A, Cantucci B, Coltorti M, 2017, Learning from soil gas change and isotopic signatures during 2012 Emilia seismic sequence, Sci Rep, 7(1): 14187.
Sherman L S, Blum J D, Johnson K P, et al, 2010, Mass-independent fractionation of mercury isotopes in Arctic snow driven by sunlight, Nat Geosci, 3(3): 173-177.
Skelton A, Andrén M, Kristmannsdóttir H, et al, 2014, Changes in groundwater chemistry before two consecutive earthquakes in Iceland, Nat Geosci, 7(10): 752-756.
Skelton A, Liljedahl-Claesson L, Wästeby N, et al, 2019, Hydrochemical changes before and after earthquakes based on long-term measurements of multiple parameters at two sites in northern iceland-a review, J Geophys Res Solid Earth, 124(3): 2702-2720.
Soldati G, Cannelli V, Piersanti A, 2020, Monitoring soil radon during the 2016-2017 central Italy sequence in light of seismicity, Sci Rep, 10(1): 13137.
Streig A R, Weldon II R J, Biasi G, et al, 2020, New insights into paleoseismic age models on the northern san andreas fault: charcoal inbuilt ages and updated earthquake correlations, Bull Seismol Soc Am, 110(3): 1077-1089.
Sukanya S, Noble J, Joseph S, 2022, Application of radon(222Rn)as an environmental tracer in hydrogeological and geological investigations: an overview, Chemosphere, 303: 135141.
Tardani D, Reich M, Roulleau E, et al, 2016, Exploring the structural controls on helium, nitrogen and carbon isotope signatures in hydrothermal fluids along an intra-arc fault system, Geochim Cosmochim Acta, 184: 193-211.
Umeda K, Asamori K, Kusano T, 2013, Release of mantle and crustal helium from a fault following an inland earthquake, Appl Geochem, 37: 134-141.
Umeda K, Asamori K, Makuuchi A, et al, 2015, Triggering of earthquake swarms following the 2011 Tohoku megathrust earthquake, J Geophys Res Solid Earth, 120(4): 2279-2291.
Uzelli T, Baba A, Mungan G G, et al, 2017, Conceptual model of the Gülbahçe geothermal system, Western Anatolia, Turkey: Based on structural and hydrogeochemical data, Geothermics, 68: 67-85.
Wakita H, Nakamura Y, Kita I, et al, 1980, Hydrogen release: new indicator of fault activity, Science, 210(4466): 188-190.
Wang C Y, Manga M, 2015, New streams and springs after the 2014 MW6.0 South Napa earthquake, Nat Commun, 6: 7597.
Wang G C, Zhang Z C, Wang M, et al, 2005, Implications of ground water chemistry and flow patterns for earthquake studies, Ground Water, 43(4): 478-484.
Wästeby N, Skelton A, Tollefsen E, et al, 2014, Hydrochemical monitoring, petrological observation, and geochemical modeling of fault healing after an earthquake, J. Geophys. Res. Solid Earth, 119(7): 5727-5740.
Xu S, Zheng G D, Nakai S I, et al, 2013, Hydrothermal He and CO2 at Wudalianchi intra-plate volcano, NE China, J Asian Earth Sci, 62: 526-530.
Yasuoka Y, Shinogi M, 1997, Anomaly in atmospheric radon concentration: a possible precursor of the 1995 Kobe, Japan, earthquake, Health Phys, 72(5): 759-761.
Yuce G, Fu C C, D′alessandro W, et al, 2017, Geochemical characteristics of soil radon and carbon dioxide within the Dead Sea Fault and Karasu Fault in the Amik Basin(Hatay), Turkey, Chem Geol, 469: 129-146.
Zhang L, Guo L S, Zhou X C, et al, 2021, Temporal variations in stable isotopes and synchronous earthquake-related changes in hot springs, J Hydrol, 599: 126316.
Zhang L, Yang D X, Liu Y W, et al, 2014, Impact of impoundment on groundwater seepage in the Three Gorges Dam in China based on CFCs and stable isotopes, Environ Earth Sci, 72(11): 4491-4500.
Zhao W B, Guo Z F, Liu J Q, et al, 2021, Fluxes and genesis of carbon dioxide emissions from Cenozoic volcanic fields in NE China, Acta Petrol Sin, 37(4): 1255-1269.
Zhou X C, Liu L, Chen Z, et al, 2017, Gas geochemistry of the hot spring in the Litang fault zone, Southeast Tibetan Plateau, Appl Geochem, 79: 17-26.
Zhou X C, Yan Y C, Fang W Y, et al, 2021, Short-term seismic precursor anomalies of hydrogen concentration in Luojishan hot spring bubbling gas, eastern Tibetan Plateau, Front Earth Sci, 8: 586279.