Determination of the regional deformation rates of Shanghai and Kashima VLBI stations based on ITRF97

The vertical deformation rates (VDRs) and horizontal deformation rates (HDRs) of Shanghai VLBI station in China and Kashima and Kashima34 VLBI stations in Japan were re-analysed using the baseline length change rates from Shanghai to 13 global VLBI stations, and from Kashima to 27 stations and from Kashima34 to 12 stations, based on the NASA VLBI global solution glb1123 (Ma, 1999). The velocity vectors of the global VLBI stations were referred to the ITRF97 reference frame, and the Eulerian vectors of different models of plate motion were used for comparative solutions. The VDR of Shanghai station is estimated to be −1.91±0.56 mm/yr, and those of Kashima and Kashima34 stations, −3.72±0.74mm/yr and −8.81±0.84mm/yr, respectively. The difference between the last two was verified by further analysis. Similar estimates were also made for the Kokee, Kauai and MK_VLBA VLBI stations in mid-Pacific.     

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.     

Morphological characters of the two-ribbon flare of 1981 May 13 and their explanation

We give a summary of the morphology of the two-ribbon flare of 1981 May 13. One striking feature is that the Ha flare began at about 0338' UT and the double-ribbon structure was formed about 0346, $\text{before}$ the impulsive phase of the radio 3 cm burst at 04 11 UT. The 3 cm radio burst flux beginning at 03 33 UT showed only slow, stepwise increases lasting half an hour until the impulsive phase and this type of increase is usually regarded to be a typical thermal process. Each step in the radio flux corresponded to a variation in the Ha flare, showing that the radio and Hα emissions during this period came from the same thermal source. In this paper, we explain this behaviour in terms of Hyder's model: we think that the magnetic trough supporting the solar prominence rose for some reason, causing the prominence matter (the dark filaments) to fall along the magnetic lines and to hit the chromosphere and trigger off the flare. We give rough estimates of the energy density, the height of prominence and the infall matter at the different radio increments. We also give a qualitative explanation for the appearances of the single-peak structure in the radio burst at 0411 and the covering of the sunspot shortly after at 04 13 and propose several possible mechanisms.     

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年代中国石油生产前景作了初步展望。     

Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

Sedimentary cycles of trace elements in Salawusu River Valley since 150 ka BP

The paper makes some analyses on 11 trace elements in the Milanggouwan stratigraphical section in the Salawusu River valley, which is regarded as a prototype geology-palaeoclimate record since 150 ka BP. The results show that the content and variation of trace elements has experienced remarkably regular changes in the pace with coarse and fine sedimentary cycles of palaeo-aeolian sands to its overlying fluvio-lacustrine facies or/and palaeosols. The trace elements with chemical properties of relatively active (V, Sr, Cu, Ni, As) and relatively stable (P, Pb, Rb, Mn, Nb, Zr) are a manifestation of the corresponding 27 changeable cycles between peak and valley values, appearing a multi-fiuctuational process line of relative gathering and migration since then. The low numerical value distribution of these two types of trace elements in the aeolian sand facies represents erosion and accumulation under wind force during the cold-dry climate. Whereas their enrichments in both fluvio-lacustrine facies and palaeosols are related to the valley’s special low-lying physiognomic position between the Ordos Plateau and the Loess Plateau under the warm and humid climate conditions. The above relatively migrated and gathered change of the trace elements is the result of 27 climatic cycles of cold-dry and warm-humid, which is probably caused by repeated alternations of winter monsoon and summer monsoon in the Mu Us Sandy Land influenced by the climate vicissitudes in northern hemisphere during glacial and interglacial periods since 150 ka BP.     

辽西中生代火山岩岩浆成因——稀土元素配分及同位素定量模式

本文通过对辽西及邻区中生代盆地火山岩稀土元素配分特征的研究及模拟运算，探讨了其岩浆成因及过程。结果显示南北岩区岩浆成因不同，南岩区岩源为MORB地幔和长期富集型幔熔融混合而成，北岩区为MORB和近期亏损型地幔熔融混合而成。     

草簇状石墨烯片限域Pd纳米粒子及其催化应用(英文)

通过氨基离子液体改性石墨烯,并将其固载于堇青石表面,作为负载型Pd催化剂的载体.所制备的Pd催化剂经加氢老化后,表面石墨烯呈草簇状结构,将Pd纳米粒子限域于片层内,有效防止了Pd的流失和团聚.在重要的工业反应对羧基苯甲醛(4-CBA)加氢中,此结构催化剂与传统的钯碳催化剂相比,表现出很好的稳定性