利用地震矩张量反演鲜水河断裂带现今运动学特征

本文探讨了利用地震矩反演断裂形变带运动学参数的基本理论和方法，将其初步应用于鲜水河断裂形变带变形分析和运动机制的研究。结果表明，鲜水河断裂带呈现出走向拉伸、倾向压缩的形变格局，由地震矩反演的断裂带剪切形变速率（１０．９ｍｍ／ａ）与用地质学估算方法（１７ｍｍ／ａ）和现今地壳形变测量（８ｍｍ／ａ）的结果相当。同时，反演出的应变主方向能解释鲜水河断裂现今活动分段性特征以及多种滑动方式共存的现状，从而证明该方法是目前研究区域运动学问题切实有效的手段之一。     

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.     

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

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

RETRACTION: Numerical simulation of a wind-sand flow under different conditions of atmospheric stability

A theoretical model for wind‐sand flow is developed by considering the coupling between wind flow and sand particle motion, the latter subject to the Magnus effect, under different atmospheric stability conditions. Using this model, the characteristics of the wind‐sand flow are discussed in detail. The results show that the atmospheric stability and the Magnus effect both have a strong influence on wind profiles and on the trajectories of sand particles. This approach produces results with characteristics that differ from those previously reported; the latter only applying to atmospheric conditions of neutral stability. The saltating sand reaches a greater height under non‐neutral stability than under neutral stability, while the maximum horizontal distance is greater under unstable conditions and is smaller under stable conditions than under conditions of neutral stability.     

Zircon U–Pb ages, Hf and O isotopes constrain the crustal architecture of the ultrahigh-pressure Dabie orogen in China

The crustal structure of the Dabie orogen was reconstructed by a combined study of U–Pb ages, Hf and O isotope compositions of zircons from granitic gneiss from North Dabie, the largest lithotectonic unit in the orogen. The results were deciphered from metamorphic history to protolith origin with respect to continental subduction and exhumation. Zircon U–Pb dating provides consistent ages of 751 ± 7 Ma for protolith crystallization, and two group ages of 213 ± 4 to 245 ± 17 Ma and 126 ± 4 to 131 ± 36 Ma for regional metamorphism. Majority of zircon Hf isotope analyses displays negative ε_Hf(t) values of − 5.1 to − 2.9 with crust Hf model ages of 1.84 to 1.99 Ga, indicating protolith origin from reworking of middle Paleoproterozoic crust. The remaining analyses exhibit positive ε_Hf(t) values of 5.3 to 14.5 with mantle Hf model ages of 0.74 to 1.11 Ga, suggesting prompt reworking of Late Mesoproterozoic to Early Neoproterozoic juvenile crust. Zircon O isotope analyses yield δ¹⁸O values of − 3.26 to 2.79‰, indicating differential involvement of meteoric water in protolith magma by remelting of hydrothermally altered low δ¹⁸O rocks. North Dabie shares the same age of Neoproterozoic low δ¹⁸O protolith with Central Dabie experiencing the Triassic UHP metamorphism, but it was significantly reworked at Early Cretaceous in association with contemporaneous magma emplacement. The Rodinia breakup at about 750 Ma would lead to not only the reworking of juvenile crust in an active rift zone for bimodal protolith of Central Dabie, but also reworking of ancient crust in an arc-continent collision zone for the North Dabie protolith. The spatial difference in the metamorphic age (Triassic vs. Cretaceous) between the northern and southern parts of North Dabie suggests intra-crustal detachment during the continental subduction. Furthermore, the Dabie orogen would have a three-layer structure prior to the Early Cretaceous magmatism: Central Dabie in the upper, North Dabie in the middle, and the source region of Cretaceous magmas in the lower.     

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.     

Influence of Structural Non-Stationarity of Surface Roughness on Morphological Characterization and Mechanical Deformation of Rock Joints

Summary Structural non-stationarity of surface roughness affects accurate morphological characterization as well as mechanical behaviour of rock joints at the laboratory scale using samples with a size below the stationarity threshold. In this paper, the effect of structural non-stationarity of surface roughness is investigated by studying the scale dependence of surface roughness and mechanical behaviour of rock joints. The results show that the structural non-stationarity mainly affects the accurate characterization of the surface roughness of the fracture samples. It also controls the amount and location of the contact areas during shear tests, which in turn affects the mechanical properties and asperity degradation of the samples. It is concluded that for accurate determination of the morphological and mechanical properties of rock joints at laboratory and field scales, samples with size equal to or larger than the stationarity threshold are required. Author’s address: Nader Fardin, Rock Mechanics Group, Department of Mining Engineering, Faculty of Engineering, University of Tehran, P.O. Box: 11365/4563, Tehran, Iran     

Nonlinear Analysis of Local Site Effects on Seismic Ground Response in the Bam Earthquake

The influence of local geologic and soil conditions on the intensity of ground shaking is addressed in this study. The amplification of the ground motion due to local site effects resulted in severe damage to dwellings in the Bam area during the 2003 Bam Earthquake. A unique set of strong motion acceleration recordings was obtained at the Bam accelerograph station. Although the highest peak ground acceleration recorded was the vertical component (nearly 1 g), the longitudinal component (fault-parallel motion) clearly had the largest maximum velocity as well as maximum ground displacement. Subsurface geotechnical and geophysical (down-hole) data in two different sites have been obtained and used to estimate the local site condition on earthquake ground motion in the area. The ground response analyses have been conducted considering the nonlinear behavior of the soil deposits using both equivalent linear and nonlinear approaches. The fully nonlinear method embodied in FLAC was used to evaluate the nonlinear soil properties on earthquake wave propagation through the soil layer, and compare with the response from the equivalent linear approach. It is shown that thick alluvium deposits amplified the ground motion and resulted in significant damage in residential buildings in the earthquake stricken region. The comparison of results indicated similar response spectra of the motions for both equivalent and nonlinear analyses, showing peaks in the period range of 0.3–1.5 s. However, the amplification levels of nonlinear analysis were less than the equivalent linear method especially in long periods. The observed response spectra are shown to be above the NEHRP building code design requirements, especially at high frequencies.     

Interference of Multiple Strip Footings on Sand Using Small Scale Model Tests

By using small scale model tests, the interference effect on the vertical load-deformation behavior of a number of equally spaced strip footings, placed on the surface of dry sand, was investigated. At any stage, all the footings were assumed to (i) carry exactly equal magnitude of load, and (ii) settle to the same extent. No tilt of the footing was permitted. The effect of clear spacing (s) among footings on the results was explored. A new experimental setup was proposed in which only one footing needs to be employed rather than a number of footings. The bearing capacity increases continuously with decrease in spacing among the footings. The interference effect becomes further prominent with increase in soil friction angle. In contrast to an increase in the bearing capacity, with decrease in spacing of footings, an increase in the footing settlement associated with the ultimate state of shear failure was observed. The present experimental observations were similar to those predicted by the available theory, based on the method of characteristics. As compared to the theory, the present experimental data, however, indicates much greater effect of interference especially for larger spacing among footings.