Variations of apparent resistivity before the 2022 MS6.9 Menyuan earthquake in Qinghai Province revealed by the virtual fault dislocation model
-
摘要:
选取2022年1月8日青海门源MS6.9地震震中400 km范围内四个地电阻率观测台站的观测精度高、具有稳态年变、震前无显著干扰的地电阻率$ {\rho }_{{\mathrm{s}}} $观测数据，结合ERA5同化数据集中的多层土壤温度和土壤水分含量，在利用多项式拟合获取各台站（或测道）地电阻率正常年动态的基础上，分析了门源地震前地电阻率的异常变化。结果显示：金银滩台EW测道、武威台NS测道和山丹台EW测道、N45°W测道震前存在超阈值的异常变化，并呈现各向异性特征。基于断层虚位错模式分析了地电阻率异常变化与孕震过程之间的联系，结果表明：金银滩台震前处于压缩区并受到NNE方向的挤压，与主压应变近似正交的EW测道于震前10个月出现负异常；同样位于压缩区的武威台，受到了ENE向的挤压，NS测道的地电阻率在孕震早期（震前13个月）以负异常为主，孕震中晚期（震前3个月）出现了正异常；山丹台，位于膨胀区，受到近似NS向的拉张，与主张应变平行的NS测道未发现异常，但EW测道震前一年地电阻率出现正异常，N45°W测道的地电阻率也在震前半年左右出现超阈值并呈正异常。此外，金银滩台、山丹台和武威台距离门源地震震中的距离分别为92 km，113 km和139 km，相应的地电阻率异常最大变化幅值分别为−3.0σ，2.2σ和−2.1σ。此外，门源地震前地电阻率异常变化的时空特征与岩石实验结果及理论模型一致，也符合震源区应力应变积累程度较高、向外围方向逐渐衰减的分布特征。由此推断，2022年门源MS6.9地震前地电阻率的时空变化可能与区域介质变形及应力变化有关。
Abstract:An earthquake with MS6.9 struck the Menyuan County, Qinghai Province, northwestern China on 8 January 2022. There were four apparent resistivity observation stations with anomalies in an area within 400 km from the epicenter. Using the observed data with high quality and stable annual variation as well as the soil temperature and soil water provided by ERA5 assimilation datasets, a polynomial fitting was performed so as to subtract the normal annual dynamics and to detect the anomalies of the apparent resistivity prior to the earthquake. The results indicate that the above-threshold changes and anisotropic characteristics appeared before the earthquake in the EW channel at the Jinyintan station, NS channel at the Wuwei station, and the EW and N45°W channels at the Shandan station. A virtual fault dislocation model was used to examine the relationship between apparent resistivity changes and seismogenic process of the Menyuan earthquake. In the areas with compression along NNE direction, the decrease changes in the EW channel, which is approximately orthogonal to the direction of the principal compressive strain, has been observed at the Jinyintan station since ten months before the earthquake, the negative anomaly in the NS channel in the early stages of the earthquake preparation （13 months before the main shock） appeared at the Wuwei station, while it turned to an increase change three months before the earthquake. At the Shandan station, which is located in the relative extensional area, no anomalies were detected in the NS channel parallel to the principle tensile strain. However, an increase change was observed in the EW channel one year prior to the earthquake, and an increase change was also recorded in the N45°W channel of Shandan station half-year before the earthquake. Furthermore, at the Jinyintan station, Shandan station, and Wuwei station with epicentral distance of 92 km, 113 km, and 139 km, the corresponding maximum variation of anomalies are −3.0σ, 2.2σ and −2.1σ respectively. In addition, Variations of apparent resistivity before the 2022 MS6.9 Menyuan earthquake were consistent with the results from rock experiment and theoretical models. Moreover, the spatio-temporal characteristics of the variation of apparent resistivity before the MS6.9 Menyuan earthquake are likely consistent with the stress accumulation in the source region, as well as with the characteristics of stress accumulation by a high degree in the epicenter and gradual attenuation towards the periphery. Therefore, it can be deduced that the spatio-temporal variation of apparent resistivity before the MS6.9 Menyuan earthquake may be related to regional medium deformation and stress change.
-
引言
海啸是联合国“国际减灾十年”所确定的全球重大自然灾害之一，它严重地威胁着沿海地区的安全。2011年3月11日5时46分24秒（UTC），在日本东北部海域发生MW9.0地震，该地震导致大约长450 km的日本海沟发生破裂，同时也引发了一场毁灭性的海啸，给太平洋沿岸地区造成了巨大的灾难（Fujii et al，2011）。
海啸是一种海洋重力波，其在传播过程中由于受到水深变化、海底地形以及近岸地形的反射、衍射等的影响，其传播路径随时间变化，表现为海啸能量的各向异性扩散（Titov et al，2005；Fine et al，2013；Shevchenko et al，2014）。由于海啸现象较为复杂，对海啸的形成机制、传播特征以及海啸源的定位等仍然存在诸多争议（Titov et al，2005；Fujii，Satake，2007；Fujii et al，2011）。早期研究对海平面的测量主要是沿岸记录（Pugh，Woodworth，2014），但是后续研究显示，在沿海地区，非线性、摩擦、混叠、仪器噪声、温度变化和风浪的影响会使海啸信号本身失真，这使得海啸监测更加复杂化，面临的挑战也更艰巨（Mungov et al，2013；Pugh，Woodworth，2014）。
海啸与地震有着密不可分的关系。据不完全统计，全球90%左右的海啸事件是由海底地震引发（叶琳等，2005）。通过海啸波形分析，能够推算引发海啸的地震震源参数（Fujii，Satake，2007；Satake et al，2013）。此外，海洋也是地球背景自由振荡、微地动和高频P波噪声的来源（Longuet-Higgins，Ursell，1948；Nishida et al，2000；Ekström，2001；Watada，Masters，2001；Rhie，Romanowicz，2004；Bromirski et al，2005；Xiao et al，2018），这表明海啸也会引起地震背景噪声的变化，通过揭示这种关联有望从新的角度来加强对海啸的监测。
地震背景噪声按周期从长到短可以分为地球背景自由振荡、微地动和更高频的地震噪声，目前对地震背景噪声的研究集中在反映大气-海洋-地球之间大尺度耦合的地球背景自由振荡（2—7 mHz）和微地动（0.05—0.5 Hz）上，并建立了较为系统的理论来阐述这些噪声的激发机制（Rhie，Romanowicz，2004；Rhie，Romanowicz，2006）。实际上，地球背景自由振荡和微地动都是地震台站所记录到的地面运动，由于这些运动是由天然地震以外的“震源”引起的，在以研究天然地震为主的记录上表现为需要被去除的背景噪声，因此通常将其统称为地震背景噪声.
地震观测研究显示，在没有地震的影响下仍然存在长周期的地面振动，即地球背景自由振荡（2—7 mHz） （Benioff et al，1959；Gilbert，MacDonald，1959；Kobayashi，Nishida，1998；Suda et al，1998），其来源最初被认为是大气的干扰（Kobayashi，Nishida，1998；Nishida，Kobayashi，1999；Tanimoto，Um，1999；Fukao et al，2002），后来关于地球背景自由振荡的季节性变化以及海底压力变化谱与地震噪声谱之间的相似性等方面的研究结果表明，海洋是地球背景自由振荡的来源（Nishida et al，2000；Ekström，2001；Watada，Masters，2001；Rhie，Romanowicz，2004；Webb，2007），Tanimoto （2005）基于此建立了海洋作为地球背景自由振荡噪声源的模型。
而这里提到的微地动并非是微小地震，而是海洋活动通过海陆耦合引发的连续地震波形记录上频率处于0.05—0.5 Hz的噪声（Longuet-Higgins，1950；Hasselmann，1963；Cicconi et al，1983；Bromirski et al，2005）。微地动一般分为两个频带：0.05—0.08 Hz的单频微地动和0.1—0.5 Hz的双频微地动，而且双频微地动的峰值有时会分裂成两个，其中峰值出现在0.1—0.15 Hz频带的被称为长周期双频微地动，峰值大约出现在0.18—0.4 Hz频带的被称为短周期双频微地动（Tanimoto，2007；Bromirski，Gerstoft，2009；Zhang et al，2010a；Sun et al，2013；Koper，Burlacu，2015）。一般认为单频微地动通过海浪的破裂及其对海岸的直接影响而在浅水中激发；双频微地动则由两列传播方向相反且频率相同或相近的波列，通过非线性波-波相互作用所激发，其频率是单频微地动的两倍，且其能量一般比单频微地动大得多（Longuet-Higgins，1950）。长周期双频微地动可能的源区有海岸（Bromirski et al，2005；Stehly et al，2006；Tanimoto et al，2006；Bromirski，2009）或开放大洋（Obrebski et al，2012），短周期双频微地动可能是由台站附近的重力波非线性作用所激发（Bromirski et al，2005），但也有可能是不同源区激发的特定频率微地动的叠加造成了微地动的分裂（Koper，Burlacu，2015），目前微地动分裂的机理仍不清楚。
除了地球背景自由振荡和微地动，连续地震波形记录上有时候会记录到更高频段（0.4—6 Hz）的高频噪声，其在地震背景噪声谱上的部分被称为火鲁（Holu）谱，这可能是由局部海洋风激发的海浪与海底作用而产生（McCreery et al，1993；Bromirski et al，2005；Sun et al，2013）。目前已有部分研究识别出火鲁谱内0.6—2.0 Hz频段的P波噪声（Gerstoft et al，2008；Zhang et al，2009，2010a，b），但高频P波地震噪声的产生机理尚不确定，Zhang等（2009，2010a）和Sun等（2013）认为这种噪声是由本地近海风通过风生海浪所激发。
目前地震背景噪声的源定位于海洋，其形成机制归因于海浪之间以及海浪与固体地球之间的相互作用，而海啸事件作为海浪运动的一种极端形式，是对整个海洋产生影响的大事件，其对连续地震波形记录上的背景噪声具有不可忽视的影响。但是，当前对海啸事件在连续地震波形记录的特征却知之甚少。关于海啸的地震学特征及其对地球背景自由振荡、微地动以及高频噪声的影响都是亟需研究的重要内容。
鉴于此，本研究以2011年3月11日日本东海岸MW9.0强震引发的海啸为例，详细阐述该海啸事件对上述三类地震背景噪声的影响，试图从连续地震波形记录中拾取海啸，以期为海啸预警、海啸物理学、背景噪声的应用以及海洋-固体地球系统耦合机制研究提供参考。
1. 数据和方法
1.1 数据
本研究从美国地震学研究联合会（Incorporated Research Institutions for Seismology，缩写为IRIS）网站https://ds.iris.edu上下载了宽频地震台站TA.O02D，TA.M31A和TA.TPFO从2011年3月10日至2011年3月15日的连续地震波形记录，台站位置分布如图1所示。台站选择基于以下考虑：选择距离震源较远的台站，利用海啸波与地震波之间的传播速度差，通过较远距离的传播，得到与主震激发的地震波分离较好的海啸波所激发的地震信号。此外，本研究还收集了地震台站附近深海评估和海啸预报系统（Deep-Ocean Assessment and Reporting of Tsunamis，简写为DART） DART46411的海底压力数据（图1），基于此得到海啸首波到达台站较精确的时刻。研究中还需要考虑潮汐的影响，潮汐预测数据源于美国海洋及大气管理局（National Oceanic and Atmospheric Administration，缩写为NOAA）网站https://www.noaa.gov。
![]()
图 1 本文所用的地震台站、DART站、潮汐计以及2011年3月10日至3月15日期间全球M≥5.0地震震中位置的空间分布图（圆圈大小表示地震震级的相对大小）Figure 1. The spatial distributions of seismic stations，DART station and tide meters used in this study as well as the global M≥5.0 earthquake epicenter locations from March 10 to 15，2011 （The circle size indicates the relative magnitude of the earthquake）1.2 方法
1.2.1 时频分析
为了研究海啸在连续地震波形记录上的特征，本文应用时频分析方法来分析连续地震波形记录。对于每个连续地震波形记录，首先进行预处理，包括去均值、去线性趋势、去仪器响应，并进行0.004—4 Hz的带通滤波（考虑到本研究中宽频带地震仪的最长有效周期为240 s）（图2a）；然后将连续地震波形记录分解为连续半小时的时间序列，对其进行快速傅里叶变换，以获得它们各自的频谱；最后按时间顺序排列这些频谱，得到连续地震波形记录的时频图，如图2b所示。
![]()
图 2 2011年3月10日至15日TA.O02D台站垂直分量的时频分析结果（a） 波形图；（b） 时频图；（c） 不同频段的半小时平均振幅A曲线；（d） 微地动频带内不同频段的半小时平均振幅A曲线。图（c）和（d）中红色圆圈表示M≥5地震的发生时刻和震级大小。地震目录引自美国国家地震信息中心（USGS，2020）Figure 2. Time-frequency analysis results of the vertical component of the station TA.O02D from March 10 to 15，2011（a） Waveform diagram；（b） Time-frequency diagram；（c） Half-hour average amplitude curves in different frequency bands；（d） Half-hour average amplitude curves of different frequency bands in microseisms. In Figs. （c） and （d），red circles indicate the time and magnitude of M≥5 earthquakes. Earthquake catalogue comes from USGS （2020）由于此次日本MW9.0地震的震级较大，所产生地震波的频率范围较宽，低频成分不少，可能存在与海啸事件产生的信号重叠的频带，于是本文以多种处理方式尝试压制主震强振幅的影响。例如，采用绝对均值归一化方法对连续地震波形记录进行归一化处理（Bensen et al，2007），以减弱地震发生时刻较强的地震能量，但出现了过校正的现象，并未取得预想的效果；又如，从实际观测的连续地震波形记录中减去正演所得的此次地震的理论连续地震波形记录，以消除大地震的影响，但是由于IRIS提供的理论合成连续地震波形记录与实际观测的连续地震波形记录的匹配度不够，也未能取得预想的效果。地震台站能够记录到由主震直接激发的地震波，以及由主震激发的海啸波传播到地震台站附近海域后与海岸作用所激发的地震波。由于本文所选取的地震台站和主震分别位于太平洋两岸，考虑到地震波与海啸波传播速度的差异，我们能够使两种不同激发机制所产生的地震信号在连续地震波形记录的时域上拉开距离。地震激发的地震波记录通常小于一个小时，在连续地震波形记录的时频分析图中仅占1到2个数据点，表现为尖峰（Sun et al，2013）。但此次引起海啸的日本大地震的矩震级高达MW9.0，震源持续时间超过三分钟（Ye et al，2013），直至震后近一天还能观测到该地震激发的绕地球表面传播的全球面波，不过后续的全球面波的到时有规律可循，且幅度显著降低，因此，本文选取距离震中较远的地震台站，以期在时域上能凸显海啸的特征，降低主震所激发的地震信号的影响。为便于比较，我们还在时频图上绘制了M≥5.0地震的发生时间，如图2c和2d所示。
1.2.2 极化分析
为了研究海啸激发地震波的源区，本文采取极化分析方法对地震台站接收到的三分量记录在不同频段进行极化分析，以此来判断该频段的能量从哪个方位到达地震台站，从而对不同频带的源区方位有个大致了解。与此同时，使用主极化方向方法来研究地震台站接收到的不同频带能量的到达方位（Rhie，Romanowicz，2006；Zhang et al，2010a）。在台站数目较少时，该方法比较有效（Sun et al，2013）。
应用主极化方向方法获取能量到达方位的具体过程如下：① 针对每个台站的水平分量记录中正北和正东分量的记录，从0°到360°以1°为间隔进行旋转；② 对所有旋转后的连续地震波形记录进行快速傅里叶变换得到其频谱；③ 计算不同旋转角度下目标频带内连续地震波形记录的平均振幅；④ 在目标频带内比较不同方位角的平均振幅大小，确定目标频带地震背景噪声的主要极化方向。与时频方法类似，我们将连续地震波形记录分解成连续半小时的时间序列，然后对每半小时时间序列应用主极化方向方法，获得目标频带随时间变化的主极化方向。
2. 结果
2.1 海啸对不同频段地球背景噪声的影响
我们首先对TA.O02D台站从2011年3月10日至15日海啸发生前后5天BHZ分量的波形连续记录进行时频分析，以判断海啸对地球背景自由振荡、微地动和高频噪声的影响（图2a和2b），并提取了地震频谱在地球背景自由振荡、微地动和高频噪声三个频带的平均振幅随时间变化的曲线（图2c和2d）。
图2b中最为显著的是0.1—0.25 Hz频段，该频段代表微地动水平状高能量条带，其次是与MW9.0地震相关的竖直状高能量条带。大地震发生时，由于地震的能量极大且激发的地震波频率成分丰富，在0—1.5 Hz频带内的地震波能量都呈突增状态。另外，在时频图上还可观测到不同时间出现的包含较多频率成分的数个竖直状高能量条带，尤其是在地震发生后的几个小时内，这些条带较为密集且能量水平较高，主要对应于该地震的余震。随着时间的推移，代表微地动的水平能量条带的高频成分能量逐渐增加，这可能反映了海啸作用于海底或海岸所导致的海洋整体活动性的增加。
由图2c可以看出，地震发生后，地球背景自由振荡、微地动和高频噪声的振幅均有一定程度的陡增，且随着时间的推移，振幅逐渐降低，其中：地球背景自由振荡（0.004—0.007 Hz）的衰减速率最小，震后近四天内逐渐衰减至震前的背景水平；微地动频段的能量衰减得相对快一点，大约在震后24小时振幅衰减至背景水平，而在震后48小时左右微地动频段噪声水平又逐渐上升；在1.3—1.5 Hz频带内的高频噪声，随着时间的推移振幅快速地降至背景水平。以上观测表明，海啸对地球背景自由振荡频带的影响最为显著，持续时间也最长，而对高频噪声频带的影响最小。这也表明海啸波激发的地震背景噪声的频率较低，从而对频率较低的地球背景自由振荡的影响较大且作用时间较长。而海啸波对微地动的影响主要体现在海啸在整个海洋传播开来并开始与海底和海岸作用，因此海啸对微地动的影响会延后发生，在地震造成微地动能量陡增并衰减至背景水平后，由海啸波作用于海底和海岸所引发的微地动能量开始逐渐增加并持续。
为了研究海啸对不同类型微地动的影响，我们给出了海啸期间不同频带微地动的平均振幅谱随时间的详细变化，如图2d所示，可见：在地震发生时刻，微地动在各频段均出现能量的激增，但是能量主要集中在0.15 Hz以下的长周期双频微地动和单频微地动频段，尤其是单频微地动频段的能量更高；随着震后时间的推移，各频段的能量均呈衰减状态，但衰减的速率有所差异。短周期双频微地动在地震发生时刻振幅突增，随后又快速衰减至背景水平并保持，未体现出海啸波的传播和到达，表明该频带主要受地震而非海啸的影响。相比之下，长周期双频微地动和单频微地动的能量增幅更大，衰减至背景水平所用时间更长，恢复至背景水平后，能量又开始逐渐增加并持续，这主要由延后到达的海啸波作用于海底和海岸而引发。与长周期、短周期双频微地动的平均振幅曲线相比，单频微地动出现了许多尖峰，也印证了单频微地动和双频微地动激发机制的差异：双频微地动的激发机制为两列传播方向相反而频率相同或相近的波，通过非线性波-波相互作用所激发，而单频微地动的激发机制为海浪与海底的直接作用，其激发条件较双频微地动更容易实现，从而表现为对地震较敏感的响应，平均振幅曲线上出现较多与地震对应的尖峰，也不排除其它来源所引起的尖峰，如沉积物中的剪切波共振也会产生时频谱峰值（McCreery et al，1993）。
此外，我们还比较了海啸期间水平地震分量与垂直地震分量的时频谱。以TA.O02D台站记录的三分量记录为例，图3给出了2011年3月10日至15日TA.O02D台站记录的5天时长的水平和垂直分量在不同频带的平均振幅随时间变化的曲线。可见：在地球背景自由振荡频带（0.004—0.007 Hz）内，垂直分量的能量在地震发生前低于水平分量的能量，在地震发生时刻，垂直分量的能量突增，三个分量的能量水平基本相当，震后垂直分量的能量水平逐渐下降，至震后8小时又低于水平分量（图3a）；而微地动（包括微地动中的各频带成分）和高频噪声的水平分量与垂直分量的平均振幅曲线基本重合，能量水平和变化趋势都极为相似（图3b-f）。这表明：在高频噪声和微地动频段，海洋与陆地的耦合在垂向和水平方向上激发的地震波能量相当；对更低频段的地球背景噪声而言，海陆耦合和大地震的发生均能激发这一频段的能量，但是海陆耦合激发的能量在水平方向更强，而此次MW9.0地震激发的能量在垂直方向更强（图3a）。
![]()
图 3 2011年3月10日至15日TA.O02D台站三分量地震连续波形记录在不同频带的半小时平均振幅A的变化曲线（a） 地球背景自由振荡频带（0.004—0.007 Hz）；（b） 单频微地动频带（0.05—0.08 Hz）；（c） 长周期双频微地动频带（0.1—0.15 Hz）；（d） 短周期双频微地动频带（0.18—0.4 Hz）；（e） 微地动频带（0.05—0.5 Hz）；（f） 高频噪声频带（1.3—1.5 Hz）。红色圆圈的意思与图2相同。Figure 3. The variation curves of half-hour average amplitude A of the three-component seismic continuous waveform records for the station TA.O02D for different frequency bands from March 10 to 15，2011（a） The Earth’s background free oscillation frenquency band （0.004−0.007 Hz）；（b） Single-frequency microseism frenquency band （0.05−0.08 Hz）；（c） Long-period double-frequency microseism frenquency band （0.1−0.15 Hz）；（d） Short-period double-frequency microseism frenquency band （0.18−0.4 Hz）；（e） Microseism frenquency band （0.05−0.5 Hz）；（f） The high-frequency noise （1.3−1.5 Hz）. The meanings of the red circles are the same as in Fig. 22.2 海啸期间地球背景噪声的空间特征
前面部分中，我们已经讨论了海啸在地球背景噪声上的时频特征，下面将使用主极化分析方法来研究海啸激发的地球背景噪声的空间特征。尽管由主极化分析方法得到的噪声源方向具有180°的简并性，但是由于地球背景自由振荡和微地动的激发源在海洋，而且对于北美地区的微地动频段，观测一般以瑞雷波为主导（Toksöz，Lacoss，1968；Lacoss et al，1969），因此这里选择指向开放水域的方向为主极化方向（Rhie，Romanowicz，2004；Bromirski et al，2005；Shapiro et al，2006）。
图4为日本MW9.0地震发生前后近岸台站TA.O02D和TA.PTFO附近海域的DART46411台站记录的海底压力变化，可见大约震后7小时海啸波到达该台站。根据地震发生后DART46411台站接收到的海啸预警信号为依据（图4）选取六个时长为半小时的时窗，图5—10分别展示了海啸作用期间所选取三个台站的连续地震波形记录在这六个时窗的主极化分析结果。六个时窗分别为：① 主震到达各地震台站后半小时（图5）；② 海啸首波到达DART46411台站附近海岸后的连续两个半小时（图6和图7）；③ 海啸到达后，随机选取的三个时窗，以显示海啸的主要能量来源方向随时间的变化（图8—10），其中为了对比地震对极化分析的影响，图9给出了针对无M≥5.0地震发生的时窗所进行的极化分析。在每个时窗内，针对所选取台站的地震数据在四个频段内进行极化分析：① 地球背景自由振荡频带0.004—0.007 Hz （图5—10中a1，b1，c1）；② 单频微地动频带0.05—0.08 Hz （图5—10中a2，b2，c2）；③ 短周期双频微地动频带0.18—0.4 Hz （图5—10中a3，b3，c3）；④ 高频噪声频带1.3—1.5 Hz （图5—10中a4，b4，c4）。针对上述六个时窗内四个频段的噪声极化分析显示，噪声的来源均为全方位，但均存在一个主导方向，且主导方向随时间和记录点位置的变化而呈现不同的特征。为了便于研究地震对主极化分析的影响，在每个主极化方向图中标出该时窗内所发生地震到达对应台站的大圆路径方位。
![]()
图 4 日本MW9.0地震前后DART46411台站记录到的2011年3月10日至15日的海底压力值箭头指示台站记录的海啸频段的海浪实测与预测高度残差幅值超过3 cm阈值进入快速汇报模式（提供海啸预警信息）的时刻Figure 4. The seafloor pressure values recorded by the station DART46411 during the 2011 Japanese tsunami from March 10 to 15，2011The arrows indicate the moments when the residuals between measured and predicted waves heightin the tsunami band recorded by the station exceeds the 3 cm threshold and enters the rapid reporting mode （which provides tsunami warning information）![]()
图 5 台站TA.O02D （a）、TA.TPFO （b）和TA.M31A （c）的连续地震波形记录在3月11日6:00—6:30时窗内的极化分析结果红色直线表示在该时间窗口内各频段内连续地震波形记录的主极化方向，蓝色实线表示各频段内连续地震波形记录在各方位上的振幅。第一列中，红色星形表示对应时间窗口内M≥5.0余震相对地震台站的大圆路径方位，绿色曲线表示该时间窗口内连续地震波形记录的R-T水平面内的质点振动图，橙色曲线表示R-Z竖直面内的质点振动图，蓝色虚线表示质点振动图中的R分量方向，红色虚线对R-T和R-Z质点振动图来说分别代表T分量和垂直分量Z的方向Figure 5. Polarization analysis results of the continuous seismic waveforms of the stations TA.O02D （a）， TA.TPFO （b） and TA.M31A （c） recorded from 6:00 to 6:30 on March 11Red lines indicate the predominant polarization direction of the continuous seismic waveforms recorded in each frequency band，and blue solid lines indicate the amplitude of the continuous seismic waveform recorded in each azimuth in each frequency band. In the first column subfigures，the red stars indicate the position of the great circle path of the M≥5 aftershocks relative to the seismic stations in the corresponding time window，green curves indicates the particle vibration diagram in the R-T horizontal plane，orange curves indicate the particle vibration diagram in the R-Z vertical plane. Blue dotted lines indicate the direction of the component R in the particle vibration diagrams，red dotted lines represent the direction of the component T and the vertical component Z of the R-T and R-Z particle vibration diagrams，respectively![]()
图 6 台站TA.O02D （a），TA.TPFO （b）和TA.M31A （c）的连续地震波形记录在时间窗3月11日15:00—15:30内的极化分析结果（图注说明同图5）Figure 6. Polarization analysis results of the continuous seismic waveforms of the stations TA.O02D （a），TA.TPFO （b） and TA.M31A （c） recorded from 15:00 to 15:30 on March 11 （The illustrations are the same as Fig. 5）![]()
图 7 台站TA.O02D （a），TA.TPFO （b）和TA.M31A （c）的连续地震波形记录在时间窗3月11日15:30—16:00内的极化分析结果（图注说明同图5）Figure 7. Polarization analysis results of the continuous seismic waveform of the stations TA.O02D （a），TA.TPFO （b） and TA.M31A （c） recorded from 15:30 to 16:00 on March 11，2011 （The illustrations are the same as Fig.5）![]()
![]()
图 9 台站TA.O02D （a），TA.TPFO （b）和TA.M31A （c）的连续地震波形记录在时间窗3月11日22:00—22:30内的极化分析结果（该时窗内无M≥5.0地震发生，图注说明同图5）Figure 9. Polarization analysis results of the continuous seismic waveforms of the stations TA.O02D （a），TA.TPFO （b） and TA.M31A （c） recorded from 22:00 to 22:30 on March 11 （Note：No M≥5.0 earthquake occurred in the time window，and the illustrations are the same as Fig. 5）![]()
图 10 台站TA.O02D （a），TA.TPFO （b）和TA.M31A （c）的连续地震波形记录在时间窗3月12日15:00—15:30内的极化分析结果（图注说明同图5）Figure 10. Polarization analysis results of the continuous seismic waveforms of the stations TA.O02D （a），TA.TPFO （b） and TA.M31A （c） recorded from 15:00 to 15:30 on March 12，2011 （Illustrations are the same as Fig.5）由图5可见，在日本MW9.0地震发生并到达地震台站后的半小时内，沿岸台站TA.O02D，TA.TPFO和内陆台站TA.M31A在地球背景自由振荡频段的主导极化方向均为NE−SW，并未指向主震地震波到达的大圆路径方向。结合三个台站在该时窗内地球背景自由振荡频段的径向分量R和垂向分量Z平面内以及径向分量R和切向分量T平面内的质点振动曲线可知，在水平面内垂直于径向的切向上的质点的振动能量较强，反映了该地震在这一频段主要激发的地震波是沿切向极化的勒夫波，其极化方向正好与主震地震波到达的大圆路径垂直（图5a1，b1，c1），而在微地动频段，沿岸的台站TA.O02D和TA.TPFO的主极化方向的一致性略差，且二者均偏离主震地震波到达的大圆路径方向：两个台站的单频微地动主极化方向分别指向WNW和WSW （图5a2和b2），双频微地动主极化方向分别指向NW和NS （图5a3和b3），而高频噪声的主极化方向分别指向WNW和NNW （图5a4和b4）。内陆台站TA.M31A的主极化方向在单频微地动、双频微地动和高频噪声频段的能量来源方向的一致性较好，均指向NW−SE，即主震的大圆路径方向（图5c2，c3，c4）。总的来说，内陆台站各频带的主极化方向较沿岸台站一致性更好，这主要是由于与近岸台站相比，内陆台站受地震影响的同时受海洋活动影响较小。而近岸台站由于离海洋激发源较近，而且单频微地动和双频微地动的激发机制不同，源区位置不一定相同，造成近岸两个台站微地动频段的极化方向比较分散。
在日本MW9.0地震发生约7小时后海啸波到达沿岸台站附近的DART46411台站（图4），图6和图7分别给出了海啸波到达后该地震台站记录到的两个连续的半小时地震数据的主极化分析结果。可见：在海啸波到达后第一个半小时即2011年3月11日15:00—15:30，各频段的能量来源方向表现出较好的一致性，台站TA.O02D，TA.PTFO和TA.M31A的主极化方向均为NW （图6），基于此综合三个地震台站的能量来源方向可知，主要能量源区的位置指向北太平洋东北海域；而半小时后，即2011年3月11日15:30—16:00，三个台站的主极化方向分别为SE、NW和SW （图7），综合三个地震台站的能量来源方向可知，主要能量源区的位置指向紧邻沿岸台站的北太平洋东部海域。这与海岸附近的潮汐计的海啸波观测相对应：位于北美西北海岸亚历山大港口（Port Alexander）的潮汐计于2011年3月11日14:09监测到海啸首波，而位于北美西海岸洪堡湾（Humboldt Bay）的潮汐计于2011年3月11日15:41监测到海啸首波（Heidarzadeh，Satake，2013），本文主极化分析所得的噪声源区方位与其一一对应，表明海啸波在地震台站在这些时窗观测到的所有频段的噪声源中占主导地位。需要指出的是，内陆台站TA.M31A在地球背景自由振荡频段测到的主极化方向为ENE−WSW （图6c1），与日本MW9.0地震发生后半小时地球背景自由振荡频段的主极化方向一致（图5c1），而同一时段近岸台站TA.O02D和TA.TPFO在地球自由背景振荡频段的主极化方向指向海啸到达的区域（图6a1和图6b1），这主要是因为近岸台站更容易受到海啸激发的地球背景自由振荡的影响，而内陆台站的主极化方向主要还是受主震激发的环型地球自由振荡面波的影响（图5和图6c1）。这一点可以由内陆台站TA.M31A的质点振动图加以佐证（图6c1），该时间窗口内沿切向T的质点振动略强于指向海啸波到达台站近岸的NW方向，说明该台站由于位于离岸较远的内陆，海啸波激发的地震波的影响小于主震激发的勒夫波的影响，其在地球背景自由振荡频段的主极化方向大致与大圆路径方向垂直，保持着日本MW9.0地震后半小时相同的主极化方向（图5和图6c1）。
为了研究噪声源区位置随海啸波传播的变化，随机选取海啸首波到达后4小时、7小时和12小时后半小时时长的连续地震波形记录，分析各频段能量的主极化方向，结果如图8—10所示。可见在各个时窗内，三个地震台站在各频段的主极化方向均较一致：在海啸首波到达4小时后，即2011年3月11日19:00—19:30，台站TA.O02D，TA.PTFO和TA.M31A的主极化方向分别为近EW，WNW和NNW （图8），其主极化方向指向北太平洋东北海域，说明该海域为这个时窗内各频段的主要能量来源方向；在海啸首波到达7小时后，即2011年3月11日22:00—22:30，三个台站的主极化方向分别为近NW，NW和WNW （图9），其主极化方向指向北太平洋东北海域，说明该海域为这个时窗内各频段的主要能量来源方向；在海啸首波到达24小时后，即2011年3月12日15:00—15:30，三个台站的主极化方向分别为近NS，NW和SW （图10），其主极化方向又指向北太平洋东海域，说明该海域为这个时窗内各频段主要能量的来源方向。结果表明，海啸与海底作用并激发地震噪声的源区位置并非固定不变，而是随时间变化的。这与海啸波本身的传播特性相关，即受海底地形、水深等因素的影响，海啸波的多次反射和衍射导致了海啸波能量集中区所处位置在海域呈现出随时间而不断变化的情况（Munk，1963；Shevchenko et al，2014）。
需要指出的是，在海啸到达DART46411台站后的各个时窗内（图6—8，10），发生的地震均为MW9.0地震的余震，其发生位置接近于主震位置，且震级较小，为M5.0左右。而极化分析得到的以上时窗各频段能量源区的主极化方向各不相同，并不指向地震方位，且三个台站的主极化方向交会于海啸到达的海域。另外，在无M≥5.0地震发生的时窗内（图9），主极化分析结果展现出与存在M≥5.0地震发生的时窗一致的结果，如图8b （有地震发生）与图9b （无地震发生）。综上推断，在海啸波到达近岸后，各频段的主极化方向由海啸主导。
3. 讨论与结论
本文利用北美沿岸和内陆的地震台站记录的连续数据，分析了2011年日本大地震激发的海啸对背景噪声的影响，并结合DART系统记录的海底压力记录和潮汐预测数据，分析了海啸期间在各频段上背景噪声的时空分布特征，主要结果如下：
1） 在此次大海啸期间，高频噪声（1.3—1.5 Hz）受海啸影响最小；海啸对地球背景自由振荡（0.004—0.007 Hz）的影响最大，影响的持续时间也最长；海啸对微地动频带（0.05—0.5 Hz）的影响较小，其中对短周期双频微地动（0.18—0.4 Hz）的影响最小，对单频微地动（0.05—0.08 Hz）和长周期双频微地动（0.1—0.15 Hz）的影响略大，持续时间略长。
2） 海啸到达近岸时，对附近台站各频段的噪声均有影响，是各频段噪声的主控能量来源，且在后续过程中各噪声主控能量来源的位置会随时间变化。这与海啸的传播特性有关，海啸在传播过程中受水深、海底地形以及近岸地形的反射、衍射等的影响，海啸能量聚集区域随时间而变化，并非均匀传播至海岸。
由于海啸现象较为复杂，目前的研究对于海啸的激发机制、传播机制、海啸源的定位等的认识仍具有局限性，并且由地震激发的海啸在时间和空间上均存在不易分离的重叠问题，这些因素严重制约了对海啸特征的进一步研究，这也是下一步工作的重点。
-
![]()
![]()
图 2 山丹台（左）和金银滩台（右）降雨量与表层土壤水分含量（a）、气温与表层土壤温度（b）的相关性
Figure 2. Correlations between precipitation and surface soil water content （a） and air temperature andsurface soil temperature （b） for the stations Shandan （left panels） and Jinyintan （right panels）
![]()
图 3 定西台NS测道（上）和EW测道（下）土壤水分含量和土壤温度拟合获取的地电阻率时间序列（a）及残差曲线（b）
Figure 3. Time series curves of apparent resistivity （a） and residual variance curves （b） at the NS channel （upper） and EW channel （lower） of Dingxi station obtained by polynomial fitting soil water content and soil temperature
![]()
图 4 山丹台、金银滩台和武威台土壤水分含量和土壤温度拟合获取的地电阻率时间序列（左）及残差（右）
（a） 金银滩台EW测道；（b） 山丹台NS测道；（c） 山丹台EW测道；（d） 山丹台N45°W测道；（e） 武威台NS测道；（f） 武威台EW测道
Figure 4. Time series curves of apparent resistivity （left panels） and residual variance curves （right panels） at Shandan，Jinyintan and Wuwei stations obtained by a polynomial fitting of soil water content and soil temperature
（a） The EW channel of Jinyintan station；（b） The NS channel of Shandan station；（c） The EW channel of Shandan station；（d） The N45°W channel of Shandan station；（e） The NS channel of Wuwei station；（f） The EW channel of Wuwei station
![]()
图 5 基于断层虚位错模式计算的门源MS6.9地震前震中区域的面应变及主应变
图中红色为面应变膨胀区，蓝色为面应变压缩区，白色箭头为主张应变，黑色箭头为主压应变
Figure 5. Distribution of surface strain and principal strain calculated using the virtual fault dislocation model in the epicenter before the Menyuan MS6.9 earthquake
The region in red color is the relative expansion area of surfacestrain，and the blue is the compression area of surface strain.The white arrows are the principal tensile strain，andthe black arrows are the principal compressive strain
表 1 2022年门源MS6.9地震周边400 km范围内地电阻率台站的基础信息
Table 1 Basic information of the apparent resistivity stations within 400 km to the epicenter of the 2022 Menyuan MS6.9 earthquake
台站 测道 震中距/km 数据质量 历史数据干扰因素 观测仪器 年变动态 备注 嘉峪关 N50°E 347 差 降雨、线路漏电、公路 ZD8M 无规律 2021年8—9月漏电干扰 N45°W 中 降雨、线路漏电、公路 夏高冬低 山丹 NS 113 优 线路漏电 ZD8B 夏高冬低 EW 优 线路漏电 夏低冬高 N45°W 优 线路漏电 夏低冬高 白水河 NS 345 优 大风 ZD8M 夏低冬高 2021年8月—2022年1月电极故障 EW 差 大风 夏高冬低 金银滩 NS 92 差 大风 ZD8M 无规律 EW 优 大风 夏低冬高 拦隆口 NS 114 差 大风、灌溉 ZD8M 夏低冬高 EW 差 大风、灌溉 夏高冬低 兰州 NS 296 差 地铁、塑料大棚 ZD8M 无规律 EW 差 地铁、塑料大棚 无规律 武威 NS 139 优 仪器故障 ZD8MI 夏低冬高 EW 优 仪器故障 夏高冬低 临夏 NS 310 中 漏电、塑料大棚、金属网 ZD8M 夏低冬高 2020—2021年测区存在施工建设 EW 差 漏电、塑料大棚、金属网 夏低冬高 定西 NS 388 优 漏电、灌溉 ZD8M 夏高冬低 EW 优 漏电、灌溉 夏高冬低 注：蓝色字为最终参与地电阻率异常分析的台站及测道。 
下载: 导出CSV
表 2 地电阻率与土壤水分含量和土壤温度的相关系数
Table 2 Correlative coefficient between apparent resistivity and soil water content and soil temperature
台站 测道 不同深度土壤水分含量与地电阻率的相关系数 不同深度土壤温度与地电阻率的相关系数 0—7 cm 7—28 cm 28—100 cm 100—289 cm 0—7 cm 0—28 cm 28—100 cm 100—289 cm 金银滩 EW −0.32 −0.39 0.08 −0.16 −0.63 −0.68 −0.78 −0.75 定西 NS 0.47 0.42 0.09 −0.25 0.62 0.67 0.77 0.80 EW 0.51 0.43 −0.02 −0.45 0.61 0.65 0.76 0.82 武威 NS 0.12 0.54 −0.33 0.29 −0.35 −0.30 −0.16 0.19 EW −0.03 −0.36 0.75 −0.79 0.01 0.00 −0.05 −0.12 山丹 NS 0.43 0.10 −0.27 −0.30 0.82 0.82 0.75 0.37 EW −0.32 −0.41 −0.34 −0.82 −0.21 −0.24 −0.27 −0.24 N45°W −0.65 −0.70 0.04 −0.22 −0.56 −0.66 −0.80 −0.87 注：蓝色数字为与地电阻率相关系数较高并参与拟合的土壤水分含量和土壤温度。
下载: 导出CSV
表 3 2022年门源MS6.9地震的震源参数（引自潘家伟等，2022）
Table 3 The source parameters of the Menyuan MS6.9 earthquake in 2022 （after Pan et al，2022）
断层中心位置 长度/km 宽度/km 中心深度/km 滑动量/cm 节面Ⅰ 东经/° 北纬/° 滑动角/° 走向/° 倾角/° 101.26 37.77 31 16 4 300 21 284 82
下载: 导出CSV
-
邓起东,张培震,冉勇康,杨晓平,闵伟,楚全芝. 2002. 中国活动构造基本特征[J]. 中国科学(D辑),32(12):1020–1030. Deng Q D,Zhang P Z,Ran Y K,Yang X P,Min W,Chu Q Z. 2003. Basic characteristics of active tectonics of China[J]. Science in China:Series D,46(4):356–372. doi: 10.1360/03yd9032
杜学彬,李宁,叶青,马占虎,闫睿. 2007. 强地震附近视电阻率各向异性变化的原因[J]. 地球物理学报,50(6):1802–1810. doi: 10.3321/j.issn:0001-5733.2007.06.021 Du X B,Li N,Ye Q,Ma Z H,Yan R. 2007. A possible reason for the anisotropic changes in apparent resistivity near the focal region of strong earthquake[J]. Chinese Journal of Geophysics,50(6):1802–1810 (in Chinese).
杜学彬. 2010. 在地震预报中的两类视电阻率变化[J]. 中国科学(地球科学),40(10):1321–1330. Du X B. 2010. Two types of changes in apparent resistivity in earthquake prediction[J]. Science China Earth Sciences,54(1):145–156.
杜学彬,孙君嵩,陈军营. 2017. 地震预测中的地电阻率数据处理方法[J]. 地震学报,39(4):531–548. doi: 10.11939/jass.2017.04.008 Du X B,Sun J S,Chen J Y. 2017. Processing methods for the observation data of apparent resistivity in earthquake prediction[J]. Acta Seismologica Sinica,39(4):531–548 (in Chinese).
冯万鹏,何骁慧,张逸鹏,房立华,Sergey S,张培震. 2023. 2022年青海门源MW6.6地震的发震断层及孕震构造模式[J]. 科学通报,68(2/3):254–270. Feng W P,He X H,Zhang Y P,Fang L H,Sergey S,Zhang P Z. 2023. Seismic faults of the 2022 MW6.6 Menyuan,Qinghai earthquake and their implication for the regional seismogenic structures[J]. Chinese Science Bulletin,68(2/3):254–270 (in Chinese).
高曙德,汤吉,杜学彬,刘小凤,苏永刚,陈彦平,狄国荣,梅东林,詹艳,王立凤. 2010. 汶川8.0级地震前后电磁场的变化特征[J]. 地球物理学报,53(3):512–525. Gao S D,Tang J,Du X B,Liu X F,Su Y G,Chen Y P,Di G R,Mei D L,Zhan Y,Wang L F. 2010. The change characteristics of electromagnetic field before to after Wenchuan MS8.0 earthquake[J]. Chinese Journal of Geophysics,53(3):512–525 (in Chinese).
桂燮泰,关华平,戴经安. 1989. 唐山、松潘地震前视电阻率短临异常图象重现性[J]. 西北地震学报,11(4):71–75. Gui X T,Guan H P,Dai J A. 1989. The short-term and immediate anomalous pattern recurrences of the apparent resistivity before the Tangshan and Songpan earthquakes of 1976[J]. Northwestern Seismological Journal,11(4):71–75 (in Chinese).
李新,黄春林. 2004. 数据同化:一种集成多源地理空间数据的新思路[J]. 科学导报,22(2):13–16. Li X,Huang C L. 2004. Date assimilation:A new means for multi-source geospatial data integration[J]. Science &Technology Review,22(12):13–16 (in Chinese).
李新艳,曾宪伟,卢军,马禾青,崔瑾,卫定军. 2022. 基于断层虚位错模式分析2015年内蒙古阿左旗MS5.8地震前地电阻率变化[J]. 地震,42(2):89–99. doi: 10.12196/j.issn.1000-3274.2022.02.007 Li X Y,Zeng X W,Lu J,Ma H Q,Cui J,Wei D J. 2022. Anomalies of apparent resistivity before the 2015 Alxa Zuoqi MS5.8 earthquake in Inner Mongolia of China based on fault virtual dislocation model[J]. Earthquake,42(2):89–99 (in Chinese).
潘家伟,李海兵,Chevalier M L,刘栋梁,李超,刘富财,吴琼,卢海建,焦利青. 2022. 2022年青海门源MS6.9地震地表破裂带及发震构造研究[J]. 地质学报,96(1):215–231. doi: 10.3969/j.issn.0001-5717.2022.01.018 Pan J W,Li H B,Chevalier M L,Liu D L,Li C,Liu F C,Wu Q,Lu H J,Jiao L Q. 2022. Coseismic surface rupture and seismogenic structure of the 2022 MS6.9 Menyuan earthquake,Qinghai Province,China[J]. Acta Geologica Sinica,96(1):215–231 (in Chinese).
钱复业,赵玉林. 1980. 地震前地电阻率变化十例[J]. 地震学报,2(2):186–197. Qian F Y,Zhao Y L. 1980. Ten examples of changes in earth-resistivity prior to strong earthquakes[J]. Acta Seismologica Sinica,2(2):186–197 (in Chinese).
钱复业,赵玉林,于谋明,王志贤,刘小伟,常思敏. 1982. 地震前地电阻率的异常变化[J]. 中国科学(B辑),12(9):831–839. Qian F Y,Zhao Y L,Yu M M,Wang Z X,Liu X W,Chang S M. 1982. The anomalous changes of apparent resistivity before earthquakes[J]. Science in China:Series B,12(9):831–839 (in Chinese).
钱复业,赵玉林,许同春. 1987. 地电阻率季节干扰变化分析[J]. 地震学报,9(3):289–302. Qian F Y,Zhao Y L,Xu T C. 1987. Analyses of the seasonal variation of disturbance in apparent resistivity[J]. Acta Seismologica Sinica,9(3):289–302 (in Chinese).
钱家栋,陈有发,金安忠. 1985. 地电阻率法在地震预报中的应用[M]. 北京:地震出版社:196−205. Qian J D,Chen Y F,Jin A Z. 1985. Application of Earth-Resistivity Method in Earthquake Prediction[M]. Beijing:Seismological Press:196−205 (in Chinese).
石富强,张国强,方炜,邵辉成,张国苓. 2014. 陕西周至地电台地电阻率年变特征分析[J]. 地震学报,36(6):1113–1123. Shi F Q,Zhang G Q,Fang W,Shao H C,Zhang G L. 2014. Annual variation characteristics of apparent resistivity at Zhouzhi geoelectric station,Shaanxi Province[J]. Acta Seismologica Sinica,36(6):1113–1123 (in Chinese).
解滔,卢军. 2016. 地表固定干扰源影响下地电阻率观测随时间变化特征分析[J]. 地震地质,38(4):922–936. doi: 10.3969/j.issn.0253-4967.2016.04.010 Xie T,Lu J. 2016. Apparent resistivity temporal variation characteristics affected by the fixed disturbance source on surface of measuring area[J]. Seismology and Geology,38(4):922–936 (in Chinese).
解滔,于晨,王亚丽,李美,卢军. 2020. 基于断层虚位错模式讨论2008年汶川MS8.0地震前视电阻率变化[J]. 中国地震,36(3):492–501. doi: 10.3969/j.issn.1001-4683.2020.03.012 Xie T,Yu C,Wang Y L,Li M,Lu J. 2020. Apparent resistivity variations before 2008 Wenchuan MS8.0 earthquake based on fault virtual dislocation model[J]. Earthquake Research in China,36(3):492–501 (in Chinese).
解滔,薛艳,卢军. 2022a. 中国MS≥7.0地震前视电阻率变化及其可能原因[J]. 地球物理学报,65(8):3064–3077. Xie T,Xue Y,Lu J. 2022a. Changes in apparent resistivity and its possible reasons before earthquakes of MS≥7.0 in China[J]. Chinese Journal of Geophysics,65(8):3064–3077 (in Chinese).
解滔,卢军,杜学彬. 2022b. 自适应变化幅度方法提取直流视电阻率中短期异常[J]. 中国地震,38(1):52–60. Xie T,Lu J,Du X B. 2022b. An adaptive variation amplitude method for extracting medium- and short-term anomalies of apparent resistivity[J]. Earthquake Research in China,38(1):52–60 (in Chinese).
解滔,于晨,王亚丽,李美,王中平,姚丽,卢军. 2022c. 2013年岷县-漳县MS6.6地震前通渭台的视电阻率变化[J]. 地震地质,44(3):701–717. Xie T,Yu C,Wang Y L,Li M,Wang Z P,Yao L,Lu J. 2022c. Apparent resistivity variation of Tongwei seismic station before the Minxian-Zhangxian MS6.6 earthquake in 2013[J]. Seismology and Geology,44(3):701–717 (in Chinese).
汪志亮,郑大林,余素荣. 2002. 地震地电阻率前兆异常现象[M]. 北京:地震出版社:45−49. Wang Z L,Zheng D L,Yu S R. 2002. Geoelectric Resistivity Precursor Anomalies of Earthquake[M]. Beijing:Seismological Press:45−49 (in Chinese).
许英才,郭祥云,冯丽丽. 2022. 2022年1月8日青海门源MS6.9地震序列重定位和震源机制解研究[J]. 地震学报,44(2):195–210. doi: 10.11939/jass.20220008 Xu Y C,Guo X Y,Feng L L. 2022. Relocation and focal mechanism solutions of the MS6.9 Menyuan earthquake sequence on January 8,2022 in Qinghai Province[J]. Acta Seismologica Sinica,44(2):195–210 (in Chinese).
杨龙翔,刘学谦,孙召华. 2020. 金属大棚对周口台地电阻率干扰定量分析[J]. 地震工程学报,42(2):447–452. doi: 10.3969/j.issn.1000-0844.2020.02.447 Yang L X,Liu X Q,Sun Z H. 2020. Quantitative analysis of the disturbance to earth resistivity caused by a metal greenhouse at Zhoukou station[J]. China Earthquake Engineering Journal,42(2):447–452 (in Chinese).
姚赛赛,高曙德,陈雪梅,张晓阳,杨超,醴武权. 2023. 2022年青海门源MS6.9地震前地电阻率变化分析[J]. 大地测量与地球动力学,43(2):141–147. Yao S S,Gao S D,Chen X M,Zhang X Y,Yang C,Li W Q. 2023. Analysis of earth resistivity changes before the 2022 Qinghai Menyuan MS6.9 earthquake[J]. Journal of Geodesy and Geodynamics,43(2):141–147 (in Chinese).
张国苓,王庆林,乔子云,贾立峰,张素欣. 2017. 地表点电流对地电阻率观测的影响[J]. 中国地震,33(1):122–128. doi: 10.3969/j.issn.1001-4683.2017.01.012 Zhang G L,Wang Q L,Qaio Z Y,Jia L F,Zhang S X. 2017. Influence from the surface disturbance current on geoelectrical resistivity observation[J]. Earthquake Research in China,33(1):122–128 (in Chinese).
张金铸,陆阳泉. 1983. 不同三轴应力条件下岩石电阻率变化的试验研究[J]. 地震学报,5(4):440–445. Zhang J Z,Lu Y Q. 1983. An experimental study on the variation of rock resistivity under triaxially different stresses[J]. Acta Seismologica Sinica,5(4):440–445 (in Chinese).
张学民,李美,关华平. 2009. 汶川8.0级地震前的地电阻率异常分析[J]. 地震,29(1):108–115. doi: 10.3969/j.issn.1000-3274.2009.01.014 Zhang X M,Li M,Guan H P. 2009. Anomaly analysis of earth resistivity observations before the Wenchuan earthquake[J]. Earthquake,29(1):108–115 (in Chinese).
赵和云,张文孝. 1985. 银川台地电阻率趋势变化及年变的分析[J]. 地震,(1):32–38. Zhao H Y,Zhang W X. 1985. Analysis of potential and annual variations of earth resistivity at Yinchuan station[J]. Earthquake,(1):32–38 (in Chinese).
赵玉林,钱复业,杨体成,刘建毅. 1983. 原地电阻率变化的实验[J]. 地震学报,5(2):217–225. Zhao Y L,Qian F Y,Yang T C,Liu J Y. 1983. Experiments in situ of electrical resistivity changes[J]. Acta Seismologica Sinica,5(2):217–225 (in Chinese).
赵玉林,卢军,李正南,钱复业,张洪魁. 1996. 唐山地震应变-电阻率前兆及虚错动模式[J]. 地震学报,18(1):78–82. Zhao Y L,Lu J,Li Z N,Qian F Y,Zhang H K. 1996. The precursors of strain-resistivity and virtual-dislocation model before Tangshan earthquake[J]. Acta Seismologica Sinica,18(1):78–82 (in Chinese).
Brace W F,Orange A S,Madden T R. 1965. The effect of pressure on the electrical resistivity of water-saturated crystalline rocks[J]. J Geophys Res,70(22):5669–5678. doi: 10.1029/JZ070i022p05669
Dobrovolsky I P,Zubkov S I,Miachkin V I. 1979. Estimation of the size of earthquake preparation zones[J]. Pure Appl Geophys,117(5):1025–1044. doi: 10.1007/BF00876083
Hillel D. 1998. Environmental Soil Physics[M]. San Diego:Academic Press:771.
Lin J,Stein R S. 2004. Stress triggering in thrust and subduction earthquakes and stress interaction between the southern San Andreas and nearby thrust and strike-slip faults[J]. J Geophys Res:Solid Earth,109(B2):B02303.
Lu J,Qian F Y,Zhao Y L. 1999. Sensitivity analysis of the Schlumberger monitoring array:Application to changes of resistivity prior to the 1976 earthquake in Tangshan,China[J]. Tectonophysics,307(3/4):397–405.
Lu J,Xue S Z,Qian F Y,Zhao Y L,Guan H P,Mao X J,Ruan A G,Yu S R,Xiao W J. 2004. Unexpected changes in resistivity monitoring for earthquakes of the Longmen Shan in Sichuan,China,with a fixed Schlumberger sounding array[J]. Phys Earth Planet Inter,145(1/2/3/4):87–97.
Lu J,Xie T,Li M,Wang Y L,Ren Y X,Gao S D,Wang L W,Zhao J L. 2016. Monitoring shallow resistivity changes prior to the 12 May 2008 M8.0 Wenchuan earthquake on the Longmen Shan tectonic zone,China[J]. Tectonophysics,675:244–257. doi: 10.1016/j.tecto.2016.03.006
Oki T,Kanae S. 2006. Global hydrological cycles and world water resources[J]. Science,313(5790):1068–1072. doi: 10.1126/science.1128845
Thomas J F, Christopher H S. 1971. Mechanism of underthrusting in southwest Japan: A model of convergent plate interactions[J]. J Geophys Res, 76（29）: 7260–7292. Thomas J F,Christopher H S. 1971. Mechanism of underthrusting in southwest Japan:A model of convergent plate interactions[J]. J Geophys Res,76(29):7260–7292.
Wan Y G,Shen Z K,Bürgmann R,Sun J B,Wang M. 2017. Fault geometry and slip distribution of the 2008 MW7.9 Wenchuan,China earthquake,inferred from GPS and InSAR measurements[J]. Geophys J Int,208(2):748–766. doi: 10.1093/gji/ggw421
Zhu Q,Nie X F,Zhou X B,Liao K H,Li H P. 2014. Soil moisture response to rainfall at different topographic positions along a mixed land-use hillslope[J]. Catena,119:61–70. doi: 10.1016/j.catena.2014.03.010
-
期刊类型引用(0)
其他类型引用(2)
计量
- 文章访问数: 112
- HTML全文浏览量: 35
- PDF下载量: 64
- 被引次数: 2
