The Wenchuan earthquake of 12 May 2008 is the most destructive earthquake in China in the past 30 years in terms of property damage and human losses. In order to understand the earthquake process and the geo-morpholog...The Wenchuan earthquake of 12 May 2008 is the most destructive earthquake in China in the past 30 years in terms of property damage and human losses. In order to understand the earthquake process and the geo-morphological factors affecting the seismic hazard, we simulated the strong ground mo-tion caused by the earthquake, incorporating three-dimensional (3D) earth structure, finite-fault rupture, and realistic surface topography. The simulated ground motions reveal that the fault rupture and basin structure control the overall pattern of the peak ground shaking. Large peak ground velocity (PGV) is distributed in two narrow areas: one with the largest PGV values is above the hanging wall of the fault and attributed to the locations of fault asperities and rupture directivity; the other is along the north-western margin of the Sichuan Basin and caused by both the directivity of fault rupture and the ampli-fication in the thick sediment basin. Rough topography above the rupture fault causes wave scattering, resulting in significantly larger peak ground motion on the apex of topographic relief than in the valley. Topography and scattering also reduce the wave energy in the forward direction of fault rupture but increase the PGV in other parts of the basin. These results suggest the need for a localized hazard as-sessment in places of rough topography that takes the topographic effects into account. Finally, had the earthquake started at the northeast end of the fault zone and ruptured to the southwest, Chengdu would have suffered a much stronger shaking than it experienced on 12 May, 2008.展开更多
This paper reports a series solution of wave functions for two-dimensional scattering and diffraction of plane SH waves induced by a symmetrical V-shaped canyon with different shape ratios. A half-space with a symmetr...This paper reports a series solution of wave functions for two-dimensional scattering and diffraction of plane SH waves induced by a symmetrical V-shaped canyon with different shape ratios. A half-space with a symmetrical V-shaped canyon is divided into two sub-regions by using a circular-arc auxiliary boundary. The two sub-regions are represented by global and local cylindrical coordinate systems, respectively. In each coordinate system, the wave field satisfying the Helmholtz equation is represented by the separation of variables method, in terms of the series of both Bessel functions and Hankel functions with unknown complex coefficients. Then, the two wave fields are described in the local coordinate system using the Graf addition theorem. Finally, the unknown coefficients are sought by satisfying the continuity conditions of the auxiliary boundary. To consider the phase characteristics of the wave scattering, a parametric analysis is carried out in the time domain by assuming an incident signal of the Ricker type. Surface and subsurface transient responses demonstrate the characteristics and mechanisms of wave propagating and scattering.展开更多
A linear response history analysis method is used to determine the influence of three factors:geometric incoherency,wave-passage,and local site characteristics on the response of lnulti-support structures subjected to...A linear response history analysis method is used to determine the influence of three factors:geometric incoherency,wave-passage,and local site characteristics on the response of lnulti-support structures subjected to differential ground motions.A one-span frame and a reduced model of a 24-span bridge,located in Las Vegas,Nevada are studied,in which the influence of each of the three factors and their combinations are analyzed.It is revealed that the incoherency of earthquake ground motion can have a dramatic influence on structural response by modifying the dynamics response to uniform excitation and inducing pseudo-static response,which does not exist in structures subjected to uniform excitation.The total response when all three sources of ground motion incoherency are included is generally larger than that of uniform excitation.展开更多
The scattering and reflection of SH waves by a slope on an elastic wedged space is investigated. A series solution is obtained by using the wave function expansion method. The slope on a wedged space is divided into t...The scattering and reflection of SH waves by a slope on an elastic wedged space is investigated. A series solution is obtained by using the wave function expansion method. The slope on a wedged space is divided into two subregions by an artificial, auxiliary circular arc. The wave fields with unknown complex coefficients within each sub-region are derived. Applying Graf addition theorem, the scattered waves in the sub-regions are expressed in a global coordinate system. Fourier transform is adopted to derive a consistent form of standing waves in the inner region using the orthogonality of the cosine functions. The boundary-valued problem is solved by stress and displacement continuity along the artificial, auxiliary arc to obtain the unknown complex coefficients. Parametric studies are next performed to investigate how the topography from the slope on the wedged space will affect the scattering and diffraction, and hence the amplification and de-amplification of the SH waves. Numerical results show that the surface motions on the slope of the wedged space is influenced greatly by the topography. Amplification of the surface motions near the slope vertex is significant. The corresponding phases along the wedged space surfaces are consistent with the direction that the SH waves are propagating.展开更多
High-quality datasets are critical for the development of advanced machine-learning algorithms in seismology.Here,we present an earthquake dataset based on the ChinArray Phase I records(X1).ChinArray Phase I was deplo...High-quality datasets are critical for the development of advanced machine-learning algorithms in seismology.Here,we present an earthquake dataset based on the ChinArray Phase I records(X1).ChinArray Phase I was deployed in the southern north-south seismic zone(20°N-32°N,95°E-110°E)in 2011-2013 using 355 portable broadband seismic stations.CREDIT-X1local,the first release of the ChinArray Reference Earthquake Dataset for Innovative Techniques(CREDIT),includes comprehensive information for the 105,455 local events that occurred in the southern north-south seismic zone during array observation,incorporating them into a single HDF5 file.Original 100-Hz sampled three-component waveforms are organized by event for stations within epicenter distances of 1,000 km,and records of≥200 s are included for each waveform.Two types of phase labels are provided.The first includes manually picked labels for 5,999 events with magnitudes≥2.0,providing 66,507 Pg,42,310 Sg,12,823 Pn,and 546 Sn phases.The second contains automatically labeled phases for 105,442 events with magnitudes of−1.6 to 7.6.These phases were picked using a recurrent neural network phase picker and screened using the corresponding travel time curves,resulting in 1,179,808 Pg,884,281 Sg,176,089 Pn,and 22,986 Sn phases.Additionally,first-motion polarities are included for 31,273 Pg phases.The event and station locations are provided,so that deep learning networks for both conventional phase picking and phase association can be trained and validated.The CREDIT-X1local dataset is the first million-scale dataset constructed from a dense seismic array,which is designed to support various multi-station deep-learning methods,high-precision focal mechanism inversion,and seismic tomography studies.Additionally,owing to the high seismicity in the southern north-south seismic zone in China,this dataset has great potential for future scientific discoveries.展开更多
In this paper, mantle circulation flow, continental drift, earthquake origin and other mechanical principles are examined as they apply to earthquake engineering, seismology and dynamics of fluid saturated porous medi...In this paper, mantle circulation flow, continental drift, earthquake origin and other mechanical principles are examined as they apply to earthquake engineering, seismology and dynamics of fluid saturated porous medium. The relationship of mantle flow to earthquakes is examined and clarified, and a new model, different from Haskell’s, is proposed for the earthquake mechanism. The proposed new model is based on the discovery that two pairs of jump stress and jump velocity will start to act from the fault plane. Records obtained directly from recent earthquakes nearby and right on the fault break show a very large velocity impulse, which verify, indirectly, the new mechanism proposed by the author. Further, at least two physical parameters that characterize the seismic intensity must be specified, because according to the discontinuous (jump) wave theory, at the earthquake source, the stress jump and the velocity jump of particle motion should act simultaneously when a sudden break occurs. The third key parameter is shown to be the break (fracture) propagation speed together with the break plane area. This parameter influences the form of the unloading time function at the source. The maximum seismic stress in and displacement of a building are estimated for two unfavorable combinations of the building and its base ground in terms of their relative rigidity. Finally, it is shown that Biot’s theory of wave propagation in fluid saturated porous media is valid only when fluid flow cannot occur.展开更多
The planetary geometries were studied at the moments of 92 strong earthquakes with a magnitude of more than 8 on the Richter scale (R8+) for the period from 1900 to 2011. Three main planetary schemes were specified, n...The planetary geometries were studied at the moments of 92 strong earthquakes with a magnitude of more than 8 on the Richter scale (R8+) for the period from 1900 to 2011. Three main planetary schemes were specified, namely the schemes of the generalized Archimedes lever (gAL) and the Kepler conjunction (gKc), as well as a new geometry of the triangles of the remote signal catcher (cRS). It was discovered 22 gAL, 42 gKc and 28 cRS geometries, which are 23.9%, 45.7%, and 30.4%, respectively, from the total 92 studying cases. It was found that in some earthquakes;the planetary geometries are absolutely identical, which indicates the universality of the mechanism that caused the earthquake. The triggered effect does not depend on the distance between the planets and the mass of the planets, so the mechanism is identified as an inertial gravitational interaction. The triggered effect increases with the multiplicity of the ratio of distances between planets, as well as with pairwise planetary parallelism, which probably indicates about the wave nature of inertial effects. The triggering effect increases with increasing multiplicity of the ratio of distances between planets, as well as with pairwise planetary parallelism, which probably indicates about the wave nature of inertial effects. According to the Archimedes’ lever principle, the Third Law of Motion is changed by adding a few words. It is assumed that inertia is a special case of gravity, namely gravitational self-induction, which really depends, like any self-induction, only on the geometry of the task.展开更多
A number of aftershocks of the May 10th 1997, Zirkuh (Ghaen-Birjand) destructive earthquake have been used to investigate the anisotropy in the upper crust by observing shear wave splitting. Particle motion diagram an...A number of aftershocks of the May 10th 1997, Zirkuh (Ghaen-Birjand) destructive earthquake have been used to investigate the anisotropy in the upper crust by observing shear wave splitting. Particle motion diagram and aspect ratio methods were used as two different approaches to obtain splitting parameters. Clear shear wave splitting was observed on the records of the selected aftershocks, indicating that the media in the region was highly anisotropic. By using particle motion method, the direction of fast shear wave was found 22°N±19°E, while the delay time between the fast and slow shear waves was obtained to be (65±16) ms. By aspect ratio method, the direction of fast shear wave was determined to be 35°N±18°E and the delay time between fast and slow shear waves was found to be (49±10) ms. For a simple horizontal layer with a thickness about 5 km and uniformly distributed anisotropy, a stress aligned cracks model was used and the result was interpreted in terms of vertical aligned cracks in the direction of N22°E, having a density about 0.01. It is assumed that cracks are fluid-filled since they are located in the upper crust. Finally, by using Hudson cracks model for three crack densities 0.005, 0.01, 0.03, the velocity curves of shear wave were plotted as a function of angle between the symmetrical axis of cracks and the azimuth of source to receiver. It was concluded that when shear wave was polarized parallel to the crack surface, the velocity was uniform, but the velocity curve varied clearly if shear wave was polarized perpendicular to the crack surface.展开更多
To reveal the seismogenic mechanism of the Luding earthquake, we employed the 118 China Seismic Network stations to collect the P-wave polarity data from each station, which was then used in the P-wave first motion ap...To reveal the seismogenic mechanism of the Luding earthquake, we employed the 118 China Seismic Network stations to collect the P-wave polarity data from each station, which was then used in the P-wave first motion approach to calculate the focal mechanism solution of the M6.8 Luding earthquake that occurred on September 5,2022. We have also studied the loading effect of tectonic stress on the Luding earthquake fault based on the stress field data for the research area. The results indicate that this earthquake was a strike-slip type, the nodal plane I:strike 167°, dip angle 78°, slip angle 2°;Nodal plane II: strike 77°, dip angle 88°, slip angle 168°. The two fault planes’ instability coefficients of the Luding earthquake are examined considering the region’s background stress field’s condition. The nodal plane I in the Moho circle is discovered to practically coincide with the Coulomb failure line and the tangent point of the Moho circle, indicating that this nodal plane has a high instability coefficient compared to the nodal plane II. The conclusion is that the nodal plane I has a higher likelihood of being the seismogenic fault plane, which is congruent with the seismogenic fault plane suggested by the aftershock distribution, the earthquake radiation energy distribution of a single station, and seismic intensity distribution.The Luding earthquake’s focal mechanism is highly like the theoretical focal mechanism of the fault situated at the location where the Coulomb failure line intersects the Mohr circle, demonstrating that background stress is what caused the earthquake. The substantial fault instability and similarity between the solved and theoretical focal mechanisms make it easier to comprehend the loading effect of tectonic stress on the Luding earthquake fault.展开更多
基金the U.S. National Science Foundation (Grant Nos. EAR 0738779 and OCE 0727919)the National Basic Research Program of China (Grant No. 2004CB418404)partially by the National Nature Science Foundation of China (Grant No. 40521002)
文摘The Wenchuan earthquake of 12 May 2008 is the most destructive earthquake in China in the past 30 years in terms of property damage and human losses. In order to understand the earthquake process and the geo-morphological factors affecting the seismic hazard, we simulated the strong ground mo-tion caused by the earthquake, incorporating three-dimensional (3D) earth structure, finite-fault rupture, and realistic surface topography. The simulated ground motions reveal that the fault rupture and basin structure control the overall pattern of the peak ground shaking. Large peak ground velocity (PGV) is distributed in two narrow areas: one with the largest PGV values is above the hanging wall of the fault and attributed to the locations of fault asperities and rupture directivity; the other is along the north-western margin of the Sichuan Basin and caused by both the directivity of fault rupture and the ampli-fication in the thick sediment basin. Rough topography above the rupture fault causes wave scattering, resulting in significantly larger peak ground motion on the apex of topographic relief than in the valley. Topography and scattering also reduce the wave energy in the forward direction of fault rupture but increase the PGV in other parts of the basin. These results suggest the need for a localized hazard as-sessment in places of rough topography that takes the topographic effects into account. Finally, had the earthquake started at the northeast end of the fault zone and ruptured to the southwest, Chengdu would have suffered a much stronger shaking than it experienced on 12 May, 2008.
基金National Natural Science Foundation of China Under Grant No.51278382
文摘This paper reports a series solution of wave functions for two-dimensional scattering and diffraction of plane SH waves induced by a symmetrical V-shaped canyon with different shape ratios. A half-space with a symmetrical V-shaped canyon is divided into two sub-regions by using a circular-arc auxiliary boundary. The two sub-regions are represented by global and local cylindrical coordinate systems, respectively. In each coordinate system, the wave field satisfying the Helmholtz equation is represented by the separation of variables method, in terms of the series of both Bessel functions and Hankel functions with unknown complex coefficients. Then, the two wave fields are described in the local coordinate system using the Graf addition theorem. Finally, the unknown coefficients are sought by satisfying the continuity conditions of the auxiliary boundary. To consider the phase characteristics of the wave scattering, a parametric analysis is carried out in the time domain by assuming an incident signal of the Ricker type. Surface and subsurface transient responses demonstrate the characteristics and mechanisms of wave propagating and scattering.
基金the China Scholarship Council and the Teaching and Research Award Program for Outstanding Young Teachers (TRAPOYT) in Higher Education Institutions of MOE,PRC.
文摘A linear response history analysis method is used to determine the influence of three factors:geometric incoherency,wave-passage,and local site characteristics on the response of lnulti-support structures subjected to differential ground motions.A one-span frame and a reduced model of a 24-span bridge,located in Las Vegas,Nevada are studied,in which the influence of each of the three factors and their combinations are analyzed.It is revealed that the incoherency of earthquake ground motion can have a dramatic influence on structural response by modifying the dynamics response to uniform excitation and inducing pseudo-static response,which does not exist in structures subjected to uniform excitation.The total response when all three sources of ground motion incoherency are included is generally larger than that of uniform excitation.
文摘The scattering and reflection of SH waves by a slope on an elastic wedged space is investigated. A series solution is obtained by using the wave function expansion method. The slope on a wedged space is divided into two subregions by an artificial, auxiliary circular arc. The wave fields with unknown complex coefficients within each sub-region are derived. Applying Graf addition theorem, the scattered waves in the sub-regions are expressed in a global coordinate system. Fourier transform is adopted to derive a consistent form of standing waves in the inner region using the orthogonality of the cosine functions. The boundary-valued problem is solved by stress and displacement continuity along the artificial, auxiliary arc to obtain the unknown complex coefficients. Parametric studies are next performed to investigate how the topography from the slope on the wedged space will affect the scattering and diffraction, and hence the amplification and de-amplification of the SH waves. Numerical results show that the surface motions on the slope of the wedged space is influenced greatly by the topography. Amplification of the surface motions near the slope vertex is significant. The corresponding phases along the wedged space surfaces are consistent with the direction that the SH waves are propagating.
基金funded by the National Key R&D Program of China (No. 2021YFC3000702)the Special Fund of the Institute of Geophysics, China Earthquake Administration (No. DQJB20B15)+2 种基金the National Natural Science Foundation of China youth Grant (No. 41804059)the Joint Funds of the National Natural Science Foundation of China (No. U223920029)the Science for Earthquake Resilience of China Earthquake Administration (No. XH211103)
文摘High-quality datasets are critical for the development of advanced machine-learning algorithms in seismology.Here,we present an earthquake dataset based on the ChinArray Phase I records(X1).ChinArray Phase I was deployed in the southern north-south seismic zone(20°N-32°N,95°E-110°E)in 2011-2013 using 355 portable broadband seismic stations.CREDIT-X1local,the first release of the ChinArray Reference Earthquake Dataset for Innovative Techniques(CREDIT),includes comprehensive information for the 105,455 local events that occurred in the southern north-south seismic zone during array observation,incorporating them into a single HDF5 file.Original 100-Hz sampled three-component waveforms are organized by event for stations within epicenter distances of 1,000 km,and records of≥200 s are included for each waveform.Two types of phase labels are provided.The first includes manually picked labels for 5,999 events with magnitudes≥2.0,providing 66,507 Pg,42,310 Sg,12,823 Pn,and 546 Sn phases.The second contains automatically labeled phases for 105,442 events with magnitudes of−1.6 to 7.6.These phases were picked using a recurrent neural network phase picker and screened using the corresponding travel time curves,resulting in 1,179,808 Pg,884,281 Sg,176,089 Pn,and 22,986 Sn phases.Additionally,first-motion polarities are included for 31,273 Pg phases.The event and station locations are provided,so that deep learning networks for both conventional phase picking and phase association can be trained and validated.The CREDIT-X1local dataset is the first million-scale dataset constructed from a dense seismic array,which is designed to support various multi-station deep-learning methods,high-precision focal mechanism inversion,and seismic tomography studies.Additionally,owing to the high seismicity in the southern north-south seismic zone in China,this dataset has great potential for future scientific discoveries.
文摘In this paper, mantle circulation flow, continental drift, earthquake origin and other mechanical principles are examined as they apply to earthquake engineering, seismology and dynamics of fluid saturated porous medium. The relationship of mantle flow to earthquakes is examined and clarified, and a new model, different from Haskell’s, is proposed for the earthquake mechanism. The proposed new model is based on the discovery that two pairs of jump stress and jump velocity will start to act from the fault plane. Records obtained directly from recent earthquakes nearby and right on the fault break show a very large velocity impulse, which verify, indirectly, the new mechanism proposed by the author. Further, at least two physical parameters that characterize the seismic intensity must be specified, because according to the discontinuous (jump) wave theory, at the earthquake source, the stress jump and the velocity jump of particle motion should act simultaneously when a sudden break occurs. The third key parameter is shown to be the break (fracture) propagation speed together with the break plane area. This parameter influences the form of the unloading time function at the source. The maximum seismic stress in and displacement of a building are estimated for two unfavorable combinations of the building and its base ground in terms of their relative rigidity. Finally, it is shown that Biot’s theory of wave propagation in fluid saturated porous media is valid only when fluid flow cannot occur.
文摘The planetary geometries were studied at the moments of 92 strong earthquakes with a magnitude of more than 8 on the Richter scale (R8+) for the period from 1900 to 2011. Three main planetary schemes were specified, namely the schemes of the generalized Archimedes lever (gAL) and the Kepler conjunction (gKc), as well as a new geometry of the triangles of the remote signal catcher (cRS). It was discovered 22 gAL, 42 gKc and 28 cRS geometries, which are 23.9%, 45.7%, and 30.4%, respectively, from the total 92 studying cases. It was found that in some earthquakes;the planetary geometries are absolutely identical, which indicates the universality of the mechanism that caused the earthquake. The triggered effect does not depend on the distance between the planets and the mass of the planets, so the mechanism is identified as an inertial gravitational interaction. The triggered effect increases with the multiplicity of the ratio of distances between planets, as well as with pairwise planetary parallelism, which probably indicates about the wave nature of inertial effects. The triggering effect increases with increasing multiplicity of the ratio of distances between planets, as well as with pairwise planetary parallelism, which probably indicates about the wave nature of inertial effects. According to the Archimedes’ lever principle, the Third Law of Motion is changed by adding a few words. It is assumed that inertia is a special case of gravity, namely gravitational self-induction, which really depends, like any self-induction, only on the geometry of the task.
文摘A number of aftershocks of the May 10th 1997, Zirkuh (Ghaen-Birjand) destructive earthquake have been used to investigate the anisotropy in the upper crust by observing shear wave splitting. Particle motion diagram and aspect ratio methods were used as two different approaches to obtain splitting parameters. Clear shear wave splitting was observed on the records of the selected aftershocks, indicating that the media in the region was highly anisotropic. By using particle motion method, the direction of fast shear wave was found 22°N±19°E, while the delay time between the fast and slow shear waves was obtained to be (65±16) ms. By aspect ratio method, the direction of fast shear wave was determined to be 35°N±18°E and the delay time between fast and slow shear waves was found to be (49±10) ms. For a simple horizontal layer with a thickness about 5 km and uniformly distributed anisotropy, a stress aligned cracks model was used and the result was interpreted in terms of vertical aligned cracks in the direction of N22°E, having a density about 0.01. It is assumed that cracks are fluid-filled since they are located in the upper crust. Finally, by using Hudson cracks model for three crack densities 0.005, 0.01, 0.03, the velocity curves of shear wave were plotted as a function of angle between the symmetrical axis of cracks and the azimuth of source to receiver. It was concluded that when shear wave was polarized parallel to the crack surface, the velocity was uniform, but the velocity curve varied clearly if shear wave was polarized perpendicular to the crack surface.
基金supported by the Special Found of the Institute of Geophysics, China Earthquake Administration (DQJB22B18)
文摘To reveal the seismogenic mechanism of the Luding earthquake, we employed the 118 China Seismic Network stations to collect the P-wave polarity data from each station, which was then used in the P-wave first motion approach to calculate the focal mechanism solution of the M6.8 Luding earthquake that occurred on September 5,2022. We have also studied the loading effect of tectonic stress on the Luding earthquake fault based on the stress field data for the research area. The results indicate that this earthquake was a strike-slip type, the nodal plane I:strike 167°, dip angle 78°, slip angle 2°;Nodal plane II: strike 77°, dip angle 88°, slip angle 168°. The two fault planes’ instability coefficients of the Luding earthquake are examined considering the region’s background stress field’s condition. The nodal plane I in the Moho circle is discovered to practically coincide with the Coulomb failure line and the tangent point of the Moho circle, indicating that this nodal plane has a high instability coefficient compared to the nodal plane II. The conclusion is that the nodal plane I has a higher likelihood of being the seismogenic fault plane, which is congruent with the seismogenic fault plane suggested by the aftershock distribution, the earthquake radiation energy distribution of a single station, and seismic intensity distribution.The Luding earthquake’s focal mechanism is highly like the theoretical focal mechanism of the fault situated at the location where the Coulomb failure line intersects the Mohr circle, demonstrating that background stress is what caused the earthquake. The substantial fault instability and similarity between the solved and theoretical focal mechanisms make it easier to comprehend the loading effect of tectonic stress on the Luding earthquake fault.
