Three-dimensional P-wave velocity structure derived from local earthquakes at the Katmai group of volcanoes,Alaska

The three-dimensional P-wave velocity structure beneath the Katmai group of volcanoes is determined by inversion of more than 10,000 rays from over 1000 earthquakes recorded on a local 18 station short-period network between September 1996 and May 2001. The inversion is well constrained from sea level to about 6 km below sea level and encompasses all of the Katmai volcanoes; Martin, Mageik, Trident, Griggs, Novarupta, Snowy, and Katmai caldera. The inversion reduced the average RMS travel-time error from 0.22 s for locations from the standard one-dimensional model to 0.13 s for the best three-dimensional model. The final model, from the 6th inversion step, reveals a prominent low velocity zone (3.6–5.0 km/s) centered at Katmai Pass and extending from Mageik to Trident volcanoes. The anomaly has values about 20–25% slower than velocities outboard of the region (5.0–6.5 km/s). Moderately low velocities (4.5–6.0 km/s) are observed along the volcanic axis between Martin and Katmai Caldera. Griggs volcano, located about 10 km behind (northwest of) the volcanic axis, has unremarkable velocities (5.0–5.7 km/s) compared to non-volcanic regions. The highest velocities are observed between Snowy and Griggs volcanoes (5.5–6.5 km/s). Relocated hypocenters for the best 3-D model are shifted significantly relative to the standard model with clusters of seismicity at Martin volcano shifting systematically deeper by about 1 km to depths of 0 to 4 km below sea level. Hypocenters for the Katmai Caldera are more tightly clustered, relocating beneath the 1912 scarp walls. The relocated hypocenters allow us to compare spatial frequency-size distributions (b-values) using one-dimensional and three-dimensional models. We find that the distribution of b is significantly changed for Martin volcano, which was characterized by variable values (0.8 < b < 2.0) with standard locations and more uniform values (0.8 < b < 1.2) after relocation. Other seismic clusters at Mageik (1.2 < b < 2.2), Trident (0.5 < b < 1.5) and Katmai Caldera (0.8 < b < 1.8) had stable b-values indicating the robustness of the observations. The strong high b-value region at Mageik volcano is mainly associated with an earthquake swarm in October, 1996 that possibly indicates a shallow intrusion or influx of gas. The new velocity and spatial b-value results, in conjunction with prior gravity (Bouguer anomalies up to − 40 mgal) and interferometry (several cm uplift) data, provide strong evidence in favor of partially molten rock at shallow depths beneath the Mageik–Katmai–Novarupta region. Moderately low velocities beneath Martin and Katmai suggest that old, mostly solidified intrusions exist beneath these volcanoes. Higher relative velocities beneath the Griggs and Snowy vents suggest that no magma is resident in the shallow crust beneath these volcanoes.     

3-D isotropic and anisotropic tomography of P-wave travel times from the Anhui airgun experiment in the Yangtze River

The Middle-Lower Yangtze River is a typical transition region between the nearly NW-oriented Tethys and NE-trending Pacific tectonic regimes.Structures of different periods and directions overlap strongly during these processes.The NE-trending Yangtze River compound structural belt and NW-trending Tongling-Hangzhou struc-tural belt both control the magmatic activities and distribu-tions of the metallogenic belts in the area.Here,we obtain 3-D high-resolution isotropic and azimuthally anisotropic velocity structures at depths of 1-10 km using the first arrivals from airgun sources.The velocity maps correspond well with the tectonic structures,with high-velocity anomalies distributed in ore-concentrated districts and low-velocity anomalies distributed along the Yangtze River.The fast directions are generally consistent with the fault strike,indicating that the azimuthal anisotropy is mainly dominated by the fault and fracture trends in the upper crust.The complicated fast directions near the Luzong and Tongling ore deposits reveal complex deformations in the upper crust,which are mainly caused by the intersection of the Yangtze River compound and Tongling-Hangzhou structural belts.The magma intrusion beneath the two ore deposits(Luzong and Tongling)are connected at depths of 5-10 km.     

三维复杂速度模型的交切法地震定位

地震定位是地震监测与减灾研究重要基础.基于均匀或横向均匀介质模型,利用震源轨迹确定震源位置的交切法具有稳健和效率高的优点,但定位精度较低,特别是震源深度.为提高震源定位精度,我们提出适用于三维复杂速度模型的地震定位交切法.将地壳速度模型由均匀或横向均匀介质模型扩展为三维复杂速度模型;均匀或横向均匀介质模型对应的原假设为球面或双曲面的震源轨迹通过最小走时树射线追踪技术予以确定.确定震源位置的震源轨迹以到时差作为约束条件;将震源定位于震源轨迹交汇最密集的点处,即总的到时差残差(RDT)最小的点处.定位结果的不确定性可通过RDT值较小节点的空间分布予以定性表示.考察了准确速度模型、扰动速度模型、扰动观测到时及地震在台网外等4种情况下改进方法的地震定位效果,结果表明改进的交切法可用于三维复杂速度模型的地震定位;综合利用P波与S波的到时差信息,可明显改善震源位置约束;使用多条震源轨迹进行定位,有助于减少由随机因素导致的定位误差.     

巧家地震和鲁甸地震震源区周边P波速度结构

本文利用2010年1月至2020年6月巧家地震和鲁甸地震震源区周围发生地震事件的走时观测资料,应用双差层析成像方法获得了2014年MS6.5鲁甸地震和2020年5月18日巧家MS5.0地震周边区域中上地壳的P波速度结构.成像结果显示:整个研究区域的速度结构存在很强的非均匀性:在鲁甸地震震源区附近,浅部存在速度高达6.4...     

Ground motion estimation in Delhi from postulated regional and local earthquakes

Ground motions are estimated at 55 sites in Delhi, the capital of India from four postulated earthquakes (three regional M _w?=?7.5, 8.0, and 8.5 and one local). The procedure consists of (1) synthesis of ground motion at a hard reference site (NDI) and (2) estimation of ground motion at other sites in the city via known transfer functions and application of the random vibration theory. This work provides a more extensive coverage than earlier studies (e.g., Singh et al., Bull Seism Soc Am 92:555–569, 2002; Bansal et al., J Seismol 13:89–105, 2009). The Indian code response spectra corresponding to Delhi (zone IV) are found to be conservative at hard soil sites for all postulated earthquakes but found to be deficient for M _w?=?8.0 and 8.5 earthquakes at soft soil sites. Spectral acceleration maps at four different natural periods are strongly influenced by the shallow geological and soil conditions. Three pockets of high acceleration values are seen. These pockets seem to coincide with the contacts of (a) Aravalli quartzite and recent Yamuna alluvium (towards the East), (b) Aravalli quartzite and older quaternary alluvium (towards the South), and (c) older quaternary alluvium and recent Yamuna alluvium (towards the North).     

利用区域地震体波走时评价华北地区三维地壳速度结构模型

本研究利用国家地震台网131个地震台站2009—2016年记录的1 749次近震的初至P波和S波走时数据,与采用快速行进方法正演计算得到的华北地区4个三维地壳速度模型对应的走时数据进行对比,通过统计分析的方法,评价这4个速度模型与真实地下结构的近似程度。结果表明:4个速度模型在大范围内存在较高的一致性,在整个研究区内（111°E—119.5°E,37°N—42°N）,Shen等的模型（简称“S模型”）相对优于Fang等的模型（简称“F模型”）和Duan等的模型（简称“D模型”）,Laske等的Crust1.0模型（简称“C模型”）相对较差。我们认为该结果与上述几个模型所使用的数据及其分辨率有关。对于研究区域内的构造单元,D模型在燕山褶皱带西南部、太行山山前构造带西北部和沧县隆起区表现较好,F模型在太行山隆起区中部、沧县隆起北部、黄骅凹陷区和燕山褶皱带表现较好,S模型在西部地块、山西凹陷区、太行山山前构造带和冀中凹陷区表现较好,C模型无明显连片表现较好区域。      

One-dimensional velocity model of the Middle Kura Depresion from local earthquakes data of Azerbaijan

We present the method for determining the velocity model of the Earth’s crust and the parameters of earthquakes in the Middle Kura Depression from the data of network telemetry in Azerbaijan. Application of this method allowed us to recalculate the main parameters of the hypocenters of the earthquake, to compute the corrections to the arrival times of P and S waves at the observation station, and to significantly improve the accuracy in determining the coordinates of the earthquakes. The model was constructed using the VELEST program, which calculates one-dimensional minimal velocity models from the travel times of seismic waves.     

Imaging 3-D crustal P-wave velocity structure of western Yunnan with bulletin data

Western Yunnan is a region with intensive tectonic activity and serious earthquake risk. It is of significant importance to study three dimensional crustal structure of this region to understand the tectonic setting and disaster mechanism. Densification and digitalization of seismic networks in this region provides an opportunity to study the velocity structure with bulletin data. In this study, we collect P-wave data of 10 403 regional earthquakes recorded by 79 seismic stations from January 2008 to December 2010. In addition to first arrivals data (Pg with epicentral distance less than 200 km and Pn), the Pg (or P) data with epicentral distance more than 200 km are also considered as later direct arrivals in the tomographic inversion. We also compare the quantity and the quality of the seismic data before 2010 and after 2010. The test results show that adding the follow-up Pg phase can effectively improve the inversion ability of crustal imaging, and quantity and the data quality are significantly improved since 2010. The tomographic results show that: (1) The Honghe fault zone, which is the major fault systems in this region, may cut through the entire crust, and the velocity contrasts between two sides at lower crust beneath the Honghe fault are estimated at higher than 10%, while the velocity difference below Nujiang fault zone extends only in the upper crust; (2) Most of the earthquakes in the region occurred at the interface of high-velocity media and low-velocity media, i.e., the areas with high velocity gradient, which has been validated in other areas.     

海南岛及邻区地壳三维P波速度结构

利用海南岛及邻区1999——2005年9个地震数字化台站记录到的3 500余条区域初至波到时资料,确定了海南岛及邻区地壳三维P波速度结构. 结果表明,该区域浅部上地壳的波速异常与地表地质构造有一定的相关性,即王五——文教断裂带以北地区波速较低,而以南地区为相对高的波速异常区. 这与琼北地区为凹陷、并伴有多期火山活动和较高的地热背景,而琼中南地区为较稳定的隆起等构造特征相一致. 在其中下地壳,琼东地区相对琼西地区波速较低,可能暗示了地幔热物质上涌等动力学过程. 波速异常分布模式还暗示出北西向断裂带(如铺前——清澜断裂带)较浅,而近东西向的王五——文教断裂带较深,有可能延伸至莫霍面或更深.      

Computation of travel times and rays in a medium of two-dimensional velocity distribution

A ray-tracing computer program is described for a two-dimensional velocity distribution defined by the local wave velocitiesv _{i, j} in points at arbitrary depthv _{i, j} below the surface points with the horizontal coordinatesx _i . The velocity variation is assumed to be linear in the triangles formed by three neighbouring points. Travel times and rays are then calculated after the exact analytic formulae for any position of the source within the model. No assumptions other than of a piecewise linear velocity structure are made. A first-order discontinuity can be approximated by a thin layer with a strong velocity gradient and refracted waves or wide-angle reflections obtained in this way. As an example,P-wave rays were computed for section No. 05 of the Alpine Longitudinal Profile. The model includes a low-velocity channel which is cut off on the eastern side, first-order discontinuities and a sediment basin.Paper presented at the ESC-Workshop Meeting Seismic Waves in Laterally Inhomogeneous Media, Liblice, 1978.     

芦山震区地壳三维P波速度精细结构及地震重定位研究

联合芦山地震序列5285个地震的50711条P波初至绝对到时数据及7294691条高质量的相对到时数据,利用双差地震层析成像方法联合反演了芦山震源区高分辨率的三维P波速度精细结构及5115个地震震源参数.反演结果表明,芦山主震震中为30.28°N,103.98°E,震源深度为16.38 km,主震南西段余震扩展长度约23 km,余震前缘倾角较和缓,主震北东段余震扩展长度约12 km,余震前缘呈铲形,倾角较陡.芦山震源区P波三维速度结构表现出明显的横向不均匀性,近地表处的P波速度异常与地形起伏及地质构造密切相关:宝兴杂岩对应明显的高速异常,此异常由地表延伸到地下15 km深度附近,而中新生代岩石表现为低速异常;大兴附近区域亦显示出小范围的大幅度高速异常,宝兴高速异常与大兴高速异常在10 km深度附近相连,进而增加了芦山震源区的高低速异常对比幅度.在芦山主震的南西、北东两段速度结构存在着较大差异,芦山主震在水平向位于宝兴及大兴高速异常所包围的低速异常的前缘.主震南西段余震主要发生在倾向北西的高低速异常转换带上并靠近低速一侧,其下盘为低速异常,上盘为高速异常.而芦山主震北东段的余震主要分布在宝兴高速体与大兴高速体之间,主发震层向北西倾斜,主发震层上方的宝兴高速异常下边界出现一条南东倾向的反冲地震带,两地震带呈"y"型分布.     

A test of new Earth models against PcP and PmKP travel times

PcP and PmKP travel times are computed for three simple or parametric Earth models, based on free-oscillation and travel-time data B₁, PEM-A and HB₁ and compared with PcP and PmKP travel times from different sources. This comparison is made only for the region above and below the core-mantle boundary and is of interest because of the current search for a standard Earth model. The comparison shows that only model B₁ does not need a correction for its PcP travel times. For the PmKP travel times for the three models, corrections of the form Δt = a + bm were obtained. The models need the following corrections for b: ?1.3 for B₁, 2.8 for HB₁ and 0.6 for PEM-A. The corrections a are shown to be equal to the observed corrections for PcP at large epicentral distances. The inversions of free-oscillation data to obtain Earth models are most successful when body-wave phases that interact with the core are included.     

A numerical study of 1-D nonlinear P-wave propagation in solid

IntroductionBecauseoftheextensivedistributionofruptures,micro-cracksandcrystalfracturesintheearth,therelationshipsbetweenthestressandstrainarenolongerlinear,infact,theyarenonlinear.Inordertoinvestigateandusethenonlinearcharacteristicsofsolidmediumintheearth,weshouldconsidertheinfluenceofnonlinearresponseduringthecourseofseismicmodelingandinversion.Thisisoneoftheimportantstudyfieldsthathavebeenpaidgreatattentionsinthere-centyearsintheworld(Minster,etal,1991;ZHANG,TENG,1993).Thenonlinearchar…     

天山东北部地震的重新定位和一维地壳速度模型的改善

根据布设在乌鲁木齐市活断层探测区内的流动宽频带地震台阵,结合区域地震台网的走时数据,利用3种不同的定位方法对新疆天山东北部地区(E85°30′～ 88°30′,N43°00′～44°40′) 2004年8月至2005年8月发生的599个地震进行了重新定位.通过比较不同方法的结果合理性,确定了适合于当地震源精定位的程序,并改善了一维地壳速度模型.结果表明:联合使用流动地震台阵和区域台网的资料,显著提高了研究区的地震定位能力,精定位后震中分布图像更加集中,展示出了天山东北部地区更为明显的与活动构造相关的条带状地震活动分布图像,除了一些与已知断层相关的地震事件外,还发现一些有待证实的活断层.     

Regional minimum 1-D P-wave velocity model for a new seismicity catalogue with precise and consistent earthquake locations in southern Iran

For faster and more robust ray tracing in 1-D velocity models and also due to the lack of reliable 3-D models, most seismological centers use 1-D models for routine earthquake locations. In this study, as solution to the coupled hypocenter-velocity problem, we compute a regional P-wave velocity model for southern Iran that can be used for routine earthquake location and also a reference initial model for 3-D seismic tomography. The inversion process was based on travel time data from local earthquakes paired reports obtained by merging the catalogues of Iranian Seismic Center (IRSC, 6422 events) and by the Broadband Iranian National Seismic Network (BIN, 4333 events) for southern Iran in the period 2006 through July 2017. After cleaning the data set from large individual reading errors and by identifying event reports from both networks belonging to same earthquake (a process called event pairing), we obtained a data set of 1115 well-locatable events with a total number of 24,606 P-wave observations. This data set was used to calculate a regional minimum 1-D model for southern Iran as result of an extensive model search by trial-and-error process including several dozens of inversions. Significantly different from previous models, we find a smoothly increasing P-velocity by depth with velocities of 5.8 km/s at shallow and velocities of 6.4 km/s at deepest crustal levels. For well-locatable events, location uncertainties are estimated in the order of ±?3 km for epicenter and double this uncertainty for hypocentral depth. The use of the minimum 1-D model with appropriate station delays in routine hypocenter location processing will yield a high-quality seismic catalogue with consistent uncertainty estimates across the region and it will also allow detection of outlier observations. Based on the two catalogues by IRSC and BIN and using the minimum 1-D model and station delays for all stations in the region, we established a new combined earthquake catalogue for southern Iran. While the general distribution of the seismicity corresponds well with that of the two individual catalogues by IRSC and BIN, the new catalogue significantly enhances the correlation of seismicity with the regional fault systems within and between the major crustal blocks that as an assembly build this continental region. Furthermore, the unified seismic catalogue and the minimum 1-D model resulting from this study provide important ingredients for seismic hazard studies.     

2-D P-wave velocity structure of lithosphere in the North China tectonic zone: Constraints from the Yancheng-Baotou deep seismic profile

We obtained the 2-D P-wave velocity structure of the lithosphere in the eastern North China Craton, Shanxi fault subsidence zone, and Yinchuan-Hetao fault subsidence zone by ray tracking technology based on six groups of clearly identified crustal phases and one group of lithospheric interface reflection phases from seismic recording sections of 21 shots along the 1300-km-long Yancheng-Baotou deep seismic wide-angle reflection/refraction profile. The results indicate significant differences between the lithospheric structure east and west of the Taihang Mountains, which is a gravity-gradient zone as well as a zone of abrupt change in lithospheric thickness and a separation zone of different rock components. East of the Taihang Mountains, the Mesozoic and Cenozoic lithospheric structure of the North China Craton has undergone strong reformation and destruction, resulting in the lithosphere thickness decreasing to 70–80 km. The North China Basin has a very thick Cenozoic sedimentary cover and the deepest point of crystalline basement is about 7.0 km, with the crustal thickness decreasing to about 31.0 km. The crystalline basement of the Luxi uplift zone is relatively shallow with a depth of 1.0–2.0 km and crustal thickness of 33.0–35.0 km. The Subei Basin has a thicker Cenozoic sedimentary cover and the bottom of its crystalline basement is at about 5.0–6.0 km with a crustal thickness of 31.0–32.0 km. The Tanlu fault is a deep fracture which cuts the lithosphere with a significant velocity structure difference on either side of the fault. The Tanlu fault plays an important role in the lithospheric destruction in the eastern part of the North China Craton. West of the Taihang Mountains, the crustal thickness increases significantly. The crust thickness beneath the Shanxi fault depression zone is about 46 km, and there is a low-velocity structure with a velocity of less than 6.1 km s?? in the upper part of the middle crust. Combined with other geophysical study results, our data shows that the lithospheric destruction at the Shaanxi-Shanxi fault depression zone and the Yinchuan-Hetao rift surrounding the Ordos block is non-uniform. The lithosphere thickness is about 80–90 km in the Datong-Baotou area, 75–137 km at the Dingxiang-Shenmu region, and about 80–120 km in the Anyang-Yichuan area. The non-uniform lithospheric destruction may be related to the ancient tectonic zone surrounding the Ordos block. This zone experienced multi-period tectonic events in the long-term process of its tectonic evolution and was repeatedly transformed and weakened. The weakening level is related to the interactions with the Ordos block. The continental collision between the Cenozoic India and Eurasia plates and N-E thrusting by the Qinghai Tibet Plateau block is causing further reformation and reduction of the lithosphere.     

2-D P-wave velocity structure in the mideast segment of Zhangjiakou-Bohai tectonic zone: Anxin-Xianghe-Kuancheng DSS profile

Introduction The Capital-circle area is the center of politics, economy and culture of China, and manystrong earthquakes occurred here in history. The Zhangjiakou-Bohai tectonic zone located in thearea is an important active zone. In the period of 500 years between AD 1481~1981, 18 earth-quakes with magnitude more than 6 occurred in the area and 11 of them lay in this active tectoniczone. During 1992~2001, continuous GPS observation found sinistral relative movement betweenYanshan Mo…     

2-D P-wave velocity structure in the mideast segment of Zhangjiakou-Bohai tectonic zone: Anxin-Xianghe-Kuancheng DSS profile

In order to get the 3-D fine velocity structure in the Capital-circle area of China, 6 explosions, ranging from 1800 to 2500 kg, were conducted and recorded by an array of 240 seismographs. A reflection/refraction survey was carried out along the profile extending from Anxin county, Hebei Province northeastward to Yanshan Mountains, crossing the Zhangjiakou-Bohai tectonic zone. The 2-D velocity structure of P wave was imaging along the profile. The results show that abnormality exists in the deep structure of the Zhangjiakou-Bohai tectonic zone: The basement is significantly depressed, the interfaces and Moho are uplifted, and a strong velocity gradient layer is existed above the Moho that may be dislocated by deep fault. The crust of Huabei basin is thin and low velocity body exists in the crust. The Yanshan Mountains’ crust is thick, the layers in the crust are quite clear and the velocity in the layer is homogeneous. Huabei basin differs from Yanshan Mountains in structure. Foundation item: Preliminary study of State Key Basic Research Project (2001CCA02100) and National Natural Science Foundation of China (40234038). Contribution No.RCEG200311, Geophysical Exploration Center, China Earthquake Administration.