With co-seismic surface rupture slip displacements provided by the field observation for the 2001 MS8.1 West Kunlun Mountain Pass earthquake, this paper estimates the rupture speed on the main faulting segment with a ...With co-seismic surface rupture slip displacements provided by the field observation for the 2001 MS8.1 West Kunlun Mountain Pass earthquake, this paper estimates the rupture speed on the main faulting segment with a long straight fault trace on the surface based on a simple slip-weakening rupture model, in which the frictional overshoot or undershoot are involved in consideration of energy partition during the earthquake faulting. In contrast to the study of Bouchon and Vallée, in which the rupture propagation along the main fault could exceed the local shear-wave speed, perhaps reach the P-wave speed on a certain section of fault, our results show that, under a slip-weakening assumption combined with a frictional undershoot (partial stress drop model), average rupture speed should be equal to or less than the Rayleigh wave speed with a high seismic radiation efficiency, which is consistent with the result derived by waveform inversion and the result estimated from source stress field. Associated with the surface rupture mechanism, such as partial stress drop (frictional undershoot) associated with the apparent stress, an alternative rupture mechanism based on the slip-weakening model has also been discussed.展开更多
There are overshoot and undershoot phenomenon and end swing phenomenon in the cubic spline fitting in Hil- bert-Huang transform. The two problems influence data quality of the empirical mode decomposition seriously. T...There are overshoot and undershoot phenomenon and end swing phenomenon in the cubic spline fitting in Hil- bert-Huang transform. The two problems influence data quality of the empirical mode decomposition seriously. The cubic spline fitting has been analysed, and the causes of producing the overshoot and undershoot phenomenon and the end swing phenomenon have been pointed out in this paper. Two new methods of cubic spline fitting and sine spline fitting and the new technique of handling the end points of the original data curve can completely re- move the overshoot and undershoot phenomenon and the end swing phenomenon on the condition of unchanging original data, and have the advantages of the continuous fitting functions and its continuous one order derivative, the simple and convenient calculations, the small calculation amount and the easy work on it.展开更多
To suppress the overshoots and undershoots in the envelope fitting for empirical mode decomposition (EMD), an alternative cubic spline interpolation method without overshooting and undershooting is proposed. On the ...To suppress the overshoots and undershoots in the envelope fitting for empirical mode decomposition (EMD), an alternative cubic spline interpolation method without overshooting and undershooting is proposed. On the basis of the derived slope constraints of knots of a non-overshooting and non-undershooting cubic interpolant, together with "not-a-knot" conditions the cubic spline interpolants are constructed by replacing the requirement for equal second order derivatives at every knot with Brodlie' s derivative formula. Analysis and simulation experiments show that this approach can effectively avoid generating new extrema, shifting or exaggerating the existing ones in a signal, and thus significantly improve the decomposition performance of EMD.展开更多
The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially i...The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.展开更多
基金National Natural Science Foundation of China (40574022)the Central Public Affair Program Award for the Specific Academic Research (ZDJ2007-1)One-Hundred Individual Program of Chinese Academy of Sciences (99T3004WAZ)
文摘With co-seismic surface rupture slip displacements provided by the field observation for the 2001 MS8.1 West Kunlun Mountain Pass earthquake, this paper estimates the rupture speed on the main faulting segment with a long straight fault trace on the surface based on a simple slip-weakening rupture model, in which the frictional overshoot or undershoot are involved in consideration of energy partition during the earthquake faulting. In contrast to the study of Bouchon and Vallée, in which the rupture propagation along the main fault could exceed the local shear-wave speed, perhaps reach the P-wave speed on a certain section of fault, our results show that, under a slip-weakening assumption combined with a frictional undershoot (partial stress drop model), average rupture speed should be equal to or less than the Rayleigh wave speed with a high seismic radiation efficiency, which is consistent with the result derived by waveform inversion and the result estimated from source stress field. Associated with the surface rupture mechanism, such as partial stress drop (frictional undershoot) associated with the apparent stress, an alternative rupture mechanism based on the slip-weakening model has also been discussed.
基金The Foundation Research and Development Programs of China (2004CB418404).
文摘There are overshoot and undershoot phenomenon and end swing phenomenon in the cubic spline fitting in Hil- bert-Huang transform. The two problems influence data quality of the empirical mode decomposition seriously. The cubic spline fitting has been analysed, and the causes of producing the overshoot and undershoot phenomenon and the end swing phenomenon have been pointed out in this paper. Two new methods of cubic spline fitting and sine spline fitting and the new technique of handling the end points of the original data curve can completely re- move the overshoot and undershoot phenomenon and the end swing phenomenon on the condition of unchanging original data, and have the advantages of the continuous fitting functions and its continuous one order derivative, the simple and convenient calculations, the small calculation amount and the easy work on it.
基金the Ministerial Level Advanced Research Foundation (445030705QB0301)
文摘To suppress the overshoots and undershoots in the envelope fitting for empirical mode decomposition (EMD), an alternative cubic spline interpolation method without overshooting and undershooting is proposed. On the basis of the derived slope constraints of knots of a non-overshooting and non-undershooting cubic interpolant, together with "not-a-knot" conditions the cubic spline interpolants are constructed by replacing the requirement for equal second order derivatives at every knot with Brodlie' s derivative formula. Analysis and simulation experiments show that this approach can effectively avoid generating new extrema, shifting or exaggerating the existing ones in a signal, and thus significantly improve the decomposition performance of EMD.
基金the National Supercomputer Center in Tianjin for their patient assistance in providing the compilation environment.We thank the editor,Huajian Yao,for handling the manuscript and Mingming Li and another anonymous reviewer for their constructive comments.The research leading to these results has received funding from National Natural Science Foundation of China projects(Grant Nos.92355302 and 42121005)Taishan Scholar projects(Grant No.tspd20210305)others(Grant Nos.XDB0710000,L2324203,XK2023DXC001,LSKJ202204400,and ZR2021ZD09).
文摘The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.