A 3D crustal model was constructed using a combination of cutting-edge techniques,which were integrated to provide a density model for Egypt and address the sporadic distribution of seismic data.These techniques inclu...A 3D crustal model was constructed using a combination of cutting-edge techniques,which were integrated to provide a density model for Egypt and address the sporadic distribution of seismic data.These techniques include obtaining gravity data from the Gravity Field and Steady-State Ocean Circulation Explorer(GOCE),creating seismic profiles,analyzing the receiver functions of seismic data,obtaining information from boreholes,and providing geological interpretations.GOCE satellite gravity data were processed to construct a preliminary model based on nonlinear inversions of the data.A regional crustal thickness model was developed using receiver functions,seismic refraction profiles,and geological insights.The inverted model was validated using borehole data and compared with seismic estimates.The model exhibited strong consistency and revealed a correlation between crustal thickness,geology,and tectonics of Egypt.It showed that the shallowest depths of the Moho are located in the north along the Mediterranean Sea and in the eastern part along the Red Sea,reflecting an oceanic plate with a thin,high-density crust.The deepest Moho depths are located in the southwestern part of Egypt,Red Sea coastal mountains,and Sinai Peninsula.The obtained 3D model of crustal thickness provided finely detailed Moho depth estimates that aligned closely with geology and tectonic characteristics of Egypt,contributing valuable insights into the subsurface structure and tectonic processes of region.展开更多
Knowledge regarding earthquake hazards and seismicity is crucial for crisis management, and the occurrence of foreshocks, seismic activity patterns, and spatiotemporal variations in seismic activity have been studied....Knowledge regarding earthquake hazards and seismicity is crucial for crisis management, and the occurrence of foreshocks, seismic activity patterns, and spatiotemporal variations in seismic activity have been studied. Furthermore, the estimation of the region-time-length (RTL) parameter has been proposed to detect seismic quiescence before the occurrence of a large earthquake. In addition, the time-to-failure method has been used to estimate the time occurrence of large earthquakes. Hence, in this study, to gain deeper insight into seismic activity in the southern Zagros region, we utilized the RTL algorithm to identify the quiescence and activation phases leading to the Fin doublet earthquakes. Temporal variations in the RTL parameter showed two significant anomalies. One corresponded to the occurrence time of the first earthquake (2017-12-12);the other anomaly was associated with the occurrence time of the second event (2021-11-14). Based on a negative value of the RTL parameter observed in the vicinity of the Fin epicenters (2021), seismic quiescence (a decrease in seismicity compared to the preceding background rate) was identified. The spatial distribution of the RTL prognostic parameters confirms the appearance of seismic quiescence surrounding the epicenter of the Fin doublet earthquakes (2021). The time-to-failure method was designed using precursory events that describe the acceleration of the seismic energy release before the mainshock. Using the time-to-failure method for the earthquake catalog, it was possible to estimate both the magnitude and time of failure of the Fin doublet. Hence, the time-tofailure technique can be a useful supplementary method to the RTL algorithm for determining the characteristics of impending earthquakes.展开更多
Seismic records produced by different seismic sources vary.In this study,we compared the waveform records and time-frequency characteristics of tectonic earthquakes,artificial explosions,and mine collapses in China’s...Seismic records produced by different seismic sources vary.In this study,we compared the waveform records and time-frequency characteristics of tectonic earthquakes,artificial explosions,and mine collapses in China’s Capital Region.The results show that tectonic earthquakes are characterized by stronger S-wave energy than P-wave energy,obvious high-frequency components,and wide frequency bands of P and S waves.Artificial explosions are characterized by greater P-wave amplitude than S-wave amplitude and near-station surface wave development.Mine collapses are characterized by lower overall frequency,more obvious surface waves,and longer duration.We extracted quantitative discriminants based on the analysis of different event records,with 31 feature values in 7 categories(P/S maximum amplitude ratio,high/low frequency energy ratio,P/S spectral ratio,corner frequency,duration,the second-order moment of spectrum,and energy strongest point).A comparison of the ability of these feature values to recognize distinct events showed that the 6-17 Hz P/S spectral ratio was able to completely distinguish artificial explosions from the other two types of events.The S-wave corner frequency performed relatively well in identifying all three types of events,with an accuracy of over 90%.Additionally,a support vector machine was used to comprehensively distinguish multiple features,with an accuracy for all three types of events reaching up to 100%.展开更多
A two and a half dimensional(2.5D)multi-domain indirect boundary element method(IBEM)is developed to study the wave scattering of obliquely incident P-,SV-and SH-waves by a hill-valley staggered topography in a multi-...A two and a half dimensional(2.5D)multi-domain indirect boundary element method(IBEM)is developed to study the wave scattering of obliquely incident P-,SV-and SH-waves by a hill-valley staggered topography in a multi-layered half-space.The IBEM algorithm includes using 2.5D full-space and half-space Green’s functions to construct scattered fields in decomposed closed and opened half-space regions,respectively,and using the dynamic stiffness method to solve the free fields.All regions are finally integrated by introducing the compatibility conditions to obtain the total wave fields.The proposed 2.5D IBEM has the flexibility in dealing with complex boundaries by directly applying the fictitious loads on the regions’boundaries,with a less storage requirement compared to the full 3D models.Besides,by combining the specific advantages of the two kinds of Green’s functions,the method is well suitable for handling coupled topographies with high accuracy.The method is validated by comparison with published results for a single valley as well as a single hill topography.The effects of height-to-width ratio of hill and layering on dynamic responses are further parametrically investigated by numerical implementations in frequency domain.Results show that the interaction between valley and hills can lead to a more significant amplification within the valley region,and dynamic responses are deeply influenced by the height-to-width of hill,simultaneously depending on incident angle and frequency.Besides,the site effects become more complex when the stratification feature is taken into account.展开更多
During the ruptures of an earthquake, the strain energy, AE, will be transferred into, at least, three parts, i.e., the seismic radiation energy （Es）, fracture energy （Eg）, and frictional energy （Ef）, that is, A...During the ruptures of an earthquake, the strain energy, AE, will be transferred into, at least, three parts, i.e., the seismic radiation energy （Es）, fracture energy （Eg）, and frictional energy （Ef）, that is, AE = Es ＋ Eg ＋ El. Friction, which is represented by a velocity- and state-de- pendent friction law by some researchers, controls the three parts. One of the main parameters of the law is the char- acteristic slip displacement, De. It is significant and nec- essary to evaluate the reliable value of Dc from observed and inverted seismic data. Since Dc controls the radiation efficiency, ηR = Es/（Es ＋ Eg）, the value of qR is a good constraint of estimating Dc. Integrating observed data and inverted results of source parameters from recorded seis- mograms, the values of Es and Eg of an earthquake can be measured, thus leading to the value of ηR. The constraint used to estimate the reliable value of Dc will be described in this work. An example of estimates of Dc based on the observed and inverted values of source parameters of the September 20, 1999 Ms 7.6 Chi-Chi （Ji-Ji）, Taiwan region, earthquake will be presented.展开更多
The topography of the inner core is crucial to understand its growth process and interaction with the geodynamo. With the accuracy of teleseismic waveform doublets in determining the travel-time shifts between PKPcd a...The topography of the inner core is crucial to understand its growth process and interaction with the geodynamo. With the accuracy of teleseismic waveform doublets in determining the travel-time shifts between PKPcd and PKPdf inner-core phases, we examined the temporal change of the inner-core boundary sampled by new earthquake doublets that occurred in the Western Pacific and those previously found from the South Sandwich Islands. The receiver stations are those within the distance range of 128°–142° from the hypocenters of the waveform doublets. Our results suggest that temporal changes in PKPcd-PKPdf differential times are very subtle except some isolated regions under central America and Africa. The observations may indicate localized topography of the inner core and/or the inner-core boundary as transient slurry in isolated regions.展开更多
The occurrence of earthquakes in oceanic tren- ches can pose a tsunami threat to lives and properties in active seismic zones. Therefore, the knowledge of focal depth, magnitude, and time distribution of earthquakes a...The occurrence of earthquakes in oceanic tren- ches can pose a tsunami threat to lives and properties in active seismic zones. Therefore, the knowledge of focal depth, magnitude, and time distribution of earthquakes along the trenches is needed to investigate the future occurrence of earthquakes in the zones. The oceanic trenches studied, were located from the seismicity map on： latitude ＋51 ° to＋53° and longitude -160° to 176° （Aleutian Trench）, latitude ＋40° to ＋53° and longitude ＋148° to＋165° （Japan Trench）, and latitude -75° to -64° and longitude -15° to ＋30° （Peru-Chile Trench）. The following features of seis- mic events were considered： magnitude distribution, focal depth distribution, and time distribution of earthquake. The results obtained in each trench revealed that the earthquakes increased with time in all the regions. This implies that the lithospheric layer is becoming more unstable. Thus, tectonic stress accumulation is increasing with time. The rate of increase in earthquakes at the Peru-Chile Trench is higher than that of the Japan Trench and the Aleutian Trench. This implies that the convergence of lithospheric plates is higher in the Peru-Chile Trench. Deep earthquakes were observed across all the trenches. The shallow earthquakes were more prominent than intermediate and deep earthquakes in all thetrenches. The seismic events in the trenches are mostly of magnitude range 3.0-4.9. This magnitude range may indi- cate the genesis of mild to moderate tsunamis in the trench zone in near future once sufficient slip would occur with displacement of water column.展开更多
This article presents an application of a procedure to modify the intensity distribution by assessing the reliability. There are two potential possibilities that may influence the intensity distribution： （1） For th...This article presents an application of a procedure to modify the intensity distribution by assessing the reliability. There are two potential possibilities that may influence the intensity distribution： （1） For the interpolation error, we generate a measured grid across the calculation region. When the point to station spacing is 〈5 km, we consider the results precise; however, some points have less precision because these are farther from the corresponding stations. When the spacing is between 5 and 50 km, we consider the results imprecise and define a reliability factor that correlates with the distance. （2） Some records may have errors that result from local site conditions, equipment problems, or some disturbance such as lightning stroke, which will lead to some grid points having an incorrect intensity. We regress the attenuation relation for sites with abnormal intensities and consider the results to be accurate when the standard deviation （STD） is 〈σ and inaccurate when the STD is 〉 2σ. We then define a reliability factor to correlate with STD between σ and 2σ, such that the intensity distribution is in accord with both wave propagation theory and the investigation intensity.展开更多
We conducted comprehensive receiver function analyses for a large amount of high-quality broadband teleseismic waveforms data recorded at 19 China National Digital Seismic Network （CNDSN） stations deployed in Northe...We conducted comprehensive receiver function analyses for a large amount of high-quality broadband teleseismic waveforms data recorded at 19 China National Digital Seismic Network （CNDSN） stations deployed in Northeast China. An advanced H-κ domain search method was adopted to accurately estimate the crustal thickness and vp/vs ratio. The crust has an average thickness of about 34.4 km. The thinnest crust occurs in the central region of Northeast China, while the thickest crust is beneath the Yanshan belt. The vp/vs ratio is relatively uniform with an average of about 1.733. The highest vp/vs ratio is found beneath the Changbaishan, likely associated with its volcanic activities. We found significant lateral heterogeneity beneath three stations CN2, MDJ, and MIH located along the Suolon suture from the back-zimuthal dependence of Moho depth. The velocity modeling from receiver functions indicated complicated Earth structure beneath these stations with large crust-mantle transition zone, noticeable velocity jump in upper mantle, and low velocity zone in middle crust. Dipping velocity interface in the crust with strike approximately parallel to the Suolon suture and down-dip to the south or southeast might explain the observed lateral heterogeneity.展开更多
To better understand repeatability of strong earthquakes in previously ruptured zones during one seismogenic period, we studied the rupture zones of the doublet of M6 earthquakes in Zhongba region of southcentral Tibe...To better understand repeatability of strong earthquakes in previously ruptured zones during one seismogenic period, we studied the rupture zones of the doublet of M6 earthquakes in Zhongba region of southcentral Tibet, China, in 11 July 2004 and 7 April 2005, respectively. We focused on the overlapping degree of two strong quakes’ aftershock areas one week after the mainshocks by using the SQH station in China Seismic Network and a 68-stations temporary broadband seismic array, a part of the international HI-CLIMB project. About 115 local earthquakes were recorded in one week after the mainquakes, and we located these earthquakes by master event relative location （MERL） method. We also used this method to relocate 31 other M3.7 ＋ earthquakes from 1 July 2004 to 1 July 2005. Meanwhile, we studied two mainshocks’ coseismic ruptures with satellite interferometric synthetic aperture radar （InSAR）. Our results show that the ruptured zones of the two earthquakes do not overlapp substantially, either from early aftershock data or from InSAR inversions.展开更多
We collect 1974 broad-band velocity records of 94 earthquakes (ML=2.84.9, △=13462 km) from seven stations of the Fujian Seismic Network from March 1999 to March 2007. Using real-time simulation, we obtain the corresp...We collect 1974 broad-band velocity records of 94 earthquakes (ML=2.84.9, △=13462 km) from seven stations of the Fujian Seismic Network from March 1999 to March 2007. Using real-time simulation, we obtain the corresponding acceleration and then adopt different models to analyze the seismic data. As a result, a new attenuation relationship between PGA and PGV of the small and moderate earthquakes on bedrock site in Fujian region is established. The Yongchun earthquake occurred recently verifies the attenuation relationship well. This paper provides a new approach for studying the ground motion attenuation relationship using velocity records.展开更多
The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The co...The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.展开更多
A new concept is suggested on tectonomagnetic research about the noise in simultaneous geomagnetic difference data caused by the effect of Sq local-time variation, together with the method of theoretical calculation. ...A new concept is suggested on tectonomagnetic research about the noise in simultaneous geomagnetic difference data caused by the effect of Sq local-time variation, together with the method of theoretical calculation. The level of the noise and its contribution to the total noises of the differences data are analyzed. The result indicates that the noise increases linearly with the increase of the distance between the two stations in the range of 40° longitude-difference, and its increasing rate is about 0.4 nT/（°）at latitude 40°N. The example calculated at a pair of sites with longitude-difference 0.357°, shows that the noise is about one fifth of the total noises of the difference data on geomagnetic quiet-day.展开更多
On the basis of the results from spatial-temporal scanning of moderate seismic activities prior to 63 strong earth- quakes occurred since 1966, we investigate in this paper the characteristics of anomalous evolution o...On the basis of the results from spatial-temporal scanning of moderate seismic activities prior to 63 strong earth- quakes occurred since 1966, we investigate in this paper the characteristics of anomalous evolution of moderate seismic activities prior to strong earthquakes with different focal mechanisms in different tectonic blocks, as well as the correlation between principal seismic activities and nodal planes of focal mechanisms. The research results indicate that the moderate seismic activities prior to strong earthquakes display two significant stages of stress concentration and stress weakening, which are consistent to the evolvement of local tectonic stress field. The first stage, defined as medium-term stage, is characterized by seismic belts, seismic gaps and centralized activities; The second stage, defined as short-term stage, is predominated by continuous quiescence. For the earthquakes occurred on the strike-slip and normal faults, the two stages are relatively evident; while for the earthquakes occurred on the thrust faults, there are not characteristic images and seismic quiescence in both stages. The spatial distribution of moderate shock belts before strong earthquakes show a good coherence with that of the two nodal planes of focal mechanisms, while the spatial distribution of seismic gaps enclosed by moderate shocks display a poor coherence with the nodal planes.展开更多
A network communication and recording system based on China-made ARCA SOC and embedded Linux operating system is introduced in this paper. It supports TCP/IP network communication protocol and mass storage medium. It ...A network communication and recording system based on China-made ARCA SOC and embedded Linux operating system is introduced in this paper. It supports TCP/IP network communication protocol and mass storage medium. It has strong points of self-monitor, low power consumption, high timing accuracy, high reliability of operation, etc. It can serve up to 20 centers real-time waveform data at the same time. It meets not only the requirements of physical networking observation, but also virtual networking observation based on Internet in which real-time data transmission is required. Its ability of field recording also meets the requirements of portable seismic observation, strong motion observation and seismic exploration observation, etc.展开更多
According to the rupture dynamics of earthquakes, variations of the apparent stress and the difference between the static stress drop and the dynamic stress drop during the rupture of earthquakes are analyzed for the ...According to the rupture dynamics of earthquakes, variations of the apparent stress and the difference between the static stress drop and the dynamic stress drop during the rupture of earthquakes are analyzed for the July 20, 1995 ML=4.1 Shacheng, Hebei, China, earthquake sequence. Results obtained show that the apparent stress for main-shock is about 5 MPa, and the average apparent stress for aftershocks 0.047 MPa. During the rupture of the main-shock, the dynamic stress drop is approximately 1.6 times greater than the static stress drop with the difference of nearly 2.7 MPa. The dynamic stress drop is less than the static stress drop for all aftershocks with the average difference of -0.75 MPa. Therefore, when the mainshock occurs the final stress on the focal fault is higher than the dynamic frictional stress, corresponding to that the fault is abruptly locked. When the aftershocks occur the final stress on the focal fault is lower than the dynamic frictional stress, corresponding to that the fault overshoots. It can be seen from the above results that there could be some differences in the physic processes between the mainshock and the aftershocks.展开更多
According to the deficiency of the strain accumulating and releasing curves and the previous models, the strain-accumulating rate of the strain accumulating and releasing model has been deduced based on the G-R relati...According to the deficiency of the strain accumulating and releasing curves and the previous models, the strain-accumulating rate of the strain accumulating and releasing model has been deduced based on the G-R relation and the empirical formula between energy release and earthquake magnitude, where the strain-accumulating rate is relative independent of the strain-releasing rate. Five typical areas in Chinese mainland are selected on the basis of the hypothesis on active tectonic block, and small earthquakes from 1970 are imported to calculate the annual strain-accumulating rates considering the completeness of historical seismic data. Having introduced the strain-accumulating rates into the amended model, present strain phases are got. According to the present stages in their own cycles, the future earthquake tendency of each sub-region is discussed.展开更多
The data of short-period (1~18 s) surface waves recorded by 23 stations belonging to the digital seismic network of Yunnan Province of China are used in this paper. From these data, the dispersion curves of phase velo...The data of short-period (1~18 s) surface waves recorded by 23 stations belonging to the digital seismic network of Yunnan Province of China are used in this paper. From these data, the dispersion curves of phase velocities of the fundamental mode Rayleigh wave along 209 paths are determined by using the two-station narrowband filtering cross-correlation method. Adopting tomography method, the distribution maps of phase velocities at various peri-ods in Yunnan region are inverted. The maps of phase velocities on profiles along 24N, 25N, 26N, 27N and 100.5E and the distribution maps of phase velocities at 3 periods in the study region are given. The results show that the phase velocity distribution in Yunnan region has strong variations in horizontal direction, and the phase velocity distribution in short-period range is closely related to the thickness of sedimentary layers in the shallow crust. The phase velocity in southern part of the Sichuan-Yunnan rhombic block encircled by the Honghe fault and Xiaojiang fault is obviously lower than that in surrounding areas. The epicentral locations of strong earthquakes in Yunnan region are mainly distributed in transitional zones between low and high phase velocities.展开更多
The Biot and Squirt-flow are the two most important mechanisms of fluid flow in the porous medium with fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, dispersion and ...The Biot and Squirt-flow are the two most important mechanisms of fluid flow in the porous medium with fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, dispersion and attenuation of elastic waves in the porous medium are widely investigated in recent years. However, we have not read any reports on numerical simulation based on the BISQ equation. In this paper, following the BISQ equation, elastic wave propagation in the transversely isotropic porous medium filled with fluids is simulated by the stag-gered grid method for different frequency and phase boundary cases and the two-layer medium. And propagating characteristics of seismic and acoustic waves and various phenomena occured in the propagating process are in-vestigated when the two mechanisms are considered simultaneously.展开更多
We developed a GA-BP algorithm by combining the genetic algorithm (GA) with the back propagation (BP) algorithm and established a genetic BP neural network. We also applied the BP neural network based on the BP algori...We developed a GA-BP algorithm by combining the genetic algorithm (GA) with the back propagation (BP) algorithm and established a genetic BP neural network. We also applied the BP neural network based on the BP algorithm and the genetic BP neural network based on the GA-BP algorithm to discriminate earthquakes and explosions. The obtained result shows that the discriminating performance of the genetic BP network is slightly better than that of the BP network.展开更多
文摘A 3D crustal model was constructed using a combination of cutting-edge techniques,which were integrated to provide a density model for Egypt and address the sporadic distribution of seismic data.These techniques include obtaining gravity data from the Gravity Field and Steady-State Ocean Circulation Explorer(GOCE),creating seismic profiles,analyzing the receiver functions of seismic data,obtaining information from boreholes,and providing geological interpretations.GOCE satellite gravity data were processed to construct a preliminary model based on nonlinear inversions of the data.A regional crustal thickness model was developed using receiver functions,seismic refraction profiles,and geological insights.The inverted model was validated using borehole data and compared with seismic estimates.The model exhibited strong consistency and revealed a correlation between crustal thickness,geology,and tectonics of Egypt.It showed that the shallowest depths of the Moho are located in the north along the Mediterranean Sea and in the eastern part along the Red Sea,reflecting an oceanic plate with a thin,high-density crust.The deepest Moho depths are located in the southwestern part of Egypt,Red Sea coastal mountains,and Sinai Peninsula.The obtained 3D model of crustal thickness provided finely detailed Moho depth estimates that aligned closely with geology and tectonic characteristics of Egypt,contributing valuable insights into the subsurface structure and tectonic processes of region.
文摘Knowledge regarding earthquake hazards and seismicity is crucial for crisis management, and the occurrence of foreshocks, seismic activity patterns, and spatiotemporal variations in seismic activity have been studied. Furthermore, the estimation of the region-time-length (RTL) parameter has been proposed to detect seismic quiescence before the occurrence of a large earthquake. In addition, the time-to-failure method has been used to estimate the time occurrence of large earthquakes. Hence, in this study, to gain deeper insight into seismic activity in the southern Zagros region, we utilized the RTL algorithm to identify the quiescence and activation phases leading to the Fin doublet earthquakes. Temporal variations in the RTL parameter showed two significant anomalies. One corresponded to the occurrence time of the first earthquake (2017-12-12);the other anomaly was associated with the occurrence time of the second event (2021-11-14). Based on a negative value of the RTL parameter observed in the vicinity of the Fin epicenters (2021), seismic quiescence (a decrease in seismicity compared to the preceding background rate) was identified. The spatial distribution of the RTL prognostic parameters confirms the appearance of seismic quiescence surrounding the epicenter of the Fin doublet earthquakes (2021). The time-to-failure method was designed using precursory events that describe the acceleration of the seismic energy release before the mainshock. Using the time-to-failure method for the earthquake catalog, it was possible to estimate both the magnitude and time of failure of the Fin doublet. Hence, the time-tofailure technique can be a useful supplementary method to the RTL algorithm for determining the characteristics of impending earthquakes.
文摘Seismic records produced by different seismic sources vary.In this study,we compared the waveform records and time-frequency characteristics of tectonic earthquakes,artificial explosions,and mine collapses in China’s Capital Region.The results show that tectonic earthquakes are characterized by stronger S-wave energy than P-wave energy,obvious high-frequency components,and wide frequency bands of P and S waves.Artificial explosions are characterized by greater P-wave amplitude than S-wave amplitude and near-station surface wave development.Mine collapses are characterized by lower overall frequency,more obvious surface waves,and longer duration.We extracted quantitative discriminants based on the analysis of different event records,with 31 feature values in 7 categories(P/S maximum amplitude ratio,high/low frequency energy ratio,P/S spectral ratio,corner frequency,duration,the second-order moment of spectrum,and energy strongest point).A comparison of the ability of these feature values to recognize distinct events showed that the 6-17 Hz P/S spectral ratio was able to completely distinguish artificial explosions from the other two types of events.The S-wave corner frequency performed relatively well in identifying all three types of events,with an accuracy of over 90%.Additionally,a support vector machine was used to comprehensively distinguish multiple features,with an accuracy for all three types of events reaching up to 100%.
文摘A two and a half dimensional(2.5D)multi-domain indirect boundary element method(IBEM)is developed to study the wave scattering of obliquely incident P-,SV-and SH-waves by a hill-valley staggered topography in a multi-layered half-space.The IBEM algorithm includes using 2.5D full-space and half-space Green’s functions to construct scattered fields in decomposed closed and opened half-space regions,respectively,and using the dynamic stiffness method to solve the free fields.All regions are finally integrated by introducing the compatibility conditions to obtain the total wave fields.The proposed 2.5D IBEM has the flexibility in dealing with complex boundaries by directly applying the fictitious loads on the regions’boundaries,with a less storage requirement compared to the full 3D models.Besides,by combining the specific advantages of the two kinds of Green’s functions,the method is well suitable for handling coupled topographies with high accuracy.The method is validated by comparison with published results for a single valley as well as a single hill topography.The effects of height-to-width ratio of hill and layering on dynamic responses are further parametrically investigated by numerical implementations in frequency domain.Results show that the interaction between valley and hills can lead to a more significant amplification within the valley region,and dynamic responses are deeply influenced by the height-to-width of hill,simultaneously depending on incident angle and frequency.Besides,the site effects become more complex when the stratification feature is taken into account.
基金financially supported by Academia Sinica and Ministry of Science and Technology under Grand Nos.of MOST 103-2116-M-001-010 and MOST 104-2116-M001-007
文摘During the ruptures of an earthquake, the strain energy, AE, will be transferred into, at least, three parts, i.e., the seismic radiation energy （Es）, fracture energy （Eg）, and frictional energy （Ef）, that is, AE = Es ＋ Eg ＋ El. Friction, which is represented by a velocity- and state-de- pendent friction law by some researchers, controls the three parts. One of the main parameters of the law is the char- acteristic slip displacement, De. It is significant and nec- essary to evaluate the reliable value of Dc from observed and inverted seismic data. Since Dc controls the radiation efficiency, ηR = Es/（Es ＋ Eg）, the value of qR is a good constraint of estimating Dc. Integrating observed data and inverted results of source parameters from recorded seis- mograms, the values of Es and Eg of an earthquake can be measured, thus leading to the value of ηR. The constraint used to estimate the reliable value of Dc will be described in this work. An example of estimates of Dc based on the observed and inverted values of source parameters of the September 20, 1999 Ms 7.6 Chi-Chi （Ji-Ji）, Taiwan region, earthquake will be presented.
基金supported by the Natural Science Foundation of China(41330209)the Ministry of Science and Technology of China(2014CB845901)the China Earthquake Administration
文摘The topography of the inner core is crucial to understand its growth process and interaction with the geodynamo. With the accuracy of teleseismic waveform doublets in determining the travel-time shifts between PKPcd and PKPdf inner-core phases, we examined the temporal change of the inner-core boundary sampled by new earthquake doublets that occurred in the Western Pacific and those previously found from the South Sandwich Islands. The receiver stations are those within the distance range of 128°–142° from the hypocenters of the waveform doublets. Our results suggest that temporal changes in PKPcd-PKPdf differential times are very subtle except some isolated regions under central America and Africa. The observations may indicate localized topography of the inner core and/or the inner-core boundary as transient slurry in isolated regions.
文摘The occurrence of earthquakes in oceanic tren- ches can pose a tsunami threat to lives and properties in active seismic zones. Therefore, the knowledge of focal depth, magnitude, and time distribution of earthquakes along the trenches is needed to investigate the future occurrence of earthquakes in the zones. The oceanic trenches studied, were located from the seismicity map on： latitude ＋51 ° to＋53° and longitude -160° to 176° （Aleutian Trench）, latitude ＋40° to ＋53° and longitude ＋148° to＋165° （Japan Trench）, and latitude -75° to -64° and longitude -15° to ＋30° （Peru-Chile Trench）. The following features of seis- mic events were considered： magnitude distribution, focal depth distribution, and time distribution of earthquake. The results obtained in each trench revealed that the earthquakes increased with time in all the regions. This implies that the lithospheric layer is becoming more unstable. Thus, tectonic stress accumulation is increasing with time. The rate of increase in earthquakes at the Peru-Chile Trench is higher than that of the Japan Trench and the Aleutian Trench. This implies that the convergence of lithospheric plates is higher in the Peru-Chile Trench. Deep earthquakes were observed across all the trenches. The shallow earthquakes were more prominent than intermediate and deep earthquakes in all thetrenches. The seismic events in the trenches are mostly of magnitude range 3.0-4.9. This magnitude range may indi- cate the genesis of mild to moderate tsunamis in the trench zone in near future once sufficient slip would occur with displacement of water column.
文摘This article presents an application of a procedure to modify the intensity distribution by assessing the reliability. There are two potential possibilities that may influence the intensity distribution： （1） For the interpolation error, we generate a measured grid across the calculation region. When the point to station spacing is 〈5 km, we consider the results precise; however, some points have less precision because these are farther from the corresponding stations. When the spacing is between 5 and 50 km, we consider the results imprecise and define a reliability factor that correlates with the distance. （2） Some records may have errors that result from local site conditions, equipment problems, or some disturbance such as lightning stroke, which will lead to some grid points having an incorrect intensity. We regress the attenuation relation for sites with abnormal intensities and consider the results to be accurate when the standard deviation （STD） is 〈σ and inaccurate when the STD is 〉 2σ. We then define a reliability factor to correlate with STD between σ and 2σ, such that the intensity distribution is in accord with both wave propagation theory and the investigation intensity.
文摘We conducted comprehensive receiver function analyses for a large amount of high-quality broadband teleseismic waveforms data recorded at 19 China National Digital Seismic Network （CNDSN） stations deployed in Northeast China. An advanced H-κ domain search method was adopted to accurately estimate the crustal thickness and vp/vs ratio. The crust has an average thickness of about 34.4 km. The thinnest crust occurs in the central region of Northeast China, while the thickest crust is beneath the Yanshan belt. The vp/vs ratio is relatively uniform with an average of about 1.733. The highest vp/vs ratio is found beneath the Changbaishan, likely associated with its volcanic activities. We found significant lateral heterogeneity beneath three stations CN2, MDJ, and MIH located along the Suolon suture from the back-zimuthal dependence of Moho depth. The velocity modeling from receiver functions indicated complicated Earth structure beneath these stations with large crust-mantle transition zone, noticeable velocity jump in upper mantle, and low velocity zone in middle crust. Dipping velocity interface in the crust with strike approximately parallel to the Suolon suture and down-dip to the south or southeast might explain the observed lateral heterogeneity.
基金supported jointly by Joint Seismological Science Foundation of China (Nos.200808078 and 200708035)
文摘To better understand repeatability of strong earthquakes in previously ruptured zones during one seismogenic period, we studied the rupture zones of the doublet of M6 earthquakes in Zhongba region of southcentral Tibet, China, in 11 July 2004 and 7 April 2005, respectively. We focused on the overlapping degree of two strong quakes’ aftershock areas one week after the mainshocks by using the SQH station in China Seismic Network and a 68-stations temporary broadband seismic array, a part of the international HI-CLIMB project. About 115 local earthquakes were recorded in one week after the mainquakes, and we located these earthquakes by master event relative location （MERL） method. We also used this method to relocate 31 other M3.7 ＋ earthquakes from 1 July 2004 to 1 July 2005. Meanwhile, we studied two mainshocks’ coseismic ruptures with satellite interferometric synthetic aperture radar （InSAR）. Our results show that the ruptured zones of the two earthquakes do not overlapp substantially, either from early aftershock data or from InSAR inversions.
基金Joint Seismological Science Foundation of China (105034)
文摘We collect 1974 broad-band velocity records of 94 earthquakes (ML=2.84.9, △=13462 km) from seven stations of the Fujian Seismic Network from March 1999 to March 2007. Using real-time simulation, we obtain the corresponding acceleration and then adopt different models to analyze the seismic data. As a result, a new attenuation relationship between PGA and PGV of the small and moderate earthquakes on bedrock site in Fujian region is established. The Yongchun earthquake occurred recently verifies the attenuation relationship well. This paper provides a new approach for studying the ground motion attenuation relationship using velocity records.
基金National Natural Science Foundation of China (50608024 and 50538050).
文摘The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.
基金Joint Seismological Science Foundation of China (198009).
文摘A new concept is suggested on tectonomagnetic research about the noise in simultaneous geomagnetic difference data caused by the effect of Sq local-time variation, together with the method of theoretical calculation. The level of the noise and its contribution to the total noises of the differences data are analyzed. The result indicates that the noise increases linearly with the increase of the distance between the two stations in the range of 40° longitude-difference, and its increasing rate is about 0.4 nT/（°）at latitude 40°N. The example calculated at a pair of sites with longitude-difference 0.357°, shows that the noise is about one fifth of the total noises of the difference data on geomagnetic quiet-day.
基金State Science and Techniology Key Program (2004BA601B01-01, 2004BA601B01-01-04).
文摘On the basis of the results from spatial-temporal scanning of moderate seismic activities prior to 63 strong earth- quakes occurred since 1966, we investigate in this paper the characteristics of anomalous evolution of moderate seismic activities prior to strong earthquakes with different focal mechanisms in different tectonic blocks, as well as the correlation between principal seismic activities and nodal planes of focal mechanisms. The research results indicate that the moderate seismic activities prior to strong earthquakes display two significant stages of stress concentration and stress weakening, which are consistent to the evolvement of local tectonic stress field. The first stage, defined as medium-term stage, is characterized by seismic belts, seismic gaps and centralized activities; The second stage, defined as short-term stage, is predominated by continuous quiescence. For the earthquakes occurred on the strike-slip and normal faults, the two stages are relatively evident; while for the earthquakes occurred on the thrust faults, there are not characteristic images and seismic quiescence in both stages. The spatial distribution of moderate shock belts before strong earthquakes show a good coherence with that of the two nodal planes of focal mechanisms, while the spatial distribution of seismic gaps enclosed by moderate shocks display a poor coherence with the nodal planes.
文摘A network communication and recording system based on China-made ARCA SOC and embedded Linux operating system is introduced in this paper. It supports TCP/IP network communication protocol and mass storage medium. It has strong points of self-monitor, low power consumption, high timing accuracy, high reliability of operation, etc. It can serve up to 20 centers real-time waveform data at the same time. It meets not only the requirements of physical networking observation, but also virtual networking observation based on Internet in which real-time data transmission is required. Its ability of field recording also meets the requirements of portable seismic observation, strong motion observation and seismic exploration observation, etc.
文摘According to the rupture dynamics of earthquakes, variations of the apparent stress and the difference between the static stress drop and the dynamic stress drop during the rupture of earthquakes are analyzed for the July 20, 1995 ML=4.1 Shacheng, Hebei, China, earthquake sequence. Results obtained show that the apparent stress for main-shock is about 5 MPa, and the average apparent stress for aftershocks 0.047 MPa. During the rupture of the main-shock, the dynamic stress drop is approximately 1.6 times greater than the static stress drop with the difference of nearly 2.7 MPa. The dynamic stress drop is less than the static stress drop for all aftershocks with the average difference of -0.75 MPa. Therefore, when the mainshock occurs the final stress on the focal fault is higher than the dynamic frictional stress, corresponding to that the fault is abruptly locked. When the aftershocks occur the final stress on the focal fault is lower than the dynamic frictional stress, corresponding to that the fault overshoots. It can be seen from the above results that there could be some differences in the physic processes between the mainshock and the aftershocks.
基金State Key Basic Research Development and Programming Project of China (G19980407) and Social Commonweal Research Project of the Ministry of Science and Technology (2002DIA10001).
文摘According to the deficiency of the strain accumulating and releasing curves and the previous models, the strain-accumulating rate of the strain accumulating and releasing model has been deduced based on the G-R relation and the empirical formula between energy release and earthquake magnitude, where the strain-accumulating rate is relative independent of the strain-releasing rate. Five typical areas in Chinese mainland are selected on the basis of the hypothesis on active tectonic block, and small earthquakes from 1970 are imported to calculate the annual strain-accumulating rates considering the completeness of historical seismic data. Having introduced the strain-accumulating rates into the amended model, present strain phases are got. According to the present stages in their own cycles, the future earthquake tendency of each sub-region is discussed.
基金Joint Seismological Science Foundation of China (101086) and the key project "Digital Crustal and Mantle Structure of Chinese Mainland" from China Earthquake Administration.
文摘The data of short-period (1~18 s) surface waves recorded by 23 stations belonging to the digital seismic network of Yunnan Province of China are used in this paper. From these data, the dispersion curves of phase velocities of the fundamental mode Rayleigh wave along 209 paths are determined by using the two-station narrowband filtering cross-correlation method. Adopting tomography method, the distribution maps of phase velocities at various peri-ods in Yunnan region are inverted. The maps of phase velocities on profiles along 24N, 25N, 26N, 27N and 100.5E and the distribution maps of phase velocities at 3 periods in the study region are given. The results show that the phase velocity distribution in Yunnan region has strong variations in horizontal direction, and the phase velocity distribution in short-period range is closely related to the thickness of sedimentary layers in the shallow crust. The phase velocity in southern part of the Sichuan-Yunnan rhombic block encircled by the Honghe fault and Xiaojiang fault is obviously lower than that in surrounding areas. The epicentral locations of strong earthquakes in Yunnan region are mainly distributed in transitional zones between low and high phase velocities.
基金State Natural Sciences Foundation of China (No. 40174012) and the Key Laboratory Foundation of the CNPC (No. GPKL0104).
文摘The Biot and Squirt-flow are the two most important mechanisms of fluid flow in the porous medium with fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, dispersion and attenuation of elastic waves in the porous medium are widely investigated in recent years. However, we have not read any reports on numerical simulation based on the BISQ equation. In this paper, following the BISQ equation, elastic wave propagation in the transversely isotropic porous medium filled with fluids is simulated by the stag-gered grid method for different frequency and phase boundary cases and the two-layer medium. And propagating characteristics of seismic and acoustic waves and various phenomena occured in the propagating process are in-vestigated when the two mechanisms are considered simultaneously.
文摘We developed a GA-BP algorithm by combining the genetic algorithm (GA) with the back propagation (BP) algorithm and established a genetic BP neural network. We also applied the BP neural network based on the BP algorithm and the genetic BP neural network based on the GA-BP algorithm to discriminate earthquakes and explosions. The obtained result shows that the discriminating performance of the genetic BP network is slightly better than that of the BP network.
