Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

应用GPS观测青藏高原东北缘应力场变化

采用各向同性弹性地球模型推导了地面位移场速率与地壳内任意点应力场变化的边界积分关系,同时利用青藏高原东北缘1999～2001年观测的GPS资料对观测区地壳深度为5 km和25 km的主应力和最大剪应力进行了计算分析.结果表明,青藏高原东北缘的主应力变化主要集中在祁连山断裂、海原断裂等,在1920～1954年间历史上发生过多次震级为7.0～8.5级强震的断裂附近,并具有主应力变化沿断层走向分布、最大剪应力沿断层走向交替变化等特征.     

Progress in Earthquake Science and Technology in China: Review and Prospects (Ⅰ)

In the article the author looks back the hard development course and great progress in earth quake science and technology in China during the last near a half of century and expounds the following 3 aspects: (1) The strong desire of the whole society to mitigate seismic disasters and reduce the effect of earthquakes on social-economic live is a great driving force to push forward the development of earthquake science and technology in China; (2) To better ensure people‘ s life and property, sustainable economic development, and social stability is an essential purpose to drive the development of earthquake science and technology in China; and (3) To insist on the dialectical connection of setup of technical system for seismic monitoring with the scientific research of earthquakes and to better handle the relation between crucial task, current scientif ic level, and the feasibility are the important principles to advance the earthquake science and technology in China. Some success and many setbacks in earthquake disaster mitigation consistently enrich our knowledge regarding the complexity of the conditions for earthquake occurrence and the process of earthquake preparation, promote the reconstruction and modernization of technical system for earthquake monitoring, and deepen the scientific research of earthquakes. During the last 5 years, the improvement and modernization of technical system for earthquake monitoring have clearly provided the technical support to study and practice of earthquake prediction and pre caution, give prominence to key problems and broaden the field of scientific research of earth quakes. These have enabled us to get some new recognition of the conditions for earthquake oc currence and process of earthquake preparation, characteristics of seismic disaster, and mecha nism for earthquake generation in China‘s continent. The progress we have made not only en courages us to enhance the effectiveness of earthquake disaster mitigation, but also provides a basis for accelerating further development of earthquake science and technology in China in the new century, especially in the 10th five-year plan. Based on the history reviewed, the author sets forth a general requirement for develop ment of earthquake science and technology in China and brings out 10 aspects to be stressed and strengthened at present and in the future. These are: upgrade and setup of the network of digitized seismic observation; upgrade and setup of the network for observation of seismic pre cursors; setup of the network for observation of strong motion; setup of the laboratories for ex periment on seismic regime; establishment of technical system for seismic information, emer gency command and urgent rescue; research on short-term and imminent earthquake predic tion; research on intermediate- and long-term earthquake prediction; research on attenuation of seismic ground motion, mechanism for seismic disaster, and control on seismic disaster; ba sic research fields related to seismology and geoscience. We expect that these efforts will signifi cantly elevate the level of earthquake science and technology in China to the advanced interna tional level, improve theories, techniques, and methods for earthquake precaution and predic tion, and enhance the effectiveness of earthquake disaster mitigation.     

Microboudin structures of piemontite along the Sambagawa metamorphic belt, Japan: implications for lateral variation of differential stress during regional metamorphism

Evidence is presented of a lateral variation in differential stress during metamorphism along a regional metamorphic belt on the basis of the proportion of microboudinaged piemontite grains (p) in a quartz matrix in metacherts. It is proposed that p is a practical indicator of relative differential stress. Analysis of 123 metacherts from the 800 km long Sambagawa metamorphic belt, Japan, reveals that p‐values range from < 0.01 to 0.7 in this region. Most samples from Wakayama in the mid‐belt area have p‐values of 0.4–0.6, whereas those from western Shikoku have p‐values of < 0.1. This difference cannot be explained by variations in metamorphic temperature, and is instead attributed to a regional, lateral variation in differential stress during metamorphism.     

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

Interpolating paleovegetation data with an artificial neural network approach

To drive an atmospheric general circulation model (AGCM), land surface boundary conditions like albedo and morphological roughness, which depend on the vegetation type present, have to be prescribed. For the late Quaternary there are some data available, but they are still sparse. Here an artificial neural network approach to assimilate these paleovegetation data is investigated. In contrast to a biome model the relation between climatological parameters and vegetation type is not based on biological knowledge but estimated from the available vegetation data and the AGCM climatology at the corresponding locations. For a test application, a data set for the modern vegetation reduced to the amount of data available for the Holocene climate optimum (about 6000 years B.P.) is used. From this, the neural network is able to reconstruct the complete global vegetation with a kappa value of 0.56. The most pronounced errors occur in Australia and South America in areas corresponding to large data gaps.     

Application of stress release models to historical earthquakes from North China

Univariate and multivariate stress release models are fitted to historical earthquake data from North China. It is shown that a better fit is obtained by treating separately the Eastern part of the region, including the North China Plain and Bohai Sea, and the Western part of the region, including the Ordos Plateau and its Eastern boundary. Further improvement is obtained by fitting the large events (M7.6) and smaller events in the Western region by different stress release models. The comparisons are made by computing the likelihoods of the fitted models and discounting the number of parameters used by Akaike's AIC criterion. The models are used to develop long-term risk scenarios for the East and West regions.     

Age estimates of coastal terraces in the Andaman and Nicobar Islands and their tectonic implications

The great Indian Ocean earthquake of December 26, 2004 caused significant vertical changes in its rupture zone. About 800 km of the rupture is along the Andaman and Nicobar Islands, which forms the outer arc ridge of the subduction zone. Coseismic deformation along the exposed land could be observed as uplift/subsidence. Here we analyze the morphological features along the coast of the Andaman and Nicobar Islands, in an effort to reconstruct the past tectonics, taking cues from the coseismic effects. We obtained radiocarbon dates from coastal terraces of the island belt and used them to compute uplift rates, which vary from 1.33 mm yr^− 1 in the Little Andaman to 2.80 mm yr^− 1 in South Andaman and 2.45 mm yr^− 1 in the North Andaman. Our radiocarbon dates converge on  600 yr and  1000 yr old coastal uplifts, which we attribute to the level changes due to two major previous subduction earthquakes in the region.