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.

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.

Three-Dimensional Density Distribution and Seismic Activity along the Guxiang–Tongmai Segment of the Jiali Fault,Tibet

The Guxiang-Tongmai segment of the Jiali fault is situated northeast of the Namche Barwa Syntaxis in northeastern Tibet.It is one of the most active strike-slip faults near the syntaxis and plays a pivotal role in the examination of seismic activity within the eastern Himalayan Syntaxis.New study in the research region has yielded a 1:200000 gravity dataset covering an area 1500 km^(2).Using wavelet transform multiscale decomposition,scratch analysis techniques,and 3D gravity inversion methods,gravity anomalies,fault distributions,and density structures were determined across various scales.Through the integration of our new gravity data with other geophysical and geological information,our findings demonstrate substantial variations in the overall crustal density within the region,with the fault distribution closely linked to these density fluctuations.Disparities in stratigraphic density are important causes of variations in the capacity of geological formations to endure regional tectonic stress.Earthquakes are predominantly concentrated within the density transition zone and are primarily situated in regions of elevated density.The hanging wall stress within the Guxiang-Tongmai segment of the Jiali fault exhibits a notable concentration,marked by pronounced anisotropy,and is positioned within the density differential zone,which is prone to earthquakes.

Seismic energy dispersion compensation by multi-scale morphology

Seismic energy decays while propagating subsurface, which may reduce the resolution of seismic data. This paper studies the method of seismic energy dispersion compensation which provides the basic principles for multi-scale morphology and the spectrum simulation method. These methods are applied in seismic energy compensation. First of all, the seismic data is decomposed into multiple scales and the effective frequency bandwidth is selectively broadened for some scales by using a spectrum simulation method. In this process, according to the amplitude spectrum of each scale, the best simulation range is selected to simulate the middle and low frequency components to ensure the authenticity of the simulation curve which is calculated by the median method, and the high frequency component is broadened. Finally, these scales are reconstructed with reasonable coefficients, and the compensated seismic data can be obtained. Examples are shown to illustrate the feasibility of the energy compensation method.

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 smgma 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

A three-dimensional local-scale P-velocity model down to 25 km depth around the main shock epicenter region was constructed using 83821 event-to-receiver seismic rays from 5856 aftershocks recorded by a newly deployed temporary seismic network. Checkerboard tests show that our tomographic model has lateral and vertical resolution of -2 km. The high-resolution P-velocity model revealed interesting structures in the seismogenic layer： （1） The Guanxian-Anxian fault, Yingxiu-Beichuan fault and Wenchuan-Maoxian fault of the Longmen Shan fault zone are well delineated by sharp upper crustal velocity changes; （2） The Pengguan massif has generally higher velocity than its surrounding areas, and may extend down to at least -10 km from the surface; （3） A sharp lateral velocity variation beneath the Wenchuan-Maoxian fault may indicate that the Pengguan massif＇s western boundary and/or the Wenchuan-Maoxian fault is vertical, and the hypocenter of the Wenchuan earthquake possibly located at the conjunction point of the NW dipping Yingxiu-Beichuan and Guanxian-Anxian faults, and vertical Wenchuan-Maoxian fault; （4） Vicinity along the Yingxiu- Beichuan fault is characterized by very low velocity and low seismicity at shallow depths, possibly due to high content of porosity and fractures; （5） Two blocks of low-velocity anomaly are respectively imaged in the hanging wall and foot wall of the Guanxian-Anxian fault with a -7 km offset with -5 km vertical component.

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 tomo- graphy. 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.

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 accurately 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, Changping to the ENE direction.

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.

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.

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.

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.

龙门山双复杂区表层结构调查方法研究

龙门山山前带地表地质条件复杂,浅表层速度及厚度差异大,给激发、接收和静校正等工作带来较大困难。表层调查工作直接影响到地下介质成像效果。为使该地区进行的地震勘探攻关能够获取准确的静校正量数据,给井深设计提供依据,研究适合该区域的表层调查方法,结合地质剖面、浅震、小层析、地面微测井、钻井取心、高密度电法、三分量共振表层调查方法的工作原理,分别应用不同方法做比对,分析了不同调查方法的应用效果。通过不同表层调查的比对,分析了不同方法的适用范围,为该区选择合适的表层调查方法及方法参数选取提供依据,同时为其他地表复杂区域进行表层调查提供参考。