Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

A traveltime reciprocity discrepancy in the Podvin & Lecomte time3d finite difference algorithm

We have found that the extensively used finite difference scheme time3d produces time fields which are 'asymmetric' in the sense that traveltimes computed to the right of the source are faster than traveltimes computed to the left. All finite difference schemes will, as they are approximations to the wave equation, to some extent fail to obey reciprocity perfectly. We show, however, that the errors in time3d may be significant—and unnecessarily large. An asymmetry in the point source initialization has been identified, and after correction time3d produces time fields with an improved reciprocity.     

The Earthquake Readiness Scale: The Development of a Valid and Reliable Unifactorial Measure

Preparedness is a key dependent variable in many studies examining people’s response to disasters such as earthquakes. A feature of many studies on this issue, however, is the lack of attention given to psychometric issues when constructing measures of preparedness. With regard to earthquake preparation, for example, many studies could be greatly improved by the use of a valid and reliable measure of preparedness. This research developed such a measure that assessed both low-level preparedness, such as having an emergency kit, and high-level preparedness, such as altering home structures to mitigate damage. Studies of Wellington (New Zealand) residents using two samples totalling n=652 showed that 23 items measuring these different aspects of earthquake preparation could be combined into a reliable, valid, unifactorial scale. This brief scale should have utility in multivariate studies of earthquake preparation, either as a dependent variable, where preparation is the outcome variable of primary interest, or as one of several independent variables, where preparation and other measures predict another outcome variable.     

基于有限差分法的抗滑桩计算机辅助设计

基于地基系数“m-m”法、“m-k”法、“k-k”法的原理,考虑桩顶和桩底边界条件以及桩在滑动面处位移、转角、弯矩和剪力的连续条件,可解得桩身各节点的位移和内力,提出了进行抗滑桩全桩内力计算的有限差分法。根据差分方程并用VB6.0编制了实用的计算程序,既可避免繁琐的查表计算,提高计算速度,又可提高计算精度,直观生动,真正实现了人机交互,在界面的引导下,设计人员可完成全部计算,并绘出内力图形和抗滑桩截面配筋图,使设计更方便快捷,该软件可以极大地提高生产效率,降低工程造价,从而实现抗滑桩的优化设计。最后采用上述方法对某滑坡的悬臂抗滑桩进行了设计与计算。     

Seismic Hazards Assessment for Radioactive Waste Disposal Sites in Regions of Low Seismic Activity

A comparative tectonic quiescence and lack of earthquakes make the stable centres of continents attractive for siting long-term radioactive waste storage facilities. The low rates of deformation in such regions, however, make it difficult to characterize their long-term seismotectonic behaviour, leading to uncertain estimates for the very low probability hazard estimates required by society. In an attempt to overcome the deficiency of both contemporary seismicity and paleoseismic data in central Canada, we have used earthquake histories from regions with similar seismotectonic characteristics from around the world. Substituting space for time, we estimate a long-term rate per 10⁶ km² of 0.004 magnitude ≥6 earthquakes per annum, of which 33–100% might rupture to the surface.     

Early Estimation of Seismic Hazard for Strong Earthquakes in Taiwan

A shakemap system providing rapid estimates of strong ground shaking could be useful for emergency response providers in a damaging earthquake. A hybrid procedure, which combines site-dependent ground motion prediction models and the limited observations of the Real-Time Digital stream output system (RTD system operated by Central Weather Bureau, CWB), was set up to provide a high-resolution shakemap in a near-real-time manner after damaging earthquakes in Taiwan. One of the main factors that affect the result of ground motion prediction analysis is the existence of site effects. The purpose of this paper is to investigate the local site effects and their influence in the ground shaking and then establish an early estimation procedure of potential hazard for damaging earthquakes. Based on the attenuation law, the site effects of each TSMIP station are discussed in terms of a bias function that is site and intensity-level dependent function. The standard deviation of the site-dependent ground motion prediction model can be significantly reduced. The nonlinear behavior of ground soil is automatically taken into account in the intensity-level dependent bias function. Both the PGA and the spectral acceleration are studied in this study. Based on the RTD data, event correctors are calculated and applied to precisely estimate the shakemap of damaging earthquakes for emergency response.