A note on the algorithm of symplectic integrators

For a Hamiltonian that can be separated into N+1(N\geq 2) integrable parts, four algorithms can be built for a symplectic integrator. This research compares these algorithms for the first and second order integrators. We found that they have similar local truncation errors represented by error Hamiltonian but rather different numerical stability. When the computation of the main part of the Hamiltonian, H ₀, is not expensive, we recommend to use S ^* type algorithm, which cuts the calculation of the H ₀ system into several small time steps as Malhotra(1991) did. As to the order of the N+1 parts in one step calculation, we found that from the large to small would get a slower error accumulation. This revised version was published online in July 2006 with corrections to the Cover Date.     

广西中三叠统浊流流向及坡向

通过对广西中三叠统浊积岩系中溢堤流沉积的识别,结合同一鲍马层序Tc段斜层理倾向与其底部槽模走向不一致的分析,以及横向冲刷模、U形障碍冲刷痕等冲蚀构造、大型滑动模等压刻构造的认识,筛选出宽阔水道充填沉积、扇状体及分支水道沉积中底模构造的指向数据,来指示浊流主体流向,并区分出溢堤流流向。从而确定了本区中三叠世时的浊流主体流向为NW向,平均为300°。根据滑塌沉积中滑塌包卷层轴面、滑塌面的测定,确定凤山县金牙处的古坡向平均为322°,与浊流主体流向基本一致。还认为中三叠世时右江盆地沉积物源来自桂东和桂东南蚀源区,海水由其南的越南境内进入。     

The mineral chemistry of pyroxenite xenoliths in the volcanic rocks of Hoh Xil and their significance

The pyroxenite xenoliths in the volcanic rocks of Hoh Xil consist of clinopyroxenes and orthopyroxenes. The mineral composition of these pyroxenes is similar to that of mantle xenoliths including peridotite and pyroxenite from China and abroad, and different from that of granulites. The pyroxenes formed at 1101–1400°C (averaging 1250°C) and under 30–60 kb (averaging 46 kb). We deduced that the magma was derived from the mantle at a depth of more than 150 km, which fits in with the geophysical conclusion that the low-velocity layer existed in the mantle under 150 km.

     

The role of magnetohydrodynamics in heliospheric space plasma physics research

Magnetohydrodynamics (MHD) is a fairly recent extension of the field of fluid mechanics. While much remains to be done, it has successfully been applied to the contemporary field of heliospheric space plasma research to evaluate the macroscopic picture of some vital topics via the use of conducting fluid equations and numerical modeling and simulations. Some representative examples from solar and interplanetary physics are described to demonstrate that the continuum approach to global problems (while keeping in mind the assumptions and limitations therein) can be very successful in providing insight and large scale interpretations of otherwise intractable problems in space physics.Paper dedicated to Professor Hannes Alfvén on the occasion of his 80th birthday, 30 May 1988.     

Age-metallicity relation for F-stars

We derived the age-metallicity relation for the solar neighbourhood from theuvby, H photometry of F-stars using the evolutionary models of Hejlesen (1980). It is shown that disk formation is divided into two phases according to the fast and slow collapse which accompany the initial rapid enrichment and slow increase of the metallicity of the galactic disk, respectively.Paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 September–6 October, 1984.     

Geochemistry of re and trace elements in ferromanganese nodules and their genesis

Trace and RE element geochemistry and genesis have been studied with respect to ferromanganese nodules from the sediments of the Pacific, Atlantic and Indian oceans.     

Evolution of electron and proton temperatures in a flaring loop

We have investigated numerically how a temperature difference between electrons and protons is produced in a flaring loop by adopting a one-fluid, two-temperature model instead of a single-temperature model. We have treated a case in which flare energy is released in the form of heating of electrons located in the top part of the loop.In this case, a large temperature difference (T _e/T _p 10) appears in the corona in the energy-input phase of the flare. When the material evaporated from the chromosphere fills the corona, the temperature difference in the loop begins to shrink rapidly from below. Eventually, in the loop apex, the proton temperature exceeds the electron temperature mainly due to cooling of the electrons by conduction down the loop and heating of the protons by compression of the ascending material. In the late phase of the flare (t 15 min from the flare onset), the temperature difference becomes less than 2% of the mean temperature of electrons and protons at every point in the loop.     

A preliminary stable isotope study on Mogok Ruby, Myanmar

The primary occurrence of ruby in the Mogok area, northern Myanmar is exclusively found in marble along with spinel–forsterite-bearing marble and phlogopite–graphite marble. These marble units are enclosed within banded biotite–garnet–sillimanite–oligoclase gneisses. Samples of these marbles collected for C–O stable isotope analysis show two trends of δ¹³C–δ¹⁸O variation resulting most likely from fluid–rock interactions. Ruby-bearing marble and phlogopite–graphite marble follow a trend with coupled C–O depletion, whereas spinel–forsterite-bearing marble follows a δ¹⁸O depletion trend with relatively constant δ¹³C values. Ruby formation might have resulted from CO₂-rich fluid–rock interaction, while spinel–forsterite-bearing marble was genetically related to CO₂-poor fluid–rock interaction. Both fluids may have arisen from external sources. Based on graphite Raman spectral thermometry, the estimated temperature for phlogopite–graphite marble, and probably ruby-bearing marble, was lower than 607 °C, and for spinel–forsterite-bearing marble, lower than 710 °C. Contrasting C/O diffusion between graphite/ruby/spinel/forsterite and calcite, local variations of isotopic compositions of newly formed minerals as a result of non-pervasive fluid infiltration, and open-system isotopic disturbance during cooling may have affected C-/O-isotopic fractionations between minerals. The estimated high formation temperatures for ruby and spinel/forsterite imply that the parental fluids may have been related to nearby igneous intrusions and/or metamorphic processes. Whether these two types of fluid were genetically related is unclear based on the present data.     

Preliminary results for a semi-automated quantification of site effects using geomorphometry and ASTER satellite data for Mozambique,Pakistan and Turkey

Estimation of the degree of local seismic wave amplification (site effects) requires precise information about the local site conditions. In many regions of the world, local geologic information is either sparse or is not readily available. Because of this, seismic hazard maps for countries such as Mozambique, Pakistan and Turkey are developed without consideration of site factors and, therefore, do not provide a complete assessment of future hazards. Where local geologic information is available, details on the traditional maps often lack the precision (better than 1:10,000 scale) or the level of information required for modern seismic microzonation requirements. We use high-resolution (1:50,000) satellite imagery and newly developed image analysis methods to begin addressing this problem. Our imagery, consisting of optical data and digital elevation models (DEMs), is recorded from the ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) sensor system. We apply a semi-automated, object-oriented, multi-resolution feature segmentation method to identify and extract local terrain features. Then we classify the terrain types into mountain, piedmont and basin units using geomorphometry (topographic slope) as our parameter. Next, on the basis of the site classification schemes from the Wills and Silva (1998) study and the Wills et al (2000) and Wills and Clahan (2006) maps of California, we assign the local terrain units with V _s 30 (the average seismic shear-wave velocity through the upper 30m of the subsurface) ranges for selected regions in Mozambique, Pakistan and Turkey. We find that the applicability of our site class assignments in each region is a good first-approximation for quantifying local site conditions and that additional work, such as the verification of the terrain’s compositional rigidity, is needed.     

流体包裹体成分测定的低温相变和显微拉曼光谱分析技术研究进展

论述了包裹体低温分析技术的原理以及进展。传统的流体包裹体低温分析技术是以显微冷冻测温测定无机盐体系为主,目前已发展到应用低温原位拉曼光谱技术对包裹体中的阳离子和阴离子进行定量-半定量分析。油气包裹体原位低温分析技术也取得了较大的突破,对不同成分油气包裹体低温下的相变过程取得了一定的认识,据此对油气包裹体进行分类,推测其主要成分,为油气包裹体计算提供基础参数。