基于同震库仑应力变化的汶川地震余震频次研究

在计算同震静态库仑应力变化的基础上,基于速率-状态摩擦定律及与同震静态库仑应力变化相关的余震频次估计方法 ,研究了汶川地震余震区及附近不同区域内与汶川地震同震静态库仑应力触发相关的“直接”余震频次。结果显示,主震破裂带尤其是主震破裂带南段的“直接”余震频次明显低于实际情况,其原因在于同震库仑应力变化导致的主震破裂面上应力水平的降低。研究结果还显示,汶川地震同震库仑应力变化对余震活动持续时间的影响与震级下限有关。对ML4以上余震,持续时间约为震后15~16个月;对ML3.5以上余震,接近60个月。在上述时段内,与同震库仑应力变化相关的“直接”余震占全部余震的比例为44.7%～48.6%。这表明,即使在同震库仑应力变化的“有效”作用时段内,主震破裂面上也大约有50%的余震活动不是缘于同震库仑应力变化的影响,这可能与震后余滑及粘弹松弛等时间相关因素的影响有关。     

基于同震库仑应力变化的汶川地震余震频次研究

在计算同震静态库仑应力变化的基础上,基于速率-状态摩擦定律及与同震静态库仑应力变化相关的余震频次估计方法,研究了汶川地震余震区及附近不同区域内与汶川地震同震静态库仑应力触发相关的"直接"余震频次。结果显示,主震破裂带尤其是主震破裂带南段的"直接"余震频次明显低于实际情况,其原因在于同震库仑应力变化导致的主震破裂面上应力水平的降低。研究结果还显示,汶川地震同震库仑应力变化对余震活动持续时间的影响与震级下限有关。对ML4.0以上余震,持续时间约为震后15~16个月;对ML3.5以上余震,接近60个月。在上述时段内,与同震库仑应力变化相关的"直接"余震占全部余震的比例为44.7%~48.6%。这表明,即使在同震库仑应力变化的"有效"作用时段内,主震破裂面上也大约有50%的余震活动不是缘于同震库仑应力变化的影响,这可能与震后余滑及粘弹松弛等时间相关因素的影响有关。     

1976年唐山强震群震后库仑应力演化及其与2020年古冶5.1级地震的关系

采用能够综合协调长期变形和震后短期变形的Burgers流变模型,模拟了1976年唐山强震群引起的震后形变场以及同震和震后库仑应力变化。结果显示:1976年唐山强震群中主震的两个破裂面以及滦县和宁河两次强余震均对2020年古冶5.1级地震表现为库仑应力加载。岩石圈粘弹性松弛效应引起的库仑应力变化显示,震后15年前后,库仑应力演化状态呈现显著的差异性:震后15年内,库仑应力变化剧烈;而震后15年后,库仑应力呈现缓慢的稳定增加状态。该过程与唐山强震群余震区地震活动过程相似,可能暗示1976年唐山强震群余震区应力调整过程已基本稳定。     

断层面上库仑破裂应力变化的相关研究进展

地震的孕育、发生和震后调整被认为是震源区及其邻区的应变积累释放过程.地震引起的静态库仑应力变化可有力地解释余震分布、强震序列等地震观测,并为探索地震发生机制和地震预测提供新的线索.地震学家们多采用弹性位错模型计算同震位错引起的库仑应力变化,以研究余震分布及地震间相互作用;随着流变学的发展,粘弹模型因可以很好的解释大量震后形变观测而被广泛用于断层面上的震后应力调整研究;此外,构造应力加载作用在更长时间尺度上造成断层面的应力积累,基于负位错模型计算震间库仑应力变化为又一研究内容.本文从同震、震后、震间库仑应力变化的角度提炼出断层面上库仑破裂应力变化的大量研究成果,介绍了库仑应力变化的基本计算原理及引起断层面上库仑应力变化的主要原因,论述了通过库仑应力演化来研究地震活动的方法和应用,进一步讨论断层面上库仑应力演化研究存在的问题和近期的进展.     

2015年尼泊尔M_W7.9地震对青藏高原活动断裂同震、震后应力影响

2015年尼泊尔MW7.9地震重烈度区从震中向东延伸,致灾范围包括尼泊尔、印度北部、巴基斯坦、孟加拉和中国藏南地区,其应力调整对邻区和周边活动断裂可能产生重要影响.本文基于地震应力触发理论,采用岩石圈地壳分层黏弹性位错模型,计算了尼泊尔MW7.9地震引起的周边断裂,特别是青藏高原活动断裂的同震和震后库仑应力变化.结果显示,尼泊尔地震同震效应引起大部分震区库仑应力升高,余震主要分布在最大同震滑动等值线外部库仑应力升高区域;少量余震靠近最大滑动量区域,可能该区域积累的地震能量在主震期间没有完全释放.尼泊尔地震同震库仑应力对青藏高原,特别是中尼边境区域活动断裂有一定影响.亚东—谷露地堑南段、北喜马拉雅断裂西段、当惹雍错—定日断裂和甲岗—定结断裂同震库仑应力升高,其中当惹雍错—定日断裂南端,北喜马拉雅断裂西段同震库仑应力变化峰值超过0.01 MPa;帕龙错断裂、班公错断裂、改则—洞措断裂库仑应力降低,其地震发生概率有所降低.震后应力影响方面,未来40年内黏弹性松弛作用导致北喜马拉雅断裂、改则—洞措断裂和喀喇昆仑断裂整体应力卸载;藏南一系列正断层震后应力持续上升,其中帕龙错断裂南段受到震后黏弹性库仑应力影响,由应力阴影区逐渐转化为应力增强区,当惹雍错—定日断裂南段应力进一步加强,震后40年其南端应力变化峰值达到0.1345 MPa,亚东—谷露断裂南段应力亦持续增强.藏南正断层的地震活动性值得进一步关注.     

自维持断层蠕滑过程和余震发生率的衰减特征—以汶川MW7.9地震序列为例

基于Dieterich地震活动性理论,本文推导出计算余震发生率和余震累积次数的一般表达式,其中主震后发震断层内部的剪切应力随时间的演化过程遵从Jeffreys-Lomnitz蠕变模型,且与修正Omori定律直接相关。修正Omori定律中的p值与震后断层的短时应力加卸载过程正相关。采用Rubin和Ampuero 给出的震后断层自维持蠕滑模型本文得出计算余震发生率的近似表达式,并对2008年汶川地震序列进行拟合。结果表明,p值的大小直接对应了速率-状态摩擦定律中摩擦参量b/a,而修正Omori定律中的c值则与速率-状态摩擦定律中的临界滑移D_c相关。对于汶川余震序列而言,拟合结果显示b/a约为1.13,D_c约为2—3 cm。Rubin-Ampuero震后自维持蠕滑描述了震后孕震层内部短暂的速率变化特征,是孕震断层演化过程不可缺少的环节。      

岩石圈结构对大地震震后形变的影响——以1976年唐山大地震和2001年昆仑山大地震为例

大地震发生之后通常会诱发一系列的余震序列,对比1976年M_S7.8唐山大地震和2001年M_S8.1昆仑山大地震周边区域的地震事件可以看出,唐山大地震余震活动时间要明显长于昆仑山大地震余震活动时间.余震序列往往与震后形变密切相关,而影响震后形变的因素不仅与地震发震断层和震级有关,同时与岩石圈的结构有关.考虑到唐山大地震的发震区华北地块和昆仑山大地震的发震区青藏高原有着较大的岩石圈结构差异,本文采用PSGRN/PSCMP软件计算了岩石圈分层模型的大地震同震和震后形变,分析了地壳弹性模量、弹性厚度以及黏滞性系数对同震和震后形变的影响,进而讨论了影响唐山地震和昆仑山地震余震序列差异的原因.计算结果显示,震后形变会在黏弹性效应的作用下逐渐调整,震后形变的持续时间与地壳弹性模量、地壳弹性厚度和下地壳黏滞性系数有关.上地壳和下地壳弹性模量越大,震后形变达到稳定值的时间越短,弹性模量对震后形变稳定值影响很小.地壳弹性厚度越大,震后形变达到稳定值的时间越短,当断层面底端深度小于地壳弹性厚度时,地壳弹性厚度的增加会引起震后形变稳定值的减小;下地壳厚度对震后形变达到稳定值的时间和稳定值基本无影响.下地壳黏滞性系数越大,震后形变达到稳定值的时间越长,反之亦然.结合唐山地震区的华北地块和昆仑山地震的青藏高原深部结构发现,两者之间的上地壳弹性模型差别不大,唐山地震区地壳弹性厚度略大于昆仑山地震区,但昆仑山地震区下地壳黏滞性系数明显低于唐山地震区.这些因素均决定了昆仑山地震的震后形变持续时间短（余震时间序列短）而唐山地震的震后形变持续时间长（余震时间序列长）.由此可见,岩石圈结构差异可能是导致唐山地震和昆仑山地震余震序列差异的主要因素之一.     

昆仑山口西M_S 8.1地震同震影响场的数值模拟

根据野外调查得到的地震破裂以及同震位错分布,利用三维有限元模型研究了昆仑山口西MS8.1地震的同震位移场和同震应力场的分布,讨论了震后地震活动与同震应力场变化的关系。模拟的MS8.1地震同震位移场显示,东昆仑断裂带南北两侧的同震位移衰减程度存在差异,断裂带南部的同震位移衰减大于北部同震位移的衰减;同震应力场研究结果则表明,MS8.1地震的同震最大剪切应力主要集中在东昆仑断裂带的周边地区,最大剪切应力等值线沿断裂带密集分布,形成了高梯度带。MS8.1地震的发生使青藏活动地块区的地壳应力发生了剧烈的变化,可能导致研究区的地壳应力场得到松弛。库伦破裂应力的研究表明,昆仑山地震对其后发生的5次强余震有一定的触发作用,昆仑山口西MS8.1地震可能也触发了2003年的德令哈6.6级地震     

从应力演化角度讨论龙门山区域7级强震序列对芦山地震的影响

计算包括同震静态库仑破裂应力变化以及震后黏滞松弛引起的应力变化等,可以更好地解释余震分布、地震序列等地震观测结果。在芦山地震之前,从1900年以来龙门山区域发生了4次7级以上地震,分析这几次地震的同震应力变化以及震后黏弹性松弛对芦山地震的产生的影响,芦山地震以后区域断裂带上的应力伴随强震如何演化,芦山地震与汶川地震的破裂空段呈现何种应力状态,探讨这些问题可能会为了解芦山地震震源处震前的应力状态及该区域未来地震风险评估提供一定的依据。文中采用有限元数值模拟方法,根据地质构造、速度、密度结构深部反演结果以及GPS及应力观测资料等,建立龙门山地区三维黏弹性有限元模型进行研究。模拟结果显示:龙门山断裂带南段及鲜水河断裂带的库仑应力年变化速率在研究区域中相对更高,这与研究区域的地震活动性一致。芦山地震的前4次地震,除叠溪地震外,康定、松潘、汶川等3次地震在芦山地震震中位置产生的同震库仑破裂应力变化大于0,表明这3次地震可能促进了芦山地震的发生,汶川地震的同震库仑破裂应力超过了0.01MPa,同震触发效应十分显著。震间的黏弹性松弛对芦山震源处起加载作用,从1900年以来这种持续的加载作用也超过0.01MPa,因此在模拟应力演化的时候,介质的黏弹性松弛效应不能被忽略。从库仑破裂应力的角度计算龙门山区域断层的应力演化,可以发现龙门山断裂带上汶川地震和芦山地震破裂的空段,在芦山地震之后仍然属于相对应力水平较高的区域。     

松原5.7级地震震中区土壤氢气变化特征

2018年5月28日吉林松原发生5.7级地震,为了探究地震对断层气的影响以及震后断层气变化特征,震后在震中区进行了多期土壤氢气浓度观测工作,得到了震中区氢气浓度的最大值、背景值等多项指标。结合地震前后氢气浓度数据,分析震中区震前、震后氢气浓度趋势性变化情况,结果显示:在此次地震活动中,氢气浓度变化有"震前缓慢升高—临震下降—震后迅速升高—强余震前再次骤升"的现象;氢气浓度变化与地震活动间有很好的映震关系,主震对氢气浓度的上升起主要作用,而余震则起诱发作用,余震的活动会使赋存在地下岩石裂隙的氢气释放量增大,表明氢气对余震活动的响应较灵敏;另外,震后某一时刻空气中的氢气浓度也会升高,这可能与区域断层中的氢气逸散到大气中有一定关系。     

The 20th September 1999 Chi-Chi Earthquake (Taiwan): a case of study for its aftershock seismic sequence

The purpose of this work is to highlight some methodological aspects related to the observation of possible anomalies in the temporal decay of aftershocks temporal series following a mainshock with magnitude M ≥ 7.0. In this paper we present the results for the Taiwan seismic sequence started on 20 September 1999 (M = 7.7) by tuning some seismic parameters that show considerable variations during the aftershock decay process. In here we also present the results obtained using a fractal approach for the seismic sequence. Earthquakes belong to a class of phenomena known as multifractals. In general it is important to define the fractal dimension D, but sometimes is not useful if we are describing a natural phenomenon; so it is necessary to define D ₀ called box-counting dimension and D ₂ called correlation dimension, usually D ₀ ≥ D ₂. In the elaborations of the fractal dimensions, for this sequence, we have obtained values lesser than 1, with a greater tendency of aftershocks to clusterize in time before a large aftershock. This is coherent with the possible existence of seismic anomalies, that could occur before the large aftershock. We also report the results obtained by using the delta/sigma method described firstly in [Caccamo et al., 2005] and later applied to different seismic sequence. The observed temporal series of the aftershocks per day can be considered as a sum of a deterministic and a stochastic contribution. If the decay can be modeled as a non-stationary Poissonian process, the number of aftershocks in a small time interval Δt is the mean value n(t) Δt, with a standard deviation (δ = √n(t)Δt. Investigating both aftershock behavior and a wide spectrum of parameters may find the key to explain better the mechanism of seismicity as a whole.     

2018年9月4日新疆伽师MS5.5地震序列及发震构造讨论

2018年9月4日新疆伽师发生M_S5.5地震,震中处于塔里木地块西北缘,位于1997~1998年伽师强震群震区内。此次伽师地震前发生了M_S4.7前震,截至9月30日最大余震震级为M_S4.6（M_L5.0）,初步判定为前-主-余型地震序列。序列精定位结果显示,余震沿近NE向展布,主震震源深度与1997~1998年伽师强震主震基本一致,发震断层陡立。本文从区域的构造环境、地震震源机制解和余震分布特征等方面分析认为,地震发生在伽师隐伏断裂东南端部,为1997~1998年伽师强震群震区的一次新的构造活动。序列参数、视应力等计算结果显示,伽师M_S5.5地震的预测最大余震震级与最大余震震级M_S4.6接近,表明序列最大余震已经发生。     

Seismic anomalies in the aftershock sequence of November 16, 2000, in Papua New Guinea

The study of seismic anomalies, related both to the temporal trend of aftershock sequences and to the temporal series of mainshocks, is important for an understanding of the physical processes relating to the existence and the characteristics of seismic precursors. The purpose of this work is to highlight some methodological aspects related to the observation of possible anomalies in the temporal decay of an aftershock sequence. It is realized by means of several parameters. We focused our work on an analysis of the Papua New Guinea seismic sequence that occurred on November 16, 2000. The magnitude of the mainshock is M = 8.2. The observed temporal series of shocks per day can be considered as a sum of a deterministic contribution and a stochastic contribution. If the decay can be modeled as a nonstationary Poisson process where the intensity function is equal to n(t) = K(t + c)^−p + K ₁, the number of aftershocks in a small time interval Δt is the mean value n(t)Δt, with a standard deviation σ = √n(t)Δt. We observe that there are some variations in seismicity that can be considered as seismic anomalies before the occurrence of a large aftershock. The data, checked according to completeness criteria, come from the website of the USGS NEIC data bank (). The text was submitted by the authors in English.     

Monitoring the West Bohemian earthquake swarm in 2008/2009 by a temporary small-aperture seismic array

The most recent intense earthquake swarm in West Bohemia lasted from 6 October 2008 to January 2009. Starting 12 days after the onset, the University of Potsdam monitored the swarm by a temporary small-aperture seismic array at 10 km epicentral distance. The purpose of the installation was a complete monitoring of the swarm including micro-earthquakes (M _L < 0). We identify earthquakes using a conventional short-term average/long-term average trigger combined with sliding-window frequency-wavenumber and polarisation analyses. The resulting earthquake catalogue consists of 14,530 earthquakes between 19 October 2008 and 18 March 2009 with magnitudes in the range of − 1.2 ≤ M _L ≤ 2.7. The small-aperture seismic array substantially lowers the detection threshold to about M _c = − 0.4, when compared to the regional networks operating in West Bohemia (M _c > 0.0). In the course of this work, the main temporal features (frequency–magnitude distribution, propagation of back azimuth and horizontal slowness, occurrence rate of aftershock sequences and interevent-time distribution) of the recent 2008/2009 earthquake swarm are presented and discussed. Temporal changes of the coefficient of variation (based on interevent times) suggest that the swarm earthquake activity of the 2008/2009 swarm terminates by 12 January 2009. During the main phase in our studied swarm period after 19 October, the b value of the Gutenberg–Richter relation decreases from 1.2 to 0.8. This trend is also reflected in the power-law behavior of the seismic moment release. The corresponding total seismic moment release of 1.02×10¹⁷ Nm is equivalent to M _L,max = 5.4.     

基于ETAS模型对青海地区余震预测效能评估的研究

为系统评估青海地区余震短期发生率的预测效能,以及构建适合地震活动特点的余震早期预测策略和预测指标体系,利用国际上当前较为前沿的时间序列ETAS模型和"瘦化算法"对青海地区2009年以来的8个地震的早期余震序列参数进行拟合,并利用N-test检验方法对预测结果进行回溯性的效能评估。研究表明:ETAS模型和"瘦化算法"对青海地区的余震发生率具有很好的预测能力,建议采用3天的预测时间窗,且对序列早期阶段进行应用,或可取得"最佳"的预测效果。     

The time of occurrence and the magnitude of the largest aftershock over India

Summary The time of occurrence and the magnitude of the largest aftershock in relation to the main shock have been studied for India and its neighbourhood based on the USCGS data during the years 1963–1971. It is found that the largest aftershock occurs within 2 hours after the main shock in about 50% of the cases and frequency of occurrencen(t) of the largest aftershocks decreases hyper-bolically with the intervalt after the main event and could be represented by a law of the formn(t)=At ^–h whereA andh are constants. The probability of occurrence of the largest aftershock within 2 hours of the main shock is found to be higher over island are regions of the world. The difference (M ₀–M ₁) of the magnitude of the largest aftershockM ₁ to that of the main shockM ₀ as a measure of aftershock activity does not show any marked regional variation over India and its neighbourhood, as was reported by Mogifor Japan. Examination of the values ofM ₁/M ₀ and the constantb in Gutenberg-Richter's frequency magnitude relationship reveals a range of variation in both; high values ofM ₁/M ₀ have been found to be associated with high values ofb in many tectonic earthquakes and thus not, restricted to reservoir associated seismic activity.     

Development of an aftershock occurrence model calibrated for Turkey and the resulting likelihoods

This paper presents the calibration of Omori's aftershock occurrence rate model for Turkey and the resulting likelihoods. Aftershock occurrence rate models are used for estimating the probability of an aftershock that exceeds a specific magnitude threshold within a time interval after the mainshock. Critical decisions on the post-earthquake safety of structures directly depend on the aftershock hazard estimated using the occurrence model. It is customary to calibrate models in a region-specific manner. These models depend on rate parameters(a, b, c and p) related to the seismicity characteristics of the investigated region. In this study, the available well-recorded aftershock sequences for a set of Mw ≥ 5.9 mainshock events that were observed in Turkey until 2012 are considered to develop the aftershock occurrence model. Mean estimates of the model parameters identified for Turkey are a =-1.90, b = 1.11, c = 0.05 and p = 1.20. Based on the developed model, aftershock likelihoods are computed for a range of different time intervals and mainshock magnitudes. Also, the sensitivity of aftershock probabilities to the model parameters is investigated. Aftershock occurrence probabilities estimated using the model are expected to be useful for post-earthquake safety evaluations in Turkey.     

2013年前郭MS5.8震群矩张量研究

2013年吉林前郭M_S5.8震群为爆发性震群,目前余震活动仍然在持续.基于吉林、辽宁、黑龙江、内蒙古四省地震台网记录的前郭震群波形资料,利用波形信噪比、震源类型、台站及速度模型组合的指标选择最佳的反演方案,应用矩张量的三种反演模式,对序列中5个M_S≥5.0地震进行矩张量反演研究,获得了全矩张量、偏矩张量和纯双力偶的矩张量.使用F-test对地震的三种模式的矩张量反演结果进行显著性检验来确定最佳反演模式.结果显示,5个地震的最优矩张量解均为全矩张量模式反演获得的结果,其双力偶分量仅有20%~65%,矩心深度位于地下3~4 km处,地震在Hudson震源类型图上的投影远离双力偶震源类型区域.这些结果表明,震源类型并非典型的构造地震,推断前郭地震可能是与人类活动有关的诱发地震.     

Relocation of the MS5.7 Earthquake Sequence in Songyuan, Jilin Province, 2018 and the Analysis of Its Seismogenic Structure

The 2018,Songyuan,Jilin M_S5. 7 earthquake occurred at the intersection of the FuyuZhaodong fault and the Second Songhua River fault. The moment magnitude of this earthquake is M_W5. 3,the centroid depth by the waveform fitting is 12 km,and it is a strike-slip type event. In this paper,with the seismic phase data provided by the China Earthquake Network, the double-difference location method is used to relocate the earthquake sequence,finally the relocation results of 60 earthquakes are obtained. The results show that the aftershock zone is about 4. 3km long and 3. 1km wide,which is distributed in the NE direction. The depth distribution of the seismic sequence is 9km-10 km. 1-2 days after the main shock,the aftershocks were scattered throughout the aftershock zone,and the largest aftershock occurred in the northeastern part of the aftershock zone. After 3-8 days,the aftershocks mainly occurred in the southwestern part of the aftershock zone. The profile distribution of the earthquake sequence shows that the fault plane dips to the southeast with the dip angle of about 75°. Combined with the regional tectonic setting,focal mechanism solution and intensity distribution,we conclude that the concealed fault of the Fuyu-Zhaodong fault is the seismogenic fault of the Songyuan M_S5. 7 earthquake. This paper also relocates the earthquake sequence of the previous magnitude 5. 0 earthquake in 2017. Combined with the results of the focal mechanism solution,we believe that the two earthquakes have the same seismogenic structure,and the earthquake sequence generally develops to the southwest. The historical seismic activity since 2009 shows that after the magnitude 5. 0 earthquake in 2017,the frequency and intensity of earthquakes in the earthquake zone are obviously enhanced,and attention should be paid to the development of seismic activity in the southwest direction of the earthquake zone.