1989年大同—阳高6.1级震群破裂方向讨论

本文利用余震震中分布、震源机制解、地震记录波形、等震线分布形态、地裂缝及地质构造资料,对大同一阳高震群序列中5.7、6.1、5.6级三次地震的破裂方向进行探讨。认为5.7级地震破裂方向为近东西向,6.1、5.6级地震破裂方向为北北东向,由二组破裂形成共轭剪切破裂,大同一阳高震群序列就发生在共轭剪切破裂面上。 1989年10月18日—19日在山西大同—阳高—带相继发生了5.7、6.1、5.6级地震,地震发生后笔者参加了地震现场分析预报组的工作,工作中发现5.7、6.1、5.6级地震破裂方向并不一致,有必要进行讨论。     

2021年5月21日云南漾濞Ms6.4地震:一次破裂在隐伏断层上的浅源走滑事件

2021年5月21日,在中国云南省大理州漾濞县发生了 Ms6.4地震.该地震震中位于川滇菱形块体西部边缘的维西—乔后—巍山断裂附近.作为该区域过去45年以来发生的最大地震,快速的灾害反应及震源参数反演将为防震减灾提供重要的帮助.由于未搜集到与此次地震地表破裂相关的报告,这次地震的发震断层目前仍难以确定.本文利用哨兵1A/B卫星的合成孔径雷达干涉(InSAR)数据获得了与此次地震相关的同震形变,并反演了这次事件的断层几何参数及详细的滑动分布.反演结果显示:破裂断层以右旋走滑运动为主,兼有少量的正断分量;滑动主要分布在2到12 km深度范围内,在约7 km深度处滑动达到最大,约为0.64 m.破裂断层浅部0到2 km范围内滑动的缺失揭示了一个重要的滑动亏损区,该区域很可能被震后余滑及震间浅部滑动所弥补.通过对反演结果、余震序列和局部断层构造特性的综合分析,认为此次地震可能破裂在维西—乔后—巍山断裂的隐伏分支断层或一个独立的未知的隐伏主断层上.计算的同震库仑应力变化表明维西—乔后—巍山断裂的巍山盆地段北端、维西—乔后—巍山断裂的玉狮场—乔后段最南端及洱源—弥渡断裂的中段北端具有较高的破裂风险.     

2014年云南鲁甸M_W6.1地震:一次共轭破裂地震

烈度与余震分布显示2014年云南鲁甸MW6.1(MS6.5)地震的发震构造较复杂.为深入了解鲁甸地震的发震断层与破裂特征,本文考虑了单一断层破裂和共轭断层破裂的情况,对震中距250km范围内的近震资料(宽频带资料和强震资料)和远震体波资料进行了反演,得到了鲁甸地震的破裂过程,探讨了滑动分布与余震分布之间的关系.根据反演得到的滑动分布、震源时间函数和波形拟合,认为鲁甸地震是一次在北西向主压应力与北东向主张应力的统一应力场下发生的两条共轭断层先后破裂的一次复杂地震事件.在破裂开始后0~2s,破裂主要发生在ENE—WSW向(近东西向)的断层上,随后NNW—SSE向(近南北向)断层开始破裂,释放了大部分的地震矩.由于近南北向断层南段(即震中以南)的破裂规模较大,且以左旋走滑为主,对近东西向断层的西段起到了一定程度的解锁作用,可能是震中以西无明显主震破裂但存在密集余震分布的主要原因.     

地震供能、储能与放能机制的研究

“双向应变结构”是孕育和发生地震可能的供能、储能与放能的一种机制。它既能缓慢地储存弹性势能与重力势能而孕育地震,又能突然快速地释放弹性势能与重力势能而发生地震。由于地震发生时孕震体是(原应变的)应变逆转的自破裂过程,因此当初始破裂发生后,在破裂的前沿区,地震剪切应力在原有应力集中的基础上,还会再次产生一个突然快速增强(使破裂加速)到快速减弱(使破裂减速到停止)的过程,即产生“应变逆转破裂效应”。这样,不仅会大大增加地震破裂前沿区的剪切应力,而且还会大大减小破裂中的阻力,可能使地震发生中的断层强度佯谬得到解释。“双向应变结构”破裂发生地震的具体物理机制,印证了震源体断层运动的复杂性：既产生了剪切错动的破裂运动,又产生了单向的拉张与单向的压缩运动,同时还产生了整体的膨胀与整体的压缩运动。     

2017年精河MS6.6地震的断层参数和破裂过程

利用远震资料、近场强震资料和合成孔径雷达干涉同震形变资料确定了2017年8月9日精河M_S6.6地震的断层面参数及震源破裂细节。为得到可靠的断层几何参数,发展了一套基于InSAR数据滑动分布反演的三维格点搜索流程,对本次地震断层面的走向、倾角和震源深度进行了格点搜索。结果显示,地震断层面走向为95°,倾角为47°,震源深度为14 km。基于搜索得到的断层模型进行破裂过程联合反演的结果显示:精河M_S6.6地震为一次单侧破裂事件,最大滑动量约为0.8 m,滑动区域集中在断层面上震源以西5—15 km,沿倾向15—25 km,破裂主要发生在10 km深度以下区域。断层面上的平均滑动角为106°。整个破裂过程释放的标量地震矩为3.6×1018 N·m,对应矩震级为M_W6.3。破裂过程持续约9 s,期间的破裂速度约为2.1—2.6 km/s。由于地震破裂主要集中在10 km以下,未来可能需要关注该区域0—10 km发生潜在地震的可能性。      

1978和2005年宫城近海地震的动力学震源模拟:破裂能解释

对1978和2005年发生在同一板块边界而震级不同的日本宫城近海地震建立了自发动力破裂模型,再现了由运动学震源模型所预测的滑动和破裂速度。对这两次地震的动力学震源参数进行比较,可使我们看到同一断层上不同震源过程的破裂动力学特征。特别是,我们比较了应力降和破裂能,在这两次地震中,二者可以进行稳定的评估。在整个断层面上和两次地震破裂的凹凸体上,两个模型中应力降的最大值和平均值几乎相同。另外,1978年地震的震级比2005年的大,它的破裂能值不仅是在整个断层面上,而且在多次破裂的凹凸体上也比2005年的大。这些结果与先前用地震矩所做的破裂能标度的研究结果一致。我们还将破裂能值与以前研究中所估算的其他内陆地震的能值进行了对比,我们得到的宫城近海地震的值较小。造成这种不一致的原因可能是内陆浅源地震和俯冲带深源地震在构造背景上的差异。     

2000年1月15日姚安6.5级地震较强余震的应力触发

20 0 0年 1月 15日在云南省姚安县发生了MS6 5级地震 ,震后发生了 7次M≥ 4 0级的较强余震。本文计算了姚安MS6 5级主震后 ,主震和前震在 7次强余震破裂面上诱发的库仑破裂应力变化 (ΔCFS)。结果表明 ,有5次较强余震发生在库仑破裂应力增加 (ΔCFS >0 )的地区 ,增加的范围为 10 -2 ～ 10 -1MPa。并对结果的可靠性进行了检验 ,发现在前震和主震的破裂面长宽度和位错量估计值分别取上下限的 16种组合情况下 ,5次较强余震破裂面上的库仑破裂应力变化都为正 ,只是数值在 10 -3 ～ 10 -1MPa之间变化 ,而另 2次较强余震破裂面上的库仑破裂应力变化正负值发生变化。结果表明 ,主震和前震位错产生的库仑破裂应力变化是姚安MS6 5级地震较强余震活动的重要原因。     

含坚固体试样破裂孕育过程中声发射的强度特征及其地震学意义

用类似于天然地震中尾波持续时间确定震级的方法,通过岩石力学实验研究了含坚固体试样破裂孕育过程中声发射的强度特征。结果表明:在加压初期,声发射强度相对平稳,当应力达78 % 的破裂应力时,声发射强度显著增加,当应力达92% 的破裂应力时,声发射活动强度达到极大值,比加压初期大20 % 以上;随着应力的进一步增加,声发射强度减小,达极小值,出现一个相对平静的阶段;应力再增加,声发射活动强度回升,紧接着,主破裂发生。从声发射强度的空间分布上看,从加压初期的平稳到极大值的出现,声发射都发生在坚固体的外围,随着活动强度减弱,并出现低值之后,强度回升,这个过程基本上发生于坚固体及其上方,随后主破裂发生,坚固体也局部破裂。这和澜沧耿马地震前的地震时间进程具有现象上的相似性     

2010年4月14日玉树Ms7.1地震加速度场预测

基于有限断层震源、且使用动力学拐角频率的地震动随机模拟方法预测玉树地震近断层的加速度场.首先,基于有限断层震源建模方法建立该次地震的震源模型;然后,基于上述地震动模拟方法预测玉树地震近断层191个节点的加速度时程.在此基础上,取每个结点的加速度峰值绘制该次地震的近断层加速度场.结果表明:(1)近断层加速度场主要受震源破裂过程和断层面上滑动分布的影响.断层面上凹凸体投影到地表的区域附近,加速度峰值最大,也是震害最严重的区域;(2)对于走滑地震,断层沿线附近的场地并非均会发生破裂方向性效应;发生破裂方向性效应的场地与凹凸体在断层面上的位置有关.     

1970年通海7.7级大地震强余震触发

1970年1月5日云南通海发生了M_S7.7地震, 震后发生了多次M_S＞5.0的强余震。 文中计算了1970年通海7.7级大地震后, 主震分别在5次强余震破裂面上诱发的库仑破裂应力变化(ΔCFS)。 结果表明, 有4次强余震发生在库仑破裂应力增加(ΔCFS＞0)的地区, 增加的范围为10^-2～10^-1 MPa; 有1次强余震按2种震源机制解结果给出的破裂面计算, 得到2种结果, 分别发生在库仑破裂应力变化为正和在库仑破裂应力变化为负的地区。 研究结果表明, 主震位错产生的库仑破裂应力变化可能是1970年通海7.7级大地震强余震活动的重要原因。     

系统聚类法在华北地区的地震分期和预测中的应用

本文用修改了的系统聚类法研究华北地区地震活动期,结果表明,当类间距离用离差平方和增量距离,且样品间用欧氏距离时,能较好地划分该地区的地震活动期。而样品间用相似系数时,不论类间距离是哪一种,所得结果均较差。本文还作了地震活动期的后验预测尝试,还通过逐次筛选方法对21个特征进行了选择。     

欧洲中小震基岩水平向地震动衰减关系研究

大震与中小震之间的地震动衰减规律有所不同.本文使用了132条欧洲基岩水平向记录研究中小震地震动衰减规律.震级范围在地方震级4～6级,距离在震源距70km以内.本文得到了峰值加速度和5%阻尼比、周期0.04～4.00秒之间的加速度反应谱的衰减关系.通过与欧洲和美国西部地震动衰减关系对比,分析了在地震活动性不同的地区内,大震与中小震之间存在的差异有何不同.     

两种质量矩阵在梁模态分析中差异的比较

一致质量矩阵与集中质量矩阵是结构动力分析中常用的质量矩阵形式。为分析它们对结构频率的影响,本文将梁单元的一致质量矩阵与集中质量矩阵的差别分解为4个部分,通过实例讨论了各个部分所引起的频率差异在这两种质量矩阵形式的频率总差异中所占的比重,指出弯曲转动惯性的影响及单元自由度之间的相互耦合是这两种质量矩阵在梁模态分析中表现出来的总差异中的主要原因。     

Linking ground motion measurements and macroseismic observations in France: a case study based on accelerometric and macroseismic databases

Comparison between accelerometric and macroseismic observations is made for three M _w?=?4.5 earthquakes, which occurred in north-eastern France and south-western Germany in 2003 and 2004. Scalar and spectral instrumental parameters are processed from the accelerometric data recorded by nine accelerometric stations located between 29 and 180 km from the epicentres. Macroseismic data are based on French Internet reports. In addition to the single questionnaire intensity, analysis of the internal correlation between the encoded answers highlights four predominant fields of questions bearing different physical meanings: (1) “vibratory motions of small objects”, (2) “displacement and fall of objects”, (3) “acoustic noise” and (4) “personal feelings”. Best correlations between macroseismic and instrumental observations are obtained when the macroseismic parameters are averaged over 10-km-radius circles around each station. Macroseismic intensities predicted by published peak ground velocity (PGV)–intensity relationships agree with our observed intensities, contrary to those based on peak ground acceleration (PGA). Correlation between the macroseismic and instrumental data for intensities between II and V (EMS-98) is better for PGV than for PGA. Correlation with the response spectra exhibits clear frequency dependence for all macroseismic parameters. Horizontal and vertical components are significantly correlated with the macroseismic parameters between 1 and 10 Hz, a range corresponding to both natural frequencies of most buildings and high energy content in the seismic ground motion. Between 10 and 25 Hz, a clear lack of correlation between macroseismic and instrumental observations exists. It could be due to a combination of the decrease in the energy signal above 10 Hz, a high level of anthropogenic noise and an increase in variability in soil conditions. Above 25 Hz, the correlation coefficients between the acceleration response spectra and the macroseismic parameters are close to the PGA correlation level.     

Faulting in zones of lithospheric extension: Quantitative analysis of natural and experimental data

Quantitative relationships between major fault parameters from geological observations and laboratory experiments are compared. Relationships are established between fault length, number, depth of fault penetration, amplitude of displacement and other characteristics. The width of destruction zones is estimated. Spacing between parallel faults of compatible length is evaluated. It is shown that there is a stable correlation between fault length and number, which is independent from the mode of loading the material under destruction. Destruction of the lithosphere is believed to occur according to the laws of deformation of Maxwell elasto-viscous body.     

Spatiotemporal variability of drought and the potential climatological driving factors in the Liao River basin

Spatiotemporal characteristics of drought based on the Standardized Precipitation Index (SPI) in the Liao River basin (LRB) are investigated in this study. High autocorrelation in SPI seems to lend itself to drought prediction. Drought is becoming more frequent, widespread, and severe in the LRB during the past several decades. Major factors affecting drought in this basin are analysed by investigating relationship between SPI and several circulations including western Pacific Subtropical High (WPSH), East Asian Summer Monsoon (EASM) and El Niño‐Southern Oscillation (ENSO) indices. Different correlation patterns between WPSH indices and SPI are obtained. Several significant positive correlations between the area, intensity of WPSH and SPI are observed in the west and the centre of the study area, while negative correlations are observed in the east. Reverse patterns are observed in the correlation between the ridge of westward longitude of WPSH and SPI. Corresponding lag‐correlation is dominated by positive correlations between the area, intensity of WPSH and SPI, and by negative correlation between the ridge of westward longitude of WPSH and SPI. EASM is mainly negative related with drought in the east of the LRB. Significant positive correlation between ENSO and SPI is mainly located in the east while negative correlation is located in west of the basin. Lag‐correlation (with lags of 1 to 12 months) between them is also investigated and results show that significant negative correlation is located in a broad area extending from the west to the centre of the basin, while less positive correlation is observed with the increase of lags. The possibility of employing general circulation models (GCMs) for drought prediction is discussed based on the above analyses. Copyright © 2011 John Wiley & Sons, Ltd.     

Suspended and Dissolved Forms of Nutrients: Their Relationships and Interaction in the Major Tributaries of the Caspian Sea

Normal annual data are used to evaluate the within-year variations in the concentrations of biogenic substances (organic and mineral components of N and P, dissolved organic carbon, inorganic Si) in the major tributaries of the Caspian Sea. Variations in the concentrations of biogenic substances are analyzed, the relationships between their mass concentrations in individual months are evaluated, and correlation between the organic and mineral components of N and P and their aggregated fractions are established. The obtained high-significance correlation between the concentrations of N and P components suggests that the rates of biochemical processes of biogenic substance transformation taking place in the water of tributaries are commensurable.     

Multifractal property and long-range cross-correlation behavior of particulate matters at urban traffic intersection in Shanghai

At urban road intersection, the levels of particulate matters within different size groups present multi-variable relationships and have been attracting increasing attentions. In this study, we attend to apply the recently developed multifractal detrended cross-correlation analysis method to investigate the fractal property and the cross-correlation behavior among the non-stationary time series of particulate matters. Six groups of particulate matters with different sizes are measured at a typical traffic intersection of street canyon under different seasons and weather conditions. Based on the collected database, the statistical analyses are carried out and the results indicate close relationships among these groups. Then the multifractal detrended cross-correlation analyses are performed to explore the interaction patterns among all groups, i.e., the fine particulate matters group of 0.3–0.49 cm with other groups, the coarse particulate matters group of 10 μm above with other groups, and the groups of particulate matters between the sizes of 0.5–9.99 μm respectively. In terms of the results, the multifractal property and long-range cross-correlation behavior are observed among all pairs. Comparing to coarse particulate matters, the distinct multifractal spectrum between fine particulate matters with other size groups are observed, which imply that the relevant cross-correlation behaviors are stronger than that in coarse group. It is also found that the cross-correlation behaviors between fine particulate matters with other size groups are highly dependent on the weather conditions while the cross-correlation behaviors between coarse particulate matters with others tend to more depend on the season variations. Finally, the long-range cross-correlation behaviors between them are also confirmed with randomly shuffled series of the observed particulate matters.     

TOMOGRAPHIC INVERSION OF NORMALIZED DATA: DOUBLE-TRACE TOMOGRAPHY ALGORITHMS1

Tomography is widely used in geophysics as a technique for imaging geological structures by means of data that are line integrals of physical characteristics. In some transmission measurements, due to various kinds of normalization, the measured data are related to two (the current and the reference) raypaths and can be expressed as a function of differences between line integrals. This is the case, for example, in seismo-acoustic emission measurements, when (since the exact start time is unknown) only the differences between traveltimes (differences between line integrals of the slowness) can be determined. Similarly the use of normalized Fourier amplitudes results in data dependent upon the difference between line integrals of the absorption coefficient (computed along the actual and the reference raypaths). In order to invert these data the ordinary tomography algorithms should be modified. Some generalizations are presented for series expansion tomography methods in order to make them applicable to reconstruction problems in which the input data are differences between two line integrals. The conjugate gradient and the simultaneous iterative reconstruction technique (SIRT) methods were adapted and tested. It is shown that the modified tomography algorithms are stable and sufficiently accurate for practical use. In the reconstruction of noise-free difference data, the conjugate gradient algorithm is found to be faster and more accurate while, in the case of noisy difference data, the modified SIRT algorithm is more stable and insensitive to noise.     

Relationships between moment magnitude and fault parameters: theoretical and semi-empirical relationships

Fault parameters are important in earthquake hazard analysis. In this paper, theoretical relationships between moment magnitude and fault parameters including subsurface rupture length, downdip rupture width, rupture area, and average slip over the fault surface are deduced based on seismological theory. These theoretical relationships are further simplified by applying similarity conditions and an unique form is established. Then, combining thc simplified theoretical relationships between moment magnitude and fault parameters with seismic source data selected in this study, a practical semi-empirical relationship is established. The seismic source data selected is also to used to derive empirical relationships between moment magnitude and fault parameters by the ordinary least square regression method. Comparisons between semi-empirical relationships and empirical relationships show that the former depict distribution trends of data better than the latter. It is also observed that downdip rupture widths of strike slip faults are saturated when moment magnitude is more than 7.0, but downdip rupture widths of dip slip faults are not saturated in the moment magnitude ranges of this study.