In this study,the vertical components of broadband teleseismic P wave data recorded by China Earthquake Network are used to image the rupture processes of the February 6th,2023 Turkish earthquake doublet via back proj...In this study,the vertical components of broadband teleseismic P wave data recorded by China Earthquake Network are used to image the rupture processes of the February 6th,2023 Turkish earthquake doublet via back projection analysis.Data in two frequency bands(0.5-2 Hz and 1-3 Hz)are used in the imaging processes.The results show that the rupture of the first event extends about 200 km to the northeast and about 150 km to the southwest,lasting~90 s in total.The southwestern rupture is triggered by the northeastern rupture,demonstrating a sequential bidirectional unilateral rupture pattern.The rupture of the second event extends approximately 80 km in both northeast and west directions,lasting~35 s in total and demonstrates a typical bilateral rupture feature.The cascading ruptures on both sides also reflect the occurrence of selective rupture behaviors on bifurcated faults.In addition,we observe super-shear ruptures on certain fault sections with relatively straight fault structures and sparse aftershocks.展开更多
Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a j...Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a joint inversion of receiver functions and surface wave dispersions with P-wave velocity constraints using data from the Chin Array Ⅱ temporary stations deployed across the Qinghai-Xizang Plateau. Prior to joint inversion, we applied the H-κ-c method(Li JT et al., 2019) to the receiver function data in order to correct for the back-azimuthal variations in the arrival times of Ps phases and crustal multiples caused by crustal anisotropy and dipping interfaces. High-resolution images of vS, crustal thickness, and vP/vSstructures in the Qinghai-Xizang Plateau were simultaneously derived from the joint inversion. The seismic images reveal that crustal thickness decreases outward from the Qinghai-Xizang Plateau. The stable interiors of the Ordos and Alxa blocks exhibited higher velocities and lower crustal vP/vSratios. While, lower velocities and higher vP/vSratios were observed beneath the Qilian Orogen and Songpan-Ganzi terrane(SPGZ), which are geologically active and mechanically weak, especially in the mid-lower crust.Delamination or thermal erosion of the lithosphere triggered by hot asthenospheric flow contributes to the observed uppermost mantle low-velocity zones(LVZs) in the SPGZ. The crustal thickness, vS, and vP/vSratios suggest that whole lithospheric shortening is a plausible mechanism for crustal thickening in the Qinghai-Xizang Plateau, supporting the idea of coupled lithospheric-scale deformation in this region.展开更多
Mount Semeru,an active volcano in East Java,Indonesia,erupted on December 4,2021,following extreme rainfall that caused an avalanche of hot pyroclastic flows and lava.The tropospheric conditions and dominant particle ...Mount Semeru,an active volcano in East Java,Indonesia,erupted on December 4,2021,following extreme rainfall that caused an avalanche of hot pyroclastic flows and lava.The tropospheric conditions and dominant particle components in the atmosphere can be monitored using Global Navigation Satellite System(GNSS)technology and remote sensing satellites.GNSS signal propagation delay in Precise Point Positioning(PPP)processing can be used to determine Zenith Tropospheric Delay(ZTD)and Precipitable Water Vapor(PWV)variables so that atmospheric conditions can be generated.In addition,by using remote sensing satellite data,it is possible to obtain rainfall data with high temporal resolution as well as the dominant particle and gas content values during eruptions.During the eruption period,the high value of PWV was dominated by the high intensity of precipitation during the rainy season.High rainfall before the eruption caused activity inside the mountain to increase,which occurred in avalanche type eruption.Apart from that,the atmosphere around Semeru was also dominated by SO_(2)content,which spreaded for tens of kilometers.SO_(2)content began to be detected significantly by remote sensing sensors on December 7,2021.In this study,deformation and atmospheric monitoring were also carried out using low-cost GNSS at the Semeru Monitoring Station on September 9-15,2022.The results of the ZTD and ZWD values show the dominance of the wet component,which is directly proportional to rainfall activity in this period.展开更多
Magnetosonic waves are an important medium for energy transfer in collisionless space plasma.Magnetosonic waves have been widely investigated in the upstream of the bow shock at Mars.These waves are believed to origin...Magnetosonic waves are an important medium for energy transfer in collisionless space plasma.Magnetosonic waves have been widely investigated in the upstream of the bow shock at Mars.These waves are believed to originate from pickup ions or reflected particles.By utilizing MAVEN spacecraft data,we have observed the occurrence of quasi-perpendicularly propagating magnetosonic emissions near the proton gyrofrequency in the Martian magnetotail region.These plasma waves are associated with a significant enhancement of proton and oxygen flux.The excited magnetosonic waves could possibly heat the protons through resonance and facilitate the ionospheric plasma escape.Our results could be helpful to better understand the Mars’magnetospheric dynamics and offer insights into possible energy redistribution between waves and plasma in the Martian nightside magnetosphere.展开更多
This research presents the variation of the gravity field and associated gravity field components over the continental area of Nigeria to provide data for geoscience research,geodetic and engineering works,aerodynamic...This research presents the variation of the gravity field and associated gravity field components over the continental area of Nigeria to provide data for geoscience research,geodetic and engineering works,aerodynamic studies and deep crustal inferences.Accurate positions and elevations were observed at 58 of the 59 base stations of the Primary Gravity Network of Nigeria(PGNN),whose absolute gravity values had been accurately determined.The absolute gravity values were plotted against their respective positions to reveal the distribution pattern and strength of the gravity field within the study area.Theoretical gravity values at each base station were generated using the Somigliana's equation.The free-air gravity and free-air anomaly gravity values were generated with respect to the World Geodetic System 1984(WGS84)ellipsoid using GPS-derived elevation data.Then,the perturbing potential,free-air gravity with respect to the geoid,and the indirect effects were evaluated.The average of the indirect effects was used to adjust the WGS84 gravity formula to produce a gravity formula that better approximates the geoid across the continental area of Nigeria,compatible with the heights measured relative to the geoid,which can serve as a reference for establishing a vertical height control.The Bouguer gravity and Bouguer gravity anomalies across Nigeria revealed a“trans-southern gravity high strip”interpreted to be associated with mantle upwelling.Two new major mega-lineaments related to mantle upwelling were mapped.A batholith province trending NWeSE was delineated,occurring from north central Nigeria to the north western region and containing closures of“Bouguer gravity lows”interpreted as batholiths.A separate closure of“Bouguer gravity low”was detected at Azare,north eastern Nigeria,which may be due to the presence of intrusive granitic body.It is recommended that the mantle structure beneath“the trans-southern gravity high strip”,“delineated batholith province”and“isolated gravity closures”around the northeast of Nigeria should be studied from seismic shear wave splitting analysis for better understanding of the deep lithospheric structures and moho relief.展开更多
Pn velocity lateral variation and anisotropy images were reconstructed by adding about 50 000 travel times from the regional seismic networks to the datum set of near 40 000 travel times from National Seismic Network ...Pn velocity lateral variation and anisotropy images were reconstructed by adding about 50 000 travel times from the regional seismic networks to the datum set of near 40 000 travel times from National Seismic Network of China used by WANG, et al. We discussed the relation of Pn velocity variation to Moho depth, Earths heat flow, distribution of Cenozoic volcanic rock and the result of rock experiment under high pressure and high temperature. The result of quantitative analysis indicates that Pn velocity is positively correlated with the crust thickness and negatively correlated with the Earths heat flow. Two linear regression equations, one between Pn velocity and crust thickness, and the other between Pn velocity and heat flow, were obtained. The rate of variation of Pn veloc-ity vP with pressure P, Pv/p, estimated from the velocity variation with crust thickness Hv/p, is close to the result obtained from the rock experiment under high pressure and high temperature. If the effect of crust thick-ness on Pn velocity is deducted from the velocity variation, then the low Pn velocity beneath Qinghai-Xizang pla-teau is more notable. The low Pn velocity regions well agree with the Cenozoic volcanic rock. In the several re-gions with significant anisotropy, the direction of fast Pn velocity is consistent with the orientation of maximum principal crustal compressive stress, and also with the direction of present-day crustal movement. It indicates that the fast Pn velocity direction may be related to the deformation or flow of top mantle material along the direction of maximum pressure.展开更多
Based on the discrete wavenumber method, we calculate the fields of dynamic Coulomb rupture stress changes and static stress changes caused by M6.5 earthquake in Wuding, and study their relationship with the subsequen...Based on the discrete wavenumber method, we calculate the fields of dynamic Coulomb rupture stress changes and static stress changes caused by M6.5 earthquake in Wuding, and study their relationship with the subsequent after- shocks. The results show that the spatial distribution patterns of the positive region of dynamic stress peak value and static stress peak value are similarly asymmetric, which are basically identical with distribution features of aftershock. The dynamic stress peak value and the static stress in the positive region are more than 0.1 MPa and 0.01 MPa of the triggering threshold, respectively, which indicates that the dynamic and static stresses are helpful for the occurrence of aftershock. This suggests that both influences of dynamic and static stresses should be con- sidered other than only either of them when studying aftershock triggering in near field.展开更多
Two central schemes of finite difference (FD) up to different accuracy orders of space sampling step Dx (Fourth order and Sixth order respectively) were used to study the 1-D nonlinear P-wave propagation in the nonlin...Two central schemes of finite difference (FD) up to different accuracy orders of space sampling step Dx (Fourth order and Sixth order respectively) were used to study the 1-D nonlinear P-wave propagation in the nonlinear solid media by the numerical method. Distinctly different from the case of numerical modeling of linear elastic wave, there may be several difficulties in the numerical treatment to the nonlinear partial differential equation, such as the steep gradients, shocks and unphysical oscillations. All of them are the great obstacles to the stability and conver-gence of numerical calculation. Fortunately, the comparative study on the modeling of nonlinear wave by the two FD schemes presented in the paper can provide us with an easy method to keep the stability and convergence in the calculation field when the product of the absolute value of nonlinear coefficient and the value of u/x are small enough, namely, the value of bu/x is much smaller than 1. Several results are founded in the numerical study of nonlinear P-wave propagation, such as the waveform aberration, the generation and growth of harmonic wave and the energy redistribution among different frequency components. All of them will be more violent when the initial amplitude A0 is larger or the nonlinearity of medium is stronger. Correspondingly, we have found that the nonlinear P-wave propagation velocity will change with different initial frequency f of source wave or the wave velocity c (equal to the P-wave velocity in the same medium without considering nonlinearity).展开更多
In recent twenty years, much numerical simulation work has been done on the evolution of Qinghai-Xizang (Tibetan) plateau. In this paper some principal numerical models and results are reviewed and analyzed. The earli...In recent twenty years, much numerical simulation work has been done on the evolution of Qinghai-Xizang (Tibetan) plateau. In this paper some principal numerical models and results are reviewed and analyzed. The earlier plane stress or plane strain model has much discrepancy with the actual deformation of Qinghai-Xizang plateau, such as the thickening of Tibetan crust and the lateral extrusion of Tibet along strike-slip faults. The thin viscous sheet model and the thin-plate model may simulate the change of the crustal thickness and the deformation pro-duced by gravitational force. It is suitable for studying the large-scale and long-time deformation. The influence of faults on the deformation of Tibetan plateau should be further studied.展开更多
Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by ...Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by the displacement or velocity field from the earthquake under sea floor,usually no interaction between them is consid-ered in pure liquid model.In this study,the potential flow theory and the finite element method with the interaction between liquid and solid are employed to model the dynamic processes of the earthquake and tsunami.For model-ing the earthquake,firstly the initial stress field to generate the earthquake is set up,and then the occurrence of the earthquake is simulated by suddenly reducing the elastic material parameters inside the earthquake fault.It is dif-ferent from seismic dislocation theory in which the relative slip on the fault is specified in advance.The modeling results reveal that P,SP and the surface wave can be found at the sea surface besides the tsunami wave.The surface wave arrives at the distance of 600 km from the epicenter earlier than the tsunami 48 minutes,and its maximum amplitude is 0.55 m,which is 2 times as large as that of the sea floor.Tsunami warning information can be taken from the surface wave on the sea surface,which is much earlier than that obtained from the seismograph stations on land.The tsunami speed on the open sea with 3 km depth is 175.8 m/s,which is a little greater than that pre-dicted by long wave theory,(gh)1/2=171.5 m,and its wavelength and amplitude in average are 32 km and 2 m,respectively.After the tsunami propagates to the continental shelf,its speed and wavelength is reduced,but its amplitude become greater,especially,it can elevate up to 10 m and run 55 m forward in vertical and horizontal directions at sea shore,respectively.The maximum vertical accelerations at the epicenter on the sea surface and on the earthquake fault are 5.9 m/s2 and 16.5 m/s2,respectively,the later is 2.8 times the former,and therefore,sea water is a good shock absorber.The acceleration at the sea shore is about 1/10 as large as at the epicenter.The maximum vertical velocity at the epicenter is 1.4 times that on the fault.The maximum vertical displacement at the fault is less than that at the epicenter.The difference between them is the amplitude of the tsunami at the epicenter.The time of the maximum displacement to occur on the fault is not at the beginning of the fault slipping but retards 23 s.展开更多
In large continental orogens, an important research topic is the behavior of deep crustal and upper mantle deformation, and the flow styles of ductile material. The morphology of the eastern margin of the Tibetan plat...In large continental orogens, an important research topic is the behavior of deep crustal and upper mantle deformation, and the flow styles of ductile material. The morphology of the eastern margin of the Tibetan plateau, adjacent to the Sichuan basin, is characterized by very steep relief with high mountain ranges. The crust beneath this region slows the velocities in the middle and lower crust. We have adopted a relatively dense network to inverse the detailed structure of the crust and upper mantle along the eastern margin of the Tibetan plateau and Sichuan basin, using teleseismic data via receiver function analysis. The results are in-line with the hypothesis that viscous crustal material is flowing beneath the eastern margin of the Tibetan plateau and that this process drives overlying crustal material around the strong and rigid Sichuan basin. When the viscous material hits this obstruction, flows are divided into two or more branches with different directions. The upper part of the upwelling viscous flow produces the pressure to intrude the upper crust, thereby driving uplift of mountain ranges and high peaks. In contrast, the lower part of the downwelling viscous flow produces the pressure to intrude the lower crust and upper mantle to deepen the Moho discontinuity, causing observed crustal thickening.展开更多
Compared with the laboratory acoustic measurement of sediment samples,the in-situ acoustic measurement in marine sediment is considered more accurate and reliable,because it covers all of the surrounding environment f...Compared with the laboratory acoustic measurement of sediment samples,the in-situ acoustic measurement in marine sediment is considered more accurate and reliable,because it covers all of the surrounding environment factors and avoids the disturbance during the course of sampling and transporting of sediment samples.A new multi-frequency in-situ geo-acoustic measurement system (MFISGAMS) has been developed.The system can provide acoustic velocity (compressional wave) and attenuation profiles of the uppermost 4~8 m sediment in the seafloor.It consists of 8 channels with 12 frequencies (multi-frequencies) and 0.5~2 MHz sampling rates.The data collected can be transmiuted in real-time.Associated with inclinometer and altimeter,it can provide the data for depth emendation.Acoustic velocity and attenuation data have been obtained from two in-situ experiments conducted in the Hangzhou Bay.展开更多
We present an equivalent form of the expressions first obtained by Tada(Geophys J Int 164:653–669,2006. doi:10.1111/j.1365-246 X.2006.03868.x), which represents the transient stress response of an infinite, homogeneo...We present an equivalent form of the expressions first obtained by Tada(Geophys J Int 164:653–669,2006. doi:10.1111/j.1365-246 X.2006.03868.x), which represents the transient stress response of an infinite, homogeneous and isotropic medium to a constant slip rate on a triangular fault that continues perpetually after the slip onset. Our results are simpler than Tada's, and the corresponding codes have a higher running speed.展开更多
This is the first of two papers that describes a regional tomography investigation,which combines P-wave arrival times of both regional and teleseismic earthquakes to obtain 3D mantle structures of East Asia up to 1 0...This is the first of two papers that describes a regional tomography investigation,which combines P-wave arrival times of both regional and teleseismic earthquakes to obtain 3D mantle structures of East Asia up to 1 000 km depth.The most important findings of this tomography study are reported in this paper as follows.(1) No fast P-wave velocity anomalies can be related to subducted oceanic slabs beneath the 660 km discontinuity;instead the subducted oceanic slabs become flattened and stagnant within the transition zone.(2) The high velocity anomalies in the transition zone extend up to 1 500 km to the westward of the active trenches,which is a unique feature in the worldwide subduction systems.(3) Slow P-wave velocity anomalies are visible up to ~250 km underneath most of the East Asia on the east of 115°E,similar to the area of the stagnant slabs.These observations have important implications for the geodynamic process at depths beneath the East Asia,which might in turn control the widespread Cenozoic volcanism and associated extensional tectonics seen at the Earth's surface.展开更多
Nonlinear response of the driven Duffng oscillator to periodic or quasi-periodic signals has been well studied.In this paper,we investigate the nonlinear response of the driven Duffng oscillator to non-periodic,more s...Nonlinear response of the driven Duffng oscillator to periodic or quasi-periodic signals has been well studied.In this paper,we investigate the nonlinear response of the driven Duffng oscillator to non-periodic,more specifically,chaotic time series.Through numerical simulations,we find that the driven Duffng oscillator can also show regular nonlinear response to the chaotic time series with different degree of chaos as generated by the same chaotic series generating model,and there exists a relationship between the state of the driven Duffng oscillator and the chaoticity of the input signal of the driven Duffng oscillator.One real-world and two artificial chaotic time series are used to verify the new feature of Duffng oscillator.A potential application of the new feature of Duffng oscillator is also indicated.展开更多
The firework algorithm(FWA) is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model s...The firework algorithm(FWA) is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model show that the FWA performs well in both solution quality and efficiency. We apply the FWA in this study to crustal velocity structure inversion using regional seismic waveform data of central Gansu on the northeastern margin of the Qinghai-Tibet plateau. Seismograms recorded from the moment magnitude(MW) 5.4 Minxian earthquake enable obtaining an average crustal velocity model for this region. We initially carried out a series of FWA robustness tests in regional waveform inversion at the same earthquake and station positions across the study region,inverting two velocity structure models, with and without a low-velocity crustal layer; the accuracy of our average inversion results and their standard deviations reveal the advantages of the FWA for the inversion of regional seismic waveforms. We applied the FWA across our study area using three component waveform data recorded by nine broadband permanent seismic stations with epicentral distances ranging between 146 and 437 km. These inversion results show that the average thickness of the crust in this region is 46.75 km, while thicknesses of the sedimentary layer, and the upper, middle, and lower crust are 3.15,15.69, 13.08, and 14.83 km, respectively. Results also show that the P-wave velocities of these layers and the upper mantle are 4.47, 6.07, 6.12, 6.87, and 8.18 km/s,respectively.展开更多
We presented a boundary element method using the approximate analytical Green's function given by Sánchez-Sesma et al.Coordinate transform is introduced to extend the method to deal with the model with consta...We presented a boundary element method using the approximate analytical Green's function given by Sánchez-Sesma et al.Coordinate transform is introduced to extend the method to deal with the model with constant-gradient velocity along oblique direction.The method is validated by comparing the numerical results with other independent methods.This method provides a useful tool for analyzing local site effects.We computed seismic response for two series of models.The results in both frequency and time domains are analyzed and show complex amplification patterns.The fundamental mode of resonance is dependent not only on the velocity at the free surface but also on the velocity distribution of the whole basin.For the higher modes of vibration the heterogeneous basin also has its own characteristic.展开更多
Scientists pay great attention to different-time-scale signals in the length of day(LOD)variations DLOD, which provide signatures of the Earth's interior structure, couplings among different layers, and potential ...Scientists pay great attention to different-time-scale signals in the length of day(LOD)variations DLOD, which provide signatures of the Earth's interior structure, couplings among different layers, and potential excitations of ocean and atmosphere. In this study,based on the ensemble empirical mode decomposition(EEMD), we analyzed the latest time series of DLOD data spanning from January 1962 to March 2015. We observed the signals with periods and amplitudes of about 0.5 month and 0.19 ms, 1.0 month and 0.19 ms, 0.5 yr and 0.22 ms, 1.0 yr and 0.18 ms, 2.28 yr and 0.03 ms, 5.48 yr and 0.05 ms, respectively, in coincidence with the results of predecessors. In addition, some signals that were previously not definitely observed by predecessors were detected in this study, with periods and amplitudes of 9.13 d and 0.12 ms, 13.69 yr and 0.10 ms, respectively. The mechanisms of the LOD fluctuations of these two signals are still open.展开更多
基金supported by the National Key R&D Program of China(No.2022YFF0800601)National Scientific Foundation of China(Nos.41930103 and 41774047).
文摘In this study,the vertical components of broadband teleseismic P wave data recorded by China Earthquake Network are used to image the rupture processes of the February 6th,2023 Turkish earthquake doublet via back projection analysis.Data in two frequency bands(0.5-2 Hz and 1-3 Hz)are used in the imaging processes.The results show that the rupture of the first event extends about 200 km to the northeast and about 150 km to the southwest,lasting~90 s in total.The southwestern rupture is triggered by the northeastern rupture,demonstrating a sequential bidirectional unilateral rupture pattern.The rupture of the second event extends approximately 80 km in both northeast and west directions,lasting~35 s in total and demonstrates a typical bilateral rupture feature.The cascading ruptures on both sides also reflect the occurrence of selective rupture behaviors on bifurcated faults.In addition,we observe super-shear ruptures on certain fault sections with relatively straight fault structures and sparse aftershocks.
基金supported by the Natural Science Basic Research Program of Shaanxi(No.2023-JC-QN-0306)the Special Fund of the Institute of Geophysics,China Earthquake Administration(No.DQJB21B32)the National Natural Science Foundation of China(No.42174069).
文摘Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a joint inversion of receiver functions and surface wave dispersions with P-wave velocity constraints using data from the Chin Array Ⅱ temporary stations deployed across the Qinghai-Xizang Plateau. Prior to joint inversion, we applied the H-κ-c method(Li JT et al., 2019) to the receiver function data in order to correct for the back-azimuthal variations in the arrival times of Ps phases and crustal multiples caused by crustal anisotropy and dipping interfaces. High-resolution images of vS, crustal thickness, and vP/vSstructures in the Qinghai-Xizang Plateau were simultaneously derived from the joint inversion. The seismic images reveal that crustal thickness decreases outward from the Qinghai-Xizang Plateau. The stable interiors of the Ordos and Alxa blocks exhibited higher velocities and lower crustal vP/vSratios. While, lower velocities and higher vP/vSratios were observed beneath the Qilian Orogen and Songpan-Ganzi terrane(SPGZ), which are geologically active and mechanically weak, especially in the mid-lower crust.Delamination or thermal erosion of the lithosphere triggered by hot asthenospheric flow contributes to the observed uppermost mantle low-velocity zones(LVZs) in the SPGZ. The crustal thickness, vS, and vP/vSratios suggest that whole lithospheric shortening is a plausible mechanism for crustal thickening in the Qinghai-Xizang Plateau, supporting the idea of coupled lithospheric-scale deformation in this region.
基金the Indonesian Collaborative Research(RKI)2022:598/IT2/T/HK.00.01/2022Center of Volcanology and Geological Hazard Mitigation of Indonesia for the data and financial support of this research。
文摘Mount Semeru,an active volcano in East Java,Indonesia,erupted on December 4,2021,following extreme rainfall that caused an avalanche of hot pyroclastic flows and lava.The tropospheric conditions and dominant particle components in the atmosphere can be monitored using Global Navigation Satellite System(GNSS)technology and remote sensing satellites.GNSS signal propagation delay in Precise Point Positioning(PPP)processing can be used to determine Zenith Tropospheric Delay(ZTD)and Precipitable Water Vapor(PWV)variables so that atmospheric conditions can be generated.In addition,by using remote sensing satellite data,it is possible to obtain rainfall data with high temporal resolution as well as the dominant particle and gas content values during eruptions.During the eruption period,the high value of PWV was dominated by the high intensity of precipitation during the rainy season.High rainfall before the eruption caused activity inside the mountain to increase,which occurred in avalanche type eruption.Apart from that,the atmosphere around Semeru was also dominated by SO_(2)content,which spreaded for tens of kilometers.SO_(2)content began to be detected significantly by remote sensing sensors on December 7,2021.In this study,deformation and atmospheric monitoring were also carried out using low-cost GNSS at the Semeru Monitoring Station on September 9-15,2022.The results of the ZTD and ZWD values show the dominance of the wet component,which is directly proportional to rainfall activity in this period.
基金the National Natural Science Foundation of China(42030101,42204158)Shanghai Pujiang Program(No.21PJD078)+1 种基金Shanghai Science and Technology Innovation Action Plan(No.21DZ1206100)the Fundamental Research Funds for the Central Universities.
文摘Magnetosonic waves are an important medium for energy transfer in collisionless space plasma.Magnetosonic waves have been widely investigated in the upstream of the bow shock at Mars.These waves are believed to originate from pickup ions or reflected particles.By utilizing MAVEN spacecraft data,we have observed the occurrence of quasi-perpendicularly propagating magnetosonic emissions near the proton gyrofrequency in the Martian magnetotail region.These plasma waves are associated with a significant enhancement of proton and oxygen flux.The excited magnetosonic waves could possibly heat the protons through resonance and facilitate the ionospheric plasma escape.Our results could be helpful to better understand the Mars’magnetospheric dynamics and offer insights into possible energy redistribution between waves and plasma in the Martian nightside magnetosphere.
文摘This research presents the variation of the gravity field and associated gravity field components over the continental area of Nigeria to provide data for geoscience research,geodetic and engineering works,aerodynamic studies and deep crustal inferences.Accurate positions and elevations were observed at 58 of the 59 base stations of the Primary Gravity Network of Nigeria(PGNN),whose absolute gravity values had been accurately determined.The absolute gravity values were plotted against their respective positions to reveal the distribution pattern and strength of the gravity field within the study area.Theoretical gravity values at each base station were generated using the Somigliana's equation.The free-air gravity and free-air anomaly gravity values were generated with respect to the World Geodetic System 1984(WGS84)ellipsoid using GPS-derived elevation data.Then,the perturbing potential,free-air gravity with respect to the geoid,and the indirect effects were evaluated.The average of the indirect effects was used to adjust the WGS84 gravity formula to produce a gravity formula that better approximates the geoid across the continental area of Nigeria,compatible with the heights measured relative to the geoid,which can serve as a reference for establishing a vertical height control.The Bouguer gravity and Bouguer gravity anomalies across Nigeria revealed a“trans-southern gravity high strip”interpreted to be associated with mantle upwelling.Two new major mega-lineaments related to mantle upwelling were mapped.A batholith province trending NWeSE was delineated,occurring from north central Nigeria to the north western region and containing closures of“Bouguer gravity lows”interpreted as batholiths.A separate closure of“Bouguer gravity low”was detected at Azare,north eastern Nigeria,which may be due to the presence of intrusive granitic body.It is recommended that the mantle structure beneath“the trans-southern gravity high strip”,“delineated batholith province”and“isolated gravity closures”around the northeast of Nigeria should be studied from seismic shear wave splitting analysis for better understanding of the deep lithospheric structures and moho relief.
基金State Key Basic Research Project of Development and Programming Mechanism and Prediction of Continental Strong Earthquakes (G1998040700).
文摘Pn velocity lateral variation and anisotropy images were reconstructed by adding about 50 000 travel times from the regional seismic networks to the datum set of near 40 000 travel times from National Seismic Network of China used by WANG, et al. We discussed the relation of Pn velocity variation to Moho depth, Earths heat flow, distribution of Cenozoic volcanic rock and the result of rock experiment under high pressure and high temperature. The result of quantitative analysis indicates that Pn velocity is positively correlated with the crust thickness and negatively correlated with the Earths heat flow. Two linear regression equations, one between Pn velocity and crust thickness, and the other between Pn velocity and heat flow, were obtained. The rate of variation of Pn veloc-ity vP with pressure P, Pv/p, estimated from the velocity variation with crust thickness Hv/p, is close to the result obtained from the rock experiment under high pressure and high temperature. If the effect of crust thick-ness on Pn velocity is deducted from the velocity variation, then the low Pn velocity beneath Qinghai-Xizang pla-teau is more notable. The low Pn velocity regions well agree with the Cenozoic volcanic rock. In the several re-gions with significant anisotropy, the direction of fast Pn velocity is consistent with the orientation of maximum principal crustal compressive stress, and also with the direction of present-day crustal movement. It indicates that the fast Pn velocity direction may be related to the deformation or flow of top mantle material along the direction of maximum pressure.
文摘Based on the discrete wavenumber method, we calculate the fields of dynamic Coulomb rupture stress changes and static stress changes caused by M6.5 earthquake in Wuding, and study their relationship with the subsequent after- shocks. The results show that the spatial distribution patterns of the positive region of dynamic stress peak value and static stress peak value are similarly asymmetric, which are basically identical with distribution features of aftershock. The dynamic stress peak value and the static stress in the positive region are more than 0.1 MPa and 0.01 MPa of the triggering threshold, respectively, which indicates that the dynamic and static stresses are helpful for the occurrence of aftershock. This suggests that both influences of dynamic and static stresses should be con- sidered other than only either of them when studying aftershock triggering in near field.
基金Project of Knowledge Innovation Program from Chinese Academy of Sciences (KZCX2-109).
文摘Two central schemes of finite difference (FD) up to different accuracy orders of space sampling step Dx (Fourth order and Sixth order respectively) were used to study the 1-D nonlinear P-wave propagation in the nonlinear solid media by the numerical method. Distinctly different from the case of numerical modeling of linear elastic wave, there may be several difficulties in the numerical treatment to the nonlinear partial differential equation, such as the steep gradients, shocks and unphysical oscillations. All of them are the great obstacles to the stability and conver-gence of numerical calculation. Fortunately, the comparative study on the modeling of nonlinear wave by the two FD schemes presented in the paper can provide us with an easy method to keep the stability and convergence in the calculation field when the product of the absolute value of nonlinear coefficient and the value of u/x are small enough, namely, the value of bu/x is much smaller than 1. Several results are founded in the numerical study of nonlinear P-wave propagation, such as the waveform aberration, the generation and growth of harmonic wave and the energy redistribution among different frequency components. All of them will be more violent when the initial amplitude A0 is larger or the nonlinearity of medium is stronger. Correspondingly, we have found that the nonlinear P-wave propagation velocity will change with different initial frequency f of source wave or the wave velocity c (equal to the P-wave velocity in the same medium without considering nonlinearity).
文摘In recent twenty years, much numerical simulation work has been done on the evolution of Qinghai-Xizang (Tibetan) plateau. In this paper some principal numerical models and results are reviewed and analyzed. The earlier plane stress or plane strain model has much discrepancy with the actual deformation of Qinghai-Xizang plateau, such as the thickening of Tibetan crust and the lateral extrusion of Tibet along strike-slip faults. The thin viscous sheet model and the thin-plate model may simulate the change of the crustal thickness and the deformation pro-duced by gravitational force. It is suitable for studying the large-scale and long-time deformation. The influence of faults on the deformation of Tibetan plateau should be further studied.
基金National Natural Science Foundation of China (40521002 and 40474013).
文摘Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by the displacement or velocity field from the earthquake under sea floor,usually no interaction between them is consid-ered in pure liquid model.In this study,the potential flow theory and the finite element method with the interaction between liquid and solid are employed to model the dynamic processes of the earthquake and tsunami.For model-ing the earthquake,firstly the initial stress field to generate the earthquake is set up,and then the occurrence of the earthquake is simulated by suddenly reducing the elastic material parameters inside the earthquake fault.It is dif-ferent from seismic dislocation theory in which the relative slip on the fault is specified in advance.The modeling results reveal that P,SP and the surface wave can be found at the sea surface besides the tsunami wave.The surface wave arrives at the distance of 600 km from the epicenter earlier than the tsunami 48 minutes,and its maximum amplitude is 0.55 m,which is 2 times as large as that of the sea floor.Tsunami warning information can be taken from the surface wave on the sea surface,which is much earlier than that obtained from the seismograph stations on land.The tsunami speed on the open sea with 3 km depth is 175.8 m/s,which is a little greater than that pre-dicted by long wave theory,(gh)1/2=171.5 m,and its wavelength and amplitude in average are 32 km and 2 m,respectively.After the tsunami propagates to the continental shelf,its speed and wavelength is reduced,but its amplitude become greater,especially,it can elevate up to 10 m and run 55 m forward in vertical and horizontal directions at sea shore,respectively.The maximum vertical accelerations at the epicenter on the sea surface and on the earthquake fault are 5.9 m/s2 and 16.5 m/s2,respectively,the later is 2.8 times the former,and therefore,sea water is a good shock absorber.The acceleration at the sea shore is about 1/10 as large as at the epicenter.The maximum vertical velocity at the epicenter is 1.4 times that on the fault.The maximum vertical displacement at the fault is less than that at the epicenter.The difference between them is the amplitude of the tsunami at the epicenter.The time of the maximum displacement to occur on the fault is not at the beginning of the fault slipping but retards 23 s.
基金partially supported by Grant-in-aid for Scientific Research(Kiban-B.11440134,Kiban-A 17204037) from Japan Society for the Promotion of Science and by some financial support from the Global Center of Excellence(G-COE) program of Tohoku University
基金supported by National Natural Science Foundation of China under grant Nos. 40839909,41074062, and 40674040
文摘In large continental orogens, an important research topic is the behavior of deep crustal and upper mantle deformation, and the flow styles of ductile material. The morphology of the eastern margin of the Tibetan plateau, adjacent to the Sichuan basin, is characterized by very steep relief with high mountain ranges. The crust beneath this region slows the velocities in the middle and lower crust. We have adopted a relatively dense network to inverse the detailed structure of the crust and upper mantle along the eastern margin of the Tibetan plateau and Sichuan basin, using teleseismic data via receiver function analysis. The results are in-line with the hypothesis that viscous crustal material is flowing beneath the eastern margin of the Tibetan plateau and that this process drives overlying crustal material around the strong and rigid Sichuan basin. When the viscous material hits this obstruction, flows are divided into two or more branches with different directions. The upper part of the upwelling viscous flow produces the pressure to intrude the upper crust, thereby driving uplift of mountain ranges and high peaks. In contrast, the lower part of the downwelling viscous flow produces the pressure to intrude the lower crust and upper mantle to deepen the Moho discontinuity, causing observed crustal thickening.
基金supported by the National High Technology Research and Development Program of China(863 Pro-gram, Grant No.2005AA615040)the National Natural Science Foundation of China(Grant No.40776038)+2 种基金Open Fund of Key Laboratory of Geo-detection(China University of Geosciences,Beijing),Ministry of Education (GrantNo. GDL0802)the Ocean Public Welfare Scientific Research Special Appropriation Project(Grant Nos.200805079and 200805005)the Open Fund of Key Laboratory of Submarine Geoscience,State Oceanic Administration(Grant No. KCSG0803)
文摘Compared with the laboratory acoustic measurement of sediment samples,the in-situ acoustic measurement in marine sediment is considered more accurate and reliable,because it covers all of the surrounding environment factors and avoids the disturbance during the course of sampling and transporting of sediment samples.A new multi-frequency in-situ geo-acoustic measurement system (MFISGAMS) has been developed.The system can provide acoustic velocity (compressional wave) and attenuation profiles of the uppermost 4~8 m sediment in the seafloor.It consists of 8 channels with 12 frequencies (multi-frequencies) and 0.5~2 MHz sampling rates.The data collected can be transmiuted in real-time.Associated with inclinometer and altimeter,it can provide the data for depth emendation.Acoustic velocity and attenuation data have been obtained from two in-situ experiments conducted in the Hangzhou Bay.
基金supported by the National Natural Science Foundation of China (Grant No. 41674050)MOST Grant (2012CB417301)
文摘We present an equivalent form of the expressions first obtained by Tada(Geophys J Int 164:653–669,2006. doi:10.1111/j.1365-246 X.2006.03868.x), which represents the transient stress response of an infinite, homogeneous and isotropic medium to a constant slip rate on a triangular fault that continues perpetually after the slip onset. Our results are simpler than Tada's, and the corresponding codes have a higher running speed.
基金grants(B-11440134,S-12002006)to Dapeng Zhao from the Japan Society for the Promotion of ScienceSupport for Shunping Pei came from a postdoct grant of Peking University+1 种基金supported by the Chinese Academy of Sciences(No.KZCX2-EW-QN102)the National Natural Science Foundation of China(Nos.41074041,90814002 and 40125011)
文摘This is the first of two papers that describes a regional tomography investigation,which combines P-wave arrival times of both regional and teleseismic earthquakes to obtain 3D mantle structures of East Asia up to 1 000 km depth.The most important findings of this tomography study are reported in this paper as follows.(1) No fast P-wave velocity anomalies can be related to subducted oceanic slabs beneath the 660 km discontinuity;instead the subducted oceanic slabs become flattened and stagnant within the transition zone.(2) The high velocity anomalies in the transition zone extend up to 1 500 km to the westward of the active trenches,which is a unique feature in the worldwide subduction systems.(3) Slow P-wave velocity anomalies are visible up to ~250 km underneath most of the East Asia on the east of 115°E,similar to the area of the stagnant slabs.These observations have important implications for the geodynamic process at depths beneath the East Asia,which might in turn control the widespread Cenozoic volcanism and associated extensional tectonics seen at the Earth's surface.
基金supported by the National Natural Science Foundation of China (Grant Nos 40574051 and 40774054)
文摘Nonlinear response of the driven Duffng oscillator to periodic or quasi-periodic signals has been well studied.In this paper,we investigate the nonlinear response of the driven Duffng oscillator to non-periodic,more specifically,chaotic time series.Through numerical simulations,we find that the driven Duffng oscillator can also show regular nonlinear response to the chaotic time series with different degree of chaos as generated by the same chaotic series generating model,and there exists a relationship between the state of the driven Duffng oscillator and the chaoticity of the input signal of the driven Duffng oscillator.One real-world and two artificial chaotic time series are used to verify the new feature of Duffng oscillator.A potential application of the new feature of Duffng oscillator is also indicated.
基金supported by the National Natural Science Foundation of China (No. 41174034)
文摘The firework algorithm(FWA) is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model show that the FWA performs well in both solution quality and efficiency. We apply the FWA in this study to crustal velocity structure inversion using regional seismic waveform data of central Gansu on the northeastern margin of the Qinghai-Tibet plateau. Seismograms recorded from the moment magnitude(MW) 5.4 Minxian earthquake enable obtaining an average crustal velocity model for this region. We initially carried out a series of FWA robustness tests in regional waveform inversion at the same earthquake and station positions across the study region,inverting two velocity structure models, with and without a low-velocity crustal layer; the accuracy of our average inversion results and their standard deviations reveal the advantages of the FWA for the inversion of regional seismic waveforms. We applied the FWA across our study area using three component waveform data recorded by nine broadband permanent seismic stations with epicentral distances ranging between 146 and 437 km. These inversion results show that the average thickness of the crust in this region is 46.75 km, while thicknesses of the sedimentary layer, and the upper, middle, and lower crust are 3.15,15.69, 13.08, and 14.83 km, respectively. Results also show that the P-wave velocities of these layers and the upper mantle are 4.47, 6.07, 6.12, 6.87, and 8.18 km/s,respectively.
基金supported by the National Science Foundation of China(Nos. D40444002 and D40521002)National Key Basic Research Program(No.2006CB705803)
文摘We presented a boundary element method using the approximate analytical Green's function given by Sánchez-Sesma et al.Coordinate transform is introduced to extend the method to deal with the model with constant-gradient velocity along oblique direction.The method is validated by comparing the numerical results with other independent methods.This method provides a useful tool for analyzing local site effects.We computed seismic response for two series of models.The results in both frequency and time domains are analyzed and show complex amplification patterns.The fundamental mode of resonance is dependent not only on the velocity at the free surface but also on the velocity distribution of the whole basin.For the higher modes of vibration the heterogeneous basin also has its own characteristic.
基金supported by National 973 Project China (2013CB733305)National Natural Science Foundation of China (NSFCs) (41174011,41429401,41210006,41128003,41021061)
文摘Scientists pay great attention to different-time-scale signals in the length of day(LOD)variations DLOD, which provide signatures of the Earth's interior structure, couplings among different layers, and potential excitations of ocean and atmosphere. In this study,based on the ensemble empirical mode decomposition(EEMD), we analyzed the latest time series of DLOD data spanning from January 1962 to March 2015. We observed the signals with periods and amplitudes of about 0.5 month and 0.19 ms, 1.0 month and 0.19 ms, 0.5 yr and 0.22 ms, 1.0 yr and 0.18 ms, 2.28 yr and 0.03 ms, 5.48 yr and 0.05 ms, respectively, in coincidence with the results of predecessors. In addition, some signals that were previously not definitely observed by predecessors were detected in this study, with periods and amplitudes of 9.13 d and 0.12 ms, 13.69 yr and 0.10 ms, respectively. The mechanisms of the LOD fluctuations of these two signals are still open.