The 2023 Turkey earthquake doublet: Earthquake relocation, seismic tomography, and stress field inversion

On February 6,2023,two earthquakes with magnitudes of M_(W) 7.8 and M_(W) 7.5 struck southeastern Turkey,causing significant casualties and economic losses.These seismic events occurred along the East Anatolian Fault Zone,a convergent boundary between the Arabian Plate and the Anatolian Subplate.In this study,we analyze the M_(W) 7.8 and M_(W) 7.5 earthquakes by comparing their aftershock relocations,tomographic images,and stress field inversions.The earthquakes were localized in the upper crust and exhibited steep dip angles.Furthermore,the aftershocks occurred either close to the boundaries of low and high P-wave velocity anomaly zones or within the low P-wave velocity anomaly zones.The East Anatolia Fault,associated with the M_(W) 7.8 earthquake,and the SürgüFault,related to the M_(W) 7.5 earthquake,predominantly experienced shear stress.However,their western sections experienced a combination of strike-slip and tensile stresses in addition to shear stress.The ruptures of the M_(W) 7.8 and M_(W) 7.5 earthquakes appear to have bridged a seismic gap that had seen sparse seismicity over the past 200 years prior to the 2023 Turkey earthquake sequence.

The Subduction Structure Beneath the New Britain Island Arc and the Adjacent Region from Double-Difference Tomography

We applied double-difference tomography to relocate seismic events and determine the lithospheric velocity structure beneath the New Britain Island arc and the South Bismarck Sea Basin,based on the local P wave arrival time dataset collected by the International Seismological Centre.Results of the seismic relocation and velocity inversion show that the subduction of Solomon Sea Plate along the New Britain Trench is spatially different above 150 km,and the subduction angle of the slab on the west side is higher than that on the east side.The relocated earthquakes also show that there are double seismic zones at the depths of about 30–90km beneath the New Britain Island Arc.The velocity structure shows that the dehydration of the subducting slab caused the low-velocity anomalies in mantle wedge above the slab,which are associated with the magmatic activities around the New Guinea-New Britain Island arc.Moreover,it shows that there is another low-velocity anomaly zone beneath the Bismarck mid-oceanic ridge with spatial variation.Beneath the west of the Bismarck mid-oceanic ridge,the low-velocity anomaly is weakly connected to the subducted Solomon Sea slab.Conversely,the low-velocity anomaly beneath the Manus Sea Basin is highly intertwined to the subducting slab and its mantle wedge,indicating that the subduction of the Solomon Sea Plate might be a key deep dynamic factor that drives the spreading of the Manus Sea Basin and the separation of the Bismarck Plate.

New Insights into the Mesozoic Large-scale Intra-plate Magmatism and Mineralization in South China from Seismic Tomography

A great amount of new S-wave data has been combined to image the mantle structure down to lower mantle depth near the southeastern margin of Eurasia and understand the tectonics in South China since the Mesozoic.Our results reveal a large-scale structure beneath the South China Sea,with a prominent,broad low-velocity feature of at least 1500-km wide in and below the mantle transition zone(MTZ)and a pronounced low-velocity feature of nearly 500-km wide in the lowermost mantle.Together these features may represent the head-and-tail of an upwelling mantle plume that has interacted with and deformed by the subducting plates in the region.The Andaman-Sumatra plate is seen subducting steeply without stagnation in the MTZ over the large-scale low-velocity structure,suggesting that the upwelling mantle may have exerted significant influence on the overlying plate tectonics.Our results suggest a long-term extensional regime near the southeastern margin of Eurasia,which may help to understand the formation of the large-scale intra-plate magmatism and mineralization with accumulation of huge amount of rare metal elements in South China during the Mesozoic,and the continental breakup,the southeastward extrusion of Indochina,and the opening of the South China Sea during the Cenozoic.

Double-difference tomography of P- and S-wave velocity structure beneath the western part of Java, Indonesia

West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults.In this study,double-difference tomography was used to obtain the 3D velocity tomograms of P and S waves beneath the western part of Java.To infer the geometry of the structure beneath the study area,precise earthquake hypo・center determination was first performed before tomographic imaging.For this,earthquake waveform data were extracted from the regional Meteorological,Climatological,Geophysical Agency(BMKG)network of Indonesia from South Sumatra to Central Java.The P and S arrival times for about 1,000 events in the period April 2009 to July 2016 were selected,the key features being events of magnitude>3,azimuthal gap<210°and number of phases>8.A nonlinear method using the oct-tree sampling algorithm from the NonLinLoc program was employed to determine the earthquake hypocenters.The hypocenter locations were then relocated using double-difference tomography(tomoDD).A significant reduction of travel-time(root mean square basis)and a better clustering of earthquakes were achieved which correlated well with the geological structure in West Java.Double-difference tomography was found to give a clear velocity structure,especially beneath the volcanic arc area,i.e.,under Mt Anak Krakatau,Mt Salak and the mountains complex in the southern part of West Java.Low velocity anomalies for the P and S waves as well as the vp/vs ratio below the volcanoes indicated possible partial melting of the upper mantle which ascended from the subducted slab beneath the volcanic arc.

Crustal structure beneath the Qilian Orogen Zone from multiscale seismic tomography

The Qilian Orogen Zone(QOZ), located in the north margin of the Tibetan Plateau, is the key area for understanding the deformation and dynamics process of Tibet. Numerous geological and geophysical studies have been carried out on the mechanics of the Tibetan Plateau deformation and uplift; however, the detailed structure and deformation style of the Qilian Orogen Zone have remained uncertain due to poor geophysical data coverage and limited resolution power of inversion algorithms. In this study, we analyze the P-wave velocity structure beneath the Qilian Orogen Zone, obtained by applying multi-scale seismic tomography technique to P-wave arrival time data recorded by regional seismic networks. The seismic tomography algorithm used in this study employs sparsity constraints on the wavelet representation of the velocity model via L1-norm regularization. This algorithm can deal efficiently with uneven-sampled volumes, and can obtain multi-scale images of the velocity model. Our results can be summarized as follows:(1) The crustal velocity structure is strongly inhomogeneous and consistent with the surface geological setting. Significant low-velocity anomalies exist in the crust of northeastern Tibet, and slight high-velocity anomalies exist beneath the Qaidam Basin and Alxa terrane.(2)The Qilian Orogen Zone can be divided into two main parts by the Laji Shan Faults: the northwestern part with a low-velocity feature, and the southeastern part with a high-velocity feature at the upper and middle crust.(3) Our tomographic images suggest that northwestern and southeastern Qilian Orogen Zones have undergone different tectonic processes. In the northwest Qilian Orogen Zone, the deformation and growth of the Northern Tibetan Plateau has extended to the Heli Shan and Beida Shan region by northward overthrusting at the upper crust and thickening in the lower crust. We speculate that in the southeast Qilian Orogen Zone the deformation and growth of the Northern Tibet Plateau were of strike-slip style at the upper crust; in the lower crust, the evidence suggests ductile shear extrusion style and active frontage extension to the Alxa terrane.(4) The multi-scale seismic tomography technique provides multiscale analysis and sparse constraints, which has allowed to us obtain stable, high-resolution results.

Three-dimensional seismic velocity tomography of the upper crust in Tengchong volcanic area, Yunnan Province

Based on data collected by deep seismic sounding carried out in 1999, a three-dimensional P wave velocity structure is determined with tomographic inversion. The tomographic result shows that there is a P wave low velocity zone (LVZ) in the upper crust beneath the Tengchong volcanic area. The LVZ is in the depth of 7~8 km and may be a small magma chamber or a partial melting body. The result also shows that the LVZ is in the northeastern side of the Rehai hydrothermal field, which is located in another LVZ near the surface. The shallow LVZ may represent a well-developed fracture zone. The strong hydrothermal activity in Rehai area can attribute to the existence of fractures between two LVZs. These fractures are the channels for going upwards of the deep hot fluid.

Crustal structure of the eastern Dabie orogenic belt from the seismic tomography and wide-angle reflection studies

The studies of seismic tomography and wide-angle reflection have been carried out to reveal the velocity structUrebeneath the eastern Dabie orogenic belt. The result from the seismic tomography shows the high velocity bodiesmight be positioned to a depth of only about 1 .5 km below sea level within the Dabie ultra-high pressure metamorphic (UHPM) belt; the fan-profile shows the Shuihou-Wuhe fault, the demarcation between the South Dabieand the North Dabie, slopes to the south-west at a dip angle of about 45° in the bottom of upper crust. The wideangle reflection shows the middle crustal boundaries and the complex features from the lower crust.

Location of high seismic activity zones and seismic hazard assessment in Zabrze Bielszowice coal mine using passive tomography

In the paper results of passive tomography calculations have been presented to assess rockburst hazard and locate high seismic activity zones in the vicinity of longwall 306 in Zabrze Bielszowice coal mine. The area of study was 1000 m in X direction by 900 m in Y direction. The zones of high values of P-wave propagation velocity have been found to correlate with the distribution of large seismic tremors.

Seismogenic Structure around the Epicenter of the May 12,2008 Wenchuan Earthquake from Micro-seismic Tomography

一个三维的本地规模的 P 速度模型在下面到 25 km，在主要吃惊震中区域附近的深度用从 5856 余震的地震光线由一个最新部署的暂时的地震网络记录了的 83821 event-to-receiver 被构造。棋盘测试证明我们的 tomographic 模型有 &#8764;2 km 的侧面、垂直的分辨率。高分辨率的 P 速度模型在 seismogenic 层揭示了有趣的结构：(1 ) Guanxian-Anxian 差错， Yingxiu-Beichuan 差错和 Longmen 掸人差错地区的 Wenchuan-Maoxian 差错好由突然地上面的外壳的速度变化描出；(2 ) Pengguan 山岳比它的包围区域有通常更高的速度，并且可以延长在下面到从表面的至少 &#8764;10 km；(3 ) 在 Wenchuan-Maoxian 差错下面的一个突然地侧面的速度变化可以显示 Pengguan 山岳的西方的边界或 Wenchuan-Maoxian 差错是垂直的，并且 Wenchuan 地震的低亚硫酸钠中心可能在蘸 Yingxiu-Beichuan 和 Guanxian-Anxian 差错，和垂直 Wenchuan-Maoxian 差错的 NW 的连词点定位；(4 ) 沿着 Yingxiu-Beichuan 差错的附近被很低的速度和低 seismicity 在浅深度描绘，可能由于孔和破裂的高内容；(5 ) 二块低速度的异例分别地与 &#8764;5 km 与 &#8764;7 km 偏移量在挂的墙和 Guanxian-Anxian 差错的脚墙中被想象垂直部件。

Passive seismic velocity tomography on longwall mining panel based on simultaneous iterative reconstructive technique (SIRT)

Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique (SIRT) inversion is used to deduce the stress redistribution around the longwall mining panel. The mining-induced microseismic events were recorded by mounting an array of receivers on the surface, above the active panel. After processing and filtering the seismic data, the three-dimensional tomography images of the p-wave velocity variations by SIRT passive seismic velocity tomography were provided. To display the velocity changes on coal seam level and subsequently to infer the stress redistribution, these three-dimensional tomograms into the coal seam level were sliced. In addition, the boundary element method (BEM) was used to simulate the stress redistribution. The results show that the inferred stresses from the passive seismic tomograms are conformed to numerical models and theoretical concept of the stress redistribution around the longwall panel. In velocity tomograms, the main zones of the stress redistribution around the panel, including front and side abutment pressures, and gob stress are obvious and also the movement of stress zones along the face advancement is evident. Moreover, the effect of the advance rate of the face on the stress redistribution is demonstrated in tomography images. The research result proves that the SIRT passive seismic velocity tomography has an ultimate potential for monitoring the changes of stress redistribution around the longwall mining panel continuously and subsequently to improve safety of mining operations.

Source side seismic tomography(3STomo):A novelmethod to image the subsurface structure beneath seismically active region

We propose a novel seismic tomography method,Source Side Seismic Tomography(3STomo),which is designed particularly to image the subsurface structure beneath seismically active regions.Unlike the teleseismic tomography,in which the data are relative traveltime residuals between closely spaced stations for each teleseismic event,3STomo uses relative traveltime shifts between earthquakes within the study region for each distant station.Given the relatively evener distribution of global seismic stations,this method has unique advantages for imaging the structure beneath regions that have numerous earthquakes but lack of dense seismic stations,for example,some subduction zones and spreading ridges in the ocean.In addition,3STomo has potentially better vertical resolution at shallow depths than the traditional teleseismic tomography.The effect of the inaccurate source parameters on its resolution can be minimized by using depth phases and the technique of joint source and structure inversion.Numerical experiments and application to Luzon Island,Philippines show that 3STomo can be a valuable tool to investigate the subsurface structure beneath some areas where the traditional method cannot be applied to,or at least it can be used as a complementary component of conventional teleseismic tomography to obtain better back-azimuth coverage and achieve higher resolution at shallow depths in the inversion.

Lithosphere structures dynamics in the central High Atlas(Morocco)by seismic tomography and gravimetric data

We investigate sedimentological and tectonic processes at the central High Atlas, in order to understand their functioning and to build a geodynamic model that placed the area in its current geological frame.The analysis used here is based on the numerical enhancement of a Landsat image where the main goal is to map surface sedimentary deposits throughout the central High Atlas in order to delimit the large geological structures. The sediment distribution throughout central High Atlas indicates that this one is a large tectonic subsiding basin, where the ongoing tectonic events and the geodynamical evolution remain to be explained by other prospecting techniques.3-D structure velocities obtained by local seismic tomography and enhancement techniques of gravimetric anomalies are used to explore and define deep structure beneath the central High Atlas. The goal is to establish the evolution of the deep structure related to the geodynamical processes. Modest crustal thickness variation beneath the central High Atlas(~20-~40 km) define by local tomography and gravimetric anomalies, confirms that, a major part of the lower crust is detached into the lithosphere by delamination. Gravimetric anomaly, local seismic tomography and vertical cross sections throughout the central High Atlas, suggest that the lower crust detached is related to the broken slab of remain northward subduction beneath High and Middle Atlas. Meanwhile, extrusions of heated Asthenosphere materials induce the rifting stage concomitant to tectonic subsidence of the basin.

Three-dimensional seismic attenuation tomography of Egypt

The main purpose of the paper is assessing the three-dimensional (3-D) seismic tomography beneath Egyptto reveal the laws of the tectonic activity, dynamic features of the crust and the upper mantle as well as the thermal structure. Thecoda wave attenuation (Q-1c ) was obtained using the single scattering theory for the central frequencies of interest laid between 1 and 24 Hz. A regionalization of the estimated Q-1c values was performed by means of a generalized inversion technique.The obtained spatial distribution of 3-D attenuation results revealshigh contrasts between East and West Egypt. A remarkable contrast in the attenuation levels was compared with the tectonic structures, geothermal gradient and heat flow features. The highest attenuations are concentrated in the east and north western offshore regionsat central frequency 1.5 Hzthat draw a good matching with the seismic andthermal features of Egypt. Smaller attenuation levels were detected with young sediments of the Nile Valley from South to the northern triangle of Nile Delta basin except seismic active areas. Low or normal attenuation was detected at western desert where there is a stable and simple shelf. We can conclude that the extended highest attenuation joins to the strong seismic sources and geothermal structures at lower frequency and the centralized high attenuation takes place at moderate seismic sources at a higher frequency. The 3D attenuation maps can draw not only tectonic and geothermal structures but also the general geologic structure map.

High Resolution in Seismic Refraction Tomography for Environmental Study

Seismic refraction tomography (SRT) involves more complex mathematic algorithms to fit more flexible model. In the field procedure SRT in generally needs more shot points than standard seismic refraction survey to obtain high resolution profile. In this seismic refraction study, we have used 9 shot-points for inline and 10 shot-points for offset in purpose to obtaine high resolution of seismic refraction tomography. During a recent geophysical test site, the subsurface material was mapped along survey line using seismic refraction method. Analyses of the site investigation data revealed that the studied site was made up of two layers of the subsurface. The upper layer has velocity values with range of 500 m/s to1500 m/s which can be classified as unconsolidated surface deposits and mixtures of unsaturated sands and gravels. Meanwhile the lower layer has velocity values with range of2000 m/s to5500 m/s which is classified as compacted fine’s soil due to high pressure of the overburden. Analysis of seismic refraction data demonstrated that refraction tomography software systems are able to reveal subsurface material which represented by their seismic velocity value. Furthermore, the velocity model obtained in this study is agreed with its synthetic modelling result as initial model. This validity and reasonable results was able to assist in interpretation of the seismic refraction method for the environmental study.

Investigating the Depth and the Geometry of the Quarzitic Panafrican Basement Using Near-Surface 3D Seismic Refraction Tomography: Case Study of the Locality of Bakel (Senegal)

Seismic refraction investigations have been carried out in Bakel, Eastern Senegal. The purpose was to map geometrical relationship between the existing rock types and the Panafrican quarzitic basement, which is valuable information for the project of the Bakel fluviatile port construction. Four seismic refraction profiles were acquired. The obtained data have been processed by inversion. The obtained four seismic P-wave velocity profiles have been integrated to obtain a 3D model. By comparing the outcropping geological formations with the observed seismic data at the surface, it was possible to identify the lithology corresponding to each measured range of seismic velocity for the alluvium, the weathered bed rock, and the fresh rock. The results showed that the depth of the fresh rock of the basement varies from 0 to 18 meters above the sea level, with a deepening toward the Senegal River and toward the Northern part of the studied area. The presence of alluviums and their thickness are linked to the existence of bays and gulfs. The results of this study give valuable information for the river bed dredging cost assessment prior to the port construction phase.

Soil and Subsurface Sediment Microzonation Using with Seismic Refraction Tomography for Site Assessment (Case Study: IKIA Airport, Iran)

The site effects relating to the amplification of ground motion under earthquake loading are strongly influenced by both the subsurface soil condition and the geologic structure. In this study, for site characterization at the Imam Khomeini International Airport (IKIA) area in south of Tehran, in-situ seismic refraction tomography were carried out as a part of site investigations project, in addition geologic setting, borehole drilling, ground waters information and measurements. Based on seismic refraction studies, three layers are separable which with increasing in depth the S and P wave velocity is increased and this indicates increasing in compaction of soil and geologic materials. In the second and third separated layers, the zones with low and high seismic shear wave velocity is approximately equal, and northeast and southwest of the airport site has the low velocities, in addition to containing loose soils, highly weathered stones, and low depth to groundwater. In terms of Poisson’s ratio, the most important and key installations of airport site are located in suitable positions. According to Iranian Seismic Code, most of the lands around the airport are in class 2 and 3. It seems that a fault or a discontinuity is passed from northwest to the southeast of the study area. This site, according to geological, subsurface geophysical, and geotechnical boreholes studies, is high risk-earthquake prone.

Seismic active faults in the northwestern Beijing by seismic tomography

Using the arrival times of 197 earthquakes well recorded by the Beijing Seismic Network and the China-Germany cooperative Yanqing Digital Seismic Array in the period from October of 2001 to December of 2005, we accu- rately relocate these earthquakes by the joint-inversion program for hypocentral position and 3-D velocity structure. The distribution of hypocenters shows that there are two major seismic active belts in the northwestern Beijing. The first belt stretches from Qinghe, Wenquan to the NW direction and the second extends from Nankou, Chang- ping to the ENE direction.

A Gradient Regularization Method in Crosswell Seismic Tomography

Crosswell seismic tomography can be used to study the lateral variation of reservoirs, reservoir properties and the dynamic movement of fluids. In view of the instability of crosswell seismic tomography, the gradient method was improved by introducing regularization, and a gradient regularization method is presented in this paper. This method was verified by processing numerical simulation data and physical model data.

A Method to Estimate Spatial Resolution in 2-D Seismic Surface Wave Tomographic Problems

A novel methodology to quantify the spatial resolution in 2-D seismic surface wave tomographic problems is proposed in this study. It is based on both the resolving kernels computed via full resolution matrix and the concept of Full Width at Half Maximum (FWHM) of a Gaussian function. This method allows estimating quantitatively the spatial resolution at any cell of a gridded area. It was applied in the northeastern Brazil and the estimated spatial resolution range is in agreement with all previous surface wave investigations in the South America continent.

层析+折射约束迭代反演近地表建模方法及应用探索

塔西南山前带是塔里木盆地重要的勘探接替领域之一,其地表为巨厚黄土覆盖,黄土厚度可达500 m以上,速度在纵、横向上变化快,表层结构复杂。为改善地震成像效果,针对复杂地表区表层模型精度难以满足需求的问题,开展了拟合时深曲线量板、层析反演及折射反演等方面研究工作。结果表明:多种表层建模方法综合运用是提高近地表模型精度的有效途径;黄土曲线量板可反映连续介质区近地表速度变化规律;层析+折射约束迭代反演近地表建模方法能够将模型精度误差控制在5%以内;逐级递进约束层析反演可拓展可信速度场深度约1500 m。以上成果认识,对类似探区构建高精度近地表模型及提高地震成像精度具有指导意义。