Microzonation of the city of Basel

During the past centuries, the city of Basel has suffered damage caused by earthquakes. One extraordinary event described in historical documents is the strong earthquake which occurred in 1356. The 1356 event, one of the strongest earthquakes in northwest-Europe, was obviously much stronger than the low-magnitude earthquakes observed in the area during this century. Even though the present seismicity in the Basel area is low, strong earthquakes have to be expected due to the city's geographical location close to the northern boundary of the African-European convergence zone, at the southern end of the Rhinegraben. A crucial step towards preparedness for future events and mitigation of earthquake risk involves a microzonation study of the city. The study is carried out in three steps: (1) a detailed mapping of the geology and geotechnical properties of the area, (2) measurement, interpretation and modelling of ambient noise data, and (3) numerical modelling of expected ground motions during earthquakes. A qualitative microzonation of the centre of Basel is presented, and it is discussed by comparing it to the historically reported damage of the 1356 earthquake.     

Modelling of ground acceleration in the near source range: the case of 1976, Friuli earthquake (M = 6.5), northern Italy

A mixed statistical-deterministic model of earthquake rupture is developed for evaluating the strong ground motion in the near source range (receiver distance comparable to the fault length). The source parametrization is based on the k-square model and the propagation is computed by asymptotic Green's functions. The method is applied to the case of 1976, Friuli earthquake (M = 6.5) in northern Italy which occurred on a low-dip thrusting fault. Acceleration records at 29 stations are computed for 100 simulations of rupture histories. The mean value map of peak ground accelerations shows clearly a maximum to the south due to the inner geometry and directivity of the source. The variation of the estimated PGA versus the epicentral distance is strongly dependent on azimuth and is not decreasing monotonically. The comparison of these curves with those predicted by empirical acceleration–distance relationships shows discrepancies in the near source distance range. This study shows the importance of considering the complexity of the source rupture process for strong motion estimate in the near source range.     

中国天山花岗岩类的时空分布及构造意义

天山花岗岩类具明显的时空分带特征,以古生代最为明显.它们可分为:1.中天山南缘中加里东期花岗岩带:2.中天山北缘晚加里东期花岗岩带;3.中天山南缘早华力西期花岗岩带;4.中天山北缘中华力西期花岗岩带;5.马宗山南缘晚华力西期花岗岩带;6.南天山南坡晚华力西期花岗岩带.在时空上,花岗岩类在中天山陆壳南、北缘呈交替分布的格局,显示出与天山陆壳南、北张合呼应、交替演化的和谐一致关系.中天山花岗岩类是几个地槽主旋回的产物.天山花岗岩带与蛇绿岩带成“对”分布,沿“中天山南缘深断裂”,及“中国天山主干断裂”两侧延伸,二者互相平行.天山花岗岩类是多成因的.分为洋壳重熔型、陆壳重熔型及洋壳陆壳混合重熔型花岗岩.     

有限元素法全倾角波动方程偏移

本文试图从两个方面提高波动方程偏移的效果:1.使用全方向波动方程进行偏移,使之适用于任意倾角界面的反射;2.使用有限元素法提高微分方程数值解法的精度。     

Estimation of strong ground motions from hypothetical earthquakes on the Cascadia subduction zone,Pacific Northwest

Strong ground motions are estimated for the Pacific Northwest assuming that large shallow earthquakes, similar to those experienced in southern Chile, southwestern Japan, and Colombia, may also occur on the Cascadia subduction zone. Fifty-six strong motion recordings for twenty-five subduction earthquakes ofM _s7.0 are used to estimate the response spectra that may result from earthquakesM _w<81/4. Large variations in observed ground motion levels are noted for a given site distance and earthquake magnitude. When compared with motions that have been observed in the western United States, large subduction zone earthquakes produce relatively large ground motions at surprisingly large distances. An earthquake similar to the 22 May 1960 Chilean earthquake (M _w 9.5) is the largest event that is considered to be plausible for the Cascadia subduction zone. This event has a moment which is two orders of magnitude larger than the largest earthquake for which we have strong motion records. The empirical Green's function technique is used to synthesize strong ground motions for such giant earthquakes. Observed teleseismicP-waveforms from giant earthquakes are also modeled using the empirical Green's function technique in order to constrain model parameters. The teleseismic modeling in the period range of 1.0 to 50 sec strongly suggests that fewer Green's functions should be randomly summed than is required to match the long-period moments of giant earthquakes. It appears that a large portion of the moment associated with giant earthquakes occurs at very long periods that are outside the frequency band of interest for strong ground motions. Nevertheless, the occurrence of a giant earthquake in the Pacific Northwest may produce quite strong shaking over a very large region.     

Predicted near-field ground motion for dynamic stress-drop models

We propose a finite difference method, using a hexagonal grid, to compute displacements (stresses, velocities, accelerations) in the near-field of a 2-D in-plane stress-drop crack, in both whole space (constant stress-drop) and half-space (depth-dependent stress-drop). To exercise the method, the stress field distribution is evaluated for both fundamental 2-D shear cracks, anti-plane. In order to test the method's reliability, the results are compared with some analytical and numerical solutions available in the literature (Kostrov, 1964;Virieux andMadariaga, 1982). For the in-plane source, the results emphasize that the method can resolve the stress concentration due to the rupture front from the stress peak associated with the shear wave propagating in front of the crack. Synthetic motions are computed on the fault, but also in an infinite medium and at the free surface. The rather complex waveforms generated in the near-field, even by simple sources, emphasize the contribution of all wave terms (near, intermediate and far-field) to the motion. The presence of near-field and the numerical procedure explain the significant low frequency content of the computed seismograms. The set of treated problems proves the method is stable and accurate.     

Numerical experiments in the 3/1 andv 6 overlapping resonance region

Numerical experiments of fictitious small bodies with initial eccentricities e=0.1 have been performed in the overlapping region of the 3/1 mean motion resonance and of thev ₆ secular resonance 2.48a2.52AU for different values of the initial inclination 16°i20°. An analysis for thev ₆ secular resonance shows that the topology is different from the one found outside the overlapping region: the critical argument for thev ₆ resonance in the overlapping region rotates in opposite direction as compared to the purev ₆ region. In the 3/1 resonance region the secular resonancev ₅ is dominant, and some secondary secular resonances asv ₆ –v ₁₆ andv ₅ +v ₆ are present.     

Numerical results to the Sitnikov-problem

We present numerical results of the so-called Sitnikov-problem, a special case of the three-dimensional elliptic restricted three-body problem. Here the two primaries have equal masses and the third body moves perpendicular to the plane of the primaries' orbit through their barycenter. The circular problem is integrable through elliptic integrals; the elliptic case offers a surprisingly great variety of motions which are until now not very well known. Very interesting work was done by J. Moser in connection with the original Sitnikov-paper itself, but the results are only valid for special types of orbits. As the perturbation approach needs to have small parameters in the system we took in our experiments as initial conditions for the work moderate eccentricities for the primaries' orbit (0.33e _primaries 0.66) and also a range of initial conditions for the distance of the 3 ^rd body (= the planet) from very close to the primaries orbital plane of motion up to distance 2 times the semi-major axes of their orbit. To visualize the complexity of motions we present some special orbits and show also the development of Poincaré surfaces of section with the eccentricity as a parameter. Finally a table shows the structure of phase space for these moderately chosen eccentricities.     

Near-fault strong motion complexity of the 2000 Tottori earthquake (Japan) from a broadband source asperity model

The October 6/2000 Tottori earthquake that occurred in central Japan was an intermediate size strike-slip event that produced a very large number of near field strong motion recordings. The large amount of recorded data provides a unique opportunity for investigating a source asperity model of the Tottori earthquake that, combined with a hybrid strong motion simulation technique, is able to reproduce the observed broadband frequency near-fault ground motion.

We investigated the optimum source asperity parameters of the Tottori earthquake, by applying a Genetic Algorithm (GA) inversion scheme to optimise the fitting between simulated and observed response spectra and Peak Ground Acceleration (PGA) values. We constrained the initial model of our inversion by using the heterogeneous slip distribution obtained from a kinematic inversion of the source of previous studies. We used all the observed near-fault ground motions (−100 m) from the borehole strong motion network of Japan (KiK-Net), which are little affected by surficial geology (site effects).

The calculation of broadband frequency strong ground motion (0.1–10 Hz) is achieved by applying a hybrid technique that combines a deterministic simulation of the wave propagation for the low frequencies and a semi-stochastic modelling approach for the high frequencies. For the simulation of the high frequencies, we introduce a frequency-dependent radiation pattern model that efficiently removes the dependence of the pattern coefficient on the azimuth and take-off angle as the frequency increases. The good agreement between the observed and simulated broadband ground motions shows that our inversion procedure is successful in estimating the optimum asperity parameters of the Tottori earthquake and provides a good test for the strong ground motion simulation technique.

The ratio of background stress drop to average asperity stress drop from our inversion is nearly 50%, in agreement with the theoretical asperity model of Das and Kostrov [Das, S., Kostrov, B.V., 1986. Fracture of a single asperity on a finite fault: a model for weak earthquakes? Earthquake Source Mechanics, AGU, pp. 91–96.], and an empirical ratio of asperities to rupture area [Seismol. Res. Lett. 70 (1999) 59–80.].

The simulated radiation pattern is very complex for epicentral distances within half the fault length, but it approaches the radiation of a double-couple point source for larger distances.

The rupture velocity and rise time have a significant influence on the Peak Ground Velocity (PGV) distribution around the fault. An increase in rupture velocity produces a similar effect on the ground motion as a reduction in rise time.     

Seismic hazard maps of Mexico, the Caribbean, and Central and South America

The growth of megacities in seismically active regions around the world often includes the construction of seismically unsafe buildings and infrastructures due to an insufficient knowledge of existing seismic hazard and/or economic constraints. Minimization of the loss of life, property damage, and social and economic disruption due to earthquakes depends on reliable estimates of seismic hazard. We have produced a suite of seismic hazard estimates for Mexico, the Caribbean, and Central and South America. One of the preliminary maps in this suite served as the basis for the Caribbean and Central and South America portion of the Global Seismic Hazard Map (GSHM) published in 1999, which depicted peak ground acceleration (pga) with a 10% chance of exceedance in 50 years for rock sites. Herein we present maps depicting pga and 0.2 and 1.0 s spectral accelerations (SA) with 50%, 10%, and 2% chances of exceedance in 50 years for rock sites. The seismicity catalog used in the generation of these maps adds 3 more years of data to those used to calculate the GSH Map. Different attenuation functions (consistent with those used to calculate the U.S. and Canadian maps) were used as well. These nine maps are designed to assist in global risk mitigation by providing a general seismic hazard framework and serving as a resource for any national or regional agency to help focus further detailed studies required for regional/local needs. The largest seismic hazard values in Mexico, the Caribbean, and Central and South America generally occur in areas that have been, or are likely to be, the sites of the largest plate boundary earthquakes. High hazard values occur in areas where shallow-to-intermediate seismicity occurs frequently.