云南测震台网台站场地响应与ML震级测定的讨论

利用2016年云南测震台网记录到的地脉动数据,运用H/V谱比法计算分析了47个测震台站在0.1～20.0Hz频段内的场地响应。选取2012～2016年期间云南测震台网记录到的ML≥2.5地震做单台震级偏差统计分析,挑出受场地放大作用而出现单台震级大于台网平均震级的台站做去场地放大校正,然后重新统计单台震级偏差情况。计算结果显示云南测震台网的台站受台基条件和所处位置等地质构造因素影响,均存在一定程度的场地放大作用,大部分台站场地响应的卓越频率为1～6Hz,放大倍数为2～4倍。受场地放大作用的影响,大姚、芒市等19个台站的单台ML震级大于台网平均震级。在去除场地放大作用后,台站单台震级与台网平均震级的偏差有所减小。     

四川地区介质衰减、场地响应与震级测定的讨论

利用四川区域地震台网的数字波形资料以及中国地震台网中心的有关震相数据,选取互相衔接的3段地震波几何衰减模型,利用基于遗传算法的Atkinson和Mereu方法,分别反演四川盆地、川西高原及攀枝花—西昌地区的介质品质因子。采用Brune的ω平方模型约束震源位移谱,使用Moya方法联合多台多地震数据求得各地震台站的场地响应。根据波形数据和观测报告,计算出2009年1月1日至2015年6月30日10535次2.5≤M_L≤4.9近震单台震级测定值与台网平均震级的差,绘制各次地震的单台震级偏差空间分布图像。结果显示:四川盆地介质衰减模型呈现高Q_0和低衰减系数γ的特点,攀西地区表现出低Q_0和较高衰减系数γ;岩石台基对地震波信号在不同频率点的放大作用不相同;单台震级偏差优势分布基本一致,在构造稳定的四川盆地,介质品质因子高(低衰减),地震波在该区域能有效传播,盆地表层松软沉积层对入射波有放大作用,造成该区域单台震级高异常;川西大部分地震台站台基增益较大,但单台震级明显偏低,推测与松潘-甘孜块体地壳深度软流圈的物质流动及该区域低Q_0异常(高衰减)有关;攀西地区断陷拉张、流体侵入及沉积填充等损耗地震波能量,常使该地区台站测定的川西藏东地震震级偏小,但同样台站测定的龙门山断裂带上中强地震震级常较台网平均震级更大,可能是受盆地边缘强烈不对称高幅值带的影响。从研究结果看,四川地区介质衰减对地震射线振幅的影响较地震台站的台基增益大,因此校正四川区域地震台网单台震级时,应结合介质衰减及台站场地响应的空间特性进行考虑。     

南天山西段地震尾波Q值研究

基于新疆地震台网的数字地震波形资料,采用Aki的地震尾波单次散射模型,选取了南天山西段10个台站周围60km范围内M_S≥2.0的397次地震事件计算并得到了不同台站的尾波Qc值。结果显示,南天山西段尾波Q_0值平均值在140~390之间,其中喀什台和乌恰台Q_0值最低;频率相关系数η的平均值在0.62~1.00之间,拟合平均Q值与对应频率之间的关系结果为Q=261.2±110.8f~(0.80±0.28),表明该地区是以低Q_0高η值为特征的构造活跃区。     

τ_c方法在内蒙古地区预警震级计算中的应用

地震预警技术是减轻地震灾害的有效手段之一,而在预警系统中,震级的测定是最重要,也是最困难的。考虑到地震预警震级估算方法的区域性特点,本文利用内蒙古测震台网2016～2018年M_L≥2. 0地震波形数据共120条,分别采用P波3s时间窗内垂直向位移时程和P波全波段垂直向位移时程计算特征周期参数值,得到了地震预警震级的计算公式,并对其结果进行对比分析。结果显示,采用P波全波段关系式得到的结果较采用P波3s固定时间窗计算值收敛性更好,关系式拟合相关系数更接近1,能够对预警地震震级做出较好的估计。同时结合数据处理结果,探讨了地震预警的警报快速性和准确性之间的关系。     

青海祁连MS5.2地震微震检测与目录完备性研究

利用模板匹配方法对2015年11月23日青海省祁连县M_S5.2地震进行遗漏地震检测研究,由于主震后短时间内目录中遗漏事件较多,故对主震后1天的连续波形进行检测。主震后1天内青海测震台网记录到的余震个数(包括单台)共62个,选取主震后M_L1.0以上余震30个作为模板事件,通过匹配滤波的方式扫描出遗漏地震31个,约为台网目录给出的0.5倍。基于包络差峰值振幅与震级的线性关系估测检测事件的震级参数,最后将检测后的余震目录与台网余震目录在主震后1天内的最小完备震级进行对比分析,结果发现检测后最小完备震级从M_L1.2降到了M_L0.7,得到青海测震台网在祁连地区最小完整性震级为M_L0.7。     

四川测震台网ML震级特征研究

基于2013—2022年四川测震台网记录到的四川及邻区2.0≤M_L≤5.5的36 693次地震事件观测资料,利用震级残差统计方法得到59个区域台站的单台M_L震级与台网平均M_L震级的偏差、各单台记录地震的平均偏差和标准偏差。结果显示：震级偏差统计直方图基本呈正态分布,相对集中于-0.5～0.5;单台M_L震级偏差沿龙门山断裂带两侧分区特征明显,呈现“东高西低”,即川西高原地区台站M_L震级偏小及四川盆地台站M_L震级偏大特征;从M_L震级随震中距的变化曲线来看,当震中距<150 km时,单台测定的M_L震级较台网平均震级偏小,当震中距≥200 km时,单台测定的M_L震级较台网平均震级偏大。     

利用S波高频衰减对天山中东段地区地壳Q值成像

利用天山中东段地区25个数字地震台站记录到的2009~2014年底5076个近震波形资料以及新疆地震台网提供的观测报告,采用遗传算法对S波位移谱的高频衰减进行拟合,得到19140条衰减算子t~*数据。根据这些t*数据,反演得到该区地壳Q_s值分布。结果表明,天山中东段地区平均Q_0值为520,其Q_s值分布及其所揭示的衰减变化特征与研究区的地表构造明显相关。天山南北两侧的山区盆地交汇部位Q_s值较低,而高Q_s值较为集中地分布于天山造山带内部。1900年以来天山中东段地区6级以上强震大多位于低Q_s值区域,该区24个高热流点也大多位于上述天山南北两侧的低Q_s值区域,热流值与衰减值呈负相关;此外,研究区的速度结构与衰减结构也呈一定的正相关,反映了二维衰减结构特征与速度结构、二维密度结构的一致性。     

内蒙古中西部地区台站震级偏差和近震量规函数研究

采用2008年至2015年11月内蒙古中西部和邻近区域43个台站所记录的内蒙古中西部地区2327次地震事件,统计各台站单台震级与台网平均震级的偏差、平均偏差和标准偏差。震级偏差随着震中距的变化呈递增趋势,在震中距范围0~150 km内,偏差值小于0;在震中距范围150~500 km内,偏差值大于0。使用台站震级偏差结果对台站记录震级修正后,统计台站记录震级随震中距的变化,对量规函数进行修正。     

2017年九寨沟MS7.0地震序列活动特征分析

中强地震序列的主震发生后,短时间内受台站距震中较远、尾波干扰和波形重叠等因素的影响,往往会遗漏大量的地震,而地震目录的完整性会直接影响到震后趋势判定和余震序列特征分析的科学性和可靠性。本文利用基于GPU加速的模板匹配方法对2017年8月1～12日的连续波形进行扫描计算,检测九寨沟M_S7.0地震前后遗漏的地震事件,选取台网目录中信噪比较高的1033个地震事件作为模板,在主震前7天至震后5天期间识别出4854个检测地震事件,为台网可定位目录的3.3倍,除去对台网单台地震事件的修正外,还检测到1797个遗漏地震事件,将完备震级从1.6级降低到1.4级。基于补充了遗漏地震的完整地震目录,对2017年8月8日九寨沟M_S7.0地震序列活动特征进行分析。结果表明,前震序列在主震前短时间内出现了地震活动的密集增强,b值也显示为低值状态,可能是深部断层发生破裂之前的加速蠕动的结果。随着时间的推移,余震序列的完备震级逐渐下降并趋于稳定,b值存在缓慢升高的趋势,未来较长时期内余震序列仍将处于持续衰减的状态。     

陕西地区近震震级的量规函数研究

选取陕西数字测震台网2002~2014年所记录到的ML≥2.0且记录台站数大于4的2 806次地震事件,采用震级残差统计方法,得到各单台震级相对于台网震级Mi的偏差εi,εi离散度较小,大体呈正态分布。分析了单台震级与台网震级的平均偏差ΔMi、标准偏差εi,只有3个台(HUAX、HZHT、TOCH)的平均偏差较大,这可能是由于场地响应放大或减小作用导致的。单台震级偏差值随震中距的变化表明,当震中距Δ≤140 km及210Δ≤250 km时,全国量规函数稳定性弱,进而修正得到陕西地区的区域量规函数。     

Use of mode-converted waves in marine seismic data to investigate the lithology of the sub-bottom sediments in Isfjorden,Svalbard

The compression wavefield is efficiently converted to shear-wave energy at post-critical angles in areas of high impedance contrast at the sea floor. We have analysed mode-converted shear waves in a data set acquired with a hybrid marine/land geometry in Isfjorden, Svalbard. Through a kinematic 2D ray-tracing modellingV _p/V_s ratios for part of the uppermost 5km of the crust are obtained. Low values (V _p /V _s =1.65) are tentatively associated with the section of Devonian sandstones which appears to attain a minimum thickness of 1.5km below 3 km depth about 10km west of Kapp Thorden.     

Moment Magnitude and Its Calculation

In this paper, we give a brief introduction to the proposal and development history of the earthquake magnitude concept. Moment magnitude M_W is the best physical quantity for measuring earthquakes. Compared with other magnitude scales used traditionally, moment magnitude is not saturated for all earthquakes, regardless of big and small earthquakes, deep and shallow earthquakes, far field and near field seismic data, geodetic and geological data, moment magnitude can be measured, and can be connected with well-known magnitude scales such as surface wave magnitude M_S. Moment magnitude is a uniform magnitude scale, which is suitable for statistics with wide magnitude range. Moment magnitude is the preferred magnitude selected by the International Seismological community, and it is preferred by the departments responsible for publishing seismic information to the public.Moment magnitude is a uniform magnitude scale, which is suitable for statistics with wide magnitude range. Moment magnitude is a preferred magnitude for international seismology, it is preferred by the agency responsible for providing information about earthquakes to the public. We provide all formulas used in the calculation of moment magnitude, and the calculation steps in detail. We also analyzed some problems and rules to solve these problems by using different formulas and numerical value calculation steps.     

The Key Technical Points and Main Characteristics of the New National Standard of Magnitude

The New Magnitude National Standard of General Rules for Earthquake Magnitude( GB17740-2017) is the state mandatory standard. It was released on May 12,2017,by the General Administration of Quality Supervision,Inspection and Quarantine of the Peoples Republic of China and the Standardization Administration of the Peoples Republic of China. This paper introduces the necessity of revising the national standard of magnitude,and the main contents,technical points and primary features of the new national standard of magnitude,so that it can be applied better in practice.     

甘肃近期几次中强地震震源机制解

12002年12月14日玉门MS5.9地震发震时间:21∶27′27″;经纬度:39.8°N,97.3°E;深度:33km;震级:MS5.9。收集了全省及邻近省份共同31个台站初动后作出震源机制解(表1,图1)。表1　玉门MS5.9地震震源参数节面A节面BP轴T轴N轴矛盾比/%走向倾角滑动角走向倾角滑动角方位角仰角方位角仰角方位角仰角147°26°35°24°75°111°277°27°140°55°18°21°0.09722003年10月25日民乐-山丹MS6.1地震发震时间:20∶41′36.4″;经纬度:38.33°N,100.95°E;深度:12km;震级:MS6.1。收集了全省21个台站初动后作出震源机制解(表2,图2)。表2　民乐-山…     

Application of Fisher method to discriminating earthquakes and explosions using criterion m b/ M S

We try to give a quantitative and global discrimination function by studying m _b/M _S data using Fisher method that is a kind of pattern recognition methods. The reliability of the function is also analyzed. The results show that this criterion works well and has a global feature, which can be used as first-level filtering criterions in event identification. The quantitative and linear discrimination function makes it possible to identify events automatically and achieve the goal to react the events quickly and effectively. Contribution No.05FE3018, Institute of Geophysics, China Earthquake Administrtion.     

Vp and Vp/Vs structures in the crust and upper mantle of the Taiwan region, China

A tomographic study of the V _p and V _p/V _s structures in the crust and upper mantle beneath the Taiwan region of China is conducted by simultaneous inversion of P and S arrival times. Compared with the previous tomographic results, the spherical finite difference technique is suitable for the strong heterogeneous velocity structure, and may improve the accuracy in the travel time and three-dimensional ray tracing calculations. The V _p and V _p/V _s structures derived from joint inversion and the relocated earthquakes can provide better constraints for analyzing the lateral heterogeneity and deep tectonic characters in the crust and upper mantle. Our tomographic results reveal significant relations between the seismic wavespeed structure and the tectonic characters. In the shallow depth, sedimentary basins and orogen show distinct wavespeed anomalies, with low V _p, high V _p/V _s in basins and high V _p, low V _p/V _s in orogen. As the suture zone of Eurasian Plate and Philippine Sea Plate, Longitudinal Valley is characterized by a significant high V _p/V _s anomaly extending to the middle-lower crust and upper mantle, which reflects the impact of rock cracking, partial melting, and the presence of fluids. In the northeast Taiwan, the V _p, V _p/V _s anomalies and relocated earthquakes depict the subducting Philippine Sea Plate under the Eurasian Plate. The high V _p of oceanic plate and the low V _p, high V _p/V _s atop the subducted oceanic plate extend to 80 km depth. Along the east-west profiles, the thickness of crust reaches 60 km at the east of Central Range with eastward dipping trend, which reveals the eastward subduction of the thickened and deformed crust of the Eurasian continental plate. Supported by Knowledge Innovation Program of the Chinese Academy of Sciences (Grant No. KZCX3-SW-234-2), National Basic Research Program of China (Grant No. 2007CB411701), National High Technology Research and Development Program of China (Grant No. 2006AA09A101-0201) and National Natural Science Foundation of China (Grant Nos. 40804016, 40704013)     

Am improved source mechanism for the 1935 Timiskaming,Quebec earthquake from regional waveforms

The Timiskaming earthquake, which occurred near the Quebec-Ontario border at the northwest end of the Western Quebec seismic zone in 1935, is one of the five largest instrumentally recorded southeastern Canadian earthquakes. Previous studies of this earthquake concentrated on modeling teismograms recorded at regional distances, a better constrained focal mechanism is obtained. The waveforms indicate thrust faulting on a moderately dipping northwest striking plane at a depth of 10 km. TheM _w of 6.1 determined in this study is in good agreement with previous magnitude estimates (m _b 6.1,M _s 6.0, andm _bLg 6.2–6.3). The focal mechanism is similar to those of many recent small to moderate earthquakes in the region, and the inferred (from theP axis) acting stress of northeast compression is consistent with the overall eastern North American stress field. The Lake Timiskaming Rift Valley in which the earthquake occurred, comprises several northwest striking faults consistent with the strike of the 1935 event. Thus, the 1935 earthquake appears to be a result of faulting on the reactivated Timiskaming graben.     

青海、新疆和西藏地区地震活动平静特征参数Wq值强震预测模式及其应用

以青海、新疆、西藏地区为研究区域,系统研究总结了该区域2010—2020年5.0级以上地震前W_q值的时空演化特征,并以2020年新疆于田6.4级和西藏尼玛6.6级地震为实例进行阐述。得出以下三点认识:①地震一般发生在W_q值异常扩展时段或扩展—减小时段,6.0级以上W_q值异常的震例中,约81%发生于异常出现后的9个月内;②地震一般位于W_q值异常面积扩展区或扩展—减小恢复区附近,近70%的震例发生在W_q值异常区内;③建立了震级与W_q值异常区面积间的正相关统计模型,二者的相关系数(R)为0.85,为预测青海、新疆和西藏地区地震的强度提供了定量关系。青海、新疆和西藏地区W_q值方法对6.0～6.9级地震的预测效果(W_q值异常的地震报准率为61%)优于5.0～5.9级地震(W_q值异常的地震报准率为26%),为我国地球物理观测程度较低地区开展强震中短期(1年内)预测提供了参考依据。     

Instantaneous global ozone balance including observed nitrogen dioxide

The catalytic destruction of stratospheric ozone by the oxides of nitrogen is believed to be an important part of the global ozone balance. The lack of sufficient measurements of NO _x concentrations has impeded efforts to quantify this process. Recent measurements of stratospheric nitrogen dioxide from ground-based stations as well as aircraft and balloons have provided a first approximation to a global distribution of NO₂ vertical columns at sunset. These observed vertical columns have been translated into time-dependent vertical NO₂ profiles by means of a one-dimensional atmospheric photochemical model. Using recent observations of air temperature and ozone along with this information, the independent instantaneous (one second) rates of ozone production from oxygen photolysis P(O₃), of ozone destruction from pure oxygen species (Chapman reactions) L(O _x ), and of ozone destruction by nitrogen oxides L(NO _x ) were estimated over the three-dimensional atmosphere. These quantities are displayed as zonal average contour maps, summed over various latitude zones, summed over various altitude bands, and integrated globally between 15 and 45 km. Although the global summation between 15 and 45 km by no means tells the complete story, these numbers are of some interest, and the relative values are: P(O₃), 100; L(O _x ), 15; L(NO _x ), 45±15. It is to be emphasized that this relative NO _x contribution to the integrated ozone balance is not a measure of the sensitivity of ozone to possible perturbations of stratospheric NO _x ; recent model results must be examined for current estimates of this sensitivity.