Global acoustic simulations are significant in revealing the internal and physical structure of the Earth. However, due to the limited flexibility of grids and the difficulties in handling boundaries, the traditional finite-difference method (FDM) is usually less used in global simulations. Nevertheless, these issues can be well resolved by employing a multi-block structured grid to discretize circular regions. In this paper, we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method (CGFDM) to solve the acoustic wave equation within this circular domain. By appropriately stretching the grid, the interconnections between each grid block are sufficiently smooth for stable information exchange. To verify the efficacy of this method, we conducted three numerical experiments, by comparing results with alternative approaches. Our test results demonstrate good agreement between our findings and the reference solutions. Since the proposed algorithm can effectively solve wave propagation problems in circular regions, it can contribute to 2D global simulation, particularly in interpreting the Earth’s interior.
{"title":"Simulation of acoustic waves in 2D circular regions using the multi-block structured grid finite-difference method","authors":"Changjiang Fang , Zhenguo Zhang , Hengxin Ren , Heng Zhang , Wenqiang Wang , Zhaoyang Lu , Qiwei Zhan , Wei Zhang , Qinghua Huang , Xiaofei Chen","doi":"10.1016/j.eqs.2025.08.005","DOIUrl":"10.1016/j.eqs.2025.08.005","url":null,"abstract":"<div><div>Global acoustic simulations are significant in revealing the internal and physical structure of the Earth. However, due to the limited flexibility of grids and the difficulties in handling boundaries, the traditional finite-difference method (FDM) is usually less used in global simulations. Nevertheless, these issues can be well resolved by employing a multi-block structured grid to discretize circular regions. In this paper, we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method (CGFDM) to solve the acoustic wave equation within this circular domain. By appropriately stretching the grid, the interconnections between each grid block are sufficiently smooth for stable information exchange. To verify the efficacy of this method, we conducted three numerical experiments, by comparing results with alternative approaches. Our test results demonstrate good agreement between our findings and the reference solutions. Since the proposed algorithm can effectively solve wave propagation problems in circular regions, it can contribute to 2D global simulation, particularly in interpreting the Earth’s interior.</div></div>","PeriodicalId":46333,"journal":{"name":"Earthquake Science","volume":"38 6","pages":"Pages 531-544"},"PeriodicalIF":4.1,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145719030","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-01DOI: 10.1016/j.eqs.2025.08.001
Ao Li , Zhengyuan Jia , Guoming Jiang , Guibin Zhang , Dapeng Zhao
The Cathaysia block in Southeast China has undergone complex compression, collision and tectonic evolution processes, which have led to the formation of rich granites and polymetallic ores. However, its history of crustal evolution and metallogenic mechanisms remain unclear. To explore the deep geodynamic processes in the central Cathaysia block, we deployed a dense broadband seismic array with a station interval of ∼5 km from July 2017 to August 2020. This seismic array, composed of 80 portable stations, traverses the Nanling and Wuyi metallogenic belts and crosses four deep faults. We employ the teleseismic receiver-function method to estimate the crustal thickness and vP/vS ratio beneath the seismic array. Our results reveal a distinct negative correlation between the crustal thickness and the vP/vS ratio within the crust, which might be closely associated with crustal deformation. The average crustal thickness is ∼30 km, indicating that the crust was thinned under an extensional environment. However, the crust beneath the Wuyi belt is slightly thicker and has a lower vP/vS ratio than its adjacent areas, which may reflect thickening of the upper crust due to its ductility and lithospheric delamination.
{"title":"Crustal structure of the central Cathaysia block in South China derived from receiver functions","authors":"Ao Li , Zhengyuan Jia , Guoming Jiang , Guibin Zhang , Dapeng Zhao","doi":"10.1016/j.eqs.2025.08.001","DOIUrl":"10.1016/j.eqs.2025.08.001","url":null,"abstract":"<div><div>The Cathaysia block in Southeast China has undergone complex compression, collision and tectonic evolution processes, which have led to the formation of rich granites and polymetallic ores. However, its history of crustal evolution and metallogenic mechanisms remain unclear. To explore the deep geodynamic processes in the central Cathaysia block, we deployed a dense broadband seismic array with a station interval of ∼5 km from July 2017 to August 2020. This seismic array, composed of 80 portable stations, traverses the Nanling and Wuyi metallogenic belts and crosses four deep faults. We employ the teleseismic receiver-function method to estimate the crustal thickness and <em>v</em><sub>P</sub>/<em>v</em><sub>S</sub> ratio beneath the seismic array. Our results reveal a distinct negative correlation between the crustal thickness and the <em>v</em><sub>P</sub>/<em>v</em><sub>S</sub> ratio within the crust, which might be closely associated with crustal deformation. The average crustal thickness is ∼30 km, indicating that the crust was thinned under an extensional environment. However, the crust beneath the Wuyi belt is slightly thicker and has a lower <em>v</em><sub>P</sub>/<em>v</em><sub>S</sub> ratio than its adjacent areas, which may reflect thickening of the upper crust due to its ductility and lithospheric delamination.</div></div>","PeriodicalId":46333,"journal":{"name":"Earthquake Science","volume":"38 6","pages":"Pages 545-563"},"PeriodicalIF":4.1,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145719031","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-01DOI: 10.1016/j.eqs.2025.08.002
Jiang Cheng, Hua Pan, Changlong Li, Yao-Hu Zhang
In this study, a seismic hazard algorithm was developed, coupled with a multi-source model. The algorithm can calculate the seismic hazard of 2D planar sources and 3D fault sources. A point source model is used to calculate the 2D planar potential seismic source by discretizing the potential seismic source into grid points which represent seismic sources. In contrast, a rupture surface model is used to calculate the 3D fault source so as to account for the influence of rupture scale in large earthquakes. The corresponding computational software was developed in Java and was tested on a study area 39.2°N–40.4°N and 116.6°E–118.7°E; the Tangshan and Xiadian faults were used to construct the fault source. Compared with the hazard results obtained using a point source model to calculate a 2D planar potential seismic source, the hazard results will be higher when a rupture surface model is used to calculate a 3D fault source. Moreover, a lower exceedance probability corresponds to a greater contribution rate of the fault source to the hazard and a greater seismic hazard closer to the fault. The algorithm can be used for seismic hazard analysis and can serve as a reference for the development of next-generation methods for zoning areas based on seismic ground motion.
{"title":"A fault source-based algorithm for probabilistic analysis of seismic hazard incorporating the rupture scale of large earthquakes","authors":"Jiang Cheng, Hua Pan, Changlong Li, Yao-Hu Zhang","doi":"10.1016/j.eqs.2025.08.002","DOIUrl":"10.1016/j.eqs.2025.08.002","url":null,"abstract":"<div><div>In this study, a seismic hazard algorithm was developed, coupled with a multi-source model. The algorithm can calculate the seismic hazard of 2D planar sources and 3D fault sources. A point source model is used to calculate the 2D planar potential seismic source by discretizing the potential seismic source into grid points which represent seismic sources. In contrast, a rupture surface model is used to calculate the 3D fault source so as to account for the influence of rupture scale in large earthquakes. The corresponding computational software was developed in Java and was tested on a study area 39.2°N–40.4°N and 116.6°E–118.7°E; the Tangshan and Xiadian faults were used to construct the fault source. Compared with the hazard results obtained using a point source model to calculate a 2D planar potential seismic source, the hazard results will be higher when a rupture surface model is used to calculate a 3D fault source. Moreover, a lower exceedance probability corresponds to a greater contribution rate of the fault source to the hazard and a greater seismic hazard closer to the fault. The algorithm can be used for seismic hazard analysis and can serve as a reference for the development of next-generation methods for zoning areas based on seismic ground motion.</div></div>","PeriodicalId":46333,"journal":{"name":"Earthquake Science","volume":"38 6","pages":"Pages 575-589"},"PeriodicalIF":4.1,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145719032","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-01DOI: 10.1016/j.eqs.2025.02.002
Yixiao Zhang, Ruifeng Liu, Zibo Wang, Zan Li
Radiated seismic energy is a quantitative characteristic of an earthquake that depends not only on the initial and final stresses, but also on the rupture history, and reflects the dynamic characteristics of the source. The seismic wave energy radiated per unit of seismic moment, or energy-moment ratio, represents the balance between the stress drop, source rupture velocity, and variation in the shear wave velocity. Earthquakes with a high energy-moment ratio thus release strain energy more rapidly. The accurate and rapid determination of radiated seismic energy and the energy-moment ratio play an important role in seismic hazard assessment, quantitative earthquake research, and engineering seismology research. In this study, waveform data from the Global Seismographic Network were used to measure the dynamic source parameters of an earthquake that occurred on January 7, 2025, in Dingri, Xizang, with the radiated energy, energy-moment ratio, slowness parameter, and apparent stress investigated. Static source parameters such as the seismic moment and moment magnitude were also determined. The main results were as follows: (1) the radiated energy of the earthquake was 9.73×1014 J, corresponding to an energy magnitude ME of 7.1, with a source rupture time of 24 s; (2) the focal mechanism was normal faulting, with the seismic moment of 4.98×1019 N·m corresponding to a moment magnitude MW of 7.1. Nodal plane I was focused at 191°/32°/−67° while plane II was at 344°/60°/−104°, and the centroid depth was 12.3 km; (3) the energy-moment ratio of the earthquake was 1.95×10–5, the slowness parameter was −4.71, and the apparent stress was 0.59 MPa. The energy-moment ratio was thus higher than the average for normal fault earthquakes on the Chinese mainland. In conclusion, the results indicated that the 2025 earthquake was a normal fault earthquake with relatively high energy release efficiency and significant potential for damage to local buildings and the infrastructure, as verified by the severe damage to ground structures and significant casualties nearby.
{"title":"Determination of radiated energy and energy-moment ratio for the 2025 Dingri, Xizang M6.8 earthquake","authors":"Yixiao Zhang, Ruifeng Liu, Zibo Wang, Zan Li","doi":"10.1016/j.eqs.2025.02.002","DOIUrl":"10.1016/j.eqs.2025.02.002","url":null,"abstract":"<div><div>Radiated seismic energy is a quantitative characteristic of an earthquake that depends not only on the initial and final stresses, but also on the rupture history, and reflects the dynamic characteristics of the source. The seismic wave energy radiated per unit of seismic moment, or energy-moment ratio, represents the balance between the stress drop, source rupture velocity, and variation in the shear wave velocity. Earthquakes with a high energy-moment ratio thus release strain energy more rapidly. The accurate and rapid determination of radiated seismic energy and the energy-moment ratio play an important role in seismic hazard assessment, quantitative earthquake research, and engineering seismology research. In this study, waveform data from the Global Seismographic Network were used to measure the dynamic source parameters of an earthquake that occurred on January 7, 2025, in Dingri, Xizang, with the radiated energy, energy-moment ratio, slowness parameter, and apparent stress investigated. Static source parameters such as the seismic moment and moment magnitude were also determined. The main results were as follows: (1) the radiated energy of the earthquake was 9.73×10<sup>14</sup> J, corresponding to an energy magnitude <em>M</em><sub>E</sub> of 7.1, with a source rupture time of 24 s; (2) the focal mechanism was normal faulting, with the seismic moment of 4.98×10<sup>19</sup> N·m corresponding to a moment magnitude <em>M</em><sub>W</sub> of 7.1. Nodal plane I was focused at 191°/32°/−67° while plane II was at 344°/60°/−104°, and the centroid depth was 12.3 km; (3) the energy-moment ratio of the earthquake was 1.95×10<sup>–5</sup>, the slowness parameter was −4.71, and the apparent stress was 0.59 MPa. The energy-moment ratio was thus higher than the average for normal fault earthquakes on the Chinese mainland. In conclusion, the results indicated that the 2025 earthquake was a normal fault earthquake with relatively high energy release efficiency and significant potential for damage to local buildings and the infrastructure, as verified by the severe damage to ground structures and significant casualties nearby.</div></div>","PeriodicalId":46333,"journal":{"name":"Earthquake Science","volume":"38 6","pages":"Pages 564-574"},"PeriodicalIF":4.1,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145719033","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-01DOI: 10.1016/j.eqs.2025.08.004
Zeyan Zhao , Han Yue , Jian Zhang , Yan Hu , Yijian Zhou , Jing Liu-Zeng , Kang Wang
The 2015 Gorkha (MW=7.8) earthquake ruptured the downdip portion of the Main Himalayan Thrust. Afterslip following this event provides valuable insights into the frictional properties on the thrust interface, yet its amplitude and distribution remain controversial. In this study, we incorporate long-term GNSS and InSAR data and correct for the viscoelastic relaxation simulated using a regional 3-D viscoelastic model. We adopt the corrected data in a novel inversion algorithm and resolve two spatially separated afterslip processes with different decay times: fast afterslip near the bottom of the coseismic rupture, possibly stopped the mainshock and triggered the MW7.3 aftershock 17 days later, and slow afterslip extending further downdip. By comparing the afterslip and aftershock patterns, we identify distinct partitioning of seismic and aseismic slip behaviors at the bottom of the seismogenic zone, which reflects local heterogeneities in frictional properties at the transition depths.
{"title":"Co-existing fast and slow afterslip processes following the 2015 Gorkha (MW7.8) earthquake resolved by full time-series inversion","authors":"Zeyan Zhao , Han Yue , Jian Zhang , Yan Hu , Yijian Zhou , Jing Liu-Zeng , Kang Wang","doi":"10.1016/j.eqs.2025.08.004","DOIUrl":"10.1016/j.eqs.2025.08.004","url":null,"abstract":"<div><div>The 2015 Gorkha (<em>M</em><sub>W</sub>=7.8) earthquake ruptured the downdip portion of the Main Himalayan Thrust. Afterslip following this event provides valuable insights into the frictional properties on the thrust interface, yet its amplitude and distribution remain controversial. In this study, we incorporate long-term GNSS and InSAR data and correct for the viscoelastic relaxation simulated using a regional 3-D viscoelastic model. We adopt the corrected data in a novel inversion algorithm and resolve two spatially separated afterslip processes with different decay times: fast afterslip near the bottom of the coseismic rupture, possibly stopped the mainshock and triggered the <em>M</em><sub>W</sub>7.3 aftershock 17 days later, and slow afterslip extending further downdip. By comparing the afterslip and aftershock patterns, we identify distinct partitioning of seismic and aseismic slip behaviors at the bottom of the seismogenic zone, which reflects local heterogeneities in frictional properties at the transition depths.</div></div>","PeriodicalId":46333,"journal":{"name":"Earthquake Science","volume":"38 6","pages":"Pages 485-503"},"PeriodicalIF":4.1,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145719035","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-01DOI: 10.1016/j.eqs.2025.08.003
Yi Wang , Baichen Wu , Yuqing Zheng , Yan Luo , Xiaohui He , Manzhong Qin
High resolution imaging of the fault zone structure is crucial to understanding the characteristics of strong earthquake activity and the deep seismogenic environment. In seismological studies, the fault zone is generally considered to be a low velocity zone with host rock on both sides. In order to determine the main parameters of fault zone, such as physical properties and interface characteristics, many efforts have been made. However, many key fault parameters still lack constraints, such as the depth extent, width and dip angle of the low velocity zone. With the advancement of the large-N array techniques in recent years, seismologists have collected high-quality data with larger apertures and denser arrays for better analysis of fault zone structures. These array data have also facilitated the development of new seismic imaging techniques. In this paper, a new waveform inversion method for fault zone parameters based on generalized teleseismic waveforms is proposed. Generalized teleseismic event is defined as the local seismic signal whose epicentral distance is greater than 7−10 times the aperture of the array. In order to efficiently simulate high frequency wavefield propagation from long distance local earthquakes, a hybrid modeling approach is proposed, which greatly reduces the computational cost for teleseismic waveform inversion. We apply the proposed new inversion method to a dense array data across an inactive fault in the Qilian Mountains, Gansu Province. As an active-source analogue of generalized teleseismic, the recorded waveforms of a 270-meter-long linear array are clearly excited by an airgun source 1.8 km away. Setting cross-correlation travel time of direct P wave as the misfit function, we perform waveform inversion for the main structural parameters of the fault zone through grid search strategy. The new method is particularly suitable for imaging fault zones with limited local seismicity.
{"title":"Waveform inversion of the fault zone structure based on generalized teleseismic wave records","authors":"Yi Wang , Baichen Wu , Yuqing Zheng , Yan Luo , Xiaohui He , Manzhong Qin","doi":"10.1016/j.eqs.2025.08.003","DOIUrl":"10.1016/j.eqs.2025.08.003","url":null,"abstract":"<div><div>High resolution imaging of the fault zone structure is crucial to understanding the characteristics of strong earthquake activity and the deep seismogenic environment. In seismological studies, the fault zone is generally considered to be a low velocity zone with host rock on both sides. In order to determine the main parameters of fault zone, such as physical properties and interface characteristics, many efforts have been made. However, many key fault parameters still lack constraints, such as the depth extent, width and dip angle of the low velocity zone. With the advancement of the large-<em>N</em> array techniques in recent years, seismologists have collected high-quality data with larger apertures and denser arrays for better analysis of fault zone structures. These array data have also facilitated the development of new seismic imaging techniques. In this paper, a new waveform inversion method for fault zone parameters based on generalized teleseismic waveforms is proposed. Generalized teleseismic event is defined as the local seismic signal whose epicentral distance is greater than 7−10 times the aperture of the array. In order to efficiently simulate high frequency wavefield propagation from long distance local earthquakes, a hybrid modeling approach is proposed, which greatly reduces the computational cost for teleseismic waveform inversion. We apply the proposed new inversion method to a dense array data across an inactive fault in the Qilian Mountains, Gansu Province. As an active-source analogue of generalized teleseismic, the recorded waveforms of a 270-meter-long linear array are clearly excited by an airgun source 1.8 km away. Setting cross-correlation travel time of direct P wave as the misfit function, we perform waveform inversion for the main structural parameters of the fault zone through grid search strategy. The new method is particularly suitable for imaging fault zones with limited local seismicity.</div></div>","PeriodicalId":46333,"journal":{"name":"Earthquake Science","volume":"38 6","pages":"Pages 504-530"},"PeriodicalIF":4.1,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145719029","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-01DOI: 10.1016/j.eqs.2025.08.006
Jiaxin Liu , Shunping Pei , Yicun Guo
The 1975 Haicheng earthquake represents the first successful prediction of a major earthquake in China and is the only earthquake forecast officially recognized by the United Nations. Understanding the crustal structure in the Haicheng region is critical for determining the seismogenic mechanisms of large earthquakes. Pg-wave tomography of the Haicheng region was used to obtain the upper crustal structure at depths of 5–10 km, in which lateral velocity variations indicate differences in tectonic activity within the seismogenic layer. A dataset comprising 62,610 Pg-wave arrival times was used to obtain high-resolution seismic velocity and anisotropy images of the upper crust in the Haicheng region. The tomography results indicate that a distinct high-velocity anomaly is located in the region that produced the Haicheng and Xiuyan earthquakes, as well as a few small earthquakes at the southern end of the Jinzhou fault. This suggests that a high-velocity asperity beneath the Haichenghe fault was able to accumulate stress due to long-term tectonic loading, eventually producing the Haicheng earthquake. The seismogenesis of Haicheng earthquake can also be used to explain other large earthquakes in the slowly deforming eastern region of China.
{"title":"High-resolution tomography of P-wave velocity structures in the Haicheng region: Implications for the seismogenesis of the 1975 MS7.3 Haicheng earthquake","authors":"Jiaxin Liu , Shunping Pei , Yicun Guo","doi":"10.1016/j.eqs.2025.08.006","DOIUrl":"10.1016/j.eqs.2025.08.006","url":null,"abstract":"<div><div>The 1975 Haicheng earthquake represents the first successful prediction of a major earthquake in China and is the only earthquake forecast officially recognized by the United Nations. Understanding the crustal structure in the Haicheng region is critical for determining the seismogenic mechanisms of large earthquakes. Pg-wave tomography of the Haicheng region was used to obtain the upper crustal structure at depths of 5–10 km, in which lateral velocity variations indicate differences in tectonic activity within the seismogenic layer. A dataset comprising 62,610 Pg-wave arrival times was used to obtain high-resolution seismic velocity and anisotropy images of the upper crust in the Haicheng region. The tomography results indicate that a distinct high-velocity anomaly is located in the region that produced the Haicheng and Xiuyan earthquakes, as well as a few small earthquakes at the southern end of the Jinzhou fault. This suggests that a high-velocity asperity beneath the Haichenghe fault was able to accumulate stress due to long-term tectonic loading, eventually producing the Haicheng earthquake. The seismogenesis of Haicheng earthquake can also be used to explain other large earthquakes in the slowly deforming eastern region of China.</div></div>","PeriodicalId":46333,"journal":{"name":"Earthquake Science","volume":"38 6","pages":"Pages 590-600"},"PeriodicalIF":4.1,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145719034","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-09-02DOI: 10.1016/j.eqs.2025.06.006
Shuguang Wang , Shuoxian Ning , Zhixiang Yao , Jiaqi Li , Wanbo Xiao , Tianfan Yan , Feng Xu
The InSight mission has obtained seismic data from Mars, offering new insights into the planet’s internal structure and seismic activity. However, the raw data released to the public contain various sources of noise, such as ticks and glitches, which hamper further seismological studies. This paper presents step-by-step processing of InSight’s Very Broad Band seismic data, focusing on the suppression and removal of non-seismic noise. The processing stages include tick noise removal, glitch signal suppression, multicomponent synchronization, instrument response correction, and rotation of orthogonal components. The processed datasets and associated codes are openly accessible and will support ongoing efforts to explore the geophysical properties of Mars and contribute to the broader field of planetary seismology.
{"title":"Basic processing of the InSight seismic data from Mars for further seismological research","authors":"Shuguang Wang , Shuoxian Ning , Zhixiang Yao , Jiaqi Li , Wanbo Xiao , Tianfan Yan , Feng Xu","doi":"10.1016/j.eqs.2025.06.006","DOIUrl":"10.1016/j.eqs.2025.06.006","url":null,"abstract":"<div><div>The InSight mission has obtained seismic data from Mars, offering new insights into the planet’s internal structure and seismic activity. However, the raw data released to the public contain various sources of noise, such as ticks and glitches, which hamper further seismological studies. This paper presents step-by-step processing of InSight’s Very Broad Band seismic data, focusing on the suppression and removal of non-seismic noise. The processing stages include tick noise removal, glitch signal suppression, multicomponent synchronization, instrument response correction, and rotation of orthogonal components. The processed datasets and associated codes are openly accessible and will support ongoing efforts to explore the geophysical properties of Mars and contribute to the broader field of planetary seismology.</div></div>","PeriodicalId":46333,"journal":{"name":"Earthquake Science","volume":"38 5","pages":"Pages 450-460"},"PeriodicalIF":4.1,"publicationDate":"2025-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144933774","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-09-02DOI: 10.1016/j.eqs.2025.06.002
Shanshan Jia , Laiyu Lu , Yutao Shi , Pingping Wu , Lijun Chang
The Sichuan-Yunnan Block is located on the southeastern margin of the Qinghai-Xizang Plateau and has frequent seismic activity on the western border, posing a potential threat to human society and economic development. Therefore, it is important to understand its geological evolution, assess earthquake risks, and formulate scientific and reasonable disaster prevention and mitigation strategies. Using 23 months of continuous ambient noise records from 81 seismic stations, we obtained 1248 phase-velocity dispersion curves of the fundamental Rayleigh wave at 5–50 s. The three-dimensional (3D) S-wave velocity structure in the northwestern Sichuan-Yunnan Block was obtained by pure-path and depth inversion. The results show that three low-velocity anomalous bands were distributed nearly north-to-south (N-S) at depths of 10–35 km. The overall shape of the low-velocity channel gradually shifted from southeast to southwest because of the influence of the Panzhihua high-velocity blocks. The low-velocity strip consists of three branches, with the first branch extending southwest from the northern part of the Lancangjiang Fault. The second branch is distributed in the N-S direction and is blocked by two high-velocity bodies near the Longpan-Qiaohou and Honghe faults. The third branch crosses the research area from N-S and gradually extends from southeast to southwest and from shallow to deep. The three low-velocity anomaly distribution areas are likely the most severely deformed areas of the collision between the Qinghai-Xizang Plateau and Yangtze Block. The results provide a more detailed understanding of the deep structure of the western boundary of the Sichuan-Yunnan Block crustal low-velocity anomalies and reliable geophysical evidence for the morphology and continuity of crustal flows.
{"title":"High-resolution 3D S-wave velocity structure in northwestern Sichuan-Yunnan Block derived from ambient noise tomography","authors":"Shanshan Jia , Laiyu Lu , Yutao Shi , Pingping Wu , Lijun Chang","doi":"10.1016/j.eqs.2025.06.002","DOIUrl":"10.1016/j.eqs.2025.06.002","url":null,"abstract":"<div><div>The Sichuan-Yunnan Block is located on the southeastern margin of the Qinghai-Xizang Plateau and has frequent seismic activity on the western border, posing a potential threat to human society and economic development. Therefore, it is important to understand its geological evolution, assess earthquake risks, and formulate scientific and reasonable disaster prevention and mitigation strategies. Using 23 months of continuous ambient noise records from 81 seismic stations, we obtained 1248 phase-velocity dispersion curves of the fundamental Rayleigh wave at 5–50 s. The three-dimensional (3D) S-wave velocity structure in the northwestern Sichuan-Yunnan Block was obtained by pure-path and depth inversion. The results show that three low-velocity anomalous bands were distributed nearly north-to-south (N-S) at depths of 10–35 km. The overall shape of the low-velocity channel gradually shifted from southeast to southwest because of the influence of the Panzhihua high-velocity blocks. The low-velocity strip consists of three branches, with the first branch extending southwest from the northern part of the Lancangjiang Fault. The second branch is distributed in the N-S direction and is blocked by two high-velocity bodies near the Longpan-Qiaohou and Honghe faults. The third branch crosses the research area from N-S and gradually extends from southeast to southwest and from shallow to deep. The three low-velocity anomaly distribution areas are likely the most severely deformed areas of the collision between the Qinghai-Xizang Plateau and Yangtze Block. The results provide a more detailed understanding of the deep structure of the western boundary of the Sichuan-Yunnan Block crustal low-velocity anomalies and reliable geophysical evidence for the morphology and continuity of crustal flows.</div></div>","PeriodicalId":46333,"journal":{"name":"Earthquake Science","volume":"38 5","pages":"Pages 408-426"},"PeriodicalIF":4.1,"publicationDate":"2025-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144933776","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-09-02DOI: 10.1016/j.eqs.2025.06.004
Lei Li , Jiacheng Zhang , Yuyang Tan , Ling Peng , Junlun Li , Jincheng Xu , Jianxin Liu
Seismic source locations can characterize the spatial and temporal distributions of seismic sources, and can provide important basic data for earthquake disaster monitoring, fault activity characterization, and fracture growth interpretation. Waveform stacking-based location methods invert the source locations by focusing the source energy with multichannel waveforms, and these methods exhibit a high level of automation and noise-resistance. Taking the cross-correlation stacking (CCS) method as an example, this work attempts to study the influential factors of waveform stacking-based methods, and introduces a comprehensive performance evaluation scheme based on multiple parameters and indicators. The waveform data are from field monitoring of induced microseismicity in the Changning region (southern Sichuan Basin of China). Synthetic and field data tests reveal the impacts of three categories of factors on waveform stacking-based location: velocity model, monitoring array, and waveform complexity. The location performance is evaluated and further improved in terms of the source imaging resolution and location error. Denser array monitoring contributes to better constraining source depth and location reliability, but the combined impact of multiple factors, such as velocity model uncertainty and multiple seismic phases, increases the complexity of locating field microseismic events. Finally, the aspects of location uncertainty, phase detection, and artificial intelligence-based location are discussed.
{"title":"Performance evaluation of the waveform stacking-based microseismic location method in the southern Sichuan Basin of China","authors":"Lei Li , Jiacheng Zhang , Yuyang Tan , Ling Peng , Junlun Li , Jincheng Xu , Jianxin Liu","doi":"10.1016/j.eqs.2025.06.004","DOIUrl":"10.1016/j.eqs.2025.06.004","url":null,"abstract":"<div><div>Seismic source locations can characterize the spatial and temporal distributions of seismic sources, and can provide important basic data for earthquake disaster monitoring, fault activity characterization, and fracture growth interpretation. Waveform stacking-based location methods invert the source locations by focusing the source energy with multichannel waveforms, and these methods exhibit a high level of automation and noise-resistance. Taking the cross-correlation stacking (CCS) method as an example, this work attempts to study the influential factors of waveform stacking-based methods, and introduces a comprehensive performance evaluation scheme based on multiple parameters and indicators. The waveform data are from field monitoring of induced microseismicity in the Changning region (southern Sichuan Basin of China). Synthetic and field data tests reveal the impacts of three categories of factors on waveform stacking-based location: velocity model, monitoring array, and waveform complexity. The location performance is evaluated and further improved in terms of the source imaging resolution and location error. Denser array monitoring contributes to better constraining source depth and location reliability, but the combined impact of multiple factors, such as velocity model uncertainty and multiple seismic phases, increases the complexity of locating field microseismic events. Finally, the aspects of location uncertainty, phase detection, and artificial intelligence-based location are discussed.</div></div>","PeriodicalId":46333,"journal":{"name":"Earthquake Science","volume":"38 5","pages":"Pages 427-440"},"PeriodicalIF":4.1,"publicationDate":"2025-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144933777","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
扫码关注我们
求助内容：
应助结果提醒方式：
京公网安备 11010802042870号