Since the non-unique solution exists in the inversion for finite-fault rupture history, the random weighting method hasbeen used to estimate error of the inversion results in this paper. The resolution distributions o...Since the non-unique solution exists in the inversion for finite-fault rupture history, the random weighting method hasbeen used to estimate error of the inversion results in this paper. The resolution distributions of slip amplitude, rake,rupture time and rise time on the finite fault were deduced quantitatively by model calculation. By using the randomweighting method, the inversion results of Taiwan Strait earthquake and Myanmar-China boundal earthquake showthat the parameters related to the rupture centers of two events have the highest resolution, and the solutinn are the mostreliable(otherwise the resolution of the slip amplitudes and rise time on the finite-fault boundary is low.展开更多
This paper derives from the representation theory the formula for calculating the radiation excited by heterogeneous fault rupture based on box-like discretization scheme. Preliminary validation indicates that our alg...This paper derives from the representation theory the formula for calculating the radiation excited by heterogeneous fault rupture based on box-like discretization scheme. Preliminary validation indicates that our algorithm has very high computation precision and efficiency; therefore, it is a very practical tool to investigate strong ground motion problems. Additionally, the equations given in this study can also be used to invert the fault rupture process.展开更多
A numerical model for generating normal fault earthquakes in the Qinghai-Xizang Plateau′S upper crust is constructed with 3-D elasto-viscous finite element method. Based on the numerical simulation calculation,some c...A numerical model for generating normal fault earthquakes in the Qinghai-Xizang Plateau′S upper crust is constructed with 3-D elasto-viscous finite element method. Based on the numerical simulation calculation,some conclusions were got:If the effective viscosity of the upper crust material is less than that of lower strata of the crust in the Qinghai-Xizang Plateau, even under the strong push of India continent,the stress state of the upper crust can still be extensional in south part of the Qinghai-Xizang Plateau.Numerical simulations show that the stress state changes with the depth of the lithosphere,from extensional stress state in upper crust to compressive in the lower part.Extensional stress state may exist mainly in the upper crust of the south part of the Qinghai-Xizang Plateau.展开更多
We use interferometric synthetic aperture radar (InSAR) and broadband seismic waveform data to estimate a source model of the 11th July, 2004 M W 6.2 Zhongba earthquake, Tibet of China. This event occurred within th...We use interferometric synthetic aperture radar (InSAR) and broadband seismic waveform data to estimate a source model of the 11th July, 2004 M W 6.2 Zhongba earthquake, Tibet of China. This event occurred within the seismically active zone of southwestern Tibetan Plateau where the east-west extension of the upper crust is observed. Because of limitations in one pair of InSAR data available, there are trade-offs among centroid depth, rupture area and amount of slip. Available seismic data tightly constrain the focal mechanism and centroid depth of the earthquake but not the horizontal location. Together, two complementary data sets can be used to identify the actual fault plane, better constrain the slip model and event location. We first use regional seismic waveform to estimate point source mechanism, then InSAR data is used to obtain better location. Finally, a joint inversion of teleseismic P-waves and InSAR data is performed to obtain a distributed model. Our preferred point source mechanism indicates a seismic moment of ~2.2×10 18 N·m (~M W 6.2), a fault plane solution of 171° (342 ° )/42 ° (48 ° )/-83 ° (-97 ° ), corresponding to strike/dip/rake, and a depth of 11 km. The fault plane with strike of 171 ? and dip of 42° is identified as the ruptured fault with the aid of InSAR data. The preferred source model features compact area of slips between depth of 5–11 km and 10 km along strike with maximum slip amplitude of about 1.5 m.展开更多
The Wenchuan Ms8.0 earthquake and the Lushan M7.0 earthquake occurred in the north and south segments of the Longmenshan nappe tectonic belt, respectively. Based on the focal mechanism and finite fault model of the We...The Wenchuan Ms8.0 earthquake and the Lushan M7.0 earthquake occurred in the north and south segments of the Longmenshan nappe tectonic belt, respectively. Based on the focal mechanism and finite fault model of the Wenchuan Ms8.0 earthquake, we calculated the coulomb failure stress change. The inverted coulomb stress changes based on the Nishimura and Chenji models both show that the Lushan MT. 0 earth- quake occurred in the increased area of coulomb failure stress induced by the Wenchuan Ms8. 0 earthquake. The coulomb failure stress increased by approximately 0. 135 - 0. 152 bar in the source of the Lushan M7.0 earthquake, which is far more than the stress triggering threshold. Therefore, the Lushan M7.0 earthquake was most likely triggered by the coulomb failure stress change.展开更多
The 2016 A/w7.8 Kaikoura(New Zealand)earthquake was the most complex event ever instrumentally recorded and geologically investigated,as it ruptured on more than 12 fault segments of various geometries.To study the ma...The 2016 A/w7.8 Kaikoura(New Zealand)earthquake was the most complex event ever instrumentally recorded and geologically investigated,as it ruptured on more than 12 fault segments of various geometries.To study the mainshock rupture characteristics,geodetic methods like InSAR and GPS play an essential role in providing satisfactory spatial resolution.However,early strong aftershocks may cause extra ground deformation which bias the mainshock rupture inversion result.In this paper,we will focus on studying the Mw 6.3 aftershock,which is the only A/6+thrust slip aftershock that occurred only 30 minutes after the Kaikoura mainshock.We will relocate the hypocenter of this event using the hypo 2000 method,make the finite fault model(FFM)inversion for the detailed rupture processes and calculate the synthetic surface displacement to compare with the observed GPS data and figure out its influence on the mainshock study.Although we are not able to resolve the real ruptured fault of this event because of limited observation data,we infer that it is a west-ward dipping event of oblique slip mechanism,consistent with the subfault geometries of the Kaikoura mainshock.According to the inverted FFM,this event can generate 10-20 cm ground surface displacement and affect the ground displacement observation at nearby GPS stations.展开更多
Inversion for the seismic fault rupture history is an important way to study the nature of the earthquake source. Inthis paper, we have selected two Taiwan earthquakes that occurred closely in time and located in the ...Inversion for the seismic fault rupture history is an important way to study the nature of the earthquake source. Inthis paper, we have selected two Taiwan earthquakes that occurred closely in time and located in the same region,inversed the distribution of the slip amplitudes, rakes, risetimes and the rupture times on the fault planes by usingGDSN broad-band and long-period records and the adaptive hybrid global search algorithm, and compared the twoevents. The slip rate of every subfault calculated provides information about the distribution of tectonic stress andfault strength. To the former event (Ms=6.0), the maximum slip amplitude 2.4 m and the minimum risetime 1.2 sare both located at the hypocentre. The latter earthquake (Ms=6.6) consisted of two subevents and the second source has 4 s delay. The maximum slip amplitUde 0.9 m located near hypocentre is corresponding to the minimumrisetime l.4 s, and the corresponding maximum slip rate 0.7 m.s~-1 is similar to the peak value of other large sliprate areas. We consider that the latter event has more complicated temporal-spatial distribution than the former.展开更多
Based on dislocation theory of 0kada, we adopted a finite-element fault model inverted by Gavin Hayes from seismic data for the 2011 Japan Mw9.0 earthquake, and obtained the corresponding surface displacement and stra...Based on dislocation theory of 0kada, we adopted a finite-element fault model inverted by Gavin Hayes from seismic data for the 2011 Japan Mw9.0 earthquake, and obtained the corresponding surface displacement and strain fields. The calculated displacement field is consistent with the observed GPS results in the trend of changes. Also the surface displacement and strain fields both show large variations in space.展开更多
The 2018 M_(W)6.7 Iburi earthquake shocked the eastern Iburi region to the west of the Hidaka Collision Zone in Hokkaido,which is a destructive inland earthquake.We resolved the kinematic rupture process of the event ...The 2018 M_(W)6.7 Iburi earthquake shocked the eastern Iburi region to the west of the Hidaka Collision Zone in Hokkaido,which is a destructive inland earthquake.We resolved the kinematic rupture process of the event by combining strong motions(SM)and synthetic aperture radar(SAR)images in a joint inversion.The results reveal that the duration of the whole rupture is about 17s,yielding a total seismic moment of 1.4×10^(19)N·m(M_(W)=6.7).The main slip area is located at a depth of approximately 24 km with a peak slip of~0.8m above the hypocenter.The comparison with the regional velocity model shows the earthquake was initiated in the upper mantle,while the majority of slips are located in the lower crust,which is an“aseismic”domain in the typical sandwich model.The location of the major slip area is consistent with a high-conductivity volume.We proposed a mechanism of low frictional property(<0.3)produced by high pore pressure to explain the abnormal high dip angle and centroid depth located in the ductile lower-crust.Aftershocks are distributed in areas where the Coulomb frictional stress increases due to co-seismic displacement with a mechanism conjugating to the mainshock.展开更多
On December 18,2023,the MS 6.2 Jishishan earthquake occurred in the northeastern region of the Qinghai-Xizang Plateau,causing heavy casualties and property damage in Gansu and Qinghai Provinces.In this study,we integr...On December 18,2023,the MS 6.2 Jishishan earthquake occurred in the northeastern region of the Qinghai-Xizang Plateau,causing heavy casualties and property damage in Gansu and Qinghai Provinces.In this study,we integrate space imaging geodesy,finite fault inversion,and back-projection methods to decipher its rupture property,including fault geometry,coseismic slip distribution,rupture direction,and propagation speed.The results reveal that the seismogenic fault dips to the southwest at an angle of 29.The major slip asperity is dominated by reverse slip and is concentrated within a depth range of 7-16 km,which explains the significant uplift near the epicenter observed by both the Sentinel-1 ascending and descending InSAR data.Moreover,the teleseismic array waveforms indicate a northwest propagating rupture with an overall slow rupture velocity of~1.91 km/s(AK array)or 1.01 km/s(AU array).展开更多
文摘Since the non-unique solution exists in the inversion for finite-fault rupture history, the random weighting method hasbeen used to estimate error of the inversion results in this paper. The resolution distributions of slip amplitude, rake,rupture time and rise time on the finite fault were deduced quantitatively by model calculation. By using the randomweighting method, the inversion results of Taiwan Strait earthquake and Myanmar-China boundal earthquake showthat the parameters related to the rupture centers of two events have the highest resolution, and the solutinn are the mostreliable(otherwise the resolution of the slip amplitudes and rise time on the finite-fault boundary is low.
基金National Natural Science Foundation of China (40474011 and 40521002).
文摘This paper derives from the representation theory the formula for calculating the radiation excited by heterogeneous fault rupture based on box-like discretization scheme. Preliminary validation indicates that our algorithm has very high computation precision and efficiency; therefore, it is a very practical tool to investigate strong ground motion problems. Additionally, the equations given in this study can also be used to invert the fault rupture process.
文摘A numerical model for generating normal fault earthquakes in the Qinghai-Xizang Plateau′S upper crust is constructed with 3-D elasto-viscous finite element method. Based on the numerical simulation calculation,some conclusions were got:If the effective viscosity of the upper crust material is less than that of lower strata of the crust in the Qinghai-Xizang Plateau, even under the strong push of India continent,the stress state of the upper crust can still be extensional in south part of the Qinghai-Xizang Plateau.Numerical simulations show that the stress state changes with the depth of the lithosphere,from extensional stress state in upper crust to compressive in the lower part.Extensional stress state may exist mainly in the upper crust of the south part of the Qinghai-Xizang Plateau.
基金study was supported jointly by National Natural Science Foundation of China (Nos.40821160549 and 41074032)CAS Knowledge Innovation Program (No. KZCX2-YW-116-1)Joint Seismological Science Fundation of China(Nos.200808078 and 200708035)
文摘We use interferometric synthetic aperture radar (InSAR) and broadband seismic waveform data to estimate a source model of the 11th July, 2004 M W 6.2 Zhongba earthquake, Tibet of China. This event occurred within the seismically active zone of southwestern Tibetan Plateau where the east-west extension of the upper crust is observed. Because of limitations in one pair of InSAR data available, there are trade-offs among centroid depth, rupture area and amount of slip. Available seismic data tightly constrain the focal mechanism and centroid depth of the earthquake but not the horizontal location. Together, two complementary data sets can be used to identify the actual fault plane, better constrain the slip model and event location. We first use regional seismic waveform to estimate point source mechanism, then InSAR data is used to obtain better location. Finally, a joint inversion of teleseismic P-waves and InSAR data is performed to obtain a distributed model. Our preferred point source mechanism indicates a seismic moment of ~2.2×10 18 N·m (~M W 6.2), a fault plane solution of 171° (342 ° )/42 ° (48 ° )/-83 ° (-97 ° ), corresponding to strike/dip/rake, and a depth of 11 km. The fault plane with strike of 171 ? and dip of 42° is identified as the ruptured fault with the aid of InSAR data. The preferred source model features compact area of slips between depth of 5–11 km and 10 km along strike with maximum slip amplitude of about 1.5 m.
基金supported by the Director Foundation of Institute of SeismologyChina Earthquake Administration(201056088)
文摘The Wenchuan Ms8.0 earthquake and the Lushan M7.0 earthquake occurred in the north and south segments of the Longmenshan nappe tectonic belt, respectively. Based on the focal mechanism and finite fault model of the Wenchuan Ms8.0 earthquake, we calculated the coulomb failure stress change. The inverted coulomb stress changes based on the Nishimura and Chenji models both show that the Lushan MT. 0 earth- quake occurred in the increased area of coulomb failure stress induced by the Wenchuan Ms8. 0 earthquake. The coulomb failure stress increased by approximately 0. 135 - 0. 152 bar in the source of the Lushan M7.0 earthquake, which is far more than the stress triggering threshold. Therefore, the Lushan M7.0 earthquake was most likely triggered by the coulomb failure stress change.
基金the New Zealand GeoNet (http://geonet.org. nz) projectfinancially supported by National Natural Science Foundation of China (Nos. 41590854 and 41461164003)
文摘The 2016 A/w7.8 Kaikoura(New Zealand)earthquake was the most complex event ever instrumentally recorded and geologically investigated,as it ruptured on more than 12 fault segments of various geometries.To study the mainshock rupture characteristics,geodetic methods like InSAR and GPS play an essential role in providing satisfactory spatial resolution.However,early strong aftershocks may cause extra ground deformation which bias the mainshock rupture inversion result.In this paper,we will focus on studying the Mw 6.3 aftershock,which is the only A/6+thrust slip aftershock that occurred only 30 minutes after the Kaikoura mainshock.We will relocate the hypocenter of this event using the hypo 2000 method,make the finite fault model(FFM)inversion for the detailed rupture processes and calculate the synthetic surface displacement to compare with the observed GPS data and figure out its influence on the mainshock study.Although we are not able to resolve the real ruptured fault of this event because of limited observation data,we infer that it is a west-ward dipping event of oblique slip mechanism,consistent with the subfault geometries of the Kaikoura mainshock.According to the inverted FFM,this event can generate 10-20 cm ground surface displacement and affect the ground displacement observation at nearby GPS stations.
文摘Inversion for the seismic fault rupture history is an important way to study the nature of the earthquake source. Inthis paper, we have selected two Taiwan earthquakes that occurred closely in time and located in the same region,inversed the distribution of the slip amplitudes, rakes, risetimes and the rupture times on the fault planes by usingGDSN broad-band and long-period records and the adaptive hybrid global search algorithm, and compared the twoevents. The slip rate of every subfault calculated provides information about the distribution of tectonic stress andfault strength. To the former event (Ms=6.0), the maximum slip amplitude 2.4 m and the minimum risetime 1.2 sare both located at the hypocentre. The latter earthquake (Ms=6.6) consisted of two subevents and the second source has 4 s delay. The maximum slip amplitUde 0.9 m located near hypocentre is corresponding to the minimumrisetime l.4 s, and the corresponding maximum slip rate 0.7 m.s~-1 is similar to the peak value of other large sliprate areas. We consider that the latter event has more complicated temporal-spatial distribution than the former.
基金supported by the Director Foundation of Institute of Seismology,CEA( IS201056088)
文摘Based on dislocation theory of 0kada, we adopted a finite-element fault model inverted by Gavin Hayes from seismic data for the 2011 Japan Mw9.0 earthquake, and obtained the corresponding surface displacement and strain fields. The calculated displacement field is consistent with the observed GPS results in the trend of changes. Also the surface displacement and strain fields both show large variations in space.
基金This work is supported by the National Key R&D Program of China(No.2018YFC1504203)the National Natural Science Foundation of China(No.42021003).
文摘The 2018 M_(W)6.7 Iburi earthquake shocked the eastern Iburi region to the west of the Hidaka Collision Zone in Hokkaido,which is a destructive inland earthquake.We resolved the kinematic rupture process of the event by combining strong motions(SM)and synthetic aperture radar(SAR)images in a joint inversion.The results reveal that the duration of the whole rupture is about 17s,yielding a total seismic moment of 1.4×10^(19)N·m(M_(W)=6.7).The main slip area is located at a depth of approximately 24 km with a peak slip of~0.8m above the hypocenter.The comparison with the regional velocity model shows the earthquake was initiated in the upper mantle,while the majority of slips are located in the lower crust,which is an“aseismic”domain in the typical sandwich model.The location of the major slip area is consistent with a high-conductivity volume.We proposed a mechanism of low frictional property(<0.3)produced by high pore pressure to explain the abnormal high dip angle and centroid depth located in the ductile lower-crust.Aftershocks are distributed in areas where the Coulomb frictional stress increases due to co-seismic displacement with a mechanism conjugating to the mainshock.
基金supported by the Open Fund of Hubei Luojia Laboratory(230100015)the Strategic Priority Research Program of the Chinese Academy of Sciences(XDB41000000)the Knowledge Innovation Program of Wuhan-Shuguang Project(2023010201020281).
文摘On December 18,2023,the MS 6.2 Jishishan earthquake occurred in the northeastern region of the Qinghai-Xizang Plateau,causing heavy casualties and property damage in Gansu and Qinghai Provinces.In this study,we integrate space imaging geodesy,finite fault inversion,and back-projection methods to decipher its rupture property,including fault geometry,coseismic slip distribution,rupture direction,and propagation speed.The results reveal that the seismogenic fault dips to the southwest at an angle of 29.The major slip asperity is dominated by reverse slip and is concentrated within a depth range of 7-16 km,which explains the significant uplift near the epicenter observed by both the Sentinel-1 ascending and descending InSAR data.Moreover,the teleseismic array waveforms indicate a northwest propagating rupture with an overall slow rupture velocity of~1.91 km/s(AK array)or 1.01 km/s(AU array).