Seismic rays in an ellipsoidal earth model: theoretical analysis and estimation of deviations due to the ellipticity

Summary. A first-order form of the Euler's equations for rays in an ellipsoidal model of the Earth is obtained. The conditions affecting the velocity law for a monotonic increase, with respect to the arc length, in the angular distance to the epicentre, and in the angle of incidence, are the same in the ellipsoidal and spherical models. It is therefore possible to trace rays and to compute travel times directly in an ellipsoidal earth as in the spherical model. Thus comparison with the rays of the same coordinates in a spherical earth provides an estimate of the various deviations of these rays due to the Earth's flattening, and the corresponding travel-time differences, for mantle P -waves and for shallow earthquakes. All these deviations are functions both of the latitude and of the epicentral distance. The difference in the distance to the Earth's centre at points with the same geocentric latitude on rays in the ellipsoidal and in the spherical model may reach several kilometres. Directly related to the deformation of the isovelocity surfaces, this difference is the only cause of significant perturbation in travel times. Other differences, such as that corresponding to the ray torsion, are of the first order in ellipticity, and may exceed 1 km. They induce only small differences in travel time, less than 0.01s. Thus, we show that the ellipticity correction obtained by Jeffreys (1935) and Bullen (1937) by a perturbational method can be recovered by a direct evaluation of the travel times in an ellipsoidal model of the Earth. Moreover, as stated by Dziewonski & Gilbert (1976), we verify the non-dependence of this correction on the choice of the velocity law.     

首都圈地区现今断层活动的数学模拟及地震危险性讨论

以首都圈地区现今活动断层上近２０年的位移测量资料为依据，用二维线弹性有限元对该区断层的活动特征进行了拟合，并结合有关资料讨论了该区的地震危险性。结果表明：１９７７－１９８６年间该区主压应力优势方位与华北较一致，约为Ｎ４５°－８０°Ｅ；１９８６－１９９０年其主压应力方向向北偏转，约为Ｎ５°－６０°Ｅ；１９９０年至现在其主压应力又向东偏转，角度大于第一阶段，约为Ｎ８０°－９５°Ｅ。该区的张家口－延庆一带近期有发生中强震的可能；丰镇－阳高－大同地区和凉城－古营盘地区应力较高，也应引起注意。     

Numerical modeling of rainstorm-induced shallow landslides in saturated and unsaturated soils

For the assessment of shallow landslides triggered by rainfall, the physically based model coupling the infinite slope stability analysis with the hydrological modeling in nearly saturated soil has commonly been used due to its simplicity. However, in that model the rainfall infiltration in unsaturated soil could not be reliably simulated because a linear diffusion-type Richards’ equation rather than the complete Richards’ equation was used. In addition, the effect of matric suction on the shear strength of soil was not actually considered. Therefore, except the shallow landslide in saturated soil due to groundwater table rise, the shallow landslide induced by the loss in unsaturated shear strength due to the dissipation of matric suction could not be reliably assessed. In this study, a physically based model capable of assessing shallow landslides in variably saturated soils is developed by adopting the complete Richards’ equation with the effect of slope angle in the rainfall infiltration modeling and using the extended Mohr–Coulomb failure criterion to describe the unsaturated shear strength in the soil failure modeling. The influence of rainfall intensity and duration on shallow landslide is investigated using the developed model. The result shows that the rainfall intensity and duration seem to have similar influence on shallow landslides respectively triggered by the increase of positive pore water pressure in saturated soil and induced by the dissipation of matric suction in unsaturated soil. The rainfall duration threshold decreases with the increase in rainfall intensity, but remains constant for large rainfall intensity.     

Seismic Hazard Analysis — Quo vadis?

The paper is dedicated to the review of methods of seismic hazard analysis currently in use, analyzing the strengths and weaknesses of different approaches. The review is performed from the perspective of a user of the results of seismic hazard analysis for different applications such as the design of critical and general (non-critical) civil infrastructures, technical and financial risk analysis. A set of criteria is developed for and applied to an objective assessment of the capabilities of different analysis methods. It is demonstrated that traditional probabilistic seismic hazard analysis (PSHA) methods have significant deficiencies, thus limiting their practical applications. These deficiencies have their roots in the use of inadequate probabilistic models and insufficient understanding of modern concepts of risk analysis, as have been revealed in some recent large scale studies. These deficiencies result in the lack of ability of a correct treatment of dependencies between physical parameters and finally, in an incorrect treatment of uncertainties. As a consequence, results of PSHA studies have been found to be unrealistic in comparison with empirical information from the real world. The attempt to compensate these problems by a systematic use of expert elicitation has, so far, not resulted in any improvement of the situation. It is also shown that scenario-earthquakes developed by disaggregation from the results of a traditional PSHA may not be conservative with respect to energy conservation and should not be used for the design of critical infrastructures without validation. Because the assessment of technical as well as of financial risks associated with potential damages of earthquakes need a risk analysis, current method is based on a probabilistic approach with its unsolved deficiencies.

Traditional deterministic or scenario-based seismic hazard analysis methods provide a reliable and in general robust design basis for applications such as the design of critical infrastructures, especially with systematic sensitivity analyses based on validated phenomenological models. Deterministic seismic hazard analysis incorporates uncertainties in the safety factors. These factors are derived from experience as well as from expert judgment. Deterministic methods associated with high safety factors may lead to too conservative results, especially if applied for generally short-lived civil structures. Scenarios used in deterministic seismic hazard analysis have a clear physical basis. They are related to seismic sources discovered by geological, geomorphologic, geodetic and seismological investigations or derived from historical references. Scenario-based methods can be expanded for risk analysis applications with an extended data analysis providing the frequency of seismic events. Such an extension provides a better informed risk model that is suitable for risk-informed decision making.     

三峡库区生基包滑坡监测及成果分析

生基包滑坡监测属于三峡库区奉节县三期地质灾害监测预警项目之一,该滑坡位于长江左岸,临近人口稠密的安坪乡集镇,地理位置重要。三峡水库175 m蓄水后,其变形破坏特征有何表现？对航道安全运营有无潜在威胁？是否会对滑坡体上的重要建筑及村民生产生活构成危害？针对这些问题,首先分析了滑坡的工程地质特征及主要的影响因素;其次,确立以4种监测手段为主、人工巡查为辅的监测方案;通过对大地变形GPS、深部位移、滑坡推力等几种监测方法的运用及对其成果进行分析研究,以实例说明其在滑坡监测中的应用;再次,结合宏观人工巡查进行对照分析,以充分说明大地变形、深部位移和滑坡推力监测在实际运用中的可行性;最后,根据监测结论提出对生基包滑坡防治的建议。     

滑坡体三维地质建模与可视化分析

针对滑坡地质研究的自身特点,提出了面向滑坡地质体三维建模的NURBS-BRep混合数据结构和地质结构单元实体构造技术.通过对滑坡基础地质数据的预处理、滑坡数字地形和滑动面的三维建模、地质结构面的系统构造、地质结构体的生成和显示,形成了一套完整的滑坡三维地质模型的构建方法.将该方法应用于某水库滑坡,建立了相应的三维滑坡地质模型,并基于此模型进行了三维剖切分析、数字钻孔、等值线生成、滑块自动剖分、滑坡失稳可视化动态模拟和滑坡体方量精确计算等一系列实用的可视化分析,为滑坡稳定性的准确计算和客观评价提供了有力的支持.     

采用钻孔倾斜仪监测滑坡深部水平位移方法和仪器的研究

介绍了目前滑坡深部水平位移钻孔倾斜仪人工监测方法存在的问题,钻孔深部位移自动化监测通信方法的选择。地质滑坡灾害监测中广泛使用钻孔倾斜仪,介绍了CZ-2型钻孔倾斜仪的构成及现场测试方法。     

基于贡献权重叠加法的滑坡风险区划

贡献权重叠加法是对滑坡本底因子和承灾体因子在滑坡发育中的贡献率进行统计后,通过贡献率均值化、规一化处理,利用权重转换模型计算出每一个因子内部的权重——自权重和因子相互之间的权重——互权重。将滑坡因子、承灾体因子和两者的自权重、互权重分别相乘叠加,得到滑坡危险度和易损度区划结果,再将2项结果相乘得出滑坡风险区划结果。     

地震折射波的虚拟射线理论及其应用

本文提出了一种地震折射液的虚拟射线理论,根据该理论,只要已知地表层的速度,即可直接由折射波信息提取地震参数,从而可实现折射界面的反演。 文中通过实际介质模型的计算机实验结果,验证了该理论的正确性。     

Deformation Patterns of an Accretionary Wedge in the Transition Zone from Subduction to Collision Offshore Southwestern Taiwan

Swath bathymetry data and seismic reflection profiles have been used to investigate details of the deformation pattern in the area offshore southwestern Taiwan where the Luzon subduction complex encroaches on the passive Chinese continental margin. Distinctive fold-and-thrust structures of the convergent zone and horst-and-graben structures of the passive margin are separated by a deformation front that extends NNW-ward from the eastern edge of the Manila Trench to the foot of the continental slope. This deformation front gradually turns into a NNE–SSW trending direction across the continental slope and the Kaoping Shelf, and connects to the frontal thrusts of the mountain belt on land Taiwan. However, the complex Penghu submarine canyon system blurs the exact location of the deformation front and nature of many morphotectonic features offshore SW Taiwan. We suggest that the deformation front offshore SW Taiwan does not appear as a simple structural line, but is characterized by a series of N–S trending folds and thrusts that terminate sequentially in an en-echelon pattern across the passive Chinese continental slope. A number of NE–SW trending lineaments cut across the fold-and-thrust structures of the frontal accretionary wedge and exhibit prominent dextral displacement indicative of the lateral expulsion of SW Taiwan. One of the prominent lineaments, named the Yung-An lineament, forms the southeastern boundary of the upper part of the Penghu submarine canyon, and has conspicuous influence over the drainage pattern of the canyon