The Xianshuihe fault is a major tectonic boundary between the Sichuan-Yunnan rhombic and Bayanhar blocks in Southwest China.With an average left-lateral strike-slip movement of 10–15 mm/yr,it is a fast-moving strike-...The Xianshuihe fault is a major tectonic boundary between the Sichuan-Yunnan rhombic and Bayanhar blocks in Southwest China.With an average left-lateral strike-slip movement of 10–15 mm/yr,it is a fast-moving strike-slip continental fault.On September 5,2022,the Ms6.8 Luding earthquake occurred along the Moxi segment of the Xianshuihe fault,reaching a maximum intensity of IX and resulting in a significant number of casualties and severe property damage.This earthquake broke the long-standing seismic quiescence of the Xianshuihe fault,which lasted for more than 40 years,and was followed by a significant number of aftershocks.An outstanding question is how the behavior of the Xianshuihe fault and major earthquakes changed following this mainshock.In this study,we examined the changes in regional seismicity following the Luding earthquake and identified the potential for future strong earthquakes along the Xianshuihe fault.We used a finite element numerical method to simulate the environment of the seismogenic fault and its adjacent areas.In addition,we used the coseismic slip model of the Luding earthquake with the split-node method to calculate how the stress and strain fields in the surrounding area were affected by the2022 mainshock.Coulomb stress changes were resolved in the main faults,and the seismicity of adjacent faults was analyzed in conjunction with the observed seismic data.The results indicate that regional tectonic movement primarily occurred to the southeast along the Moxi segment.The stress field is approximately north-south in tension and east-west in compression.Variation in the stress field in the epicentral region of the Luding earthquake exceeded 1 MPa.The maximum displacement of the coseismic deformation field was concentrated between Moxi town and Tuanjie village,and the Coulomb stress of the fault zone in this region experienced the largest decrease.However,the b-value of the Gutenberg-Richter magnitude-frequency relationship at the epicenter and the surrounding area exhibited an abnormal 展开更多
The subsurface fluid injection can cause pressure increase within faults,leading to earthquake occurrences.However,the factors controlling earthquake rupture due to pressure perturbation remain poorly understood.To re...The subsurface fluid injection can cause pressure increase within faults,leading to earthquake occurrences.However,the factors controlling earthquake rupture due to pressure perturbation remain poorly understood.To resolve this problem,we simulate the physical processes of earthquake nucleation and rupture on strike-slip faults perturbated by pressure migration based on the slip-weakening law.Multiple kinds of factors,including background stress,fluid injection rates,the area of the pressurized region,fault geometry,and fault friction coefficients,are considered in our simulations.Our simulation results reveal that the ratio of shear stress to normal stress rather than their absolute values controls the rupture behavior.With the large stress ratios,high injection rates,and large pressurized areas,earthquakes are prone to propagate as runaway ruptures.Additionally,faults with large aspect ratios of length to width are also favorable for causing runaway ruptures.In contrast,the factors of fault strike,dip angles and friction coefficients have minor influence on rupture behavior.展开更多
Based upon some simplified numerical models, a 2-D plain strain elastic FEM program is compiled to study the distributions of the stress fields produced by the volume change of the phase transformation from olivine to...Based upon some simplified numerical models, a 2-D plain strain elastic FEM program is compiled to study the distributions of the stress fields produced by the volume change of the phase transformation from olivine to spinel, by the volume change from temperature variation, and by density difference and boundary action in a piece of subducted slab located in transition zone of the mantle. Thermal stress could explain the fault plane solutions of deep focus earthquakes, but could not explain the distribution of deep seismicity. When large extent metastable olivine is included, the stress field produced by the density difference contradicts with the results of fault plane solutions and with the distribution of deep seismicity. Although the stress produced by volume change of the phase transformation from olivine to spinel dominates the stress state, its main direction is different from the observed results. We conclude that the deep seismicity could not be simply explained by elastic simulation.展开更多
Considering the heterogeneity of geomechanical materials, seismicity during brittle rock failure under compressive loading on the sample with an original weak zone is simulated by using rock failure process analysis c...Considering the heterogeneity of geomechanical materials, seismicity during brittle rock failure under compressive loading on the sample with an original weak zone is simulated by using rock failure process analysis code (RFPA2D). The run-through process of weak zone, the forming of new fault and associated micro-seismicities are studied. The modeling demonstrates the total process of source development of earthquake from deformation, micro-failure to collapse and the behavior of temporal-spatial distribution of micro-seismicities. The stress, strain and the temporal-spatial distribution of micro-seismicities life-likely portrayed the phenomena of localization and temporal-spatial transitions, which is similar to those observed in our real crust. Also, the results obtained in simulations are in agreement with or similar to the reported experimental observations.展开更多
The characteristics of the stress fields in deep subducting slabs are studied using viscoelastic plain strain finite element method. When introducing the new rheology structure given by Karato, et al into our computat...The characteristics of the stress fields in deep subducting slabs are studied using viscoelastic plain strain finite element method. When introducing the new rheology structure given by Karato, et al into our computation, there emerge two regions with great shear stress just below the olivine-spinel phase transition zone, which encompass the low viscosity zone below the lower tip of the metastable wedge. Further, the directions of the main compressional stress of these two regions are all along the dip direction of the slab. These are in accordance with the seismic observations that there are two deep seismic zones in a slab and the directions of the main compressional stress in these two seismic zones are along the dip direction of the slab. Smaller effective viscosity probably caused by smaller grain size in the phase transformation zone does not have great influence on the stress state in the deep part of the slab. There is the maximum of shear stress at the transition region from olivine to spinel and the direction of the main compressional stress in this region is roughly perpendicular to the trend of the phase boundary no matter whether there exists metastable wedge, which nevertheless do not correspond to some well-known seismic observations.展开更多
By dealing with strain and acoustic emission (AE) data for two parallelling faultss instability and failure with the same slip direction including asperities, the temporal-spatial evolution of strain and AE field dist...By dealing with strain and acoustic emission (AE) data for two parallelling faultss instability and failure with the same slip direction including asperities, the temporal-spatial evolution of strain and AE field distribution on the asperity of parallelling faults is analyzed. Furthermore the failure process of asperities and interaction among the asperities, i.e., positive and negative seismicity are discussed. Results show that instability and failure for the parallelling faults is a kind of negative seismicity.展开更多
By considering the heterogeneity of geomechanical materials, the source development of earthquake under compression boundary conditions is studied with a newly developed numerical method, Rock Failure Process Analysis...By considering the heterogeneity of geomechanical materials, the source development of earthquake under compression boundary conditions is studied with a newly developed numerical method, Rock Failure Process Analysis code (RFPA2D). The process of fault forming and associated micro seismicities in a rectangle area with a inclusion but without any clear structural features of original fault is modeled. The modeling demonstrates the whole process of source development of earthquake from deformation, micro failure to collapse and the behavior of temporal spatial distribution of micro seismicities. The stress, strain and the temporal spatial distribution of micro seismicities vividly depict the phenomena of localization, temporal transitions, dilatation or rise, elastic rebound and conjugate (X type) deformation zone.展开更多
基金supported by the National Key Research and Development Project of China(Grant No.2021YFC3000703-06)the Central Publicinterest Scientific Institution Basal Research Fund(Grant No.CEAIEF20230209)+1 种基金the National Natural Science Foundation of China(Grant Nos.U2239205,41725017,U2039207)partially supported by the National Key Scientific and Technological Infrastructure project“Earth System Science Numerical Simulator Facility(Earth Lab)”。
文摘The Xianshuihe fault is a major tectonic boundary between the Sichuan-Yunnan rhombic and Bayanhar blocks in Southwest China.With an average left-lateral strike-slip movement of 10–15 mm/yr,it is a fast-moving strike-slip continental fault.On September 5,2022,the Ms6.8 Luding earthquake occurred along the Moxi segment of the Xianshuihe fault,reaching a maximum intensity of IX and resulting in a significant number of casualties and severe property damage.This earthquake broke the long-standing seismic quiescence of the Xianshuihe fault,which lasted for more than 40 years,and was followed by a significant number of aftershocks.An outstanding question is how the behavior of the Xianshuihe fault and major earthquakes changed following this mainshock.In this study,we examined the changes in regional seismicity following the Luding earthquake and identified the potential for future strong earthquakes along the Xianshuihe fault.We used a finite element numerical method to simulate the environment of the seismogenic fault and its adjacent areas.In addition,we used the coseismic slip model of the Luding earthquake with the split-node method to calculate how the stress and strain fields in the surrounding area were affected by the2022 mainshock.Coulomb stress changes were resolved in the main faults,and the seismicity of adjacent faults was analyzed in conjunction with the observed seismic data.The results indicate that regional tectonic movement primarily occurred to the southeast along the Moxi segment.The stress field is approximately north-south in tension and east-west in compression.Variation in the stress field in the epicentral region of the Luding earthquake exceeded 1 MPa.The maximum displacement of the coseismic deformation field was concentrated between Moxi town and Tuanjie village,and the Coulomb stress of the fault zone in this region experienced the largest decrease.However,the b-value of the Gutenberg-Richter magnitude-frequency relationship at the epicenter and the surrounding area exhibited an abnormal
基金co-supported by the National Natural Science Foundation of China under Grants(42130101,42074007)the Fundamental Research Funds for the Central Universities under Grants 2042023kf1035。
文摘The subsurface fluid injection can cause pressure increase within faults,leading to earthquake occurrences.However,the factors controlling earthquake rupture due to pressure perturbation remain poorly understood.To resolve this problem,we simulate the physical processes of earthquake nucleation and rupture on strike-slip faults perturbated by pressure migration based on the slip-weakening law.Multiple kinds of factors,including background stress,fluid injection rates,the area of the pressurized region,fault geometry,and fault friction coefficients,are considered in our simulations.Our simulation results reveal that the ratio of shear stress to normal stress rather than their absolute values controls the rupture behavior.With the large stress ratios,high injection rates,and large pressurized areas,earthquakes are prone to propagate as runaway ruptures.Additionally,faults with large aspect ratios of length to width are also favorable for causing runaway ruptures.In contrast,the factors of fault strike,dip angles and friction coefficients have minor influence on rupture behavior.
基金Pre-elected National Important Fundamental Research Project (95-S-05) and Foundation for University Key Teacher
文摘Based upon some simplified numerical models, a 2-D plain strain elastic FEM program is compiled to study the distributions of the stress fields produced by the volume change of the phase transformation from olivine to spinel, by the volume change from temperature variation, and by density difference and boundary action in a piece of subducted slab located in transition zone of the mantle. Thermal stress could explain the fault plane solutions of deep focus earthquakes, but could not explain the distribution of deep seismicity. When large extent metastable olivine is included, the stress field produced by the density difference contradicts with the results of fault plane solutions and with the distribution of deep seismicity. Although the stress produced by volume change of the phase transformation from olivine to spinel dominates the stress state, its main direction is different from the observed results. We conclude that the deep seismicity could not be simply explained by elastic simulation.
基金The Development Program on National Key Basic Researches under the Project Mechanism and Prediction of Continental Strong Earthquakes (G19980407) State Natural Science Foundation (49974009).
文摘Considering the heterogeneity of geomechanical materials, seismicity during brittle rock failure under compressive loading on the sample with an original weak zone is simulated by using rock failure process analysis code (RFPA2D). The run-through process of weak zone, the forming of new fault and associated micro-seismicities are studied. The modeling demonstrates the total process of source development of earthquake from deformation, micro-failure to collapse and the behavior of temporal-spatial distribution of micro-seismicities. The stress, strain and the temporal-spatial distribution of micro-seismicities life-likely portrayed the phenomena of localization and temporal-spatial transitions, which is similar to those observed in our real crust. Also, the results obtained in simulations are in agreement with or similar to the reported experimental observations.
基金Pre-elected National Important Fundamental Research Project (95-S-05) and Foundation for University Key Teacher
文摘The characteristics of the stress fields in deep subducting slabs are studied using viscoelastic plain strain finite element method. When introducing the new rheology structure given by Karato, et al into our computation, there emerge two regions with great shear stress just below the olivine-spinel phase transition zone, which encompass the low viscosity zone below the lower tip of the metastable wedge. Further, the directions of the main compressional stress of these two regions are all along the dip direction of the slab. These are in accordance with the seismic observations that there are two deep seismic zones in a slab and the directions of the main compressional stress in these two seismic zones are along the dip direction of the slab. Smaller effective viscosity probably caused by smaller grain size in the phase transformation zone does not have great influence on the stress state in the deep part of the slab. There is the maximum of shear stress at the transition region from olivine to spinel and the direction of the main compressional stress in this region is roughly perpendicular to the trend of the phase boundary no matter whether there exists metastable wedge, which nevertheless do not correspond to some well-known seismic observations.
基金State Key Basic Research Development and Programming Project Mechanism and Prediction of Continental Strong Earthquake (G19980407) and Foundation of Laboratory of Tectonophysics China Seismological Bureau and State Natural Science Foundation of China (
文摘By dealing with strain and acoustic emission (AE) data for two parallelling faultss instability and failure with the same slip direction including asperities, the temporal-spatial evolution of strain and AE field distribution on the asperity of parallelling faults is analyzed. Furthermore the failure process of asperities and interaction among the asperities, i.e., positive and negative seismicity are discussed. Results show that instability and failure for the parallelling faults is a kind of negative seismicity.
文摘By considering the heterogeneity of geomechanical materials, the source development of earthquake under compression boundary conditions is studied with a newly developed numerical method, Rock Failure Process Analysis code (RFPA2D). The process of fault forming and associated micro seismicities in a rectangle area with a inclusion but without any clear structural features of original fault is modeled. The modeling demonstrates the whole process of source development of earthquake from deformation, micro failure to collapse and the behavior of temporal spatial distribution of micro seismicities. The stress, strain and the temporal spatial distribution of micro seismicities vividly depict the phenomena of localization, temporal transitions, dilatation or rise, elastic rebound and conjugate (X type) deformation zone.
