For the determination of the smoothing factor (also known as the regularization parameter) in the co-seismic slip distribution inversion, the compromise curve between the model roughness and the data fitting residual ...For the determination of the smoothing factor (also known as the regularization parameter) in the co-seismic slip distribution inversion, the compromise curve between the model roughness and the data fitting residual is generally used to determine (in order to distinguish the method proposed in this paper, the method is called “L curve” according to its shape). Based on the L-curve, the Eclectic Intersection curve as a new method is proposed to determine the smoothing factor in this paper. The results of the simulated experiment show that the inversion accuracy of the parameters of the seismic slip distribution with the smoothing factor determined by the Eclectic Intersection curve method is better than that of the L curve method. Moreover, the Eclectic Intersection curve method and the L curve method are used to determine the smoothing factor of L’Aquila earthquake and the Taiwan Meinong earthquake slip distribution inversion respectively, and the inversion results are compared and analyzed. The analysis results show that the L’Aquila and the Taiwan Meinong actual earthquake slip distribution results are in the range of other scholars at home and abroad, and compared with the L curve method, the Eclectic Intersection curve method has advantages of high computation efficiency, no need to depend on data fitting degree and more appropriate of smoothing factor and so on.展开更多
Data obtained by GRACE(Gravity Recovery and Climate Experiment) have been used to invert for the seismic source parameters of megathrust earthquakes under the assumption of either uniform slip over an entire fault or ...Data obtained by GRACE(Gravity Recovery and Climate Experiment) have been used to invert for the seismic source parameters of megathrust earthquakes under the assumption of either uniform slip over an entire fault or a point-like seismic source.Herein, we further extend the inversion of GRACE long-wavelength gravity changes to heterogeneous slip distributions during the 2011 Tohoku earthquake using three fault models:(Ⅰ) a constant-strike and constant-dip fault,(Ⅱ) a variable dip fault, and(Ⅲ) a realistically varying strike fault. By removing the post-seismic signal from the time series, and taking the effect of ocean water redistribution into account, we invert for slip models I, II, and III using co-seismic gravity changes measured by GRACE, de-striped by DDK3 decorrelation filter. The total seismic moments of our slip models, with respective values of 4.9×10^(22) Nm, 5.1×10^(22) Nm, and 5.0×10^(22) Nm, are smaller than those obtained by other studies relying on GRACE data. The resulting centroids are also located at greater depths(20 km, 19.8 km,and 17.4 km, respectively). By combining onshore GPS, GPS-Acoustic, and GRACE data, we obtain a jointly inverted slip model with a seismic moment of 4.8×10^(22) Nm, which is larger than the seismic moment obtained using only the GPS displacements. We show that the slip inverted from low degree space-borne gravimetric data, which contains information at the ocean region, is affected by the strike of the arcuate trench. The space-borne gravimetric data help us constrain the source parameters of a megathrust earthquake within the frame of heterogeneous slip models.展开更多
The equivalence of geophysical fields, the finiteness of measurements and the measurement errors make the result of geophysical inversion non-unique. For example, the measurements and inversion method used, the priori...The equivalence of geophysical fields, the finiteness of measurements and the measurement errors make the result of geophysical inversion non-unique. For example, the measurements and inversion method used, the priori rupture model determined and the slip distribution smoothing factor selected will have significant influences on the earthquake rupture slip distribution. Using different data and methods, different authors have given different rupture slip distribution models of the 2015 Mw7.9 Nepal earth- quake, with the maximum slip ranging from 3.0 m to 6.8 m. In this paper, geometry parameters of the single rectangular fault model in elastic half-space were inferred constraining with the Global Posi- tioning System （GPS） and Interferometric Synthetic Aperture Radar （InSAR） coseismic deformations and bounding the slip with approximate average value; and then, the single rectangular fault was divided into multiple sub-faults, and the final slip smoothing factor, the final slip distribution and the maximum slip were determined with the misfit-roughness tradeoff curve, the cross-validation sum of squares （CVSS） and the third-party observation data or indexes being comprehensively taken into account. The results show that, the rupture of the Nepal earthquake extended by over 100 km east by south. The maximum slip of the earthquake was about 6.5-6.7 m, and most of the slip is confined at depths of 8 -20 kin, consistent with the depth distribution of aftershocks. The method for reducing the multiplicity of solutions to rupture slip distribution in this paper was ever used in inversion of rupture slip distri- bution for the 2008 Wenchuan and 2013 Lushan earthquakes, and the third-party measurement - surface dislocation has very large effect on reducing the multiplicity of solutions to inversion of the Wenchuan earthquake. Other priori information or indicators, such as fault strike, dip, earthquake magnitude, seismic activity, Coulomb stress, and seismic period, can be used for beneficial validation of and comparison with inversion results.展开更多
We analyze co-seismic displacement field of the 26 December 2004, giant Sumatra–Andaman earthquake derived from Global Position System observations,geological vertical measurement of coral head, and pivot line observ...We analyze co-seismic displacement field of the 26 December 2004, giant Sumatra–Andaman earthquake derived from Global Position System observations,geological vertical measurement of coral head, and pivot line observed through remote sensing. Using the co-seismic displacement field and AK135 spherical layered Earth model, we invert co-seismic slip distribution along the seismic fault. We also search the best fault geometry model to fit the observed data. Assuming that the dip angle linearly increases in downward direction, the postfit residual variation of the inversed geometry model with dip angles linearly changing along fault strike are plotted. The geometry model with local minimum misfits is the one with dip angle linearly increasing along strike from 4.3oin top southernmost patch to 4.5oin top northernmost path and dip angle linearly increased. By using the fault shape and geodetic co-seismic data, we estimate the slip distribution on the curved fault. Our result shows that the earthquake ruptured ＊200-km width down to a depth of about 60 km.0.5–12.5 m of thrust slip is resolved with the largest slip centered around the central section of the rupture zone78N–108N in latitude. The estimated seismic moment is8.2 9 1022 N m, which is larger than estimation from the centroid moment magnitude（4.0 9 1022 N m）, and smaller than estimation from normal-mode oscillation data modeling（1.0 9 1023 N m）.展开更多
The authors analyze co_seismic displacement field derived from the Global Position System (GPS) observations collected before and after the 2001 Kokoxili earthquake, western China. Using the co_seismic displacement da...The authors analyze co_seismic displacement field derived from the Global Position System (GPS) observations collected before and after the 2001 Kokoxili earthquake, western China. Using the co_seismic displacement data, and constrained with surface rupture data, they invert co_seismic slip distribution along the seismic fault. Their result shows that the earthquake ruptured the upper crust down to a depth of 13.1～22km (at 70% certainty), with its optimal estimate at 16.5km. A 2～3m left_lateral strike slip is resolved between the Sun Lake segment and the west end of the main rupture zone, although surface rupture has not been observed there. The surface rupture of this earthquake is ended at the Sun Lake to the west, but left_lateral slip of 1.5～2.0m seems to exist beyond the east end of surface rupture observed from field geology. Seismic moment release estimated using GPS and surface rupture measurement is 6.0×10 20 N·m, which is in good agreement with the result obtained from seismic wave inversion.展开更多
基金National Natural Science Foundation of China(Nos.4187400141664001)+1 种基金Support Program for Outstanding Youth Talents in Jiangxi Province(No.20162BCB23050)National Key Research and Development Program(No.2016YFB0501405)。
文摘For the determination of the smoothing factor (also known as the regularization parameter) in the co-seismic slip distribution inversion, the compromise curve between the model roughness and the data fitting residual is generally used to determine (in order to distinguish the method proposed in this paper, the method is called “L curve” according to its shape). Based on the L-curve, the Eclectic Intersection curve as a new method is proposed to determine the smoothing factor in this paper. The results of the simulated experiment show that the inversion accuracy of the parameters of the seismic slip distribution with the smoothing factor determined by the Eclectic Intersection curve method is better than that of the L curve method. Moreover, the Eclectic Intersection curve method and the L curve method are used to determine the smoothing factor of L’Aquila earthquake and the Taiwan Meinong earthquake slip distribution inversion respectively, and the inversion results are compared and analyzed. The analysis results show that the L’Aquila and the Taiwan Meinong actual earthquake slip distribution results are in the range of other scholars at home and abroad, and compared with the L curve method, the Eclectic Intersection curve method has advantages of high computation efficiency, no need to depend on data fitting degree and more appropriate of smoothing factor and so on.
基金supported financially by the National Natural Science Foundation of China (No.41574021,41474059,41331066,41774088,and 41174063)a research grant from the Institute of Crustal Dynamics,China Earthquake Administration (No.ZDJ2017-23)+4 种基金the CAS/CAFEA International Partnership Program for Creative Research Teams (No.KZZD-EW-TZ-19)the Key Research Program of Frontier Sciences CAS (Chinese Academy of Sciences) (QYZDY-SSW-SYS003)the SKLGED foundation (SKLGED2014-1-1-E)the GOCE Italy Project (the Italian Space Agency and the European Space Agency Endorsement)the China Postdoctoral Science Foundation (No.133014)
文摘Data obtained by GRACE(Gravity Recovery and Climate Experiment) have been used to invert for the seismic source parameters of megathrust earthquakes under the assumption of either uniform slip over an entire fault or a point-like seismic source.Herein, we further extend the inversion of GRACE long-wavelength gravity changes to heterogeneous slip distributions during the 2011 Tohoku earthquake using three fault models:(Ⅰ) a constant-strike and constant-dip fault,(Ⅱ) a variable dip fault, and(Ⅲ) a realistically varying strike fault. By removing the post-seismic signal from the time series, and taking the effect of ocean water redistribution into account, we invert for slip models I, II, and III using co-seismic gravity changes measured by GRACE, de-striped by DDK3 decorrelation filter. The total seismic moments of our slip models, with respective values of 4.9×10^(22) Nm, 5.1×10^(22) Nm, and 5.0×10^(22) Nm, are smaller than those obtained by other studies relying on GRACE data. The resulting centroids are also located at greater depths(20 km, 19.8 km,and 17.4 km, respectively). By combining onshore GPS, GPS-Acoustic, and GRACE data, we obtain a jointly inverted slip model with a seismic moment of 4.8×10^(22) Nm, which is larger than the seismic moment obtained using only the GPS displacements. We show that the slip inverted from low degree space-borne gravimetric data, which contains information at the ocean region, is affected by the strike of the arcuate trench. The space-borne gravimetric data help us constrain the source parameters of a megathrust earthquake within the frame of heterogeneous slip models.
基金supported by the Director Foundation of Institute of Seismology,China Earthquake Adminstration(IS201506220)the National Natural Science Foundation of China(40974012,41304019)the Special Foundation for Seismic Research(201208006)
文摘The equivalence of geophysical fields, the finiteness of measurements and the measurement errors make the result of geophysical inversion non-unique. For example, the measurements and inversion method used, the priori rupture model determined and the slip distribution smoothing factor selected will have significant influences on the earthquake rupture slip distribution. Using different data and methods, different authors have given different rupture slip distribution models of the 2015 Mw7.9 Nepal earth- quake, with the maximum slip ranging from 3.0 m to 6.8 m. In this paper, geometry parameters of the single rectangular fault model in elastic half-space were inferred constraining with the Global Posi- tioning System （GPS） and Interferometric Synthetic Aperture Radar （InSAR） coseismic deformations and bounding the slip with approximate average value; and then, the single rectangular fault was divided into multiple sub-faults, and the final slip smoothing factor, the final slip distribution and the maximum slip were determined with the misfit-roughness tradeoff curve, the cross-validation sum of squares （CVSS） and the third-party observation data or indexes being comprehensively taken into account. The results show that, the rupture of the Nepal earthquake extended by over 100 km east by south. The maximum slip of the earthquake was about 6.5-6.7 m, and most of the slip is confined at depths of 8 -20 kin, consistent with the depth distribution of aftershocks. The method for reducing the multiplicity of solutions to rupture slip distribution in this paper was ever used in inversion of rupture slip distri- bution for the 2008 Wenchuan and 2013 Lushan earthquakes, and the third-party measurement - surface dislocation has very large effect on reducing the multiplicity of solutions to inversion of the Wenchuan earthquake. Other priori information or indicators, such as fault strike, dip, earthquake magnitude, seismic activity, Coulomb stress, and seismic period, can be used for beneficial validation of and comparison with inversion results.
基金supported by the Special Fund of Fundamental Scientific Research Business Expense for Higher School of Central Government(Projects for creation teams ZY20110101)NSFC 41090294talent selection and training plan project of Hebei university
文摘We analyze co-seismic displacement field of the 26 December 2004, giant Sumatra–Andaman earthquake derived from Global Position System observations,geological vertical measurement of coral head, and pivot line observed through remote sensing. Using the co-seismic displacement field and AK135 spherical layered Earth model, we invert co-seismic slip distribution along the seismic fault. We also search the best fault geometry model to fit the observed data. Assuming that the dip angle linearly increases in downward direction, the postfit residual variation of the inversed geometry model with dip angles linearly changing along fault strike are plotted. The geometry model with local minimum misfits is the one with dip angle linearly increasing along strike from 4.3oin top southernmost patch to 4.5oin top northernmost path and dip angle linearly increased. By using the fault shape and geodetic co-seismic data, we estimate the slip distribution on the curved fault. Our result shows that the earthquake ruptured ＊200-km width down to a depth of about 60 km.0.5–12.5 m of thrust slip is resolved with the largest slip centered around the central section of the rupture zone78N–108N in latitude. The estimated seismic moment is8.2 9 1022 N m, which is larger than estimation from the centroid moment magnitude（4.0 9 1022 N m）, and smaller than estimation from normal-mode oscillation data modeling（1.0 9 1023 N m）.
文摘The authors analyze co_seismic displacement field derived from the Global Position System (GPS) observations collected before and after the 2001 Kokoxili earthquake, western China. Using the co_seismic displacement data, and constrained with surface rupture data, they invert co_seismic slip distribution along the seismic fault. Their result shows that the earthquake ruptured the upper crust down to a depth of 13.1～22km (at 70% certainty), with its optimal estimate at 16.5km. A 2～3m left_lateral strike slip is resolved between the Sun Lake segment and the west end of the main rupture zone, although surface rupture has not been observed there. The surface rupture of this earthquake is ended at the Sun Lake to the west, but left_lateral slip of 1.5～2.0m seems to exist beyond the east end of surface rupture observed from field geology. Seismic moment release estimated using GPS and surface rupture measurement is 6.0×10 20 N·m, which is in good agreement with the result obtained from seismic wave inversion.
