An Early Cretaceous slab window beneath central Tibet,SW China: Evidence from OIB‐like alkaline gabbros in the Duolong area

The geodynamic mechanism of the late Early Cretaceous magmatic flare‐up in the collisional zone between the Lhasa and Qiangtang terranes in Tibet is controversial because of a scarcity of robust evidence. To address this problem, we report geochronological, geochemical and Hf isotopic data for the newly discovered Gufeng gabbros from the Duolong Cu–Au mineral district of the western Bangong–Nujiang Suture Zone (BNSZ). The gabbro samples, dated at 126.3 ± 1.8 Ma, show geochemical similarities to typical ocean island basalt (OIB) and have positive εHf(t) values of +3.3 to +6.9. The gabbros were generated by decompression melting of deep upwelling asthenosphere. This event is best explained by slab break‐off and the resultant development of a slab window beneath central Tibet.     

A multiscale work-analysis approach for geotechnical structures

This paper presents a second-order work analysis in application to geotechnical problems by using a novel effective multiscale approach. To abandon complicated equations involved in conventional phenomenological models, this multiscale approach employs a micromechanically-based formulation, in which only four parameters are involved. The multiscale approach makes it possible a coupling of the finite element method (FEM) and the micromechanically-based model. The FEM is used to solve the boundary value problem (BVP) while the micromechanically-based model is utilized at the Gauss point of the FEM. Then, the multiscale approach is used to simulate a three-dimensional triaxial test and a plain-strain footing. On the basis of the simulations, material instabilities are analyzed at both mesoscale and global scale. The second-order work criterion is then used to analyze the numerical results. It opens a road to interpret and understand the micromechanisms hiding behind the occurrence of failure in geotechnical issues.     

Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

四川石棉独立碲矿床蚀变岩的地质地球化学研究

四川石棉大水沟独立碲矿床堪称世界首例，其矿体两侧的围岩蚀变较窄，且蚀变带与矿脉及未蚀变岩呈清楚而截然的接触关系。主要围岩蚀变为白云石化、白（绢）云母化、云英岩化及电气石化等。这些蚀变作用分别发生于１７７．７～１６５．１Ｍａ和９１．７１～８０．１９Ｍａ。蚀变作用过程中ＳｉＯ２大量析出，Ａｌ２Ｏ３、ＴｉＯ２、ＣａＯ、ＭｇＯ、Ｋ２Ｏ、Ｈ２Ｏ、ＣＯ２、Ｆｅ２Ｏ３、ＭｎＯ则有不同程度地带入。与此同时，ＲＥＥ在此过程中包集体迁出，Ｔｅ、Ｂｉ、ＡＳ等元素的含量则决定于蚀变作用及其强度，而与原岩无关。     

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.     

Luminescence dating of sand–loess sequences and response of Mu Us and Otindag sand fields (north China) to climatic changes

The sand–loess transition zone in north China is sensitive to climate change, and is an ideal place to investigate past environmental changes. However, past climate change at millennial–centennial timescales in this region has not been well reconstructed because of limited numerical dating. Alternations of sandy loam soils with aeolian sand layers in the Mu Us and Otindag sand fields, which lie along the sand–loess transition zone, indicate multiple intervals of dune activity and stability. This change is probably a response to variations of the East Asian monsoon climate during the late Quaternary. The single aliquot regeneration (SAR) optically stimulated luminescence (OSL) dating protocol, which has been successfully applied to aeolian deposits worldwide, is applied to these two sand fields in this study. The OSL ages provide reliable constraints for reconstruction of past climate changes at suborbital timescale. Sections in both sand fields contain aeolian sand beds recording millennial‐scale episodes of dry climate and widespread dune activation, including episodes at about the same time as Heinrich Event 5 and the Younger Dryas in the North Atlantic region. These results demonstrate the potential of aeolian sediments in semi‐arid north China to record millennial‐scale climatic events, and also suggest that dry–wet climate variation at the desert margin in China may be linked to climatic change elsewhere in the Northern Hemisphere, through atmospheric circulation. This article was published online on 27 November 2008. An error was subsequently identified. This notice is included in the online and print versions to indicate that both have been corrected (16 December 2008). Copyright © 2008 John Wiley & Sons, Ltd.     

A THREE-DIMENSIONAL ELASTIC NESTED-GRID MESO-(β-γ)SCALE ATMOSPHERIC MODEL——PARTI:MODEL DESCRIPTION

A three-dimensional elastic nonhydrostatic mesoscale(β-γ)model with nested-grid is presented.It uses a set of fullequations in terrain-following coordinates as its basic dynamic frame,which is solved with a time-splitting algorithmfor acoustic and gravity waves.The model physical parameterization includes a K-theory subgrid eddy mixing for cloudand free atmosphere,a bulk planetary boundary layer parameterization,and three types of sofisticated cloudmicrophysics schemes with double-parameters for hail-bearing clouds,warm clouds and snowing clouds respectively.The model is designed to be used flexibly for simulations of a variety of meso-and small-scale atmospheric processes,and can be improved as a regional and local operational NWP system in future.     

中国石油生产的展顾

本文分析了中国石油生产近20年变化的原因,对一些学者认为80年代中国会成为超级产油大国,后又悲观断定80年代中国可能变成石油净进口国作了评述。同时对90年代中国石油生产前景作了初步展望。     

The evolution of quasars selected by slitless technique

In former works (Zhouet al., 1983, 1985), a quantitative method have been developed to take the selection effects in the identification of emission lines of quasars into account. It was proved that these selection effects may be the cause of the unevenness in the redshift distribution of quasars. The present work is a continuation and development of former works. We use results given by the surveys with same limit-apparent magnitude and choose the quasars whose absolute magnitudes are within a specific range as the analysing samples. Using the method given in the former papers we may find out the evolutionary parameter in an evolutionary law with form of (1+z) ^y from the best fitting between the calculative and observational redshift distribution. The result of analysis shows that the evolutionary law of quasars selected by slitless technique isp =p ₀(1 + z)^{6.5 ± 1} up toz=2.8. This result coincides with and generalizes the earlier result given by other authors.