
本期概要
1、琉球南部弧前浅层慢滑现象的大地测量证据
研究利用25年GNSS数据,在琉球南部弧前识别出2009-2011年一次相当于Mw 6.8的浅层慢滑事件(SSE),其与一次Mw 6.1板间地震时空耦合,可能受其库仑应力触发。结果表明浅部俯冲界面可积累并缓慢释放应变,强调近海沟海底大地测量对识别潜在海啸源区的重要性。
2、多段断层几何形状和复杂滑动分布的同步贝叶斯估计:应用于2024年能登半岛地震
研究提出一种高效贝叶斯反演算法,可同步估计多段断层几何形状和复杂滑动分布。该方法应用于2024年能登半岛地震,三段断层模型最优,显著提升与余震分布及形变观测的一致性,凸显多段建模对强震机制解析的关键价值。
3、基于分布式声学传感子阵列波束成形的地震定位及其对地震预警的意义
针对线性光纤定位歧义问题,提出利用多段DAS子阵列进行时域波束成形,融合P/S波反方位角与到时差实现稳健震中定位。实测小震验证其有效性,并指出该方法可作为自动震相拾取器,具良好实时地震监测和预警应用潜力,尤其利于提升预警时效性与覆盖密度。
琉球南部弧前浅层慢滑现象的大地测量证据
慢滑事件(SSE)是俯冲带中普遍存在的现象,但是,与典型的突然释放能量的地震不同,慢滑事件会在数天至数月内逐渐释放能量,并且通常无法在没有精密仪器的情况下被探测到。SSE可能会影响俯冲带中引发海啸的大型地震发生的时间和地点,但准确追踪浅层SSE的时空分布仍然具有挑战性。来自日本等国的学者利用25年的全球导航卫星系统的数据,研究了菲律宾海板块西部边界南琉球弧前可能发生浅层SSE的情况。该研究在宫古群岛附近发现了向海沟的瞬态运动,2009至2011年间累计位移约为10 mm,表明浅层SSE释放的测地矩(geodetic moment)相当于Mw6.8。有趣的是,这些慢滑事件与板块边界上发生的6.1级板间地震同时发生,且这些SSE可能发生在板间地震的上倾位置,与板间地震相关的正库仑应力变化可能促进了慢滑,表明它们之间可能存在联系。研究结果表明,浅层的南琉球板块界面可以积累弹性应变,并以无震或地震方式释放能量。因此认为,近海沟的海底大地测量对于全面表征南琉球弧前下方断层的滑动行为具有重要意义。相关研究发表在2025年7月的Geophysical Research Letters。
SSE指滑移速率低于常规地震的断层破裂幕,其在俯冲带中起到分担部分板块运动的作用。由于破裂速度较慢,SSE几乎不辐射地震波,因此需要连续的大地观测数据才能检测到。SSE通常可分为深部SSE和浅部SSE,其中,浅部SSE发生在靠近俯冲带的浅层,检测难度较大。尽管如此,浅部SSE可能反映了巨型逆冲断层最浅部分的摩擦特性,而这些区域是潜在的大规模地震和海啸的源头。虽然海底大地测量观测对于表征浅层SSE至关重要,但SSE持续时间的巨大差异给其探测带来了额外挑战。
研究人员聚焦于南琉球俯冲带,该区域的菲律宾海板块以每年84~88 mm的速度斜向俯冲至扬子地块,冲绳海槽持续的弧后裂谷作用导致南琉球弧向南移动,使南琉球边缘的缩短速率高达130 mm/yr。目前,已有该地区关于深部SSE和甚低频地震(VLFE)及低频地震(LFE)的研究,表明存在一定的滑移行为,但弧前区SSE的潜在发生可能性仍未被深入探究,且该地区尚未发现持续时间超过约1个月的SSE,所以需要进一步探索该地区的慢滑现象。
研究人员通过分析1997年至2023年间全球导航卫星系统(GNSS)的数据,探讨了位于日本琉球群岛南部的SSE。为了揭示慢滑事件的分布,研究使用了滑移反演方法,将板块几何形态离散化为多个子断层,并通过弹性半空间模型与GNSS观测到的位移数据进行反演。另外,估算了三个时期的瞬态位移:2009.5—2011.0、2009.5—2010.0和2010.5—2011.0,分别将其称为2009—2011瞬态、SSE1和SSE2。

图1 南琉球边缘的构造背景
1. 2009—2011年的瞬变
2009—2011年期间,GNSS观测结果显示,相对于扬子地块,宫古群岛存在显著的南向运动。在这一阶段,宫古群岛还表现出额外的7~10 mm的向海沟运动,相当于在1.5年间速度增加了5~7 mm/yr,占该区域长期运动速率的12%~20%,表明存在一次明显的瞬变事件。然而,以往的研究虽然探讨了琉球俯冲SSE的时空分布,但均未能完整捕捉2009—2011年的这一特殊过程。相比之下,该研究依托足够长时间的GNSS数据,并选择了合适的参考站点,因而能够更有效地识别出该时期的SSE及其整体过程。
2. 瞬变源
2009—2011年的瞬变事件可通过板块界面上的滑移来解释。分布滑移模型因GNSS站点空间覆盖不足可能出现模糊,但与均匀滑移模型存在广泛一致性,表明了模型的可靠性。分布滑移模型的滑移区域较大,最大滑移为33 mm,而均匀滑移模型的滑移区域较小,滑移为55 mm,滑移角约为110°。SSE1和SSE2的单独反演结果显示,滑移分布与整体瞬变一致,从均匀滑移模型得到的SSE1和SSE2的应力降符合典型SSE值。
虽然冲绳海槽弧后扩张可以解释瞬变事件,但它的应力场方向与2009—2011年的南北方向运动不符,且弧后地震活动未显著增加。SSE1可能与弧前下方的地震EQ1(发生在2009年的地震,Mw6.1)有关,EQ1震源机制和深度表明它是板间地震,滑移角约120°,符合反演结果。这表明,2009—2011年瞬变的源头是弧前板块界面下的SSE。
合并SSE1和SSE2的矩释放,重新估算的总测地矩相当于Mw6.84—6.87。尽管分辨率有限,反演结果显示出慢滑的横向迁移和滑移方向的变化。SSE1和SSE2相隔约6个月,类似于其他地区的SSE现象,表明SSE1可能促进了SSE2的发生。

图2 GNSS位置的时间序列(2005—2013年)和研究区域的地震活动性。(a–c)使用2005—2009年数据去趋势后的三维时间序列;(d)2009—2011年瞬态期的位移;(e)堆叠的GNSS时间序列和聚类中归一化的地震计数
3. 潜在机制
考虑到地震EQ1和SSE1的时空接近性,研究认为,EQ1可能促进了SSE1期间的无震滑移。为了评估EQ1的余滑是否能解释观察到的SSE现象,研究分析了2009—2011年瞬变期间的无震/地震矩比(ASMR)。结果表明,2009—2011年的瞬变不能仅归因于EQ1地震的余滑。然而,估算的ASMR仍然低于常规SSE的典型范围,表明2009—2011年的瞬态事件可能与板间地震EQ1相互作用,类似于秘鲁北部的地震-无震同步序列。应力变化分析显示,EQ1引发的应力变化在SSE推测区域下盘产生了最大30 kPa的剪切应力增量,这支持了EQ1促进SSE的假设。虽然这些应力变化不如潮汐诱发的应力显著,但SSE的成核阶段可能持续超过一天,使其对静态应力变化更为敏感。
另外,EQ1所在的摩擦闭锁区可能在上倾条件稳定的板块界面处形成了应力阴影,使得弹性应变得以积累。然而,EQ1可能并未有效弱化运动学闭锁的浅层界面,以触发动态破裂,最终导致无震滑移事件。通过之后的SSE,后续几年的地震活动显著减少,表明应力已得到充分释放。
EQ2(发生在2010年的地震,Mw6.3)的发生与SSE2时间不符,并且EQ2主要为板内事件。但是,由于SSE2使EQ2处的库仑应力降低了30~40 kPa,研究推断EQ2要么独立于SSE2,要么可能通过流体迁移与其相关。进一步研究这一机制需要更多的观测数据。

图3 在不同时间段反演的滑移。左图显示观测到的和模拟的水平位移,右图表示反演的滑移模型,其中分布式滑移模型以颜色编码,蓝色矩形框显示均匀滑移模型,白色矢量表示滑移方向
4. 对地震灾害的影响
研究推测的慢滑移区域与1771年Mw8.0—8.5的地震源接近,表明它们可能存在相互作用。许多慢滑移区域与高自由空气重力场重叠,这些区域被解释为俯冲的横向脊,与吕宋-冲绳断裂带平行。被俯冲脊带入的富流体物质可能降低有效法向应力,并增强板块界面上的岩性横向异质性,形成一个浅层的条件稳定区域,有利于SSE的发生。琉球群岛南部边缘的俯冲脊可能为浅层SSE和慢地震提供传导条件,尽管其中许多可能未被目前的观测网络探测到。俯冲脊和1771年地震源的空间互补性表明,推测的SSE易发区可能作为地震破裂的屏障,类似于其他俯冲带的观察结果。然而,另一些研究表明,如果这些区域的有效应力足够低,或由于热效应而增强动态减弱,地震破裂可能会突破SSE易发区。

图4 地震活动性时空分布及推断的SSE。(a)地震活动性和重新定位的低频事件,红色和蓝色虚线标记了弧前、弧内和弧后地震的边界,黑色虚线标记了低频活动带;(b)时空图,红色、蓝色和灰色圆圈分别代表弧前、弧后和弧内的地震
来源:Tang C H, Hsu Y J, Okada Y, et al. Geodetic evidence of shallow slow-slip phenomena beneath the southern Ryukyu forearc. Geophysical Research Letters, 2025, 52(13): e2025GL114742.
多段断层几何形状和复杂滑动分布的同步贝叶斯估计:应用于2024年能登半岛地震
对于地壳形变数据的反演,需要同时估计断层几何形状和滑动分布,人们尝试在贝叶斯推断框架内同时估计它们。在这些方法中,由于概率分布的设定,无法获得参数后验分布的解析表达式。蒙特卡罗方法通常用于通过多个样本来近似概率分布,然而,随着参数维数的增加,蒙特卡罗抽样的计算成本会变得很高,使得将这些方法应用于具有复杂破裂模式的大型地震的逆分析(其中待估计的参数是高维的)具有挑战性。来自日本的研究者发明了一种高效的算法,用于同时对断层几何形状和滑动分布进行贝叶斯估计,该算法对高维参数具有鲁棒性,有助于估计多段断层几何形状和复杂的滑动分布。在高维参数的数值试验中,研究人员将所开发方法与传统方法进行比较,发现估计结果的准确性和收敛性得到了提高,同时降低了计算成本。作为所开发方法的一个应用示例,研究根据2024年能登半岛地震期间观测到的地壳变形,估算了断层的几何形状和滑动分布,并考虑了单段、双段和三段断层模型。通过比较分析结果,证明三段模型最为合理,这强调了改进估算方法对于准确估算具有复杂滑动分布的多段断层几何形状的重要性。相关研究发表在2025年8月的Geophysical Journal International。
在地震机制研究中,反演分析通过大地测量数据估算断层滑移,可帮助研究库仑应力变化、地震运动、海啸等现象。对于难以通过地下结构确定断层几何形状的地震,通常采用两步法:首先,假设滑移分布均匀,通过反演数据估算断层几何形状;然后,在固定几何形状的条件下,估算滑移分布。然而,这种方法无法保证最终的断层几何形状与滑移分布一致,且均匀滑移假设可能引入偏差。因此,贝叶斯方法得到应用,它能同时估计断层几何形状和滑移分布,将两者视为随机变量,通过计算后验概率分布来优化模型。
然而,当断层破裂模式复杂且由多个片段组成时,参数维度的增加使得贝叶斯估计变得更加困难。在高维情况下,蒙特卡罗采样的计算成本极高。与此同时,合成孔径雷达(SAR)等新技术的进步促进了高密度地壳形变的探测,为大规模地震滑移的详细估算提供了宝贵数据。
研究人员开发了一种同时进行断层几何形状(θ)和滑移分布(s)的贝叶斯估计方法,并利用从全球导航卫星系统(GNSS)和SAR获得的地壳形变观测数据,对2024年能登半岛地震的断层几何形状和滑动分布进行了贝叶斯估计。在模型参数化上,断层滑移被建模为由断层几何形状和滑移分布两个主要部分组成。断层通过多个矩形断层段进行表示,每个段的几何形态由其位置、角度和尺寸来定义。概率分布的设定考虑了不同数据源的观测误差。

图1 2024年能登半岛地震中观测到的地壳变形,基于GNSS和SAR数据。(a)、(b)基于GNSS数据得到的水平和垂直位移,白色方块表示GNSS观测点;(c)基于SAR数据得到的视线位移
研究人员开发的方法具有两个特点:①采用自适应调整的反射哈密尔顿蒙特卡罗(ReHMC)方法,其结合了哈密尔顿蒙特卡罗方法和反射机制,使得采样过程更加高效,特别适用于带有约束条件的高维滑移分布问题。在s维数设置为300的情况下,进行的边际似然计算数值试验表明,即使样本量达到105,Metropolis核无法准确计算边际似然。相反,研究人员开发的核在使用样本量为103的实际观测数据进行分析时,达到了足够的精度。②对于θ的采样,采用了改进的广义过渡马尔可夫链蒙特卡罗(GTMCMC),其中高斯分布近似于后验概率分布,被当作重要性分布q(θ)来使用。与q(θ)=P(θ)(先验概率分布)的情况相比,该算法提高了θ样本分布的收敛性,减少了连续蒙特卡罗(SMC)每次迭代的计算时间,提高了采样的并行性能,并降低了总体计算成本。

图2 断层模型参数化示意图。(a–c)单段、双段和三段断层模型的断层迹线参数化;(d)沿倾角方向的平面图,显示地下几何形状的参数;(e)正滑动方向的定义,滑动方向受分析中非负约束的限制
针对2024年能登半岛地震,研究人员采用了所开发的能够表示复杂断层几何形状和滑动分布的模型。结果表明,该方法能够推导出一个合理的断层模型,并通过在贝叶斯推断框架内整合来自滑动约束和余震分布的先验信息,以及观测到的地壳形变数据,来量化其不确定性。对具有一个、两个和三个断层段的模型进行比较发现,增加断层段的数量可以提高所估计断层几何形状与余震分布之间的一致性,以及计算和观测到的地壳形变之间的一致性。此外,统计模型选择依次倾向于具有三个、两个和一个断层段的模型。虽然三个以上的断层段可能更可取,但这些发现凸显了改进估算方法(以使其能够解释多个断层段)的重要性。
1. 断层段数的统计模型选择
研究分别应用了一个、两个和三个断层段的模型来估算2024年能登半岛地震的断层几何形状和滑移分布。为了评估哪个模型最为合理,进行了模型的比较选择。
首先,通过比较每个模型的P(d)(衡量模型对观测d的解释力),来评估模型的拟合度。在该研究中,P(d)用作模型选择标准,用于衡量模型与观测数据的一致性。根据计算得到的P(d)结果,三段模型的拟合度最佳,表明该模型能够最好地描述地震发生区域的断层几何形状和滑移分布。
然后,研究使用赤池贝叶斯信息量准则(ABIC)来进一步验证每个模型的适用性。该标准考虑了模型的解释力和复杂性之间的权衡。结果显示,三段模型的ABIC值最小,说明该模型相较于其他模型具有更高的统计优越性。这表明,从统计模型选择的角度看,在2024年能登半岛地震的断层模型中纳入多个断层段是可行的。
2. 不同采样核的P(d|θ)计算准确率
研究人员评估了使用ReHMC核函数时,边际似然P(d|θ)是否能够被足够准确的计算,以及使用ReHMC核函数时计算精度的提高(相对于Metropolis核函数)。数值实验的结果表明,增加ReHMC核中的K(样本数量)和nτ(蛙跳时间积分步数),或增加Metropolis内核中的K,计算精度都会提高,但计算时间均会有所增加。对于Metropolis核,随着K的增加,误差的收敛速度较慢,需要较大的样本量才能达到较高水平准确度。相反,通过增加nτ和K,ReHMC核可以实现P(d|θ)的更准确计算。

图3 三段断层模型的地壳形变观测值与计算值的比较
3. 重要性分布q(θ)的选择对θ采样的影响
在该研究中,分别使用了两种不同的q(θ),一种是基于先验分布P(θ)的重要性分布,另一种是基于已获得样本的高斯分布。后者通过对已有样本进行均值和协方差计算来近似真实的后验分布。实验结果表明,与使用P(θ)作为重要性分布相比,使用基于样本的高斯分布能显著提高采样效率,减少计算所需的迭代次数。研究人员还分析了不同q(θ)在不同采样样本数下的表现,结果表明,采用样本均值和协方差来定义q(θ)时,可以显著改善采样的收敛速度,尤其是在处理高维参数空间时。通过优化重要性分布,马儿可夫链蒙特卡罗(MCMC)采样的效率得到了显著提高,有助于在复杂模型中得到更准确的估计。

图4 估计的θ后验分布,该分布由GTMCMC获得的样本直方图表示。(a)单段断层模型;(b)双段断层模型;(c)三段断层模型
来源:Nakao K, Ichimura T, Munekane H, et al. Simultaneous Bayesian estimation of multi-segment fault geometry and complex slip distribution: application to the 2024 Noto Peninsula earthquake. Geophysical Journal International, 2025: ggaf231.
基于分布式声学传感子阵列波束成形的地震定位及其对地震预警的意义
地震定位是地震监测和地震预警(EEW)系统的重要组成部分。新兴的分布式声学传感(DAS)技术能够以前所未有的空间分辨率在长距离上对应变率进行阵列测量。因此,它在地震定位评估方面有着巨大的应用前景。然而,许多现有的用于地震学应用的光纤都是线性几何结构,经常导致位置估计不准确。来自以色列的学者展示了在光纤几何形态不理想情况下,利用DAS原生应变率数据进行波束成形,实现震中定位的可行性。研究表明,应用多个短DAS子阵列的时域波束成形和相位到达时差可以对震源位置进行稳健的评估。光纤上多个段的波束和相位到时的聚合可以约束震中位置。研究利用在66 km光纤上记录的两个局部小震来演示该方法,并进一步展示了该方法在实时EEW中的适用性。相关研究发表在2025年7月的Seismological Research Letters。
近年来,分布式声学传感(DAS)在地球科学领域快速发展。DAS可利用光纤获取高时空分辨率的连续应变率数据,陆地和海上光纤的全球普及进一步确立了DAS作为地震研究用途的最先进测量工具的地位。DAS在地震预警(EEW)中潜力巨大,可在断层附近实现密集监测,增加预警的提前时间。但是,目前其动态范围有限,强震信号易饱和,这构成了重大挑战。
在震源定位方面,已有研究探索基于走时、振幅及波束成形的方法。其中,波束成形最具优势。它利用DAS的密集阵列特性,无需速度模型与绝对时间戳,便于系统整合,还能作为速度滤波器,突出体波、抑制面波,并能同时完成检测与定位,减少人工震相拾取的不确定性。然而,DAS波束成形仍面临多重困难:信号相干性较低,光纤位置和几何不确定性会带来偏差;其方向性响应和单分量特征限制了对不同震相的利用;多数光纤为线性布设,难以优化阵列形态,导致定位歧义。因此,亟需发展适应既有光纤网络的框架,以更好发挥DAS在地震监测与预警中的潜力。
研究人员提出了一种专为任意几何形状光纤设计的地震震中定位方案,该方案易于实时实施。研究方法基于时域波束成形,结合了光纤沿线多个段的P波和S波反方位角(BAZ)估计值,以及P波和S波的到时差(作为震中距的替代值)。最后,使用两次当地地震的高信噪比(SNR)记录验证了该方法。

图1 (a)研究区地图,绿线表示光纤,红星表示2024年3月13日3.6级地震和2024年8月13日3.1级地震的震中位置;(b)3.6级地震的归一化DAS应变率数据;(c)3.1级地震的归一化DAS应变率数据
1. 波束成形方法讨论
在光纤几何结构可能不理想的情况下,研究利用DAS原生应变率数据进行了震中地震定位。具体而言,准线性光纤很丰富,但线性阵列几何常因慢度(SLO)和反方位角(Φ)而存在较大歧义,不利于定位。DAS的两大优势(高空间采样率与长距离探测)为克服此缺陷提供了条件:高密度采样保证了在相对短距离内具备足够多的通道,从而获得较好相干性;长距离覆盖则允许将光纤划分为多个子段,从而在震源相对位置上获得多方位的测量。尽管单一光纤段的BAZ存在歧义,但多段组合后的结果显示震中定位并无明显偏差,相较于基于走时的定位方法,可显著减少结果的模糊性。
该研究方法的关键在于光纤的分段策略。研究人员采用了简化方式,即保持恒定段长(约4.55 km)又设置重叠(250个通道重叠),使得部分段呈单向几何,部分段则呈多向几何。单向段的波束在地理坐标映射中贡献较大的泛化区域,而多向段的波束则提供更集中的定位区域。虽然未明确引入波束加权,但该方法天然赋予多向段更高的权重,从而提升了最终定位精度。因此,光纤分段是提高准线性光纤定位有效性的关键步骤。
波束成形在EEW中的可靠性取决于对直接P、S波的良好记录。然而,由于这些高速体波的应变率幅度较低,其相干性可能不足,甚至低于仪器噪声,从而限制波束成形效果。此前,有学者提出通过与同址地震仪结合,将DAS数据在空间上积分为地面运动,以增强直达体波相干性。但这种方法受限较大:同址地震仪在现有商用或暗光纤中罕见,且光纤不完全笔直时的积分会引入误差。幸运的是,对于EEW相关的中强震,直达体波的记录质量预计优于研究人员所分析的小震。因此,研究认为,波束成形方法直接作用于原始DAS数据即可满足EEW需求。
此外,尽管研究采用人工切割方式获取P波和S波时窗,但从概念上展示了波束成形在时域内可以作为震相自动拾取器。由于波束成形能聚焦于有限慢度范围,本质上相当于一种速度滤波器,能够稳健地区分地震高速体波与DAS常见的低速信号,进一步凸显其在地震监测与预警中的应用潜力。

图2 2024年3月13日地震S相到达的单光纤段波束成形。(a)按标准差缩放的应变率数据;(b)光纤段的几何形状;(c)对于固定的慢度值,相似度作为时间和方位角的函数;(d)对于固定的方位角值,相似度作为时间和慢度的函数;(e)所有Φ-SLO对在t相时的相似度(分别为角坐标和径向坐标)
2. 实时应用
在实时应用方面,研究提出的波束成形定位方法相较于其他定位方法具有诸多优势,这使得它在EEW实施中颇具吸引力。对于已投入使用的EEW系统,所提出的算法可以提供分别对应于最大震级和最小震级的最小距离和最大距离。一旦获得初步的地震位置和震级估计值,就可以预测地震动强度并发出警报。
对两次地震位置反演过程的伪实时实现结果显示:在初始阶段,首批P波波束即可提供震源大致方向;当全部P波波束到达时,若地震位于光纤网络范围内,已能较好地约束震中位置;但若地震发生在网络外,仅依靠P波仍难以准确确定距离。随后,S波波束及其震相拾取的加入与已有P波拾取结合,可利用到时差估算震中距,从而大幅提高定位精度。不论“网内”(in-network,从整个光纤到目录震中位置的BAZ范围较宽)还是“网外”(off-network,从整个光纤到目录震中位置的BAZ范围较窄)事件,最终在所有S波波束到达后,定位结果均明显收敛并接近目录位置。

图3 “网内”事件的实时定位。(a)DAS应变率数据,圆圈表示基于相似度的相位拾取;(b)由P波束、S波束和S-P到达时差组成的定位得分图,白色标记表示光纤位置
要使该定位方法真正适用于EEW的实时应用,还存在两个关键问题需要解决。首先,震相自动拾取问题。文中示例的拾取窗口是人工预先设定的,用于分别包含P波和S波。若要实现实时应用,可以在每个光纤段的相干性时间序列上运行长短时均值比(STA/LTA)算法,而不是直接作用于原始数据。基于相干性的检测有三大优势:其一,波束成形能天然滤除非相干噪声;其二,受限的慢度范围可同时滤除非体波震相和局部噪声;其三,可推算出两个相近时相(同一光纤段)的慢度差异,有助于区分P波与S波。

图4 “网外”事件的实时定位。(a)DAS应变率数据,圆圈表示基于相似度的相位拾取;(b)由P波束、S波束以及S-P到达时差组成的定位得分图
其次,固有延迟问题。由于所有相似度相关计算都需在一个保守缓冲时间之后进行,而该时间等于光纤段长度乘以最大分析慢度。缩短延迟的办法包括减小段长或降低最大慢度,但这两者都有局限:降低最大慢度可能导致对S波的漏检,而S波对“网外”事件的定位至关重要;缩短段长则需谨慎,因为段长必须足够覆盖目标地震可能产生的最长波长。另一种潜在的实时策略是忽略计算出负时移的通道,具体实现仍有待未来工作的进一步探索。
来源:Ben-Zeev S, Lior I. Earthquake location with distributed acoustic sensing subarray beamforming with implications for earthquake early warning. Seismological Research Letters, 2025.
主送:中国地震局领导
编发:中国地震局地震预测研究所





